Theoretical evaluation of structural models of the S2 state in the oxygen evolving complex of Photosystem II: protonation states and magnetic interactions

Experimentally Assessing the Electronic Structure and Spin-State Energetics in MnFe Dimers Using 1s3p Resonant Inelastic X-ray Scattering

The synergistic interaction between Mn and Fe centers is investigated via a comprehensive analysis of full 1s3p resonant inelastic X-ray scattering (RIXS) planes at both the Fe and Mn K-edges in a series of homo- and heterometallic molecular systems. Deconvolution of the experimental two-dimensional 1s3p RIXS maps provides insights into the modulation of metal-ligand covalency and variations in the metal multiplet structure induced by subtle electronic structural differences imposed by the presence of the second metal. These modulations in the electronic structure are key toward understanding the reactivity of biological systems with active sites that require heterometallic centers, including MnFe purple acid phosphatases and MnFe ribonucleotide reductases. Herein, we demonstrate the capabilities of 1s3p RIXS to provide information on the excited state energetics in both element- and spin-selective fashion. The contributing excited states are identified and isolated by their multiplicity and π- and σ-contributions, building a conceptual bridge between the electronic structures of metal centers and their reactivity. The ability of the presented 1s3p RIXS methodology to address fundamental questions in transition metal catalysis reactivity is highlighted.

Closing Kok's cycle of nature's water oxidation catalysis

The Mn4CaO5(6) cluster in photosystem II catalyzes water splitting through the Si state cycle (i = 0-4). Molecular O2 is formed and the natural catalyst is reset during the final S3 → (S4) → S0 transition. Only recently experimental breakthroughs have emerged for this transition but without explicit information on the S0-state reconstitution, thus the progression after O2 release remains elusive. In this report, our molecular dynamics simulations combined with density functional calculations suggest a likely missing link for closing the cycle, i.e., restoring the first catalytic state. Specifically, the formation of closed-cubane intermediates with all hexa-coordinate Mn is observed, which would undergo proton release, water dissociation, and ligand transfer to produce the open-cubane structure of the S0 state. Thereby, we theoretically identify the previously unknown structural isomerism in the S0 state that acts as the origin of the proposed structural flexibility prevailing in the cycle, which may be functionally important for nature's water oxidation catalysis.

Water Ligands Regulate the Redox Leveling Mechanism of the Oxygen-Evolving Complex of the Photosystem II

Understanding how water ligands regulate the conformational changes and functionality of the oxygen-evolving complex (OEC) in photosystem II (PSII) throughout the catalytic cycle of oxygen evolution remains a highly intriguing and unresolved challenge. In this study, we investigate the effect of water insertion (WI) on the redox state of the OEC by using the molecular dynamics (MD) and quantum mechanics/molecular mechanics (QM/MM) hybrid methods. We find that water binding significantly reduces the free energy change for proton-coupled electron transfer (PCET) from Mn to YZ•, underscoring the important regulatory role of water binding, which is essential for enabling the OEC redox-leveling mechanism along the catalytic cycle. We propose a water binding mechanism in which WI is thermodynamically favored by the closed-cubane form of the OEC, with water delivery mediated by Ca2+ ligand exchange. Isomerization from the closed- to open-cubane conformation at three post-WI states highlights the importance of the location of the MnIII center in the OEC and the orientation of its Jahn-Teller axis to conformational changes of the OEC, which might be critical for the formation of the O-O bond. These findings reveal a complex interplay between conformational changes in the OEC and the ligand environment during the activation of the OEC by YZ•. Analogous regulatory effects due to water ligand binding are expected to be important for a wide range of catalysts activated by redox state transitions in aqueous environments.

Unveiling the Low-Lying Spin States of [Fe3S4] Clusters via the Extended Broken-Symmetry Method

Photosynthetic water splitting, when synergized with hydrogen production catalyzed by hydrogenases, emerges as a promising avenue for clean and renewable energy. However, theoretical calculations have faced challenges in elucidating the low-lying spin states of iron-sulfur clusters, which are integral components of hydrogenases. To address this challenge, we employ the Extended Broken-Symmetry method for the computation of the cubane-[Fe3S4] cluster within the [FeNi] hydrogenase enzyme. This approach rectifies the error caused by spin contamination, allowing us to obtain the magnetic exchange coupling constant and the energy level of the low-lying state. We find that the Extended Broken-Symmetry method provides more accurate results for differences in bond length and the magnetic coupling constant. This accuracy assists in reconstructing the low-spin ground state force and determining the geometric structure of the ground state. By utilizing the Extended Broken-Symmetry method, we further highlight the significance of the geometric arrangement of metal centers in the cluster's properties and gain deeper insights into the magnetic properties of transition metal iron-sulfur clusters at the reaction centers of hydrogenases. This research illuminates the untapped potential of hydrogenases and their promising role in the future of photosynthesis and sustainable energy production.

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.

