Reactivity of the Excited States of the H-Cluster of FeFe Hydrogenases

How general is the effect of the bulkiness of organic ligands on the basicity of metal-organic catalysts? H2-evolving Mo oxides/sulphides as case studies

Tailoring the activity of an organometallic catalyst usually requires a targeted ligand design. Tuning the ligand bulkiness and tuning the electronic properties are popular approaches, which are somehow interdependent because substituents of different sizes within ligands can determine inter alia the occurrence of different degrees of inductive effects. Ligand basicity, in particular, turned out to be a key property for the modulation of protonation reactions occurring in vacuo at the metals in complexes bearing organophosphorus ligands; however, when the same reactions take place in a polar organic solvent, their energetics becomes dependent on the trade-off between ligand basicity and bulkiness, with the polarity of the solvent playing a key role in this regard [Bancroft et al., Inorg. Chem., 1986, 25, 3675; Rovaletti et al., J. Phys. Org. Chem., 2018, 31, e3748]. In the present contribution, we carried out molecular dynamics and density functional theory calculations on water-soluble Mo-based catalysts for proton reduction, in order to study the energetics of protonation reactions in complexes where the incipient proton binds a catalytically active ligand (i.e., an oxide or a disulphide). We considered complexes either soaked in water or in a vacuum, and featuring N-based ancillary ligands of different bulkiness (i.e. cages constituted either by pyridine or isoquinoline moieties). Our results show that the energetics of protonation events can be affected by ancillary ligand bulkiness even when the metal center does not play the role of the H⁺ acceptor. In vacuo, protonation at the O or S atom in the α position relative to the metal in complexes featuring the bulky isoquinoline-based ligand is more favored by around 10 kcal mol^-1 when compared to the case of the pyridine-based counterparts, a difference that is almost zero when the same reactions occur in water. Such an outcome is rationalized in light of the different electrostatic properties of complexes bearing ancillary ligands of different sizes. The overall picture from theory indicates that such effects of ligand bulkiness can be relevant for the design of green chemistry catalysts that undergo protonation steps in water solutions.

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.

Accurate Molecular Geometries in Complex Excited-State Potential Energy Surfaces from Time-Dependent Density Functional Theory

The interplay of electronic excitations and structural changes in molecules impacts nonradiative decay and charge transfer in the excited state, thus influencing excited-state lifetimes and photocatalytic reaction rates in optoelectronic and energy devices. To capture such effects requires computational methods providing an accurate description of excited-state potential energy surfaces and geometries. We suggest time-dependent density functional theory using optimally tuned range-separated hybrid (OT-RSH) functionals as an accurate approach to obtain excited-state molecular geometries. We show that OT-RSH provides accurate molecular geometries in excited-state potential energy surfaces that are complex and involve an interplay of local and charge-transfer excitations, for which conventional semilocal and hybrid functionals fail. At the same time, the nonempirical OT-RSH approach maintains the high accuracy of parametrized functionals (e.g., B3LYP) for predicting excited-state geometries of small organic molecules showing valence excited states.

Ligand effects on structural, protophilic and reductive features of stannylated dinuclear iron dithiolato complexes

The synthesis and characterization of Fe₂(CO)₅(L){μ-(SCH₂)₂SnMe₂} (L = PPh₃ (2) and P(OMe)₃ (3)) derived from the parent hexacarbonyl complex Fe₂(CO)₆{μ-(SCH₂)₂}SnMe₂ (1) are reported.

Energetics for Proton Reduction in FeFe Hydrogenase

The energetics for proton reduction in FeFe-hydrogenase has been reinvestigated by theoretical modeling, in light of recent experiments. Two different mechanisms have been considered. In the first one, the bridging hydride position was blocked by the enzyme, which is the mechanism that has been supported by a recent spectroscopic study by Cramer et al. A major difficulty in the present study to agree with experimental energetics was to find the right position for the added proton in the first reduction step. It was eventually found that the best position was as a terminal hydride on the distal iron, which has not been suggested in any of the recent, experimentally based mechanisms. The lowest transition state was surprisingly found to be a bond formation between a proton on a cysteine and the terminal hydride. This type of TS is similar to the one for heterolytic H₂ cleavage in NiFe hydrogenase. The second mechanism investigated here is not supported by the present calculations or the recent experiments by Cramer et al., but was still studied as an interesting comparison. In that mechanism, the formation of the bridging hydride was allowed. The H–H formation barrier is only 3.6 kcal/mol higher than for the first mechanism, but there are severe problems concerning the motion of the protons.

