Search

Molecular modeling techniques have become indispensable in many fields of molecular sciences in which the details related to mechanisms and reactivity need to be studied at an atomistic level. This review article provides a collection of computational modeling works on a topic of enormous interest and urgent relevance: the properties of metalloenzymes involved in the degradation and valorization of natural biopolymers and synthetic plastics on the basis of both circular biofuel production and bioremediation strategies. In particular, we will focus on lytic polysaccharide monooxygenase, laccases, and various heme peroxidases involved in the processing of polysaccharides, lignins, rubbers, and some synthetic polymers. Special attention will be dedicated to the interaction between these enzymes and their substrate studied at different levels of theory, starting from classical molecular docking and molecular dynamics techniques up to techniques based on quantum chemistry.

This perspective aims at illustrating a computational viewpoint on some specific issues concerning structure–activity relationships related to [FeFe]-hydrogenase ([FeFe]-H2ase) biomimicry. Most of the research outlined herein has been addressed by means of density functional theory (DFT) based tools, in some cases in conjunction with experimental techniques. The number of computational, experimental and mixed approach studies on the “hydrogenase models” topic is extraordinarily high. Accordingly, several comprehensive reviews have been published over the last twenty years. Moreover, the number of iron compounds of increasing nuclearity that have been considered as “biomimetic models” of [FeFe]-H2ase has grown exponentially over time. Additionally, the issues regarding some mismatches existing between the natural system and related synthetic analogues, are quite variegated. As a consequence of the countless examples that could be considered as related to “hydrogenase mimicry”, the intent of the present contribution is to provide an account of a limited number of recent study cases, in which DFT has been employed to elucidate redox properties, reactivity issues of selected examples of diiron biomimicry. Herein the treated topics have been explicitly chosen to show how DFT has allowed us to rationalize and complement experimental observations, with a special focus on electrocatalysis aspects. Finally, the specific purpose of this perspective is to show how mutually intertwined are the various issues concerning a fascinating branch of biomimetic research.

Some microorganisms, like the aerobic soil bacteria, Oligotropha carboxidovorans, have the capability to oxidize the highly toxic atmospheric gas carbon monoxide (CO) into CO2 through CO dehydrogenase enzymes, whose active site contains a bimetallic MoCu center. Over the last decades, a number of experimental and theoretical investigations were devoted to understanding the mechanism of CO oxidation and, in particular, the role of a very stable thiocarbonate intermediate that may be formed during the catalytic cycle. The occurrence of such an intermediate was reported to make the CO2 release step kinetically difficult. In this work, by using an accurate QM/MM approach and energy refinement by means of the BigQM method, we were able to determine the role of such an intermediate and propose a novel mechanism for the oxidation of CO into CO2 by Mo/Cu CO dehydrogenase. Surprisingly, we found that the detachment of CO2 occurs directly from the product of the Mo=O nucleophilic attack reaction on the carbon of CO aided by the transient coordination of the active site glutamate to the Mo ion. The estimated activation barrier is in good agreement with the experimental one, while the thiocarbonate turned out to not interfere with the CO-oxidation catalytic cycle. The results highlight the importance of the environmental effects in the assembly of the molecular model and in the choice of the computational protocol. Our accurate modeling of the enzyme also allowed us to exclude the involvement of a frustrated Lewis pair in the CO-oxidation mechanism, which has recently been suggested based on an analysis of structural and electronic features of synthetic mimics of the Mo/Cu CO dehydrogenase active site.

FeIFeI Fe2(S2C3H6)(CO)6(µ-CO) (1a–CO) and its FeIFeII cationic species (2a+–CO) are the simplest model of the CO-inhibited [FeFe] hydrogenase active site, which is known to undergo CO photolysis within a temperature-dependent process whose products and mechanism are still a matter of debate. Using density functional theory (DFT) and time-dependent density functional theory (TDDFT) computations, the ground state and low-lying excited-state potential energy surfaces (PESs) of 1a–CO and 2a+–CO have been explored aimed at elucidating the dynamics of the CO photolysis yielding Fe2(S2C3H6)(CO)6 (1a) and [Fe2(S2C3H6)(CO)6]+ (2a+), two simple models of the catalytic site of the enzyme. Two main results came out from these investigations. First, a–CO and 2a+–CO are both bound with respect to any CO dissociation with the lowest free energy barriers around 10 kcal mol−1, suggesting that at least 2a+–CO may be synthesized. Second, focusing on the cationic form, we found at least two clear excited-state channels along the PESs of 2a+–CO that are unbound with respect to equatorial CO dissociation.

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.

