Mechanisms of Mn(V)-Oxo to Mn(IV)-Oxyl Conversion: From Closed-Cubane Photosystem II to Mn(V) Catalysts and the Role of the Entering Ligands

The low activation barrier for O-O coupling in the closed-cubane Oxygen-Evolving Centre (OEC) of Photosystem II (PSII) requires water coordination with the Mn4 'dangler' ion in the Mn(V)-oxo fragment. This coordination transforms the Mn(V)-oxo complex into a more reactive Mn4(IV)-oxyl species, enhancing O-O coupling. This study explains the mechanism behind the coordination and indicates that in the most stable form of the OEC, the Mn4 fragment adopts a trigonal bipyramidal geometry but needs to transition to a square pyramidal form to be activated for O-O coupling. This transition stabilizes the Mn4 dxy orbital, enabling electron transfer from the oxo ligand to the dxy orbital, converting the oxo ligand into an oxyl species. The role of the water is to coordinate with the square pyramidal structure, reducing the energy gap between the oxo and oxyl forms, thereby lowering the activation energy for O-O coupling. This mechanism applies not only to the OEC system but also to other Mn(V)-based catalysts. For other catalysts, ligands such as OH- stabilize the Mn(IV)-oxyl species better than water, improving catalyst activation for reactions like C-H bond activation. This study is the first to explain the Mn(V)-oxo to Mn(IV)-oxyl conversion, providing a new foundation for Mn-based catalyst design.

Independent Mutation of Two Bridging Carboxylate Ligands Stabilizes Alternate Conformers of the Photosynthetic O2-Evolving Mn4CaO5 Cluster in Photosystem II

The O2-evolving Mn4CaO5 cluster in photosystem II is ligated by six carboxylate residues. One of these is D170 of the D1 subunit. This carboxylate bridges between one Mn ion (Mn4) and the Ca ion. A second carboxylate ligand is D342 of the D1 subunit. This carboxylate bridges between two Mn ions (Mn1 and Mn2). D170 and D342 are located on opposite sides of the Mn4CaO5 cluster. Recently, it was shown that the D170E mutation perturbs both the intricate networks of H-bonds that surround the Mn4CaO5 cluster and the equilibrium between different conformers of the cluster in two of its lower oxidation states, S1 and S2, while still supporting O2 evolution at approximately 50% the rate of the wild type. In this study, we show that the D342E mutation produces much the same alterations to the cluster's FTIR and EPR spectra as D170E, while still supporting O2 evolution at approximately 20% the rate of the wild type. Furthermore, the double mutation, D170E + D342E, behaves similarly to the two single mutations. We conclude that D342E alters the equilibrium between different conformers of the cluster in its S1 and S2 states in the same manner as D170E and perturbs the H-bond networks in a similar fashion. This is the second identification of a Mn4CaO5 metal ligand whose mutation influences the equilibrium between the different conformers of the S1 and S2 states without eliminating O2 evolution. This finding has implications for our understanding of the mechanism of O2 formation in terms of catalytically active/inactive conformations of the Mn4CaO5 cluster in its lower oxidation states.

Correlating Structure with Spectroscopy in Ascorbate Peroxidase Compound II

Structural and spectroscopic investigations of compound II in ascorbate peroxidase (APX) have yielded conflicting conclusions regarding the protonation state of the crucial Fe(IV) intermediate. Neutron diffraction and crystallographic data support an iron(IV)-hydroxo formulation, whereas Mössbauer, X-ray absorption (XAS), and nuclear resonance vibrational spectroscopy (NRVS) studies appear consistent with an iron(IV)-oxo species. Here we examine APX with spectroscopy-oriented QM/MM calculations and extensive exploration of the conformational space for both possible formulations of compound II. We establish that irrespective of variations in the orientation of a vicinal arginine residue and potential reorganization of proximal water molecules and hydrogen bonding, the Fe-O distances for the oxo and hydroxo forms consistently fall within distinct, narrow, and nonoverlapping ranges. The accuracy of geometric parameters is validated by coupled-cluster calculations with the domain-based local pair natural orbital approach, DLPNO-CCSD(T). QM/MM calculations of spectroscopic properties are conducted for all structural variants, encompassing Mössbauer, optical, X-ray absorption, and X-ray emission spectroscopies and NRVS. All spectroscopic observations can be assigned uniquely to an Fe(IV)═O form. A terminal hydroxy group cannot be reconciled with the spectroscopic data. Under no conditions can the Fe(IV)═O distance be sufficiently elongated to approach the crystallographically reported Fe-O distance. The latter is consistent only with a hydroxo species, either Fe(IV) or Fe(III). Our findings strongly support the Fe(IV)═O formulation of APX-II and highlight unresolved discrepancies in the nature of samples used across different experimental studies.

