Search

This work originated from the need to functionalize surfactant-coated inorganic nanoparticles for biomedical applications, a process that is limited by excess unbound surfactant. These limitations are connected to the bioconjugation of targeting molecules that are often in equilibrium between the free aliquot in solution and that which binds the surface of the nanoparticles. The excess in solution can play a role in the biocompatability in vitro and in vivo of the final nanoparticles stock. For this purpose, we tested the ability of common surfactants - monothiolated polyethylene glycol and amphiphilic polymers - to colloidally stabilize nanoparticles as excess surfactant is removed and compared them to newly appearing multidentate surfactants endowed with high avidity for inorganic nanoparticles. Our results showed that monothiolated polyethylene glycol or amphiphilic polymers have an insufficient affinity to the nanoparticles and as the excess surfactant is removed the colloidal stability is lost, while multidentate high-avidity surfactants excel in the same regard, possibly allowing improvement in an array of nanoparticle applications, especially in those stated.

In endemic regions (West Africa and the Congo Basin), the genetic diversity of monkeypox virus (MPXV) is geographically structured into two major clades (clade I and II) that differ in virulence and host associations. Clade IIb is closely related to the B.1 lineage, which is dominating a worldwide outbreak initiated in 2022. Lineage B.1 has however accumulated mutations of unknown significance that most likely result from APOBEC3 (apolipoprotein B mRNA editing catalytic polypeptide-like 3) editing. We applied a population genetics-phylogenetics approach to investigate the evolution of MPXV during historical viral spread in Africa and to infer the distribution of fitness effects. We observed a high preponderance of codons evolving under strong purifying selection, particularly in viral genes involved in morphogenesis and replication or transcription. However, signals of positive selection were also detected and were enriched in genes involved in immunomodulation and/or virulence. In particular, several genes showing evidence of positive selection were found to hijack different steps of the cellular pathway that senses cytosolic DNA. Also, a few selected sites in genes that are not directly involved in immunomodulation are suggestive of antibody escape or other immuno-mediated pressures. Because orthopoxvirus host range is primarily determined by the interaction with the host immune system, we suggest that the positive selection signals represent signatures of host adaptation and contribute to the different virulence of clade I and II MPXV. We also used the calculated selection coefficients to infer the effects of mutations that define the predominant hMPXV1 lineage B.1, as well as the changes that have been accumulating during the worldwide outbreak. Results indicated that a proportion of deleterious mutations were purged from the predominant outbreak lineage, whose spread was not driven by the presence of beneficial changes. Polymorphic mutations with a predicted beneficial effect on fitness are few and have a low frequency. It remains to be determined whether they have any significance for ongoing virus evolution.

Three hexacarbonyl diiron dithiolate complexes [Fe2(CO)6(μ-(SCH2)2X)] with different substituted bridgeheads (X=CH2, CEt2, CBn2 (Bn=CH2C6H5)), have been studied under the same experimental conditions by cyclic voltammetry in dichloromethane [NBu4][PF6] 0.2 M. DFT calculations were performed to rationalize the mechanism of reduction of these compounds. The three complexes undergo a two-electron transfer whose the mechanism depends on the bulkiness of the dithiolate bridge, which involves a different timing of the structural changes (Fe−S bond cleavage, inversion of conformation and CO bridging) vs redox steps. The introduction of a bulky group in the dithiolate linker has obviously an effect on normally ordered (as for propanedithiolate (pdt)) or inverted (pdtEt2, pdtBn2) reduction potentials. Et→Bn replacement is not theoretically predicted to alter the geometry and energy of the most stable mono-reduced and bi-reduced forms but such a replacement alters the kinetics of the electron transfer vs the structural changes.

De novo metalloprotein design is a remarkable approach to shape protein scaffolds toward specific functions. Here, we report the design and characterization of Due Rame 1 (DR1), a de novo designed protein housing a di-copper site and mimicking the Type 3 (T3) copper-containing polyphenol oxidases (PPOs). To achieve this goal, we hierarchically designed the first and the second di-metal coordination spheres to engineer the di-copper site into a simple four-helix bundle scaffold. Spectroscopic, thermodynamic, and functional characterization revealed that DR1 recapitulates the T3 copper site, supporting different copper redox states, and being active in the O 2 -dependent oxidation of catechols to o -quinones. Careful design of the residues lining the substrate access site endows DR1 with substrate recognition, as revealed by Hammet analysis and computational studies on substituted catechols. This study represents a premier example in the construction of a functional T3 copper site into a designed four-helix bundle protein.