The MoCu CO dehydrogenase enzyme not only transforms CO into CO2 but it can also oxidise H2. Even if its hydrogenase activity has been known for decades, a debate is ongoing on the most plausible mode for the binding of H2 to the enzyme active site and the hydrogen oxidation mechanism. In the present work, we provide a new perspective on the MoCu-CODH hydrogenase activity by improving the in silico description of the enzyme. Energy refinement—by means of the BigQM approach—was performed on the intermediates involved in the dihydrogen oxidation catalysis reported in our previously published work (Rovaletti, et al. “Theoretical Insights into the Aerobic Hydrogenase Activity of Molybdenum–Copper CO Dehydrogenase.” Inorganics 7 (2019) 135). A suboptimal description of the H2–HN(backbone) interaction was observed when the van der Waals parameters described in previous literature for H2 were employed. Therefore, a new set of van der Waals parameters is developed here in order to better describe the hydrogen–backbone interaction. They give rise to improved binding modes of H2 in the active site of MoCu CO dehydrogenase. Implications of the resulting outcomes for a better understanding of hydrogen oxidation catalysis mechanisms are proposed and discussed.

International audience; FeFe hydrogenases catalyse H2 oxidation and production using a "H-cluster", where two Fe ions are bound by an aza-dithiolate (adt) ligand. Various hypotheses have been proposed (by us and others) to explain that the enzyme reversibly inactivates under oxidizing, anaerobic conditions: intramolecular binding of the N atom of adt, formation of the so-called Hox/inact state or non-productive binding of H2 to isomers of the H-cluster. Here we show that none of the above explains the new finding that the anaerobic, oxidative, H2-dependent reversible inactivation is strictly dependent on the presence of Cl- or Br-. We provide experimental evidence that chloride uncompetitively inhibits the enzyme: it reversibly binds to catalytic intermediates of H2 oxidation (but not to the resting "Hox" state), after which oxidation locks the active site into a stable, saturated, inactive form, the structure of which is proposed here based on DFT calculations. The halides interact with the amine group of the H-cluster but do not directly bind to iron. It should be possible to stabilize the inhibited state in amounts compatible with spectroscopic investigations to explore further this unexpected reactivity of the H-cluster of hydrogenase.

The synthesis and characterization of a hexanuclear cobalt complex 1 involving a non-heme ligand system, L1, supported on a Sn6O6 stannoxane core are reported. Complex 1 acts as a unique catalyst for dioxygen reduction, whose selectivity can be changed from a preferential 4e(-)/4H(+) dioxygen-reduction (to water) to a 2e(-)/2H(+) process (to hydrogen peroxide) only by increasing the temperature from -50 (o)C to 25 (o)C. A variety of spectroscopic methods ((119)Sn-NMR, Magnetic circular dichroism (MCD), electron paramagnetic resonance (EPR), SQUID, UV-Vis absorption, resonance Raman (rRaman), and X-ray absorption spectroscopy (XAS) ) coupled with advanced theoretical calculations has been applied for the unambiguous assignment of the geometric and electronic structure of 1. The mechanism of the O2-reduction reaction has been clarified based on kinetic studies on the overall catalytic reaction as well as each step in the catalytic cycle and by low-temperature detection of intermediates. The O2-binding to 1 results in the efficient formation of a stable end-on μ-1,2-peroxodicobalt(III) intermediate 2 at -50 (o)C, followed by a proton-coupled electron-transfer (PCET) reduction to complete the O(2)-to-2H2O cata-lytic conversion in an overall 4e(-)/4H(+) step. In contrast, at higher temperatures (> 20 oC) the constraints provided by the stannoxane core, makes 2 unstable against a preferential proton-transfer (PT) step, leading to the generation of H2O2 by a 2e(-)/2H(+) process. The present study provides deep mechanistic insight into the dioxygen reduction process that should serve as useful and broadly applicable principles for future design of more efficient catalysts in fuel cells

A bisphosphide-bridged diiron hexacarbonyl complex 3 with NEt2 groups on the phosphide bridge was synthesized to examine a new proton relay system from the NEt2 group to the bridging hydride between the two iron centers. As a precursor of the bridging moiety, peri-Et2NP-PNEt2-bridged naphthylene 5 was synthesized by the reaction of 1,8-dilithionaphthylene with two equivalents of Cl2PNEt2 followed by reductive P-P bond formation by magnesium. The reaction of the diphosphine ligand 5 with Fe2(CO)9 gave the diiron hexacarbonyl complex 3, in which the P-P bond of the ligand was cleaved to form the bisphosphide-bridge. The molecular structure of 3 indicated that the trigonal plane of the NEt2 group was forced to face the Fe-Fe bond to avoid steric congestion with the naphthylene group linking the two phosphide groups. The NEt2 group could be protonated by p-toluenesulfonic acid. Density functional theory (DFT) calculations confirmed that the proton of the N(H)Et2 group adopted a position close to the bridging hydride. The DFT results for the ferrocene analogue 1, in which the 1,8-naphthylene group of 3 was replaced with the 1,1′-ferrocenylene group, also revealed that the most stable orientation of the protonated NHEt2 group was that in the protonated 3. As a result, electrochemical proton reduction reactions using complexes 1 and 3 proceeded with similar catalytic efficiencies. Unfortunately, the catalytic efficiencies (CEs) of these complexes were much lower than those of the complexes with a proton relay system of the terminal hydrogen, indicating that the reactive properties of the bridging hydride in the present proton relay system cannot exceed those of the terminal hydride.