Chemistry of a series of heterobimetallic complexes MnIII2(CaII/SrII)X2 (X = Cl-, Br-)

This manuscript describes the syntheses, structures and magnetism of MnIII-CaII/SrII complexes which are compositionally relevant in the context of the oxygen evolving complex (OEC) of photosystem II (PS II). A series of trimetallic tetraoxo complexes containing redox-inactive metal ions CaII or SrII were synthesized using a tetranucleating ligand framework. The structural characteristics of these complexes, with the oxido ligands bridging the redox-inactive metals and the manganese centres, make them particularly relevant to biological and heterogeneous metal-oxido clusters. Electrochemical studies of these compounds show that the reduction potentials are highly dependent upon the Lewis acidity of the redox-inactive metal, identifying the chemical basis for the observed differences in electrochemistry. This correlation provides insights into the role of the CaII/SrII ion in modulating the redox potential of the OEC and of other redox-inactive ions in tuning the redox potentials of other metal-oxide electrocatalysts. Temperature dependent magnetic measurements have also been performed for the complexes.

Interplay of two low-barrier hydrogen bonds in long-distance proton-coupled electron transfer for water oxidation

D1-Tyr161 (TyrZ) forms a low-barrier H-bond with D1-His190 and functions as a redox-active group in photosystem II. When oxidized to the radical form (TyrZ-O•), it accepts an electron from the oxygen-evolving Mn4CaO5 cluster, facilitating an increase in the oxidation state (Sn; n = 0-3). In this study, we investigated the mechanism of how TyrZ-O• drives proton-coupled electron transfer during the S2 to S3 transition using a quantum mechanical/molecular mechanical approach. In response to TyrZ-O• formation and subsequent loss of the low-barrier H-bond, the ligand water molecule at the Ca2+ site (W4) reorients away from TyrZ and donates an H-bond to D1-Glu189 at Mn4 of Mn4CaO5 together with an adjacent water molecule. The H-bond donation to the Mn4CaO5 cluster triggers the release of the proton from the lowest pKa site (W1 at Mn4) along the W1…D1-Asp61 low-barrier H-bond, leading to protonation of D1-Asp61. The interplay of the two low-barrier H-bonds, involving the Ca2+ interface and forming the extended Grotthuss-like network [TyrZ…D1-His190]-[Mn4CaO5]-[W1…D1-Asp61], rather than the direct electrostatic interaction, is likely a basis of the apparent long-distance interaction (11.4 Å) between TyrZ-O• formation and D1-Asp61 protonation.

Structural and energetic insights into Mn-to-Fe substitution in the oxygen-evolving complex

Manganese (Mn) serves as the catalytic center for water splitting in photosystem II (PSII), despite the abundance of iron (Fe) on earth. As a first step toward why Mn and not Fe is employed by Nature in the water oxidation catalyst, we investigated the Fe₄CaO₅ cluster in the PSII protein environment using a quantum mechanical/molecular mechanical (QM/MM) approach, assuming an equivalence between Mn(III/IV) and Fe(II/III). Substituting Mn with Fe resulted in the protonation of μ-oxo bridges at sites O2 and O3 by Arg357 and D1-His337, respectively. While the Mn₄CaO₅ cluster exhibits distinct open- and closed-cubane S₂ conformations, the Fe₄CaO₅ cluster lacks this variability due to an equal spin distribution over sites Fe1 and Fe4. The absence of a low-barrier H-bond between a ligand water molecule (W1) and D1-Asp61 in the Fe₄CaO₅ cluster may underlie its incapability for ligand water deprotonation, highlighting the relevance of Mn in natural water splitting.

Insights into the protonation state and spin structure for the g = 2 multiline electron paramagnetic resonance signal of the oxygen-evolving complex