Nitrogen gas is a highly inert molecule and its activation under mild conditions represents a crucial goal in current research. In a recent study, the discovery of low-valence Ca(i) compounds capable of coordinating and reducing N2 was reported [B. Rösch, T. X. Gentner, J. Langer, C. Färber, J. Eyselein, L. Zhao, C. Ding, G. Frenking and S. Harder, Science, 2021, 371, 1125]. The study of low-valence alkaline earth complexes represents a new horizon in inorganic chemistry and demonstrates examples of spectacular reactivity. For example, complexes of the [BDI]2Mg2 type are selective reducing reagents in both organic and inorganic synthesis reactions. To date, however, no activity of Mg(i) complexes in the activation of the nitrogen molecule has been reported. By computational studies, in the present work, we investigated the analogies and differences of low-valence Ca(i) and Mg(i) complexes in the coordination, activation and protonation of N2. We have shown that the possibility of alkaline earth metals to employ atomic orbitals of the d type is reflected in the differences in the N2 binding energy and its coordination mode (end-on vs. side-on), as well as in the spin state of the resulting adduct (singlet vs. triplet). These divergences are finally observed in the subsequent protonation reaction, which turned out to be prohibitive in the presence of Mg.

Oxidative stress and copper dyshomeostasis are well established factors in the course of Alzheimer's disease (AD). Indeed, the pathological interaction of copper with the amyloid peptide is able to induce reactive oxygen species (ROS) production through a three-step cycle in which O2 is gradually reduced to superoxide, oxygen peroxide and finally to highly neurotoxic OH⋅ radicals. The purpose of this study has been to investigate at a computational level, by Density Functional Theory (DFT) and molecular mechanics techniques, the final two steps of O2 reduction to OH− +OH radical i.e the superoxide and hydroperoxide reduction by ascorbate, on different Cu(II) ⋅ Aβ coordination models. At variance with the first step of O2 reduction to O2.− characterized by high energy barrier, the successive O2.− and O2H− reductions are low energy barriers processes and that the OH⋅ radical formation is further favored by the possibility of Cu(II) ⋅ Aβ to form the Asp1 amino terminal radical. We discuss the results obtained from the perspective of the recent findings in the Aβ redox chemistry in AD pathology.

Coronaviruses (CoVs) have complex genomes that encode a fixed array of structural and nonstructural components, as well as a variety of accessory proteins that differ even among closely related viruses. Accessory proteins often play a role in the suppression of immune responses and may represent virulence factors. Despite their relevance for CoV phenotypic variability, information on accessory proteins is fragmentary. We applied a systematic approach based on homology detection to create a comprehensive catalogue of accessory proteins encoded by CoVs. Our analyses grouped accessory proteins into 379 orthogroups and 12 super-groups. No orthogroup was shared by the four CoV genera and very few were present in all or most viruses in the same genus, reflecting the dynamic evolution of CoV genomes. We observed differences in the distribution of accessory proteins in CoV genera. Alphacoronaviruses harboured the largest diversity of accessory open reading frames (ORFs), deltacoronaviruses the smallest. However, the average number of accessory proteins per genome was highest in betacoronaviruses. Analysis of the evolutionary history of some orthogroups indicated that the different CoV genera adopted similar evolutionary strategies. Thus, alphacoronaviruses and betacoronaviruses acquired phosphodiesterases and spike-like accessory proteins independently, whereas horizontal gene transfer from reoviruses endowed betacoronaviruses and deltacoronaviruses with fusion-associated small transmembrane (FAST) proteins. Finally, analysis of accessory ORFs in annotated CoV genomes indicated ambiguity in their naming. This complicates cross-communication among researchers and hinders automated searches of large data sets (e.g., PubMed, GenBank). We suggest that orthogroup membership is used together with a naming system to provide information on protein function.

The electrochemical reduction of complexes [Fe2 (CO)4 (κ2 -phen)(μ-xdt)] (phen=1,10-phenanthroline; xdt=pdt (1), adtiPr (2)) in MeCN-[Bu4 N][PF6 ] 0.2 m is described as a two-reduction process. DFT calculations show that 1 and its monoreduced form 1- display metal- and phenanthroline-centered frontier orbitals (LUMO and SOMO) indicating the non-innocence of the phenanthroline ligand. Two energetically close geometries were found for the doubly reduced species suggesting an intriguing influence of the phenanthroline ligand leading to the cleavage of a Fe-S bond as proposed generally for this type of complex or retaining the electron density and avoiding Fe-S cleavage. Extension of calculations to other complexes with edt, adtiPr bridge and even virtual species [Fe2 (CO)4 (κ2 -phen)(μ-adtR )] (R=CH(CF3 )2 , H) or [Fe2 (CO)4 (κ2 -phen)(μ-pdtR )] (R=CH(CF3 )2 , iPr) showed that the relative stability between both two-electron-reduced isomers depends on the nature of the bridge and the possibility to establish a remote anagostic interaction between the iron center {Fe(CO)3 } and the group carried by the bridged-head atom of the dithiolate group.