Carbon monoxide (CO) is a highly toxic gas to many living organisms. However, some microorganisms are able to use this molecule as the sole source of carbon and energy. Soil bacteria such as the aerobic Oligotropha carboxidovorans are responsible for the annual removal of about 2x108 tons of CO from the atmosphere. Detoxification through oxidation of CO to CO2 is enabled by the MoCu-dependent CO-dehydrogenase enzyme (MoCu-CODH) which-differently from other enzyme classes with similar function-retains its catalytic activity in the presence of atmospheric O2. In the last few years, targeted advancements have been described in the field of bioengineering and biomimetics, which is functional for future technological exploitation of the catalytic properties of MoCu-CODH and for the reproduction of its reactivity in synthetic complexes. Notably, a growing interest for the quantum chemical investigation of this enzyme has recently also emerged. This mini-review compiles the current knowledge of the MoCu-CODH catalytic cycle, with a specific focus on the outcomes of theoretical studies on this enzyme class. Rather controversial aspects from different theoretical studies will be highlighted, thus illustrating the challenges posed by this system as far as the application of density functional theory and hybrid quantum-classical methods are concerned.

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

The chemical groups present at the surface of graphite have been thought for a long time to be mainly responsible for its catalytic activity in the oxygen reduction reaction. Recently, it was proposed that the surface defects of graphite also significantly contribute to promote this reaction. Although the behaviour of surface defects has been reported, only few comments have been dedicated to their involvement in the mechanism and the possible intermediate species in the oxygen reduction reaction. Herein, we aim to present a more detailed explanation of the catalytic activity of graphite particles based on the structure of their defects and their size. Structural, spectroscopic and magnetic investigation (X-ray diffraction, Raman and electron spin resonance) and electrochemical measurements were performed to describe the nature of the defects and their aptitude to transfer electrons. Computational description supplied precise details of the energy of the different defects and their ability to promote the reduction, also suggesting the structure of the intermediate adduct in the oxygen reduction. The results indicated that molecular oxygen preferentially interacts with graphite defects, which involve the π-electron system and accumulation of the spin density on the edges of the grains, in particular, on the zig-zag edges present on ball-milled graphite. This promotes the reactivity of this nanomaterial. Furthermore, the activation increases by decreasing the particle size.

We present a theoretical investigation providing key insights on a long-standing controversial issue that dominated the debate on carbon monoxide oxidation by Mo–Cu CO-dehydrogenases. Previous investigations gravitate around the possible occurrence of a thiocarbonate intermediate, that was repeatedly reported to behave as a thermodynamic sink on the catalytic energy landscape. By using a hierarchy of quantum mechanical and hybrid quantum/classical models of the enzyme, we show that no such energy sink is present on the catalytic energy profile. Consequent perspectives for the definition of a novel mechanistic proposal for the enzyme-catalyzed CO-oxidation are discussed in light of the recent literature.

Catalytic H2 oxidation has been dissected by means of DFT into the key steps common to the Fe2 unit of both the [FeFe]-hydrogenase cofactor and selected biomimics. The aim was to elucidate the molecular details underlying the very different performances of the two systems. We found that the better enzyme performance is based on a single iron atom that is maintained electron-poor, favoring H2 binding, although embedded within a highly electron-rich cofactor, ensuring a facile oxidation of the Fe2 -H2 adduct. This is due to 1) CN- coordinating to both iron atoms, due to their amphipathic Lewis acid/base properties, and 2) the 4Fe4S subunit further withdrawing electrons from the Fe2 core. Preserving a moderate electron deficiency at a single iron also helps the cofactor preserve hydride affinity, which favors H2 cleavage. Such valuable characteristics allow the biocatalyst to turnover close to equilibrium conditions. All previous biomimicry has shown, in contrast, the impossibility to properly balance the two apparently contrasting aforementioned requisites, although evident progress has been made by the H2 -ase community. Disclosure of the differences identified could inspire the design of novel biomimics, for instance, reconsidering the use of CN- in the catalyst architecture. Indeed, in the presence of bases normally employed in oxidative catalysis, undesired stable protonation at coordinated CN- , which affects the opposite process (proton reduction), could be overcome.

DFT calculations were used to investigate the cycloaddition reaction of CO2 to styrene oxide for the formation of styrene carbonate. The uncatalyzed process alongside the reactions assisted by tetrabutylammonium bromide (TBAB), the novel nonsymmetrical single-center aluminum(III) salen-acac hybrid complex (salenac) (Al1cat) and the binary Al1cat/TBAB system were all investigated and for all of them the optimized structures, rate-determining steps, and lowest energy barrier reaction pathways were intercepted for both gas-phase and solvent environments. For the catalyzed systems, the reaction mechanism consists of three key elementary steps: 1) epoxy ring opening; 2) CO2 electrophilic attack and 3) intramolecular cyclization. In the presence of Al1cat, the central metal of the catalyst coordinates with an oxygen atom of the epoxide, activating it towards a nucleophilic attack by the halide. An oxyanion species is formed that affords the corresponding cyclic carbonate after reaction with CO2. Our results provide important hints on the cycloaddition of CO2 and epoxides promoted by nonsymmetrical aluminum complex containing a single metal center, and can satisfactorily explain the previous experimental observations allowing the development of more efficient catalysts for organic carbonate production.