An artificial metalloenzyme that can oxidize water photocatalytically: design, synthesis, and characterization

In nature, light-driven water oxidation (WO) catalysis is performed by photosystem II via the delicate interplay of different cofactors positioned in its protein scaffold. Artificial systems for homogeneous photocatalytic WO are based on small molecules that often have limited solubility in aqueous solutions. In this work, we alleviated this issue and present a cobalt-based WO-catalyst containing artificial metalloenzyme (ArM) that is active in light-driven, homogeneous WO catalysis in neutral-pH aqueous solutions. A haem-containing electron transfer protein, cytochrome B5 (CB5), served to host a first-row transition-metal-based WO catalyst, CoSalen (CoIISalen, where H2Salen = N,N'-bis(salicylidene)ethylenediamine), thus producing an ArM capable of driving photocatalytic WO. The CoSalen ArM formed a water-soluble pre-catalyst in the presence of [Ru(bpy)3](ClO4)2 as photosensitizer and Na2S2O8 as the sacrificial electron acceptor, with photocatalytic activity similar to that of free CoSalen. During photocatalysis, the CoSalen-protein interactions were destabilized, and the protein partially unfolded. Rather than forming tens of nanometer sized CoOx nanoparticles as free CoSalen does under photocatalytic WO conditions, the CB5 : CoSalen ArM showed limited protein cross-linking and remained soluble. We conclude that a weak, dynamic interaction between a soluble cobalt species and apoCB5 was formed, which generated a catalytically active adduct during photocatalysis. A detailed analysis was performed on protein stability and decomposition processes during the harsh oxidizing reaction conditions of WO, which will serve for the future design of WO ArMs with improved activity and stability.

Evolution in the Design of Water Oxidation Catalysts with Transition-Metals: A Perspective on Biological, Molecular, Supramolecular, and Hybrid Approaches

Increased demand for a carbon-neutral sustainable energy scheme augmented by climatic threats motivates the design and exploration of novel approaches that reserve intermittent solar energy in the form of chemical bonds in molecules and materials. In this context, inspired by biological processes, artificial photosynthesis has garnered significant attention as a promising solution to convert solar power into chemical fuels from abundantly found H2O. Among the two redox half-reactions in artificial photosynthesis, the four-electron oxidation of water according to 2H2O → O2 + 4H+ + 4e- comprises the major bottleneck and is a severe impediment toward sustainable energy production. As such, devising new catalytic platforms, with traditional concepts of molecular, materials and biological catalysis and capable of integrating the functional architectures of the natural oxygen-evolving complex in photosystem II would certainly be a value-addition toward this objective. In this review, we discuss the progress in construction of ideal water oxidation catalysts (WOCs), starting with the ingenuity of the biological design with earth-abundant transition metal ions, which then diverges into molecular, supramolecular and hybrid approaches, blurring any existing chemical or conceptual boundaries. We focus on the geometric, electronic, and mechanistic understanding of state-of-the-art homogeneous transition-metal containing molecular WOCs and summarize the limiting factors such as choice of ligands and predominance of environmentally unrewarding and expensive noble-metals, necessity of high-valency on metal, thermodynamic instability of intermediates, and reversibility of reactions that create challenges in construction of robust and efficient water oxidation catalyst. We highlight how judicious heterogenization of atom-efficient molecular WOCs in supramolecular and hybrid approaches put forth promising avenues to alleviate the existing problems in molecular catalysis, albeit retaining their fascinating intrinsic reactivities. Taken together, our overview is expected to provide guiding principles on opportunities, challenges, and crucial factors for designing novel water oxidation catalysts based on a synergy between conventional and contemporary methodologies that will incite the expansion of the domain of artificial photosynthesis.

Comprehensive Evaluation of Models for Ammonia Binding to the Oxygen Evolving Complex of Photosystem II

The identity and insertion pathway of the substrate oxygen atoms that are coupled to dioxygen by the oxygen-evolving complex (OEC) remains a central question toward understanding Nature's water oxidation mechanism. In several studies, ammonia has been used as a small "water analogue" to elucidate the pathway of substrate access to the OEC and to aid in determining which of the oxygen ligands of the tetramanganese cluster are substrates for O-O bond formation. On the basis of structural and spectroscopic investigations, five first-sphere binding modes of ammonia have been suggested, involving either substitution of an existing H2O/OH-/O2- group or addition as an extra ligand to a metal ion of the Mn4CaO5 cluster. Some of these modes, specifically the ones involving substitution, have already been subject to spectroscopy-oriented quantum chemical investigations, whereas more recent suggestions that postulate the addition of ammonia have not been examined so far with quantum chemistry for their agreement with spectroscopic data. Herein, we use a common structural framework and theoretical methodology to evaluate structural models of the OEC that represent all proposed modes of first-sphere ammonia interaction with the OEC in its S2 state. Criteria include energetic, magnetic, kinetic, and spectroscopic properties compared against available experimental EPR, ENDOR, ESEEM, and EDNMR data. Our results show that models featuring ammonia replacing one of the two terminal water ligands on Mn4 align best with experimental data, while they definitively exclude substitution of a bridging μ-oxo ligand as well as incorporation of ammonia as a sixth ligand on Mn1 or Mn4.

Nature of S-States in the Oxygen-Evolving Complex Resolved by High-Energy Resolution Fluorescence Detected X-ray Absorption Spectroscopy

Photosystem II, the water splitting enzyme of photosynthesis, utilizes the energy of sunlight to drive the four-electron oxidation of water to dioxygen at the oxygen-evolving complex (OEC). The OEC harbors a Mn4CaO5 cluster that cycles through five oxidation states Si (i = 0-4). The S3 state is the last metastable state before the O2 evolution. Its electronic structure and nature of the S2 → S3 transition are key topics of persisting controversy. Most spectroscopic studies suggest that the S3 state consists of four Mn(IV) ions, compared to the Mn(III)Mn(IV)3 of the S2 state. However, recent crystallographic data have received conflicting interpretations, suggesting either metal- or ligand-based oxidation, the latter leading to an oxyl radical or a peroxo moiety in the S3 state. Herein, we utilize high-energy resolution fluorescence detected (HERFD) X-ray absorption spectroscopy to obtain a highly resolved description of the Mn K pre-edge region for all S-states, paying special attention to use chemically unperturbed S3 state samples. In combination with quantum chemical calculations, we achieve assignment of specific spectroscopic features to geometric and electronic structures for all S-states. These data are used to confidently discriminate between the various suggestions concerning the electronic structure and the nature of oxidation events in all observable catalytic intermediates of the OEC. Our results do not support the presence of either peroxo or oxyl in the active configuration of the S3 state. This establishes Mn-centered storage of oxidative equivalents in all observable catalytic transitions and constrains the onset of the O-O bond formation until after the final light-driven oxidation event.

Mutation of a metal ligand stabilizes the high-spin form of the S2 state in the O2-producing Mn4CaO5 cluster of photosystem II

The residue D1-D170 bridges Mn4 with the Ca ion in the O₂-evolving Mn₄CaO₅ cluster of Photosystem II. Recently, the D1-D170E mutation was shown to substantially alter the S_n+1-minus-S_n FTIR difference spectra [Debus RJ (2021) Biochemistry 60:3841-3855]. The mutation was proposed to alter the equilibrium between different Jahn-Teller conformers of the S₁ state such that (i) a different S₁ state conformer is stabilized in D1-D170E than in wild-type and (ii) the S₁ to S₂ transition in D1-D170E produces a high-spin form of the S₂ state rather than the low-spin form that is produced in wild-type. In this study, we employed EPR spectroscopy to test if a high-spin form of the S₂ state is formed preferentially in D1-D170E PSII. Our data show that illumination of dark-adapted D1-D170E PSII core complexes does indeed produce a high-spin form of the S₂ state rather than the low-spin multiline form that is produced in wild-type. This observation provides further experimental support for a change in the equilibrium between S state conformers in both the S₁ and S₂ states in a site-directed mutant that retains substantial O₂ evolving activity.

Does Serial Femtosecond Crystallography Depict State-Specific Catalytic Intermediates of the Oxygen-Evolving Complex?

Recent advances in serial femtosecond crystallography (SFX) of photosystem II (PSII), enabled by X-ray free electron lasers (XFEL), provided the first geometric models of distinct intermediates in the catalytic S-state cycle of the oxygen-evolving complex (OEC). These models are obtained by flash-advancing the OEC from the dark-stable state (S₁) to more oxidized intermediates (S₂ and S₃), eventually cycling back to the most reduced S₀. However, the interpretation of these models is controversial because geometric parameters within the Mn₄CaO₅ cluster of the OEC do not exactly match those expected from coordination chemistry for the spectroscopically verified manganese oxidation states of the distinct S-state intermediates. Here we focus on the first catalytic transition, S₁ → S₂, which represents a one-electron oxidation of the OEC. Combining geometric and electronic structure criteria, including a novel effective oxidation state approach, we analyze existing 1-flash (1F) SFX-XFEL crystallographic models that should depict the S₂ state of the OEC. We show that the 1F/S₂ equivalence is not obvious, because the Mn oxidation states and total unpaired electron counts encoded in these models are not fully consistent with those of a pure S₂ state and with the nature of the S₁ → S₂ transition. Furthermore, the oxidation state definition in two-flashed (2F) structural models is practically impossible to elucidate. Our results advise caution in the extraction of electronic structure information solely from the literal interpretation of crystallographic models and call for re-evaluation of structural and mechanistic interpretations that presume exact correspondence of such models to specific catalytic intermediates of the OEC.

Predicting the oxidation states of Mn ions in the oxygen-evolving complex of photosystem II using supervised and unsupervised machine learning

Serial Femtosecond Crystallography at the X-ray Free Electron Laser (XFEL) sources enabled the imaging of the catalytic intermediates of the oxygen evolution reaction of Photosystem II (PSII). However, due to the incoherent transition of the S-states, the resolved structures are a convolution from different catalytic states. Here, we train Decision Tree Classifier and K-means clustering models on Mn compounds obtained from the Cambridge Crystallographic Database to predict the S-state of the X-ray, XFEL, and CryoEM structures by predicting the Mn's oxidation states in the oxygen-evolving complex. The model agrees mostly with the XFEL structures in the dark S₁ state. However, significant discrepancies are observed for the excited XFEL states (S₂, S_3, and S₀) and the dark states of the X-ray and CryoEM structures. Furthermore, there is a mismatch between the predicted S-states within the two monomers of the same dimer, mainly in the excited states. We validated our model against other metalloenzymes, the valence bond model and the Mn spin densities calculated using density functional theory for two of the mismatched predictions of PSII. The model suggests designing a more optimized sample delivery and illumiation systems are crucial to precisely resolve the geometry of the advanced S-states to overcome the noncoherent S-state transition. In addition, significant radiation damage is observed in X-ray and CryoEM structures, particularly at the dangler Mn center (Mn4). Our model represents a valuable tool for investigating the electronic structure of the catalytic metal cluster of PSII to understand the water splitting mechanism.

How Metal Nuclearity Impacts Electrocatalytic H2 Production in Thiocarbohydrazone-Based Complexes

Thiocarbohydrazone-based catalysts feature ligands that are potentially electrochemically active. From the synthesis point of view, these ligands can be easily tailored, opening multiple strategies for optimization, such as using different substituent groups or metal substitution. In this work, we show the possibility of a new strategy, involving the nuclearity of the system, meaning the number of metal centers. We report the synthesis and characterization of a trinuclear nickel-thiocarbohydrazone complex displaying an improved turnover rate compared with its mononuclear counterpart. We use DFT calculations to show that the mechanism involved is metal-centered, unlike the metal-assisted ligand-centered mechanism found in the mononuclear complex. Finally, we show that two possible mechanisms can be assigned to this catalyst, both involving an initial double reduction of the system.

Alternative Mechanism for O2 Formation in Natural Photosynthesis via Nucleophilic Oxo-Oxo Coupling

O2 formation in photosystem II (PSII) is a vital event on Earth, but the exact mechanism remains unclear. The presently prevailing theoretical model is "radical coupling" (RC) involving a Mn(IV)-oxyl unit in an "open-cubane" Mn4CaO6 cluster, which is supported experimentally by the S3 state of cyanobacterial PSII featuring an additional Mn-bound oxygenic ligand. However, it was recently proposed that the major structural form of the S3 state of higher plants lacks this extra ligand, and that the resulting S4 state would feature instead a penta-coordinate dangler Mn(V)=oxo, covalently linked to a "closed-cubane" Mn3CaO4 cluster. For this proposal, we explore here a large number of possible pathways of O-O bond formation and demonstrate that the "nucleophilic oxo-oxo coupling" (NOOC) between Mn(V)=oxo and μ3-oxo is the only eligible mechanism in such a system. The reaction is facilitated by a specific conformation of the cluster and concomitant water binding, which is delayed compared to the RC mechanism. An energetically feasible process is described starting from the valid S4 state through the sequential formation of peroxide and superoxide, followed by O2 release and a second water insertion. The newly found mechanism is consistent with available experimental thermodynamic and kinetic data and thus a viable alternative pathway for O2 formation in natural photosynthesis, in particular for higher plants.

Synthesis and structural, magnetic and spectroscopic characterization of iron(III) complexes with in situ formed ligands from methyl-2-pyridyl ketone transformations

Reactions of methyl-2-pyridyl ketone, pyCOMe, with FeCl₃·6H₂O in various solvents gave complexes [Fe₄Cl₆(OMe)₂(L1)₂]·0.7MeCN·0.4MeOH (1·0.7MeCN·0.4MeOH) and [Fe₃Cl₄(bicine)(L2)]·Me₂CO·0.2H₂O (2·Me₂CO·0.2H₂O). The ligands (L1)^2- = pyCO(Me)CHCOpy (in 1) and (L2)^2- = pyCO(Me)CH₂CO(OMe)py (in 2) are formed in situ, through an aldol reaction-type mechanism between the carbanion pyC(O)CH₂^- (formed by the nucleophilic attack of the MeO^- in pyCOMe) and pyCOMe which results in the formation of a new C-C bond. The intermediate compound undergoes attack in the -CH₂- or -CO- group by a MeO^- group, and the new ligands (L1)^2- and (L2)^2-, respectively, are formed. The molecular structure of 1 consists of three corner-sharing [Fe₂O₂] rhombic units in cis-arrangement. The two terminal Fe^III ions display distorted square pyramidal geometry and the two central Fe^III ions are distorted octahedral. The molecular structure of 2 consists of two corner-sharing [Fe₂O₂] rhombic units, with the two terminal Fe^III ions in distorted square pyramidal geometry and the central Fe^III in distorted octahedral. The differentiation in the coordination environment of the Fe^III ions in 1-2 is reflected in the values of the Mössbauer hyperfine parameters. In agreement with theoretical calculations, the square pyramidal sites exhibit a smaller isomer shift value in comparison to the octahedral sites. Magnetic studies indicate antiferromagnetic interactions leading to an S = 0 ground state in 1 and to an S = 5/2 ground state in 2, consistent with Electron Paramagnetic Resonance spectroscopy. Mössbauer spectra of 2 indicate the onset of relaxation effects below 80 K. At 1.5 K the spectrum of 2 consists of magnetic sextets. The determined hyperfine magnetic fields are consistent with the exchange coupling scheme imposed by the crystal structure of 2. Theoretical calculations shed light on the differences in the electronic structure between the square pyramidal and the octahedral sites.

Water oxidation in oxygenic photosynthesis studied by magnetic resonance techniques

The understanding of light-induced biological water oxidation in oxygenic photosynthesis is of great importance both for biology and (bio)technological applications. The chemically difficult multistep reaction takes place at a unique protein-bound tetra-manganese/calcium cluster in photosystem II whose structure has been elucidated by X-ray crystallography (Umena et al. Nature 2011, 473, 55). The cluster moves through several intermediate states in the catalytic cycle. A detailed understanding of these intermediates requires information about the spatial and electronic structure of the Mn₄ Ca complex; the latter is only available from spectroscopic techniques. Here, the important role of Electron Paramagnetic Resonance (EPR) and related double resonance techniques (ENDOR, EDNMR), complemented by quantum chemical calculations, is described. This has led to the elucidation of the cluster's redox and protonation states, the valence and spin states of the manganese ions and the interactions between them, and contributed substantially to the understanding of the role of the protein surrounding, as well as the binding and processing of the substrate water molecules, the O-O bond formation and dioxygen release. Based on these data, models for the water oxidation cycle are developed.

Insight into the structural dynamics of light sensitive proteins from time-resolved crystallography and quantum chemical calculations

The structural dynamics underlying molecular mechanisms of light-sensitive proteins can be studied by a variety of experimental and computational biophysical techniques. Here we review recent progress in combining time-resolved crystallography at X-ray free electron lasers and quantum chemical calculations to study structural changes in photoenzymes, photosynthetic proteins, photoreceptors, and photoswitchable fluorescent proteins following photoexcitation.

The Role of Redox-Inactive Metals in Modulating the Redox Potential of the Mn4CaO4 Model Complex

Photosynthetic oxygen-evolving center (OEC), the "engine of life", is a unique Mn₄CaO₅ cluster catalyzing the water oxidation. The role of redox-inactive component Ca²⁺, which can only be functionally replaced by Sr²⁺ in a biological environment, has been under debate for a long time. Recently, its modulating effect on the redox potential of native OEC and artificial structural OEC model complex has received great attention, and linear relationship between the potential and the Lewis acidity of the redox-inactive metal has been proposed for the MMn₃O₄ model complex. In this work, the modulating effect has been studied in detail using the Mn₄CaO₄ model complex, which is the closest structural model to OEC to date and has a similar redox potential at the S₁-S₂ transition. We found the redox-inactive metal only has a weak modulating effect on the potential, which is comparable in strength to that of the ligand environments. Meanwhile, the net charge of the complex, which could be changed along with the redox-inactive metal, has a high impact on the potential and can be unified by protonation, deprotonation, or ligand modification. Although the modulating effect of the redox-inactive metal is not very strong, the linear relationship between the potential and the Lewis acidity is still valid for Mn₄MO₄ complexes. Our results of strong modulating effects for net charge and weak modulating effects for redox-inactive metal fit with the previous experimental observations on Mn₄MO₄ (M = Ca²⁺, Y³⁺, and Gd³⁺) model complexes, and suggest that Ca²⁺ can be structurally and electrochemically replaced with other metal cations, together with proper ligand modifications.

Reversible Structural Isomerization of Nature's Water Oxidation Catalyst Prior to O-O Bond Formation

Photosynthetic water oxidation is catalyzed by a manganese–calcium oxide cluster, which experiences five “S-states” during a light-driven reaction cycle. The unique “distorted chair”-like geometry of the Mn₄CaO₅₍₆₎ cluster shows structural flexibility that has been frequently proposed to involve “open” and “closed”-cubane forms from the S₁ to S₃ states. The isomers are interconvertible in the S₁ and S₂ states, while in the S₃ state, the open-cubane structure is observed to dominate inThermosynechococcus elongatus (cyanobacteria) samples. In this work, using density functional theory calculations, we go beyond the S₃⁺Y_z state to the S₃ⁿY_z^• → S₄⁺Y_z step, and report for the first time that the reversible isomerism, which is suppressed in the S₃⁺Y_z state, is fully recovered in the ensuing S₃ⁿY_z^• state due to the proton release from a manganese-bound water ligand. The altered coordination strength of the manganese–ligand facilitates formation of the closed-cubane form, in a dynamic equilibrium with the open-cubane form. This tautomerism immediately preceding dioxygen formation may constitute the rate limiting step for O₂ formation, and exert a significant influence on the water oxidation mechanism in photosystem II.

EPR Spectroscopy of Cu(II) Complexes: Prediction of g-Tensors Using Double-Hybrid Density Functional Theory

Computational electron paramagnetic resonance (EPR) spectroscopy is an important field of applied quantum chemistry that contributes greatly to connecting spectroscopic observations with the fundamental description of electronic structure for open-shell molecules. However, not all EPR parameters can be predicted accurately and reliably for all chemical systems. Among transition metal ions, Cu(II) centers in inorganic chemistry and biology, and their associated EPR properties such as hyperfine coupling and g-tensors, pose exceptional difficulties for all levels of quantum chemistry. In the present work, we approach the problem of Cu(II) g-tensor calculations using double-hybrid density functional theory (DHDFT). Using a reference set of 18 structurally and spectroscopically characterized Cu(II) complexes, we evaluate a wide range of modern double-hybrid density functionals (DHDFs) that have not been applied previously to this problem. Our results suggest that the current generation of DHDFs consistently and systematically outperform other computational approaches. The B2GP-PLYP and PBE0-DH functionals are singled out as the best DHDFs on average for the prediction of Cu(II) g-tensors. The performance of the different functionals is discussed and suggestions are made for practical applications and future methodological developments.

Successes and challenges in multiscale modelling of artificial metalloenzymes: the case study of POP-Rh2 cyclopropanase

Molecular modelling applications in metalloenzyme design are still scarce due to a series of challenges. On top of that, the simulations of metal-mediated binding and the identification of catalytic competent geometries require both large conformational exploration and simulation of fine electronic properties. Here, we demonstrate how the incorporation of new tools in multiscale strategies, namely substrate diffusion exploration, allows taking a step further. As a showcase, the enantioselective profiles of the most outstanding variants of an artificial Rh₂-based cyclopropanase (GSH, HFF and RFY) developed by Lewis and co-workers (Nat. Commun., 2015, 6, 7789 and Nat. Chem., 2018, 10, 318-324) have been rationalized. DFT calculations on the free-cofactor-mediated process identify the carbene insertion and the cyclopropanoid formation as crucial events, the latter being the enantiodetermining step, which displays up to 8 competitive orientations easily altered by the protein environment. The key intermediates of the reaction were docked into the protein scaffold showing that some mutated residues have direct interaction with the cofactor and/or the co-substrate. These interactions take the form of a direct coordination of Rh in GSH and HFF and a strong hydrophobic patch with the carbene moiety in RFY. Posterior molecular dynamics sustain that the cofactor induces global re-arrangements of the protein. Finally, massive exploration of substrate diffusion, based on the GPathFinder approach, defines this event as the origin of the enantioselectivity in GSH and RFY. For HFF, fine molecular dockings suggest that it is likely related to local interactions upon diffusion. This work shows how modelling of long-range mutations on the catalytic profiles of metalloenzymes may be unavoidable and software simulating substrate diffusion should be applied.

Alteration of the O2-Producing Mn4Ca Cluster in Photosystem II by the Mutation of a Metal Ligand

The O₂-evolving Mn₄Ca cluster in photosystem II (PSII) is arranged as a distorted Mn₃Ca cube that is linked to a fourth Mn ion (denoted as Mn4) by two oxo bridges. The Mn4 and Ca ions are bridged by residue D1-D170. This is also the only residue known to participate in the high-affinity Mn(II) site that participates in the light-driven assembly of the Mn₄Ca cluster. In this study, we use Fourier transform infrared difference spectroscopy to characterize the impact of the D1-D170E mutation. On the basis of analyses of carboxylate and carbonyl stretching modes and the O-H stretching modes of hydrogen-bonded water molecules, we show that this mutation alters the extensive network of hydrogen bonds that surrounds the Mn₄Ca cluster in the same manner as that of many other mutations. It also alters the equilibrium between conformers of the Mn₄Ca cluster in the dark-stable S₁ state so that a high-spin form of the S₂ state is produced during the S₁-to-S₂ transition instead of the low-spin form that gives rise to the S₂ state multiline electron paramagnetic resonance signal. The mutation may also change the coordination mode of the carboxylate group at position 170 to unidentate ligation of Mn4. This is the first mutation of a metal ligand in PSII that substantially impacts the spectroscopic signatures of the Mn₄Ca cluster without substantially eliminating O₂ evolution. The results have significant implications for our understanding of the roles of alternate active/inactive conformers of the Mn₄Ca cluster in the mechanism of O₂ formation.

Redox Isomerism in the S3 State of the Oxygen-Evolving Complex Resolved by Coupled Cluster Theory

The electronic and geometric structures of the water-oxidizing complex of photosystem II in the steps of the catalytic cycle that precede dioxygen evolution remain hotly debated. Recent structural and spectroscopic investigations support contradictory redox formulations for the active-site Mn4 CaOx cofactor in the final metastable S3 state. These range from the widely accepted MnIV 4 oxo-hydroxo model, which presumes that O-O bond formation occurs in the ultimate transient intermediate (S4 ) of the catalytic cycle, to a MnIII 2 MnIV 2 peroxo model representative of the contrasting "early-onset" O-O bond formation hypothesis. Density functional theory energetics of suggested S3 redox isomers are inconclusive because of extreme functional dependence. Here, we use the power of the domain-based local pair natural orbital approach to coupled cluster theory, DLPNO-CCSD(T), to present the first correlated wave function theory calculations of relative stabilities for distinct redox-isomeric forms of the S3 state. Our results enabled us to evaluate conflicting models for the S3 state of the oxygen-evolving complex (OEC) and to quantify the accuracy of lower-level theoretical approaches. Our assessment of the relevance of distinct redox-isomeric forms for the mechanism of biological water oxidation strongly disfavors the scenario of early-onset O-O formation advanced by literal interpretations of certain crystallographic models. This work serves as a case study in the application of modern coupled cluster implementations to redox isomerism problems in oligonuclear transition metal systems.

H2 activation by hydrogenase-inspired NiFe catalyst using frustrated Lewis pair: effect of buffer and halide ion in the heterolytic H-H bond cleavage

Hydrogen is a clean fuel alternative to fossil fuels, and it is vital to develop catalysts for its efficient activation and production. We investigate the reaction mechanism of H₂ activation in an aqueous solution by the recently developed NiFe complex (Ogo et al. Sci. Adv. 2020, 6, eaaz8181) using density functional theory (DFT) calculation. Our computational results showed that H₂ is activated using frustrated Lewis pair. That is, H₂ binds to the Fe site of the NiFe complex, acting as a Lewis acid, while the added buffer, acting as Lewis base, abstracts protons to form a hydride complex. Furthermore, the higher basicity in the proton abstraction reaction characterises reaction more exergonic and lowers the reaction barrier. In addition, in the proton abstraction by the water molecule, the reaction barrier was lowered when anion such as Cl⁻ is in the vicinity of the water. Understanding the chemical species that contribute to the catalytic process in cooperation with the metal catalyst at the atomic level should help to maximise the function of the catalyst.

Hydrogen is a clean fuel alternative to fossil fuels, and it is vital to develop catalysts for its efficient activation and production.