Ion pairing in water solutions alters both the water hydrogen-bond network and ion solvation, modifying the dynamics and properties of electrolyte water solutions. Here, we report an anomalous intrinsic fluorescence of KCl aqueous solution at room temperature and show that its intensity increases with the salt concentration. From the ab initio density functional theory (DFT) and time-dependent DFT modeling, we propose that the fluorescence emission could originate from the stiffening of the hydrogen bond network in the hydration shell of solvated ion-pairs that suppresses the fast nonradiative decay and allows the slower radiative channel to become a possible decay pathway. Because computations suggest that the fluorophores are the local ion-water structures present in the prenucleation phase, this band could be the signature of the incoming salt precipitation.

The electron rich Fe2(pdt)(RNC)6 1 (pdt=CH2(CH2S−)2) has been used as starting point for a DFT multi-functional study to assess the feasibility of designing a stable inverted (or rotated) disposition of the two FeL3 pyramidal moieties in the dimetallic core, a key feature of [FeFe]-hydrogenase cofactor. The choice of 1 was motivated by the presence of a rotated form in solution, slightly less stable than the unrotated stereoisomer. Aimed to find an upgraded version of 1, featuring the rotated isomer as ground state, various combinations of factors have been tested, for their effect on the relative stability of rotated vs unrotated isomers. The general result is that combining coordination asymmetry, electron donor presence and isocyanides R substituents able to establish intramolecular interactions is effective in stabilizing the rotated isomer. Our DFT study may inspire the design of synthetic biomimetics, with improved resemblance to the natural system.

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.

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.

The unique properties of liquid water mainly arise from its hydrogen bond network. The geometry and dynamics of this network play a key role in shaping the characteristics of soft matter, from simple solutions to biosystems. Here we report an anomalous intrinsic fluorescence of HCl and NaOH aqueous solutions at room temperature that shows important differences in the excitation and emission bands between the two solutes. From ab initio time-dependent density functional theory modeling we propose that fluorescence emission could originate from hydrated ion species contained in transient cavities of the bulk solvent. These cavities, which are characterized by a stiff surface, could provide an environment that, upon trapping the excited state, suppresses the fast nonradiative decay and allows the slower radiative channel to become a possible decay pathway.

Hydrogenases are enzymes that catalyze the oxidation and production of molecular hydrogen. For about fifteen years, there have been many reports about the successful connection of these enzymes to photosensitizers with the aim of designing H2 photoproduction systems, but relatively little attention has been paid to whether and why illumination may affect the catalytic properties of the enzyme. In all hydrogenases, hydrogen activation occurs at an inorganic active site that includes at least one Fe-carbonyl motif, which may make it sensitive to irradiation. Here we review the evidence that hydrogenases are indeed photosensitive. We focus mainly on the so-called FeFe hydrogenases; their active site, called the H-cluster, consists of a [4Fe4S] cluster that is bound by a cysteine sulfur to a diiron site. The iron atoms of the binuclear cluster are coordinated by carbonyl and cyanide ligands and an azadithiolate group. We describe the effects of UV-visible light irradiation on the enzyme under cryogenic or turnover conditions and the photoreactivity of model complexes that mimic the diiron site. We emphasize the dependence of the photochemical processes on wavelength, and warn about FeFe hydrogenase photoinhibition, which should probably be considered when attempts are made to use FeFe hydrogenases for the artificial photosynthesis of solar fuels. We also underline the relevance of studies of synthetic mimics of the H-cluster for understanding at atomistic level the photochemical processes observed in the enzyme.