Artifacts from the Egyptian grave goods of Kha and Merit preserved at the Museo Egizio in Turin were studied through a combination of non-destructive and non-invasive neutron and gamma techniques (namely neutron imaging, neutron diffraction and prompt gamma activation analysis). The results provide unprecedented morphological reconstructions of the inner parts of the two alabaster and metallic vases and their isotopic and phase composition, thereby extending our knowledge of the hitherto unknown content of the vases and their functions.

Environmentally persistent free radicals (EPFRs) are toxic products deriving from incomplete combustion and are able to generate DNA damage and pulmonary dysfunction. They are formed on particulate matter through interaction with aromatic hydrocarbons, catalyzed by transition metal oxides, and produce reactive oxygen species (ROS) in aquatic media. The processes are already described for substituted aromatic molecules, for example, phenol, but not for unsubstituted aromatic systems, such as benzene. This Article reports on the reaction of benzene with molecular oxygen in the presence of CuxO/SiO2, suggesting a mechanism based on cluster and periodic computational models. The activation of O2 by interaction with silica coordinated Cu(I) centers leads to a peroxy species that yields the phenoxy radical upon reaction with benzene. Dissociation of OH• radical eventually allows for the recovery of the catalyst. The experimental characterization of the CuxO/SiO2 catalyst regarded morphology, crystal structure, copper electronic state, and crystal field around Cu(II). Electron paramagnetic resonance (EPR) spectroscopy revealed the formation of phenoxy radical entrapped in the catalyst upon reaction between benzene and CuxO/SiO2. Moreover, EPR investigation of ROS in aqueous solution evidenced the generation of OH• radicals by benzene-contacted CuxO/SiO2. All of the experimental results nicely fit the outcomes of the computational models.

The redox properties of Ni complexes bound to a new ligand, [DMC-nit] +, where a N-heterocyclic nitrenium group is anchored on a 1,4,8,11-tetraazacyclotetradecane backbone, have been examined using spectroscopic and DFT methods. Ligand-based [(DMC-nit)Ni]2+/+ reduction and metal-based [(DMC-nit)Ni]2+/3+ oxidation processes have been established for the [(DMC-nit)Ni]+/2+/3+ redox series, which represents the first examples of nitrenium nitrogen (Nnit)-bound first-row transition-metal complexes. An unprecedented bent binding mode of Nnit in [(DMC-nit)Ni]2+ is observed, which possibly results from the absence of any Nnit→Ni σ-donation. For the corresponding [(DMC-nit)Ni(F)]2+ complex, σ-donation is dominant, and hence a coplanar arrangement at Nnit is predicted by DFT. The binding of the triazolium ion to Ni enables new chemistry (formate oxidation) that is not observed in a derivative that lacks this functional group. Thus the N-heterocyclic nitrenium ligand is a potentially useful and versatile reagent in transition-metal-based catalysis. © 2013 American Chemical Society.

International audience; [FeFe] hydrogenases are H2-evolving enzymes that feature a diiron cluster in their active site (the [2Fe]H cluster). One of the iron atoms has a vacant coordination site that directly interacts with H2, thus favoring its splitting in cooperation with the secondary amine group of a neighboring, flexible azadithiolate ligand. The vacant site is also the primary target of the inhibitor O2. The [2Fe]H cluster can span various redox states. The active-ready form (Hox) attains the Fe(II)Fe(I) state. States more oxidized than Hox were shown to be inactive and/or resistant to O2. In this work, we used density functional theory to evaluate whether azadithiolate-to-iron coordination is involved in oxidative inhibition and protection against O2, a hypothesis supported by recent results on biomimetic compounds. Our study shows that Fe-N(azadithiolate) bond formation is favored for an Fe(II)Fe(II) active-site model which disregards explicit treatment of the surrounding protein matrix, in line with the case of the corresponding Fe(II)Fe(II) synthetic system. However, the study of density functional theory models with explicit inclusion of the amino acid environment around the [2Fe]H cluster indicates that the protein matrix prevents the formation of such a bond. Our results suggest that mechanisms other than the binding of the azadithiolate nitrogen protect the active site from oxygen in the so-called H ox (inact) state.

International audience; FeFe hydrogenases catalyze H-2 oxidation and formation at an inorganic active site (the "H-cluster"), which consists of a [Fe-2(CO)(3)(CN)(2)(dithiomethylamine)] subcluster covalently attached to a Fe4S4 subcluster. This active site is photosensitive: visible light has been shown to induce the release of exogenous CO (a reversible inhibitor of the enzyme), shuffle the intrinsic CO ligands, and even destroy the H-cluster. These reactions must be understood because they may negatively impact the Use of hydrogenase for the photoproduction of H-2. Here, we explore in great detail the reactivity of the excited states of the H-duster under catalytic conditions by examining, both experimentally and using TDDFT calculations, the simplest photochemical reaction: the binding and release of exogenous CO. A simple dyad model can be used to predict which excitations are active. This could be used for probing other, aspects of the photoreactivity of the H-cluster.