Molecular Photochemistry: Recent Developments in Theory

Photochemistry is a fascinating branch of chemistry that is concerned with molecules and light. However, the importance of simulating light-induced processes is reflected also in fields as diverse as biology, material science, and medicine. This Minireview highlights recent progress achieved in theoretical chemistry to calculate electronically excited states of molecules and simulate their photoinduced dynamics, with the aim of reaching experimental accuracy. We focus on emergent methods and give selected examples that illustrate the progress in recent years towards predicting complex electronic structures with strong correlation, calculations on large molecules, describing multichromophoric systems, and simulating non-adiabatic molecular dynamics over long time scales, for molecules in the gas phase or in complex biological environments.

Spin-symmetrised structures and vibrational frequencies of iron-sulfur clusters

Calculations of relaxed geometries of multi-centre transition metal compounds are routinely carried out using Broken Symmetry Density Functional Theory. The resulting low-spin open shell electronic state is described by one single Slater determinant and is affected by spin contamination. To alleviate this symmetry breaking, the Extended Broken Symmetry (EBS) approach can be applied to complexes with an arbitrary number of local high-spin metal ions. The actual symmetry is therefore reconstructed through minimization of an effective Hamiltonian leading to a relaxed geometry consistent with the magnetic couplings. In the present work we extend the approach already introduced by [Chu et al., J. Chem. Theory Comput., 2017, 13, 4675] to the calculation of vibrational frequencies. As prototypes we have considered the iron-sulfur clusters Fe₂S₂Cl₄^2- and Fe₄S₄Cl₄. We have compared the results obtained for different spin states (high spin, broken symmetry and extended broken symmetry) and by using different DFT functionals (B3LYP, OPBE, BP, M06 and B2PLYP) and a post-HF method (SCS-MP2). The data have shown that for specific vibrational modes the EBS technique produces shifts up to 40 cm^-1 with respect to the routinely used Broken Symmetry approach, indicating that the use of a consistent spin-symmetrised state is a crucial ingredient for an accurate description of vibrational properties, as certified by the comparison with the experimental data for the Fe₂S₂Cl₄^2- cluster.

Anomalous Intrinsic Fluorescence of HCl and NaOH Aqueous Solutions

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.

Engineering an [FeFe]-Hydrogenase: Do Accessory Clusters Influence O2 Resistance and Catalytic Bias?

[FeFe]-hydrogenases, HydAs, are unique biocatalysts for proton reduction to H₂. However, they suffer from a number of drawbacks for biotechnological applications: size, number and diversity of metal cofactors, oxygen sensitivity. Here we show that HydA from Megasphaera elsdenii (MeHydA) displays significant resistance to O₂. Furthermore, we produced a shorter version of this enzyme (MeH-HydA), lacking the N-terminal domain harboring the accessory FeS clusters. As shown by detailed spectroscopic and biochemical characterization, MeH-HydA displays the following interesting properties. First, a functional active site can be assembled in MeH-HydA in vitro, providing the enzyme with excellent hydrogenase activity. Second, the resistance of MeHydA to O₂ is conserved in MeH-HydA. Third, MeH-HydA is more biased toward proton reduction than MeHydA, as the result of the truncation changing the rate limiting steps in catalysis. This work shows that it is possible to engineer HydA to generate an active hydrogenase that combines the resistance of the most resistant HydAs and the simplicity of algal HydAs, containing only the H-cluster.

Electrochemical Investigations of Hydrogenases and Other Enzymes That Produce and Use Solar Fuels