Four coronaviruses (HCoV-OC43, HCoV-HKU1, HCoV-NL63, and HCoV-229E) are endemic in human populations. All these viruses are seasonal and generate short-term immunity. Like the highly pathogenic coronaviruses, the endemic coronaviruses have zoonotic origins. Thus, understanding the evolutionary dynamics of these human viruses might provide insight into the future trajectories of SARS-CoV-2 evolution. Because the zoonotic sources of HCoV-OC43 and HCoV-229E are known, we applied a population genetics–phylogenetic approach to investigate which selective events accompanied the divergence of these viruses from the animal ones. Results indicated that positive selection drove the evolution of some accessory proteins, as well as of the membrane proteins. However, the spike proteins of both viruses and the hemagglutinin-esterase (HE) of HCoV-OC43 represented the major selection targets. Specifically, for both viruses, most positively selected sites map to the receptor-binding domains (RBDs) and are polymorphic. Molecular dating for the HCoV-229E spike protein indicated that RBD Classes I, II, III, and IV emerged 3–9 years apart. However, since the appearance of Class V (with much higher binding affinity), around 25 years ago, limited genetic diversity accumulated in the RBD. These different time intervals are not fully consistent with the hypothesis that HCoV-229E spike evolution was driven by antigenic drift. An alternative, not mutually exclusive possibility is that strains with higher affinity for the cellular receptor have out-competed strains with lower affinity. The evolution of the HCoV-OC43 spike protein was also suggested to undergo antigenic drift. However, we also found abundant signals of positive selection in HE. Whereas such signals might result from antigenic drift, as well, previous data showing co-evolution of the spike protein with HE suggest that optimization for human cell infection also drove the evolution of this virus. These data provide insight into the possible trajectories of SARS-CoV-2 evolution, especially in case the virus should become endemic.

International audience; The synthesis, characterization, and protonation of [Fe2(CO)6{(μ-SCH2)2(Et)P═O}] (1) using the moderately strong acid CF3CO2H (pKaMeCN = 12.7) are reported. Digital simulations of the cyclic voltammetry of 1 in the presence of CF3CO2H and DFT calculations have allowed us to obtain a detailed mechanistic picture of the processes underlying the catalytic hydrogen evolution reaction (HER) that 1 can mediate. Moreover, DFT has shed light on the role of the P═O functionality in the whole catalytic cycle of proton reduction. The reductive behavior of 1 features a double electron transfer with potential inversion, which is associated with deep structural rearrangement of the catalyst. The double reduction appears also functional to the intramolecular proton transfer from the P═O group to the diiron core, a crucial process for the H+/H– heterocoupling yielding H2. The key intermediate for the H2 formation and release is predicted to be a 3H+/3e– species, in which P═O is perfectly poised to shuttle protons from solution to the Fe–H–Fe moiety. Therefore, the R-P═O bridgehead installed in a dithiolato linker of a diiron core proves a valid and versatile alternative to the natural nitrogen-based Fe2 strap.

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.

Understanding the evolutionary dynamics of the four human endemic coronaviruses might provide insight into the future trajectories of SARS-CoV-2 evolution. We re-assessed the timing of endemic coronavirus emergence and we show that all viruses entered human populations in a time-frame ranging from ~500 to 55 years ago. Because the three highly pathogenic coronaviruses (SARS-CoV, MERS-CoV, and SARS-CoV-2) spilled-over in tight temporal succession, the pattern of coronavirus emergence, in analogy to that of influenza pandemics, is highly irregular. To contextualize this observation in a wider perspective of viral diseases emergence since 1945, we mined epidemiology database information. After controlling for reporting bias, we find that, contrary to widespread beliefs, the occurrence of viral diseases (either zoonotic or not) has not increased over the last decades. Analysis of the recent and ongoing evolution of HCoV-229E and HCoV-OC43 indicated that positive selection most likely contributed to fine-tune the interaction with the human interferon/inflammatory response. Conversely, integration of evolutionary inference and molecular dating provided evidence that these viruses are not undergoing antigenic drift and the temporal emergence of spike protein variants is best explained by optimization of receptor binding affinity. These data provide a fresh look on viral disease emergence and on coronavirus evolution.

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.

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., Density functional theory was used to investigate Ni,Fe-containing carbon monoxide dehydrogenase enzymes. Different coordination modes of the substrate CO2 to several forms and redox states of the C-cluster—the enzyme active site—were considered. The obtained results highlight the key role of the protein environment in tuning the reactivity and the geometry of the C-cluster. This helps to uncover which stereoelectronic properties of the active site are crucial for the efficient binding and release of CO2.