High-valent copper-nitrene intermediates have long been proposed to play a role in copper-catalyzed aziridination and amination reactions. However, such intermediates have eluded detection for decades, preventing the unambiguous assignments of mechanisms. Moreover, the electronic structure of the proposed copper-nitrene intermediates has also been controversially discussed in the literature. These mechanistic questions and controversy have provided tremendous motivation to probe the accessibility and reactivity of CuIII-NR/ CuIINR species. In this paper, we report a breakthrough in this field that was achieved by trapping a transient copper-tosylnitrene species, 3-Sc, in the presence of scandium triflate. The sufficient stability of 3-Sc at -90 °C enabled its characterization with optical, resonance Raman, NMR, and X-ray absorption near-edge spectroscopies, which helped to establish its electronic structure as CuIINTs (Ts = tosyl group) and not Cu IIINTs. 3-Sc can initiate tosylamination of cyclohexane, thereby suggesting CuIINTs cores as viable reactants in oxidation catalysis. © 2012 American Chemical Society.

Amyloids result from the aggregation of several unrelated proteins, due to either specific mutations or promoting intra- or extra-cellular conditions. Structurally, they are rich in intermolecular β-sheets and are the causative agents of several diseases, both neurodegenerative and systemic. It is believed that the most toxic species are small aggregates, referred to as oligomers, rather than the final fibrillar assemblies. Their mechanisms of toxicity are mostly mediated by aberrant interactions with the cell membranes, with resulting derangement of membrane-related functions. Much effort is being put in the search for natural antiamyloid agents, and/or in the development of synthetic molecules. Actually, it is well documented that the prevention of amyloid aggregation results in several cytoprotective effects. Here, we portray the state of the art in the field. Several natural compounds are effective antiamyloid agents, notably tetracyclines and polyphenols. They are generally non-specific, as documented by their partially overlapping mechanisms, and the capability to interfere with the aggregation of several unrelated proteins. Among rationally designed molecules, we mention the prominent examples of β-breakers peptides, whole antibodies and fragments thereof, and the special case of drugs contrasting transthyretin aggregation. In this framework, we stress the pivotal role of the computational approaches. When combined with biophysical methods, in several cases they have helped clarify in detail the protein/drug modes of interaction, which make it plausible that more effective drugs will be developed in the future.

Nature's recipe: A theoretical study analyzes how the environment of the [FeFe] hydrogenase's catalytic cofactor affects its chemical properties, particularly the relative stability of complexes with bridging and terminal hydride ligands (see picture; Fe teal, S yellow, C green, N blue, O red, H gray). The results help to elucidate key rules for the design of bioinspired synthetic catalysts for H2 production. © 2009 Wiley-VCH Verlag GmbH & Co. KGaA.

The complex (S2C3H6)[Fe-2(CO)(5)P(NC4H8)(3)] (A), a recently synthesized functional model of the active site of Fe-hydrogenases, is able to electrocatalyze proton reduction, leading to molecular hydrogen evolution. Experimental results suggested that the presence of the electron donor P(NC4H8)(3) ligand in A could favor the formation of a mu-CO species similar to that observed in the enzymatic cluster. However, insight into the structural features of key catalytic intermediates deriving from reduction and protonation of A was still lacking. Here we present results obtained using density functional theory to evaluate structures, relative stabilities, and spectroscopic properties of several species relevant for the electrocatalytic H-2 evolving process. The results enabled us to unravel the structure of the mu-CO complex experimentally detected after monoelectronic reduction of A. Moreover, we show that the introduction of the large electron-donor ligand P(NC4H8)(3) in the biomimetic complex does not favour the stabilization of terminal-hydrido adducts, which are expected to be very reactive in terms of H-2 production. The comparison of our findings with previous theoretical and experimental results obtained on similar model complexes suggests that the introduction of an electron donor ligand as good as P(NC4H8)(3), but less sterically demanding, could represent a better choice to facilitate the formation of mu-CO complexes more closely resembling the structure of the enzymatic cluster. ((C) Wiley-VCH Verlaq GmbH & Co. KGaA, 69451 Weinheim, Germany, 2007).

Ad hoc quantum chemical modeling of the acetyl-CoA synthase local structure and folding allowed us to identify an unprecedented coordination mode of histidine sidechain to protein-embedded metal ions. This journal is

