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

Understanding Ligand Effects on Bielectronic Transitions: Chemo- and Electroreduction of Rhodium Bis(Diphosphine) Complexes to Low Oxidation States

Rhodium complexes in the -I and 0 oxidation states are of great potential interest in catalytic applications. In contrast to their rhodium +I congeners, however, the structural and electronic parameters governing their access and stability are far less understood. Herein, we investigate the two-electron reduction of a parameterized series of bis(diphosphine) Rh complexes [Rh(dxpy)2]NTf2 (x=P-substituent, y=alkanediyl bridging P atoms). Through (electro)reductions from the RhI parents, Rh-I d10-complexes were obtained and characterized spectroscopically, including 103Rh NMR data. The reductive steps convolute with structural rearrangements from square planar to tetrahedral coordination. We found that the extent of these reorganisations defines whether the first E0(RhI/0) and second E0(Rh0/-I) reduction potentials are normally ordered, leading to monoelectronic stepwise transitions, or inverted, giving bielectronic events. Reductionist approaches based on Hammett parameters or the P-Rh-P bite angles provide only partial correlations with the redox potentials. However, we identified the C-O stretch of analogue diphosphine complexes as an expedient computational parameter that enables these correlations through both electronic and geometric features, even in a predictive manner. Gaining control over two-electron reduction behaviors through rationalized ligand effects has potential impact beyond Rh complexes, for molecular and enzymatic metal sites commonly exhibiting bielectronic transitions.

Electrochemical reduction and protonation of a biomimetic diiron azadithiolate hexacarbonyl complex: Mechanistic insights

The electrochemical reduction and protonation of [Fe2(adtH)(CO)6] (1, adtH = SCH2N(H)CH2S) and [Fe2(pdt)(CO)6] (2, pdt = SCH2CH2CH2S) in the presence of moderately strong acid in acetonitrile was investigated by cyclic voltammetry (CV), focusing on the catalysis of hydrogen evolution reaction (HER) by a {2e-,2H+} pathway. The turnover frequencies at zero overpotential (TOF0) of the N-protonated product 1(H)+ and 2 for the HER were estimated from simulations of the catalytic CV responses at low acid concentration using a simple ECEC mechanism (two electrochemical and chemical steps). This approach confirmed that 1(H)+ is clearly a better catalyst than 2, pointing to a possible role of the protonable and biologically relevant adtH ligand in the enhancement of the catalytic performances. Density functional theory (DFT) calculations further suggested that, owing to a strong structural rearrangement in the course of the catalytic cycle, the HER catalysis by 1(H)+ only involves the iron center adjacent to the amine group in adtH and not the two iron centers as in 2. Since terminal hydride species (FeFe-H) are known to more easily undergo protonolyse to H2 than their bridging hydride isomers (Fe-H-Fe), this may explain here the enhanced activity of 1(H)+ over 2 for the HER.

Machine Learning for Efficient Prediction of Protein Redox Potential: The Flavoproteins Case

Determining the redox potentials of protein cofactors and how they are influenced by their molecular neighborhoods is essential for basic research and many biotechnological applications, from biosensors and biocatalysis to bioremediation and bioelectronics. The laborious determination of redox potential with current experimental technologies pushes forward the need for computational approaches that can reliably predict it. Although current computational approaches based on quantum and molecular mechanics are accurate, their large computational costs hinder their usage. In this work, we explored the possibility of using more efficient QSPR models based on machine learning (ML) for the prediction of protein redox potential, as an alternative to classical approaches. As a proof of concept, we focused on flavoproteins, one of the most important families of enzymes directly involved in redox processes. To train and test different ML models, we retrieved a dataset of flavoproteins with a known midpoint redox potential (Em) and 3D structure. The features of interest, accounting for both short- and long-range effects of the protein matrix on the flavin cofactor, have been automatically extracted from each protein PDB file. Our best ML model (XGB) has a performance error below 1 kcal/mol (∼36 mV), comparing favorably to more sophisticated computational approaches. We also provided indications on the features that mostly affect the Em value, and when possible, we rationalized them on the basis of previous studies.