In photosystem II (PSII), one-electron oxidation of the most stable oxidation state of the Mn4CaO5 cluster (S1) leads to formation of two distinct states, the open-cubane S2 conformation [Mn1(III)Mn2(IV)Mn3(IV)Mn4(IV)] with low spin and the closed-cubane S2 conformation [Mn1(IV)Mn2(IV)Mn3(IV)Mn4(III)] with high spin. In electron paramagnetic resonance (EPR) spectroscopy, the open-cubane S2 conformation exhibits a g = 2 multiline signal. However, its protonation state remains unclear. Here, we investigated the protonation state of the open-cubane S2 conformation by calculating exchange couplings in the presence of the PSII protein environment and simulating the pulsed electron-electron double resonance (PELDOR). When a ligand water molecule, which forms an H-bond with D1-Asp61 (W1), is deprotonated at dangling Mn4(IV), the first-exited energy (34 cm-1) in manifold spin excited states aligns with the observed value in temperature-dependent pulsed EPR analyses, and the PELDOR signal is best reproduced. Consequently, the g = 2 multiline signal observed in EPR corresponds to the open-cubane S2 conformation with the deprotonated W1 (OH-).

Molecular Structure Related to an S = 5/2 High-Spin S2 State Manganese Cluster of Photosystem II Investigated by Q-Band Pulse EPR Spectroscopy

The high-spin S₂ state of the photosynthetic oxygen-evolving cluster Mn₄CaO₅, corresponding to the g = 4.1 signal for X-band electron paramagnetic resonance (EPR), was investigated using Q-band pulsed EPR, which detected a main peak at g = 3.10 and satellite peaks at 5.25, 4.55, and 2.80. We evaluated the spin state as the zero-field splitting of D = 0.465 cm^-1 and E/D = 0.245 with S = 5/2. The temperature dependence of the T₁ relaxation time revealed that the excited-state energy was 28.7 cm^-1 higher than that of the S = 5/2 ground state. By comparing present quantum mechanical (QM) calculation models, a closed-cubane structure with the protonation state of two oxygens, W1 (= OH^-) and W2 (= H₂O), was the most probable structure for the S = 5/2 state. The three-pulse electron spin-echo envelope modulation (ESEEM) detected the nuclear signal, which was assigned to nitrogen as His332 ligated to the Mn1 ion. The obtained hyperfine constant for the nitrogen signal was significantly reduced from that in the S = 1/2 low-spin state. These results indicate that the S = 5/2 spin state arises from the closed-cubane structure.

Coordination Number in High-Spin-Low-Spin Equilibrium in Cluster Models of the S2 State of the Oxygen Evolving Complex

The S2 state of the Oxygen Evolving Complex (OEC) of Photosystem II (PSII) shows high-spin (HS) and low-spin (LS) EPR signals attributed to distinct structures based on computation. Five-coordinate MnIII centers are proposed in these species but are absent in available spectroscopic model complexes. Herein, we report the synthesis, crystal structure, electrochemistry, SQUID magnetometry, and EPR spectroscopy of a MnIIIMnIV3O4 cuboidal complex featuring five-coordinate MnIII. This cluster displays a spin ground state of S = 5/2, while conversion to a six-coordinate Mn upon treatment with water results in a spin state change to S = 1/2. These results demonstrate that coordination number, without dramatic changes within the Mn4O4 core, has a substantial effect on spectroscopy.

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.

Identification of the protonation and oxidation states of the oxygen-evolving complex in the low-dose X-ray crystal structure of photosystem II

In photosystem II (PSII), the O3 and O4 sites of the Mn₄CaO₅ cluster form hydrogen bonds with D1-His337 and a water molecule (W539), respectively. The low-dose X-ray structure shows that these hydrogen bond distances differ between the two homogeneous monomer units (A and B) [Tanaka et al., J. Am Chem. Soc. 2017, 139, 1718]. We investigated the origin of the differences using a quantum mechanical/molecular mechanical (QM/MM) approach. QM/MM calculations show that the short O4-O_W539 hydrogen bond (~2.5 Å) of the B monomer is reproduced when O4 is protonated in the S₁ state. The short O3-Nε_His337 hydrogen bond of the A monomer is due to the formation of a low-barrier hydrogen bond between O3 and doubly-protonated D1-His337 in the overreduced states (S_-1 or S_-2). It seems plausible that the oxidation state differs between the two monomer units in the crystal.

Substitution of Ca2+ and changes in the H-bond network near the oxygen-evolving complex of photosystem II

Ca²⁺, which provides binding sites for ligand water molecules W3 and W4 in the Mn₄CaO₅ cluster, is a prerequisite for O₂ evolution in photosystem II (PSII). We report structural changes in the H-bond network and the catalytic cluster itself upon the replacement of Ca²⁺ with other alkaline earth metals, using a quantum mechanical/molecular mechanical approach. The small radius of Mg²⁺ makes W3 donate an H-bond to D1-Glu189 in Mg²⁺-PSII. If an additional water molecule binds at the large surface of Ba²⁺, it donates H-bonds to D1-Glu189 and the ligand water molecule at the dangling Mn, altering the H-bond network. The potential energy profiles of the H-bond between D1-Tyr161 (TyrZ) and D1-His190 and the interconversion between the open- and closed-cubane S₂ conformations remain substantially unaltered upon the replacement of Ca²⁺. Remarkably, the O5⋯Ca²⁺ distance is shortest among all O5⋯metal distances irrespective of the radius being larger than that of Mg²⁺. Furthermore, Ca²⁺ is the only alkaline earth metal that equalizes the O5⋯metal and O2⋯metal distances and facilitates the formation of the symmetric cubane structure.