The di-iron complex Fe-2(S2C3H6)(CO)(6) (a), one of the simplest functional models of the Fe-hydrogenases active site, is able to electrocatalyze proton reduction. In the present study, the H-2 evolving path catalyzed by a has been characterized using density functional theory. It is showed that, in the early stages of the catalytic cycle, a neutral mu-H adduct is formed; monoelectron reduction and subsequent protonation can give rise to a diprotonated neutral species (a-mu H-SH), which is characterized by a mu-H group, a protonated sulfur atom, and a CO group bridging the two iron centers, in agreement with experimental IR data indicating the formation of a long-lived mu-CO species. H-2 release from a-mu H-SH, and its less stable isomer a-H-2 is kinetically unfavorable, while the corresponding monoanionic compounds (a-mu H-SH- and a-H-2(-)) are more reactive in terms of dihydrogen evolution, in agreement with experimental data. The key species involved in electrocatalysis have structural features different from the hypothetical intermediates recently proposed to be involved in the enzymatic process, an observation that is possibly correlated with the reduced catalytic efficiency of the biomimetic di-iron assembly.

The present paper reports a study on the energetics of protonation of a hydrogenase biomimetic complex, [Fe2(μ-adt)(CO)4(PMe3)2] (adt = N-benzylazadithiolate), and of its homologue featuring triphenylphosphine ligands in place of trimethylphosphines. Formation of a terminal hydride on one of the Fe centres was considered first, given the key relevance of terminal hydride species in the enzymatic mechanism. Theoretical calculations highlight that, in a vacuum, terminal protonation of the selected Fe ion in the PPh3-bearing organometallic complex is highly favoured when compared to the analogous reaction involving the PMe3-containing species, but the trend is inverted in the case of models optimized in a continuum polarizable. An unexpected parallel is thus established between relative basicities of PPh3 and PMe3 in vacuum or in solution phase [C. A. Tolman, J. Am. Chem. Soc. 1970, 92, 2953; G. M. Bancroft, Inorg. Chem. 1986, 25, 3675], and the energetics of terminal hydride formation upon protonation of [FeFe]-hydrogenase biomimetic complexes bearing such organophosphorous ligands. Bridging hydride formation was also considered in the present study: calculations showed that protonation of the PMe3-bearing organometallic complex is again strongly favoured in vacuo, as compared to the case of the PPh3-containing model. However, protonation energies become significantly smaller when solvent effects are taken into account. Such differences between protonation reactions modelled in vacuo and in the polarizable continuum are rationalized in light of the different electrostatic properties of the diiron complexes here considered. Implications for the design and modelling of biomimetic catalysts are briefly discussed in light of recent literature

Proton-coupled electron transfer (PCET) is a fundamental process at the core of oxidation-reduction reactions for energy conversion. The [FeFe]-hydrogenases catalyze the reversible activation of molecular H2 through a unique metallocofactor, the H-cluster, which is finely tuned by the surrounding protein environment to undergo fast PCET transitions. The correlation of electronic and structural transitions at the H-cluster with proton-transfer (PT) steps has not been well-resolved experimentally. Here, we explore how modification of the conserved PT network via a Cys → Ser substitution at position 169 proximal to the H-cluster of Chlamydomonas reinhardtii [FeFe]-hydrogenase (CrHydA1) affects the H-cluster using electron paramagnetic resonance (EPR) and Fourier transform infrared (FTIR) spectroscopy. Despite a substantial decrease in catalytic activity, the EPR and FTIR spectra reveal different H-cluster catalytic states under reducing and oxidizing conditions. Under H2 or sodium dithionite reductive treatments, the EPR spectra show signals that are consistent with a reduced [4Fe-4S]H(+) subcluster. The FTIR spectra showed upshifts of νCO modes to energies that are consistent with an increase in oxidation state of the [2Fe]H subcluster, which was corroborated by DFT analysis. In contrast to the case for wild-type CrHydA1, spectra associated with Hred and Hsred states are less populated in the Cys → Ser variant, demonstrating that the exchange of -SH with -OH alters how the H-cluster equilibrates among different reduced states of the catalytic cycle under steady-state conditions.

Nature is a valuable source of inspiration in the design of catalysts, and various approaches are used to elucidate the mechanism of hydrogenases, the enzymes that oxidize or produce H 2. In FeFe hydrogenases, H 2 oxidation occurs at the H-cluster, and catalysis involves H 2 binding on the vacant coordination site of an iron centre. Here, we show that the reversible oxidative inactivation of this enzyme results from the binding of H 2 to coordination positions that are normally blocked by intrinsic CO ligands. This flexibility of the coordination sphere around the reactive iron centre confers on the enzyme the ability to avoid harmful reactions under oxidizing conditions, including exposure to O 2. The versatile chemistry of the diiron cluster in the natural system might inspire the design of novel synthetic catalysts for H 2 oxidation. © 2014 Macmillan Publishers Limited.

[FeFe]-hydrogenases harbor a {2Fe3S} assembly bearing two CO and two CN- groups, a mu-CO ligand, and a vacant coordination site trans to the mu-CO group. Recent theoretical results obtained studying the isolated {2Fe3S} subsite indicated that one of the CN- ligands can easily move from the crystallographic position to the coordination site trans to the mu-CO group; such an isomerization would have a major impact on substrates and inhibitors binding regiochemistry and, consequently, on the catalytic mechanism. To shed light on this crucial issue, we have carried out hybrid QM/MM and free energy perturbation calculations on the whole enzyme, which demonstrate that the protein environment plays a crucial role and maintains the CN- group fixed in the position observed in the crystal structure; these results strongly support the hypothesis that the vacant coordination site trans to the mu-CO group has a crucial functional relevance both in the context of CO-mediated inhibition of the enzyme and in dihydrogen oxidation/evolution catalysis.