A Modular Strategy for Expanding Electron-Sink Capacity in Noncanonical Cluster Assemblies

A modular synthetic strategy is described whereby organometallic complexes exhibiting considerable electron-sink capacity may be assembled by using only a few simple molecular components. The Fe₂(PPh₂)₂(CO)₅ fragment was selected as a common electroactive component and was assembled around aromatic cores bearing one, two, or three isocyanide functional groups, with the resultant complexes possessing electron-sink capacities of two, four, and six electrons, respectively. The latter complex is noteworthy in that its electron-sink capacity was found to rival that of large multinuclear clusters (e.g., [Ni₃₂C₆(CO)₃₆]^6– and [Ni₃₈Pt₆(CO)₄₈]^6–), which are often considered as benchmarks of electron-sink behavior. Moreover, the modular assembly bearing three Fe₂(PPh₂)₂(CO)₅ fragments was observed to undergo reduction to a hexaanionic state over a potential window of about −1.4 to −2.1 V (vs Fc/Fc⁺), the relatively compressed range being attributed to potential inversions operative during the addition of the second, fourth, and sixth electrons. Such complexes may be designated noncanonical clusters because they exhibit redox properties similar to those of large multinuclear clusters yet lack the extensive network of metal–metal bonds and the condensed metallic cores that typify the latter.

By use of a modular synthetic strategy and relatively few molecular components, organometallic complexes exhibiting considerable electron-sink capacity have been characterized. Complexes bearing one, two, or three Fe₂(PPh₂)₂(CO)₅ fragments bound to aromatic isocyanide cores were found to possess electron-sink capacities of two, four, and six electrons, respectively, the latter rivaling the electron-sink capacity of large polynuclear cluster benchmarks.

Computational and Experimental Investigations of the Fe2(μ-S2)/Fe2(μ-S)2 Equilibrium

Density functional theory (DFT) calculations on Fe2S2(CO)6-2n(PMe3)2n for n = 0, 1, and 2 reveal that the most electron-rich derivatives (n = 2) exist as diferrous disulfides lacking an S-S bond. The thermal interconversion of the FeII2(S)2 and FeI2(S2) valence isomers is symmetry-forbidden. Related electron-rich diiron complexes [Fe2S2(CN)2(CO)4]2- of an uncertain structure are implicated in the biosynthesis of [FeFe]-hydrogenases. Several efforts to synthesize electron-rich derivatives of Fe2(μ-S2)(CO)6 (1) are described. First, salts of iron persulfido cyanides [Fe2(μ-S2)(CO)5(CN)]- and [Fe2(μ-S2)(CN)(CO)4(PPh3)]- were prepared by the reactions of NaN(tms)2 with 1 and Fe2(μ-S2)(CO)5(PPh3), respectively. Alternative approaches to electron-rich diiron disulfides targeted Fe2(μ-S2)(CO)4(diphosphine). Whereas the preparation of Fe2(μ-S2)(CO)4(dppbz) was straightforward, that of Fe2(μ-S2)(CO)4(dppv) required an indirect route involving the oxidation of Fe2(μ-SH)2(CO)4(dppv) (dppbz = C6H4-1,2-(PPh2)2, dppv = cis-C2H2(PPh2)2). DFT calculations indicate that the oxidation of Fe2(μ-SH)2(CO)4(dppv) produces singlet diferrous disulfide Fe2(μ-S)2(CO)4(dppv), which is sufficiently long-lived as to be trapped by ethylene. The reaction of 1 and dppv mainly afforded Fe2(μ-SCH=CHPPh2)(μ-SPPh2)(CO)5, implicating a S-centered reaction.

The Photochemistry of Fe2(S2C3H6)(CO)6(µ-CO) and Its Oxidized Form, Two Simple [FeFe]-Hydrogenase CO-Inhibited Models. A DFT and TDDFT Investigation

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.