Metamorphic oxygen-evolving molecular Ru and Ir catalysts

Today sustainable and clean energy conversion strategies are based on sunlight and the use of water as a source of protons and electrons, in a similar manner as it happens in Photosystem II. To achieve this, the charge separation state induced by light has to be capable of oxidising water by 4 protons and 4 electrons and generating molecular oxygen. This oxidation occurs by the intermediacy of a catalyst capable of finding low-energy pathways via proton-coupled electron transfer steps. The high energy involved in the thermodynamics of water oxidation reaction, coupled with its mechanistic complexity, is responsible for the difficulty of discovering efficient and oxidatively robust molecules capable of achieving such a challenging task. A significant number of Ru coordination complexes have been identified as water oxidation catalysts (WOCs) and are among the best understood from a mechanistic perspective. In this review, we describe the catalytic performance of these complexes and focus our attention on the factors that influence their performance during catalysis, especially in cases where a detailed mechanistic investigation has been carried out. The collective information extracted from all the catalysts studied allows one to identify the key features that govern the complex chemistry associated with the catalytic water oxidation reaction. This includes the stability of trans-O-Ru-O groups, the change in coordination number from CN6 to CN7 at Ru high oxidation states, the ligand flexibility, the capacity to undergo intramolecular proton transfer, the bond strain, the axial ligand substitution, and supramolecular effects. Overall, combining all this information generates a coherent view of this complex chemistry.

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.

Protonation structure of the closed-cubane conformation of the O2-evolving complex in photosystem II

In photosystem II (PSII), one-electron oxidation of the most stable state of the oxygen-evolving Mn4CaO5 cluster (S1) leads to the S2 state formation, Mn1(III)Mn2(IV)Mn3(IV)Mn4(IV) (open-cubane S2) or Mn1(IV)Mn2(IV)Mn3(IV)Mn4(III) (closed-cubane S2). In electron paramagnetic resonance (EPR) spectroscopy, the g = 4.1 signal is not observed in cyanobacterial PSII but in plant PSII, whereas the g = 4.8 signal is observed in cyanobacterial PSII and extrinsic-subunit-depleted plant PSII. Here, we investigated the closed-cubane S2 conformation, a candidate for a higher spin configuration that accounts for g > 4.1 EPR signal, considering all pairwise exchange couplings in the PSII protein environment (i.e. instead of considering only a single exchange coupling between the [Mn3(CaO4)] cubane region and the dangling Mn4 site). Only when a ligand water molecule that forms an H-bond with D1-Asp61 (W1) is deprotonated at dangling Mn4(IV), the g = 4.1 EPR spectra can be reproduced using the cyanobacterial PSII crystal structure. The closed-cubane S2 is less stable than the open-cubane S2 in cyanobacterial PSII, which may explain why the g = 4.1 EPR signal is absent in cyanobacterial PSII.

Pulsed Multifrequency Electron Paramagnetic Resonance Spectroscopy Reveals Key Branch Points for One- vs Two-Electron Reactivity in Mn/Fe Proteins

Traditionally, the ferritin-like superfamily of proteins was thought to exclusively use a diiron active site in catalyzing a diverse array of oxygen-dependent reactions. In recent years, novel redox-active cofactors featuring heterobimetallic Mn/Fe active sites have been discovered in both the radical-generating R2 subunit of class Ic (R2c) ribonucleotide reductases (RNRs) and the related R2-like ligand-binding oxidases (R2lox). However, the protein-specific factors that differentiate the radical reactivity of R2c from the C-H activation reactions of R2lox remain unknown. In this work, multifrequency pulsed electron paramagnetic resonance (EPR) spectroscopy and ligand hyperfine techniques in conjunction with broken-symmetry density functional theory calculations are used to characterize the molecular and electronic structures of two EPR-active intermediates trapped during aerobic assembly of the R2lox Mn/Fe cofactor. A Mn^III(μ-O)(μ-OH)Fe^III species is identified as the first EPR-active species and represents a common state between the two classes of redox-active Mn/Fe proteins. The species downstream from the Mn^III(μ-O)(μ-OH)Fe^III state exhibits unique EPR properties, including unprecedented spectral breadth and isotope-dependent g-tensors, which are attributed to a weakly coupled, hydrogen-bonded Mn^III(μ-OH)Fe^III species. This final intermediate precedes formation of the Mn^III/Fe^III resting state and is suggested to be relevant to understanding the endogenous reactivity of R2lox.