Rotated geometry: The first example of a dinuclear iron(I)-iron(I) complex featuring a fully rotated geometry related to the active site of [Fe-Fe] hydrogenase is reported (see figure). Copyright © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Key stereoelectronic properties of Desulfovibrio desulfuricans [FeFe]-hydrogenase (DdH) were investigated by quantum mechanical description of its complete inorganic core, which includes a Fe6S6 active site (the H-cluster), as well as two ancillary Fe4S4 assemblies (the F and F' clusters). The partially oxidized, active-ready form of DdH is able to efficiently bind dihydrogen, thus starting H2 oxidation catalysis. The calculations allow us to unambiguously assign a mixed Fe(II)Fe(I) state to the catalytic core of the activeready enzyme and show that H2 uptake exerts subtle, yet crucial influences on the redox properties of DdH. In fact, H2 binding can promote electron transfer from the H-cluster to the solvent-exposed F'-cluster, thanks to a 50% decrease of the energy gap between the HOMO (that is localized on the H-cluster) and the LUMO (which is centered on the F'-cluster). Our results also indicate that the binding of the redox partners of DdH in proximity of its F'-cluster can trigger one-electron oxidation of the H2-bound enzyme, a process that is expected to have an important role in H2 activation. Our findings are analyzed not only from a mechanistic perspective, but also in consideration of the physiological role of DdH. In fact, this enzyme is known to be able to catalyze both the oxidation and the evolution of H2, depending on the cellular metabolic requirements. Hints for the design of targeted mutations that could lead to the enhancement of the oxidizing properties of DdH are proposed and discussed. © 2011 American Chemical Society.

International audience; Carbon monoxide is often described as a competitive inhibitor of FeFe hydrogenases, and it is used for probing H-2 binding to synthetic or in silico models of the active site H-cluster. Yet it does not always behave as a simple inhibitor. Using an original approach which combines accurate electrochemical measurements and theoretical calculations, we elucidate the mechanism by which, under certain conditions, CO binding can cause permanent damage to the H-cluster. Like in the case of oxygen inhibition, the reaction with CO engages the entire H-cluster, rather than only the Fe-2 subsite.

Recent advances aimed at modeling the chemistry of the active site of [FeFe]-hydrogenases (the H-cluster, composed by a catalytic Fe 2S 2 subcluster and an Fe 4S 4 portion) have led to the synthesis of binuclear coordination compounds containing a noninnocent organophosphine ligand [2,3-bis(diphenylphosphino)maleic anhydride, bma] that is able to undergo monoelectron reduction, analogously to the tetranuclear Fe 4S 4 subcluster portion of the H-cluster. However, such a synthetic model was shown to feature negligible electronic communication between the noninnocent ligand and the remaining portion of the cluster, at variance with the enzyme active site. Here, we report a theoretical investigation that shows why the electron transfer observed in the enzyme upon protonation of the catalytic Fe 2S 2 subsite cannot take place in the bma-containing cluster. In addition, we show that targeted modifications of the bma ligand are sufficient to restore the electronic communication within the model, such that electron density can be more easily withdrawn from the noninnocent ligand, as a result of protonation of the iron centers. Similar results were also obtained with a ligand derived from cobaltocene. The relevance of our findings is discussed from the perspective of biomimetic reproduction of proton reduction to yield molecular hydrogen. © 2011 American Chemical Society.

The presence of Fe-bound cyanide ligands in the active site of the proton-reducing enzymes [FeFe]-hydrogenases has led to the hypothesis that such Brønsted-Lowry bases could be protonated during the catalytic cycle, thus implying that hydrogen isocyanide (HNC) might have a relevant role in such crucial microbial metabolic paths. We present a hybrid quantum mechanical/molecular mechanical (QM/MM) study of the energetics of CN - protonation in the enzyme, and of the effects that cyanide protonation can have on [FeFe]-hydrogenase active sites. A detailed analysis of the electronic properties of the models and of the energy profile associated with H2 evolution clearly shows that such protonation is dysfunctional for the catalytic process. However, the inclusion of the protein matrix surrounding the active site in our QM/MM models allowed us to demonstrate that the amino acid environment was finely selected through evolution, specifically to lower the Brønsted-Lowry basicity of the cyanide ligands. In fact, the conserved hydrogen-bonding network formed by these ligands and the neighboring amino acid residues is able to impede CN- protonation, as shown by the fact that the isocyanide forms of [FeFe]-hydrogenases do not correspond to stationary points on the enzyme QM/MM potential-energy surface. © 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

[FeFe]-hydrogenases are highly efficient H(2)-evolving metalloenzymes that include cyanides and carbonyls in the active site. The latter is an Fe(6)S(6) cluster (the so-called H-cluster) that can be subdivided into a binuclear portion carrying the CO and CN(-) groups and a tetranuclear subcluster. The fundamental role of cyanide ligands in increasing the basicity of the H-cluster has been highlighted previously. Here a more subtle but crucial role played by the two CN(-) ligands in the active site of [FeFe]-hydrogenases is disclosed. In fact, QM/MM calculations on all-atom models of the enzyme from Desulfovibrio desulfuricans show that the cyanide groups fine-tune the electronic and redox properties of the active site, affecting both the protonation regiochemistry and electron transfer between the two subclusters of the H-cluster. Despite the crucial role of cyanides in the protein active site, the currently available bioinspired electrocatalysts generally lack CN(-) groups in order to avoid competition between the latter and the catalytic metal centers for proton binding. In this respect, we show that a targeted inclusion of phosphine ligands in hexanuclear biomimetic clusters may restore the electronic and redox features of the wild-type H-cluster.

A DFT study of protonation thermodynamics in H(2)-evolving biomimetic catalysts related to [FeFe]-hydrogenases active site is presented here. Taking as a reference system the electrocatalytic dihydrogen evolution mechanism recently proposed for the synthetic assembly [Fe(2)(CO)(4)(kappa(2)-Ph(2)PCH(2)CH(2)PPh(2))(mu-S(CH(2))(3)S)] (a, which is able to release H(2) after having undergone monoelectron reduction steps and three sequential protonation reactions), we show how the reduction of model complexes to oxidation states lower than those observed in [FeFe]-hydrogenases cofactor leads to a protonation regiochemistry that has no counterpart in the enzymatic mechanism of H(2) production. In particular, double protonation of the metal centers turned out to be disfavored in a by up to 12.5 kcal mol(-1) with respect to alternative protonation paths; as for the regiochemistry of triple protonation, the formation of eta(2)-H(2) adducts is disfavored by at least approximately 25 kcal mol(-1). Structural analysis of the theoretical models also revealed that over-reduction of synthetic complexes, though necessary for observing H(2) evolution from the currently available biomimetic electrocatalysts, can generally impair their structural integrity. Possible approaches for the modulation of protonation regiochemistry are then proposed; in particular, it turned out that a targeted use of sigma-donating ligands showing low basicity can favor double protonation of iron centers.

Fe-only hydrogenases are enzymes that catalyze dihydrogen production or oxidation, due to the presence of an unusual Fe(6)S(6) cluster (the so-called H-cluster) in their active site, which is composed of a Fe(2)S(2) subsite, directly involved in catalysis, and a classical Fe(4)S(4) cubane cluster. Here, we present a hybrid quantum mechanical and molecular mechanical (QM/MM) investigation of the Fe-only hydrogenase from Desulfovibrio desulfuricans, in order to unravel key issues regarding the activation of the enzyme from its completely oxidized inactive state (Hoxinact) and the influence of the protein environment on the structural and catalytic properties of the H-cluster. Our results show that the Fe(2)S(2) subcluster in the Fe(II)Fe(II) redox state - which is experimentally observed for the completely oxidized form of the enzyme - binds a water molecule to one of its metal centers. The computed QM/MM energy values for water binding to the diferrous subsite are in fact over 70 kJ mol(-1); however, the affinity toward water decreases by 1 order of magnitude after a one-electron reduction of H(ox)(inact), thus leading to the release of coordinated water from the H-cluster. The investigation of a catalytic cycle of the Fe-only hydrogenase that implies formation of a terminal hydride ion and a di(thiomethyl)amine (DTMA) molecule acting as an acid/base catalyst indicates that all steps have reasonable reaction energies and that the influence of the protein on the thermodynamic profile of H(2) production catalysis is not negligible. QM/MM results show that the interactions between the Fe(2)S(2) subsite and the protein environment could give place to structural rearrangements of the H-cluster functional for catalysis, provided that the bidentate ligand that bridges the iron atoms in the binuclear subsite is actually a DTMA residue.

The double-histone fold is a rare protein fold in which two consecutive regions characterized by the typical structure of histones assemble together, thus giving a histone pseudodimer. Previously, this fold was found in a few prokaryotic histones and in the regulatory region of guanine-nucleotide exchange factors of the Sos family. Standard methods of sequence comparison did not allow us to find new proteins containing a histone pseudodimer, as previously reported (Sondermann et al. 2003). However, a deeper investigation of protein sequences showed that the two histone folds included in Sos proteins share significant sequence similarity with nucleosomal histones. On the basis of this observation, we applied a specific strategy of sequence-homology search, which led to the identification of a new group of histone pseudodimers in Cca3 and proteins similar to Cca3 (Cca3S). A homology model of the histone pseudodimer included in rat Cca3 was constructed. A subsequent structure-function relationship study revealed that the histone pseudodimers included in Cca3 and Cca3S proteins, but not those present in Sos proteins, could retain the ability of mediating protein-DNA interactions, and could consequently act as DNA-binding modules. © Springer-Verlag 2005.