Oxidative stress and metal dyshomeostasis are considered as crucial factors in the pathogenesis of Alzheimer's disease (AD). Indeed, transition metal ions such as Cu(ii) can generate Reactive Oxygen Species (ROS) via O2 Fenton-like reduction, catalyzed by Cu(ii) coordinated to the Amyloid beta (Aβ) peptide. Despite intensive effort, the mechanisms of ROS-induced molecular damage remain poorly understood. In the present paper, we investigate on the basis of molecular modelling computations the mechanism of OH radical propagation toward the Aβ peptide, starting from the end-product of OH radical generation by Cu(ii)·Aβ. We evaluate (i) the OH oxidative capacity, as well as the energetics of the possible Aβ oxidation target residues, by quantum chemistry Density Functional Theory (DFT) on coordination models of Cu(ii)/OH/Aβ and (ii) the motion of the OH˙ approaching the Aβ target residues by classical Molecular Dynamics (MD) on the full peptide Cu(ii)/OH/Aβ(1–16). The results show that the oxidative capacity of OH coordinated Cu(ii)Aβ is significantly lower than that of the free OH radical and that propagation toward Aβ Asp and His residues is favoured over Tyr residues. These results are discussed on the basis of the recent literature on in vitro Aβ metal-catalyzed oxidation and on the possible implications for the AD oxidative stress mechanism.

A few bio-inorganic systems have been studied by means of modern quantum chemistry (QM) methods, such as the Density Functional Theory (DFT) in connection with a large set of basis functions. At this level, DFT has proven to give reliable (semi-quantitative) description of several properties of transition metals complexes which are either relevant part of a natural enzyme or ad hoc built to mimic the behaviour of enzymes. The first example considered is the N,N'-ethylenebis(acetylacetonatiminato)Co(II) (Coacacen), a very famous synthetic oxygen carrier. We report results of a DFT investigation of Coacacen in its free form as well as in penta-coordinated form with a suitable axial ligand and, finally, in hexa-coordinated form bound to dioxygen. The problems related to the electronic structure of Coacacen derivatives are very old, and have been investigated (mainly experimentally) for several years. Coacacen has been the first bio-inorganic system investigated in the group of Renato Ugo in the 1970s. In addition, a short review will be presented concerning the “modern” problems of the bio-inorganics, as they are presently investigated in our research group based at the University of Milano-Bicocca (Italy): the catalytic cycle characteristic of Fe-Fe hydrogenases and the model compounds describing the interaction between copper and amyloid β-peptides.

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 Fe

The electrochemical reduction of complexes [Fe

International audience; Electrochemical oxidation of the complex [Fe2(CO)4(κ2‐dmpe)(µ‐adtBn)] (adtBn = (SCH2)2NCH2C6H5, dmpe = Me2PCH2CH2PMe2) (1) has been studied by cyclic voltammetry (CV) in acetonitrile and in dichloromethane in presence of various substrates L (L =MeCN, trimethylphosphite, isocyanide). The oxidized species, [1‐MeCN](PF6)2, [1‐(P(OMe)3)2](PF6)2 and [1‐(RNC)4](PF6)2 (R = tButyl, Xylyl), have been prepared and characterized by IR, NMR spectroscopies and, excepted [1‐MeCN](PF6)2, by X‐ray diffraction analysis. The crystallographic structures of the new FeIIFeII complexes reveal in any case that the association of one additional ligand (P(OMe)3 or RNC) occurs and, according to the nature of the substrates, further substitutions of one or three carbonyl groups, by P(OMe)3 or RNC, respectively, arise. Density functional theory (DFT) calculations have been performed for elucidating and discriminating, in each case, the mechanisms leading to the corresponding oxidized species. Moreover, the different degree of ligand substitution in the diiron core has been theoretically rationalized.

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 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

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.

The electrochemical behavior of complexes [FeMo(CO) 5 (κ 2 -dppe)(μ-pdt)] (1) and [FeMo(CO) 4 (MeCN)( κ 2 -dppe)(μ-pdt)] (2), in the absence and in the presence of acid, has been investigated. The reduction of 1 follows at slow scan rates, in CH 2 Cl 2 -[NBu 4 ][PF 6 ] and acid-free media, an EC rev E mechanism that is supported by cyclic voltammetry (CV) experiments and digital CV simulations. In MeCN-[NBu 4 ][PF 6 ], the electrochemical reduction of 1 is the same as in dichloromethane and follows an ECE mechanism at slow scan rates, but with a positive shift of the redox potentials. In contrast, the oxidation of 1 is strongly solvent-dependent. In dichloromethane, the oxidation of 1 is reversible and involves a single electron, while in acetonitrile, it is irreversible at moderate and slow scan rates (v ≤ ca. 1 V s -1 ), and some chemical reversibility is apparent at higher scan rates (v = 10 V s -1 ). Density functional theory calculations revealed that the chemical step in the EC rev E mechanism corresponds to the dissociation of one PPh 2 end of the diphosphine ligand and the transfer of the semibridging CO to the Fe atom, similarly to the mechanism observed in the FeFe analogue complex. However, in the case of 1, the subsequent coordination of the phosphine ligand to the other metal is an unfavorable process.