Requirement of Chloride for the Downhill Electron Transfer Pathway from the Water-Splitting Center in Natural Photosynthesis

In photosystem II (PSII), Cl^- is a prerequisite for the second flash-induced oxidation of the Mn₄CaO₅ cluster (the S₂ to S₃ transition). We report proton transfer from the substrate water molecule via D1-Asp61 and electron transfer via redox-active D1-Tyr161 (TyrZ) to the chlorophyll pair in Cl^--depleted PSII using a quantum mechanical/molecular mechanical approach. The low-barrier H-bond formation between the substrate water molecule and D1-Asp61 remained unaffected upon the depletion of Cl^-. However, the binding site, D2-Lys317, formed a salt bridge with D1-Asp61, leading to the inhibition of the subsequent proton transfer. Remarkably, the redox potential (E_m) of S₂/S₃ increased significantly, making electron transfer from S₂ to TyrZ energetically uphill, as observed in Ca²⁺-depleted PSII. The uphill electron transfer pathway was induced by the significant increase in E_m(S₂/S₃) caused by the loss of charge compensation for D2-Lys317 upon the depletion of Cl^-, whereas it was induced by the significant decrease in E_m(TyrZ) caused by the rearrangement of the water molecules at the Ca²⁺ binding moiety upon the depletion of Ca²⁺.

Structural and dynamic insights into Mn4Ca cluster-depleted Photosystem II

In the first steps of natural oxygenic photosynthesis, sunlight is used to oxidize water molecules to protons, electrons and molecular oxygen. This reaction takes place on the Mn₄Ca cluster located in the reaction centre of Photosystem II (PSII), where the cluster is assembled and continuously repaired through a process known as photoactivation. Understanding the molecular details of such a process has important implications in different fields, in particular inspiring synthesis and repair strategies for artificial photosynthesis devices. In this regard, a detailed structural and dynamic characterization of Photosystem II lacking a Mn₄Ca cluster, namely apo PSII, is a prerequisite for the full comprehension of the photoactivation. Recently, the structure of the apo PSII was resolved at 2.55 Å resolution [Zhang et al., eLife, 2017, 6, e26933], suggesting a pre-organized structure of the protein cavity hosting the cluster. Anyway, the question of whether these findings are a feature of the method used remains open. Here, by means of classical Molecular Dynamics simulations, we characterized the structural and dynamic features of the apo PSII for different protonation states of the cluster cavity. Albeit an overall conformational stability common to all investigated systems, we found significant deviations in the conformation of the side chains of the active site with respect to the X-ray positions. Our findings suggest that not all residues acting as Mn ligands are pre-organized prior to the Mn₄Ca formation and previous local conformational changes are required in order to bind the first Mn ion in the high-affinity binding site.

The exchange of the fast substrate water in the S2 state of photosystem II is limited by diffusion of bulk water through channels - implications for the water oxidation mechanism

The molecular oxygen we breathe is produced from water-derived oxygen species bound to the Mn4CaO5 cluster in photosystem II (PSII). Present research points to the central oxo-bridge O5 as the 'slow exchanging substrate water (Ws)', while, in the S2 state, the terminal water ligands W2 and W3 are both discussed as the 'fast exchanging substrate water (Wf)'. A critical point for the assignment of Wf is whether or not its exchange with bulk water is limited by barriers in the channels leading to the Mn4CaO5 cluster. In this study, we measured the rates of H2 16O/H2 18O substrate water exchange in the S2 and S3 states of PSII core complexes from wild-type (WT) Synechocystis sp. PCC 6803, and from two mutants, D1-D61A and D1-E189Q, that are expected to alter water access via the Cl1/O4 channels and the O1 channel, respectively. We found that the exchange rates of Wf and Ws were unaffected by the E189Q mutation (O1 channel), but strongly perturbed by the D61A mutation (Cl1/O4 channel). It is concluded that all channels have restrictions limiting the isotopic equilibration of the inner water pool near the Mn4CaO5 cluster, and that D61 participates in one such barrier. In the D61A mutant this barrier is lowered so that Wf exchange occurs more rapidly. This finding removes the main argument against Ca-bound W3 as fast substrate water in the S2 state, namely the indifference of the rate of Wf exchange towards Ca/Sr substitution.