Many enzymes that produce or transform small molecules such as O₂, H₂, and CO₂ embed inorganic cofactors based on transition metals. Their active site, where the chemical reaction occurs, is buried in and protected by the protein matrix, and connected to the solvent in several ways: chains of redox cofactors mediate long-range electron transfer; static or dynamic tunnels guide the substrate, product and inhibitors; amino acids and water molecules transfer protons. The catalytic mechanism of these enzymes is therefore delocalized over the protein and involves many different steps, some of which determine the response of the enzyme under conditions of stress (extreme redox conditions, presence of inhibitors, light), the catalytic rates in the two directions of the reaction and their ratio (the "catalytic bias"). Understanding all the steps in the catalytic cycle, including those that occur on sites of the protein that are remote from the active site, requires a combination of biochemical, structural, spectroscopic, theoretical, and kinetic methods. Here we argue that kinetics should be used to the fullest extent, by extracting quantitative information from the comparison of data and kinetic models and by exploring the combination of experimental kinetics and theoretical chemistry. In studies of these catalytic mechanisms, direct electrochemistry, the technique which we use and contribute to develop, has become unescapable. It simply consists in monitoring the changes in activity of an enzyme that is wired to an electrode by recording an electric current. We have described kinetic models that can be used to make sense of these data and to learn about various aspects of the mechanism that are difficult to probe using more conventional methods: long-range electron transfer, diffusion along gas channels, redox-driven (in)activations, active site chemistry and photoreactivity under conditions of turnover. In this Account, we highlight a few results that illustrate our approach. We describe how electrochemistry can be used to monitor substrate and inhibitor diffusion along the gas channels of hydrogenases and we discuss how the kinetics of intramolecular diffusion relates to global properties such as resistance to oxygen and catalytic bias. The kinetics and/or thermodynamics of intramolecular electron transfer may also affect the catalytic bias, the catalytic potentials on either side of the equilibrium potential, and the overpotentials for catalysis (defined as the difference between the catalytic potentials and the open circuit potential). This is understood by modeling the shape of the steady-state catalytic response of the enzyme. Other determinants of the catalytic rate, such as domain motions, have been probed by examining the transient catalytic response recorded at fast scan rates. Last, we show that combining electrochemical investigations and MD, DFT, and TD-DFT calculations is an original way of probing the reactivity of the H-cluster of hydrogenase, in particular its reactions with CO, O₂, and light. This approach contrasts with the usual strategy which aims at stabilizing species that are presumed to be catalytic intermediates, and determining their structure using spectroscopic or structural methods.

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

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

Effects of Static Correlation between Spin Centers in Multicenter Transition Metal Complexes

Multicenter transition metal complexes are the key moieties of many processes in chemistry, biochemistry, and materials science such as in the active sites of enzymes, molecular catalysts, and biological electron carriers. Their electronic structure, often characterized by high-spin-polarized metal sites, is a challenge for theoretical chemists because of their high degree of dynamical and static correlation. Static correlation is necessary both for the appropriate description of the metal-ligand bonding and for a correct description of the multideterminant character arising from the magnetic interactions between spin centers. Density functional theory (DFT) is usually applied using a single-determinant broken-symmetry state that is lacking the correct spin symmetry when the ground state has total low-spin character. To alleviate this drawback, we use the extended broken-symmetry (EBS) approach to derive approximate ground-state energies and, for the first time, forces for the correctly symmetric ground state of an arbitrary number of spin centers within the framework of the Heisenberg-Dirac-van Vleck Hamiltonian. Remarkably, the proposed procedure supplies relaxed geometries that are fully consistent with the calculated J-coupling constants. We apply the method to investigate the relaxed geometrical structure of the low-spin ground state of iron-sulfur clusters with two, three, and four iron centers. We observed significant differences in both geometrical parameters and coupling constant J between the symmetrized ground state, the high-spin, and the broken-symmetry optimized structures. These changes are often comparable with the differences observed by using different functionals, and the use of EBS always improves the description of the studied systems. It will be therefore important to include it in any DFT attempt to quantitatively describe multicenter transition metal complexes in the future.