S3 State Models of Nature's Water Oxidizing Complex: Analysis of Bonding and Magnetic Exchange Pathways, Assessment of Experimental Electron Paramagnetic Resonance Data, and Implications for the Water Oxidation Mechanism

Broken symmetry density functional theory (BS-DFT) calculations on large models of Nature's water oxidizing complex (WOC) are used to investigate the electronic structure and associated magnetic interactions of this key intermediate state. The electronic origins of the ferromagnetic and antiferromagnetic couplings between neighboring Mn ions are investigated and illustrated by using corresponding orbital transformations. Protonation of the O4 and/or O6 atoms leads to large variation in the distribution of spin around the complex with associated changes in its magnetic resonance properties. Models for Sr²⁺ exchange and methanol addition indicate minor perturbations reflected in slightly altered spin projection coefficients for the Mn₁ and Mn₂ ions. These are shown to account for the observed changes observed experimentally via electron paramagnetic resonance methods and suggest a reinterpretation of the experimental findings. By comparison with experimental determinations, we show that the spin projections and resulting calculated ⁵⁵Mn hyperfine couplings support the open cubane form of an oxo (O5)-hydroxo (O6) cluster in all cases with no need to invoke a closed cubane intermediate. The implications of these findings for the water oxidation mechanism are discussed.

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.

CaMn3 IV O4 Cubane Models of the Oxygen-Evolving Complex: Spin Ground States S<9/2 and the Effect of Oxo Protonation

We report the single crystal XRD and MicroED structure, magnetic susceptibility, and EPR data of a series of CaMn₃ ^IV O₄ and YMn₃ ^IV O₄ complexes as structural and spectroscopic models of the cuboidal subunit of the oxygen-evolving complex (OEC). The effect of changes in heterometal identity, cluster geometry, and bridging oxo protonation on the spin-state structure was investigated. In contrast to previous computational models, we show that the spin ground state of CaMn₃ ^IV O₄ complexes and variants with protonated oxo moieties need not be S=9/2. Desymmetrization of the pseudo-C₃ -symmetric Ca(Y)Mn₃ ^IV O₄ core leads to a lower S=5/2 spin ground state. The magnitude of the magnetic exchange coupling is attenuated upon oxo protonation, and an S=3/2 spin ground state is observed in CaMn₃ ^IV O₃ (OH). Our studies complement the observation that the interconversion between the low-spin and high-spin forms of the S₂ state is pH-dependent, suggesting that the (de)protonation of bridging or terminal oxygen atoms in the OEC may be connected to spin-state changes.

Bonding and Magnetic Exchange Pathways in Nature's Water-Oxidizing Complex

The nature of the bonding and magnetic exchange pathways of the water-oxidizing complex of photosystem 2 is explored using broken symmetry density functional theory. The electronic structure and superexchange pathways are illustrated and analyzed using corresponding orbitals and intrinsic bond orbitals. These demonstrate a dominating influence on the bonding and magnetic interactions by both the geometrical structure of the Mn₄CaO₅ core complex and the ionic interactions of the oxo bridges with the neighboring Ca²⁺ ion. The demonstrated ionic nature of the Ca²⁺ bonds is proposed to contribute to the stabilization of the oxygen atoms participating in O-O bond formation.

Role of redox-inactive metals in controlling the redox potential of heterometallic manganese-oxido clusters

Photosystem II (PSII) contains Ca2+, which is essential to the oxygen-evolving activity of the catalytic Mn4CaO5 complex. Replacement of Ca2+ with other redox-inactive metals results in a loss/decrease of oxygen-evolving activity. To investigate the role of Ca2+ in this catalytic reaction, we investigate artificial Mn3[M]O2 clusters redox-inactive metals  [M] ([M]  = Mg2+, Ca2+, Zn2+, Sr2+, and Y3+), which were synthesized by Tsui et al. (Nat Chem 5:293, 2013). The experimentally measured redox potentials (Em) of these clusters are best described by the energy of their highest occupied molecular orbitals. Quantum chemical calculations showed that the valence of metals predominantly affects Em(MnIII/IV), whereas the ionic radius of metals affects Em(MnIII/IV) only slightly.

Orientational Jahn-Teller Isomerism in the Dark-Stable State of Nature's Water Oxidase

The tetramanganese–calcium cluster of the oxygen‐evolving complex of photosystem II adopts electronically and magnetically distinct but interconvertible valence isomeric forms in its first light‐driven oxidized catalytic state, S₂. This bistability is implicated in gating the final catalytic states preceding O−O bond formation, but it is unknown how the biological system enables its emergence and controls its effect. Here we show that the Mn₄CaO₅ cluster in the resting (dark‐stable) S₁ state adopts orientational Jahn–Teller isomeric forms arising from a directional change in electronic configuration of the “dangler” Mn^III ion. The isomers are consistent with available structural data and explain previously unresolved electron paramagnetic resonance spectroscopic observations on the S₁ state. This unique isomerism in the resting state is shown to be the electronic origin of valence isomerism in the S₂ state, establishing a functional role of orientational Jahn–Teller isomerism unprecedented in biological or artificial catalysis.

Two Distinct Oxygen-Radical Conformations in the X-ray Free Electron Laser Structures of Photosystem II

We report the existence of two distinct oxygen-radical-containing Mn₄CaO_5/6 conformations with short O···O bonds in the crystal structures of the oxygen-evolving enzyme photosystem II (PSII), obtained using an X-ray free electron laser (XFEL). A short O···O distance of <2.3 Å between the O4 site of the Mn₄CaO₅ complex and the adjacent water molecule (W539) in the proton-conducting O4-water chain was observed in the second flash-induced (2F) XFEL structure (2F-XFEL), which may correspond to S₃. By use of a quantum mechanical/molecular mechanical approach, the OH^• formation at W539 and the short O4···O_W539 distance (<2.3 Å) were reproduced in S₂ and S₃ with reduced Mn1(III), which lacks the additional sixth water molecule O6. As the O^•- formation at O6 and the short O5···O6 distance (1.9 Å) have been reported in another 2F-XFEL structure with reduced Mn4(III), two distinct oxygen-radical conformations exist in the 2F-XFEL crystals.

A Force Field for a Manganese-Vanadium Water Oxidation Catalyst: Redox Potentials in Solution as Showcase

We present a molecular mechanics force field in AMBER format for the mixed-valence manganese vanadium oxide cluster [Mn4V4O17(OAc)3]3−—a synthetic analogue of the oxygen-evolving complex that catalyzes the water oxidation reaction in photosystem II—with parameter sets for two different oxidation states. Most force field parameters involving metal atoms have been newly parametrized and the harmonic terms refined using hybrid quantum mechanics/molecular mechanics reference simulations, although some parameters were adapted from pre-existing force fields of vanadate cages and manganese oxo dimers. The characteristic Jahn–Teller distortions of d4 MnIII ions in octahedral environments are recovered by the force field. As an application, the developed parameters have been used to calculate the redox potential of the [MnIIIMn3IV] ⇌ [Mn4IV]+e− half-reaction in acetonitrile by means of Marcus theory.

Successes, challenges, and opportunities for quantum chemistry in understanding metalloenzymes for solar fuels research

Overview of the rich and diverse contributions of quantum chemistry to understanding the structure and function of the biological archetypes for solar fuel research, photosystem II and hydrogenases.

A Theoretical Study of the Recently Suggested MnVII Mechanism for O-O Bond Formation in Photosystem II

The mechanism for water oxidation in photosystem II has been a major topic for several decades. The active catalyst has four manganese atoms connected by bridging oxo bonds, in a complex termed the oxygen-evolving complex (OEC), which also includes a calcium atom. The O-O bond of oxygen is formed after absorption of four photons in a state of the OEC termed S₄. There has been essential consensus that in the S₄ state, all manganese atoms are in the Mn(IV) oxidation state. However, recently there has been a suggestion that one of the atoms is in the Mn(VII) state. In the present computational study, the feasibility of that proposal has been investigated. It is here shown that the mechanism involving Mn(VII) has a much higher barrier for forming O₂ than the previous proposal with four Mn(IV) atoms.

pKa of the ligand water molecules in the oxygen-evolving Mn4CaO5 cluster in photosystem II

Release of the protons from the substrate water molecules is prerequisite for O₂ evolution in photosystem II (PSII). Proton-releasing water molecules with low pK_a values at the catalytic moiety can be the substrate water molecules. In some studies, one of the ligand water molecules, W2, is regarded as OH^-. However, the PSII crystal structure shows neither proton acceptor nor proton-transfer pathway for W2, which is not consistent with the assumption of W2 = OH^-. Here we report the pK_a values of the four ligand water molecules, W1 and W2 at Mn4 and W3 and W4 at Ca²⁺, of the Mn₄CaO₅ cluster. pK_a(W1) ≈ pK_a(W2) << pK_a(W3) ≈ pK_a(W4) in the Mn₄CaO₅ cluster in water. However, pK_a(W1) ≈ pK_a(D1-Asp61) << pK_a(W2) in the PSII protein environment. These results suggest that in PSII, deprotonation of W2 is energetically disfavored as far as W1 exists.

Current Understanding of the Mechanism of Water Oxidation in Photosystem II and Its Relation to XFEL Data

The investigation of water oxidation in photosynthesis has remained a central topic in biochemical research for the last few decades due to the importance of this catalytic process for technological applications. Significant progress has been made following the 2011 report of a high-resolution X-ray crystallographic structure resolving the site of catalysis, a protein-bound Mn4CaOxcomplex, which passes through ≥5 intermediate states in the water-splitting cycle. Spectroscopic techniques complemented by quantum chemical calculations aided in understanding the electronic structure of the cofactor in all (detectable) states of the enzymatic process. Together with isotope labeling, these techniques also revealed the binding of the two substrate water molecules to the cluster. These results are described in the context of recent progress using X-ray crystallography with free-electron lasers on these intermediates. The data are instrumental for developing a model for the biological water oxidation cycle.

Untangling the sequence of events during the S2 → S3 transition in photosystem II and implications for the water oxidation mechanism

In oxygenic photosynthesis, light-driven oxidation of water to molecular oxygen is carried out by the oxygen-evolving complex (OEC) in photosystem II (PS II). Recently, we reported the room-temperature structures of PS II in the four (semi)stable S-states, S1, S2, S3, and S0, showing that a water molecule is inserted during the S2 → S3 transition, as a new bridging O(H)-ligand between Mn1 and Ca. To understand the sequence of events leading to the formation of this last stable intermediate state before O2 formation, we recorded diffraction and Mn X-ray emission spectroscopy (XES) data at several time points during the S2 → S3 transition. At the electron acceptor site, changes due to the two-electron redox chemistry at the quinones, QA and QB, are observed. At the donor site, tyrosine YZ and His190 H-bonded to it move by 50 µs after the second flash, and Glu189 moves away from Ca. This is followed by Mn1 and Mn4 moving apart, and the insertion of OX(H) at the open coordination site of Mn1. This water, possibly a ligand of Ca, could be supplied via a "water wheel"-like arrangement of five waters next to the OEC that is connected by a large channel to the bulk solvent. XES spectra show that Mn oxidation (τ of ∼350 µs) during the S2 → S3 transition mirrors the appearance of OX electron density. This indicates that the oxidation state change and the insertion of water as a bridging atom between Mn1 and Ca are highly correlated.

Rationalizing the Geometries of the Water Oxidising Complex in the Atomic Resolution, Nominal S3 State Crystal Structures of Photosystem II

Three atomic resolution crystal structures of Photosystem II, in the double flashed, nominal S₃ intermediate state of its Mn₄ Ca Water Oxidising Complex (WOC), have now been presented, at 2.25, 2.35 and 2.08 Å resolution. Although very similar overall, the S₃ structures differ within the WOC catalytic site. The 2.25 Å structure contains only one oxy species (O5) in the WOC cavity, weakly associated with Mn centres, similar to that in the earlier 1.95 Å S₁ structure. The 2.35 Å structure shows two such species (O5, O6), with the Mn centres and O5 positioned as in the 2.25 Å structure and O5-O6 separation of ∼1.5 Å. In the latest S₃ variant, two oxy species are also seen (O5, Ox), with the Ox group appearing only in S₃ , closely ligating one Mn, with O5-Ox separation <2.1 Å. The O5 and O6/Ox groups were proposed to be substrate water derived species. Recently, Petrie et al. (Chem. Phys. Chem., 2017) presented large scale Quantum Chemical modelling of the 2.25 Å structure, quantitatively explaining all significant features within the WOC region. This, as in our earlier studies, assumed a 'low' Mn oxidation paradigm (mean S₁ Mn oxidation level of +3.0, Petrie et al., Angew. Chem. Int. Ed., 2015), rather than a 'high' oxidation model (mean S₁ oxidation level of +3.5). In 2018 we showed (Chem. Phys. Chem., 2018) this oxidation state assumption predicted two energetically close S₃ structural forms, one with the metal centres and O5 (as OH^- ) positioned as in the 2.25 Å structure, and the other with the metals similarly placed, but with O5 (as H₂ O) located in the O6 position of the 2.35 Å structure. The 2.35 Å two flashed structure was likely a crystal superposition of two such forms. Here we show, by similar computational analysis, that the latest 2.08 Å S₃ structure is also a likely superposition of forms, but with O5 (as OH^- ) occupying either the O5 or Ox positions in the WOC cavity. This highlights a remarkable structural 'lability' of the WOC centre in the S₃ state, which is likely catalytically relevant to its water splitting function.