Unravelling Mn4Ca cluster vibrations in the S1, S2 and S3 states of the Kok-Joliot cycle of photosystem II

Vibrational spectroscopy serves as a powerful tool for characterizing intermediate states within the Kok-Joliot cycle. In this study, we employ a QM/MM molecular dynamics framework to calculate the room temperature infrared absorption spectra of the S1, S2, and S3 states via the Fourier transform of the dipole time auto-correlation function. To better analyze the computational data and assign spectral peaks, we introduce an approach based on dipole-dipole correlation function of cluster moieties of the reaction center. Our analysis reveals variation in the infrared signature of the Mn4Ca cluster along the Kok-Joliot cycle, attributed to its increasing symmetry and rigidity resulting from the rising oxidation state of the Mn ions. Furthermore, we successfully assign the debated contributions in the frequency range around 600 cm-1. This computational methodology provides valuable insights for deciphering experimental infrared spectra and understanding the water oxidation process in both biological and artificial systems.

A Vision for the Future of Multiscale Modeling

The rise of modern computer science enabled physical chemistry to make enormous progresses in understanding and harnessing natural and artificial phenomena. Nevertheless, despite the advances achieved over past decades, computational resources are still insufficient to thoroughly simulate extended systems from first principles. Indeed, countless biological, catalytic and photophysical processes require ab initio treatments to be properly described, but the breadth of length and time scales involved makes it practically unfeasible. A way to address these issues is to couple theories and algorithms working at different scales by dividing the system into domains treated at different levels of approximation, ranging from quantum mechanics to classical molecular dynamics, even including continuum electrodynamics. This approach is known as multiscale modeling and its use over the past 60 years has led to remarkable results. Considering the rapid advances in theory, algorithm design, and computing power, we believe multiscale modeling will massively grow into a dominant research methodology in the forthcoming years. Hereby we describe the main approaches developed within its realm, highlighting their achievements and current drawbacks, eventually proposing a plausible direction for future developments considering also the emergence of new computational techniques such as machine learning and quantum computing. We then discuss how advanced multiscale modeling methods could be exploited to address critical scientific challenges, focusing on the simulation of complex light-harvesting processes, such as natural photosynthesis. While doing so, we suggest a cutting-edge computational paradigm consisting in performing simultaneous multiscale calculations on a system allowing the various domains, treated with appropriate accuracy, to move and extend while they properly interact with each other. Although this vision is very ambitious, we believe the quick development of computer science will lead to both massive improvements and widespread use of these techniques, resulting in enormous progresses in physical chemistry and, eventually, in our society.

Theoretical elucidation of the structure, bonding, and reactivity of the CaMn4Ox clusters in the whole Kok cycle for water oxidation embedded in the oxygen evolving center of photosystem II. New molecular and quantum insights into the mechanism of the O-O bond formation

This paper reviews our historical developments of broken-symmetry (BS) and beyond BS methods that are applicable for theoretical investigations of metalloenzymes such as OEC in PSII. The BS hybrid DFT (HDFT) calculations starting from high-resolution (HR) XRD structure in the most stable S1 state have been performed to elucidate structure and bonding of whole possible intermediates of the CaMn4Ox cluster (1) in the Si (i = 0 ~ 4) states of the Kok cycle. The large-scale HDFT/MM computations starting from HR XRD have been performed to elucidate biomolecular system structures which are crucial for examination of possible water inlet and proton release pathways for water oxidation in OEC of PSII. DLPNO CCSD(T0) computations have been performed for elucidation of scope and reliability of relative energies among the intermediates by HDFT. These computations combined with EXAFS, XRD, XFEL, and EPR experimental results have elucidated the structure, bonding, and reactivity of the key intermediates, which are indispensable for understanding and explanation of the mechanism of water oxidation in OEC of PSII. Interplay between theory and experiments have elucidated important roles of four degrees of freedom, spin, charge, orbital, and nuclear motion for understanding and explanation of the chemical reactivity of 1 embedded in protein matrix, indicating the participations of the Ca(H2O)n ion and tyrosine(Yz)-O radical as a one-electron acceptor for the O-O bond formation. The Ca-assisted Yz-coupled O-O bond formation mechanisms for water oxidation are consistent with recent XES and very recent time-resolved SFX XFEL and FTIR results.

How chloride functions to enable proton conduction in photosynthetic water oxidation: Time-resolved kinetics of intermediates (S-states) in vivo and bromide substitution

Chloride (Cl-) is essential for O2 evolution during photosynthetic water oxidation. Two chlorides near the water-oxidizing complex (WOC) in Photosystem II (PSII) structures from Thermosynechococcus elongatus (and T. vulcanus) have been postulated to transfer protons generated from water oxidation. We monitored four criteria: primary charge separation flash yield (P* → P+QA-), rates of water oxidation steps (S-states), rate of proton evolution, and flash O2 yield oscillations by measuring chlorophyll variable fluorescence (P* quenching), pH-sensitive dye changes, and oximetry. Br-substitution slows and destabilizes cellular growth, resulting from lower light-saturated O2 evolution rate (-20 %) and proton release (-36 % ΔpH gradient). The latter implies less ATP production. In Br- cultures, protonogenic S-state transitions (S2 → S3 → S0') slow with increasing light intensity and during O2/water exchange (S0' → S0 → S1), while the non-protonogenic S1 → S2 transition is kinetically unaffected. As flash rate increases in Cl- cultures, both rate and extent of acidification of the lumen increase, while charge recombination is suppressed relative to Br-. The Cl- advantage in rapid proton escape from the WOC to lumen is attributed to correlated ion-pair movement of H3O+Cl- in dry water channels vs. separated Br- and H+ ion movement through different regions (>200-fold difference in Bronsted acidities). By contrast, at low flash rates a previously unreported reversal occurs that favors Br- cultures for both proton evolution and less PSII charge recombination. In Br- cultures, slower proton transfer rate is attributed to stronger ion-pairing of Br- with AA residues lining the water channels. Both anions charge-neutralize protons and shepherd them to the lumen using dry aqueous channels.

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.

The electron-proton bottleneck of photosynthetic oxygen evolution

Photosynthesis fuels life on Earth by storing solar energy in chemical form. Today's oxygen-rich atmosphere has resulted from the splitting of water at the protein-bound manganese cluster of photosystem II during photosynthesis. Formation of molecular oxygen starts from a state with four accumulated electron holes, the S4 state-which was postulated half a century ago1 and remains largely uncharacterized. Here we resolve this key stage of photosynthetic O2 formation and its crucial mechanistic role. We tracked 230,000 excitation cycles of dark-adapted photosystems with microsecond infrared spectroscopy. Combining these results with computational chemistry reveals that a crucial proton vacancy is initally created through gated sidechain deprotonation. Subsequently, a reactive oxygen radical is formed in a single-electron, multi-proton transfer event. This is the slowest step in photosynthetic O2 formation, with a moderate energetic barrier and marked entropic slowdown. We identify the S4 state as the oxygen-radical state; its formation is followed by fast O-O bonding and O2 release. In conjunction with previous breakthroughs in experimental and computational investigations, a compelling atomistic picture of photosynthetic O2 formation emerges. Our results provide insights into a biological process that is likely to have occurred unchanged for the past three billion years, which we expect to support the knowledge-based design of artificial water-splitting systems.

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.

Crystal structures of photosystem II from a cyanobacterium expressing psbA2 in comparison to psbA3 reveal differences in the D1 subunit

Three psbA genes (psbA₁, psbA₂, and psbA₃) encoding the D1 subunit of photosystem II (PSII) are present in the thermophilic cyanobacterium Thermosynechococcus elongatus and are expressed differently in response to changes in the growth environment. To clarify the functional differences of the D1 protein expressed from these psbA genes, PSII dimers from two strains, each expressing only one psbA gene (psbA₂ or psbA₃), were crystallized, and we analyzed their structures at resolutions comparable to previously studied PsbA1-PSII. Our results showed that the hydrogen bond between pheophytin/D1 (Pheo_D1) and D1-130 became stronger in PsbA2- and PsbA3-PSII due to change of Gln to Glu, which partially explains the increase in the redox potential of Pheo_D1 observed in PsbA3. In PsbA2, one hydrogen bond was lost in Pheo_D1 due to the change of D1-Y147F, which may explain the decrease in stability of Pheo_D1 in PsbA2. Two water molecules in the Cl-1 channel were lost in PsbA2 due to the change of D1-P173M, leading to the narrowing of the channel, which may explain the lower efficiency of the S-state transition beyond S₂ in PsbA2-PSII. In PsbA3-PSII, a hydrogen bond between D1-Ser270 and a sulfoquinovosyl-diacylglycerol molecule near Q_B disappeared due to the change of D1-Ser270 in PsbA1 and PsbA2 to D1-Ala270. This may result in an easier exchange of bound Q_B with free plastoquinone, hence an enhancement of oxygen evolution in PsbA3-PSII due to its high Q_B exchange efficiency. These results provide a structural basis for further functional examination of the three PsbA variants.

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.

Molecular Principles of Redox-Coupled Protonation Dynamics in Photosystem II

Photosystem II (PSII) catalyzes light-driven water oxidization, releasing O₂ into the atmosphere and transferring the electrons for the synthesis of biomass. However, despite decades of structural and functional studies, the water oxidation mechanism of PSII has remained puzzling and a major challenge for modern chemical research. Here, we show that PSII catalyzes redox-triggered proton transfer between its oxygen-evolving Mn₄O₅Ca cluster and a nearby cluster of conserved buried ion-pairs, which are connected to the bulk solvent via a proton pathway. By using multi-scale quantum and classical simulations, we find that oxidation of a redox-active Tyr_z (Tyr161) lowers the reaction barrier for the water-mediated proton transfer from a Ca²⁺-bound water molecule (W3) to Asp61 via conformational changes in a nearby ion-pair (Asp61/Lys317). Deprotonation of this W3 substrate water triggers its migration toward Mn1 to a position identified in recent X-ray free-electron laser (XFEL) experiments [Ibrahim et al. Proc. Natl. Acad. Sci. USA 2020, 117, 12,624-12,635]. Further oxidation of the Mn₄O₅Ca cluster lowers the proton transfer barrier through the water ligand sphere of the Mn₄O₅Ca cluster to Asp61 via a similar ion-pair dissociation process, while the resulting Mn-bound oxo/oxyl species leads to O₂ formation by a radical coupling mechanism. The proposed redox-coupled protonation mechanism shows a striking resemblance to functional motifs in other enzymes involved in biological energy conversion, with an interplay between hydration changes, ion-pair dynamics, and electric fields that modulate the catalytic barriers.

Structural dynamics in the water and proton channels of photosystem II during the S2 to S3 transition

Light-driven oxidation of water to molecular oxygen is catalyzed by the oxygen-evolving complex (OEC) in Photosystem II (PS II). This multi-electron, multi-proton catalysis requires the transport of two water molecules to and four protons from the OEC. A high-resolution 1.89 Å structure obtained by averaging all the S states and refining the data of various time points during the S₂ to S₃ transition has provided better visualization of the potential pathways for substrate water insertion and proton release. Our results indicate that the O1 channel is the likely water intake pathway, and the Cl1 channel is the likely proton release pathway based on the structural rearrangements of water molecules and amino acid side chains along these channels. In particular in the Cl1 channel, we suggest that residue D1-E65 serves as a gate for proton transport by minimizing the back reaction. The results show that the water oxidation reaction at the OEC is well coordinated with the amino acid side chains and the H-bonding network over the entire length of the channels, which is essential in shuttling substrate waters and protons.

Mechanism of Oxygen Evolution and Mn4CaO5 Cluster Restoration in the Natural Water-Oxidizing Catalyst

Water oxidation occurring in the first steps of natural oxygenic photosynthesis is catalyzed by the pigment/protein complex Photosystem II. This process takes place on the Mn₄Ca cluster located in the core of Photosystem II and proceeds along the five steps (S₀-S₄) of the so-called Kok-Joliot cycle until the release of molecular oxygen. The catalytic cycle can therefore be started afresh through insertion of a new water molecule. Here, combining quantum mechanics/molecular mechanics simulations and minimum energy path calculations, we characterized on different spin surfaces the events occurring in the last sector of the catalytic cycle from structural, electronic, and thermodynamic points of view. We found that the process of oxygen evolution and water insertion can be described well by a two-step mechanism, with oxygen release being the rate-limiting step of the process. Moreover, our results allow us to identify the upcoming water molecule required to regenerate the initial structure of the Mn₄Ca cluster in the S₀ state. The insertion of the water molecule was found to be coupled with the transfer of a proton to a neighboring hydroxide ion, thus resulting in the reconstitution of the most widely accepted model of the S₀ state.

Proton exit pathways surrounding the oxygen evolving complex of photosystem II

Photosystem II allows water to be the primary electron source for the photosynthetic electron transfer chain. Water is oxidized to dioxygen at the Oxygen Evolving Complex (OEC), a Mn₄CaO₅ inorganic core embedded on the lumenal side of PSII. Water-filled channels surrounding the OEC must bring in substrate water molecules, remove the product protons to the lumen, and may transport the product oxygen. Three water-filled channels, denoted large, narrow, and broad, extend from the OEC towards the aqueous surface more than 15 Å away. However, the role of each pathway in the transport in and out of the OEC is yet to be established. Here, we combine Molecular Dynamics (MD), Multi Conformation Continuum Electrostatics (MCCE) and Network Analysis to compare and contrast the three potential proton transfer paths. Hydrogen bond network analysis shows that near the OEC the waters are highly interconnected with similar free energy for hydronium at all locations. The paths diverge as they move towards the lumen. The water chain in the broad channel is better connected than in the narrow and large channels, where disruptions in the network are observed approximately 10 Å from the OEC. In addition, the barrier for hydronium translocation is lower in the broad channel. Thus, a proton released from any location on the OEC can access all paths, but the likely exit to the lumen passes through PsbO via the broad channel.

Activation energies for two steps in the S2→ S3 transition of photosynthetic water oxidation from time-resolved single-frequency infrared spectroscopy

The mechanism of water oxidation by the Photosystem II (PSII) protein-cofactor complex is of high interest, but specifically, the crucial coupling of protonation dynamics to electron transfer (ET) and dioxygen chemistry remains insufficiently understood. We drove spinach-PSII membranes by nanosecond-laser flashes synchronously through the water-oxidation cycle and traced the PSII processes by time-resolved single-frequency infrared (IR) spectroscopy in the spectral range of symmetric carboxylate vibrations of protein side chains. After the collection of IR-transients from 100 ns to 1 s, we analyzed the proton-removal step in the S₂ ⇒ S₃ transition, which precedes the ET that oxidizes the Mn₄CaO_x-cluster. Around 1400 cm^-1, pronounced changes in the IR-transients reflect this pre-ET process (∼40 µs at 20 °C) and the ET step (∼300 µs at 20 °C). For transients collected at various temperatures, unconstrained multi-exponential simulations did not provide a coherent set of time constants, but constraining the ET time constants to previously determined values solved the parameter correlation problem and resulted in an exceptionally high activation energy of 540 ± 30 meV for the pre-ET step. We assign the pre-ET step to deprotonation of a group that is re-protonated by accepting a proton from the substrate-water, which binds concurrently with the ET step. The analyzed IR-transients disfavor carboxylic-acid deprotonation in the pre-ET step. Temperature-dependent amplitudes suggest thermal equilibria that determine how strongly the proton-removal step is reflected in the IR-transients. Unexpectedly, the proton-removal step is only weakly reflected in the 1400 cm^-1 transients of PSII core complexes of a thermophilic cyanobacterium (T. elongatus).

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.

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.

Pivotal role of the redox-active tyrosine in driving the water splitting catalyzed by photosystem II

TyrZ oxidation state triggers hydrogen bond modification in the water oxidation catalysis.

Water oxidation in photosystem II

Biological water oxidation, performed by a single enzyme, photosystem II, is a central research topic not only in understanding the photosynthetic apparatus but also for the development of water splitting catalysts for technological applications. Great progress has been made in this endeavor following the report of a high-resolution X-ray crystallographic structure in 2011 resolving the cofactor site (Umena et al. in Nature 473:55-60, 2011), a tetra-manganese calcium complex. The electronic properties of the protein-bound water oxidizing Mn4OxCa complex are crucial to understand its catalytic activity. These properties include: its redox state(s) which are tuned by the protein matrix, the distribution of the manganese valence and spin states and the complex interactions that exist between the four manganese ions. In this short review we describe how magnetic resonance techniques, particularly EPR, complemented by quantum chemical calculations, have played an important role in understanding the electronic structure of the cofactor. Together with isotope labeling, these techniques have also been instrumental in deciphering the binding of the two substrate water molecules to the cluster. These results are briefly described in the context of the history of biological water oxidation with special emphasis on recent work using time resolved X-ray diffraction with free electron lasers. It is shown that these data are instrumental for developing a model of the biological water oxidation cycle.

Early Binding of Substrate Oxygen Is Responsible for a Spectroscopically Distinct S2 State in Photosystem II

The biological generation of oxygen by the oxygen-evolving complex (OEC) in photosystem II (PS II) is one of nature's most important reactions. The OEC is a Mn₄Ca cluster that has multiple Mn-O-Mn and Mn-O-Ca bridges and binds four water molecules. Previously, binding of an additional oxygen was detected in the S₂ to S₃ transition. Here we demonstrate that early binding of the substrate oxygen to the five-coordinate Mn₁ center in the S₂ state is likely responsible for the S₂ high-spin EPR signal. Substrate binding in the Mn₁-OH form explains the prevalence of the high-spin S₂ state at higher pH and its low-temperature conversion into the S₃ state. The given interpretation was confirmed by X-ray absorption spectroscopic measurements, DFT, and broken symmetry DFT calculations of structures and magnetic properties. Structural, electronic, and spectroscopic properties of the high-spin S₂ state model are provided and compared with the available S₃ state models. New interpretation of the high-spin S₂ state opens opportunity for analysis of factors controlling the oxygen substrate binding in PS II.

Implications of structural heterogeneity for the electronic structure of the final oxygen-evolving intermediate in photosystem II

Heterogeneity in intermediate catalytic states of the oxygen-evolving complex (OEC) of Photosystem II is known from a wide range of experimental and theoretical data, but its potential implications for the mechanism of water oxidation remain unexplored. We delineate the consequences of structural heterogeneity for the final step of the catalytic cycle by tracing the evolution of three spectroscopically relevant and structurally distinct components of the last metastable S₃ state to the transient O₂-evolving S₄ state of the OEC. Using quantum chemical calculations, we show that each S₃ isomer leads to a different electronic structure formulation for the active S₄ state. Crucially, in addition to previously hypothesized Mn(IV)-oxyl species, we establish for the first time, how a genuine Mn(V)-oxo can be obtained in the catalytically active S₄ state: this takes the form of a five-coordinate and locally high-spin (S_Mn = 1) Mn(V) site. This formulation for the S₄ state evolves naturally from a preceding S₃-state structural intermediate that contains a quasi-trigonal-bipyramidal Mn(IV) ion. The results strongly suggest that water binding in the S₃ state is not prerequisite for reaching the oxygen-evolving S₄ state of the complex, supporting the notion that both substrates are preloaded at the beginning of the catalytic cycle. This scenario allows true four-electron metal-centered hole accumulation to precede OO bond formation and hence the latter can proceed via a genuine even-electron mechanism. This can occur as intramolecular nucleophilic coupling of two oxo units synchronously with the binding of a water substrate for the next catalytic cycle.

On the comparison between differential vibrational spectroscopy spectra and theoretical data in the carboxyl region of photosystem II

Understanding the structural modification experienced by the Mn₄ CaO₅ oxygen-evolving complex of photosystem II along the Kok-Joliot's cycle has been a challenge for both theory and experiments since many decades. In particular, differential infrared spectroscopy was extensively used to probe the surroundings of the reaction center, to catch spectral changes between different S-states along the catalytic cycle. Because of the complexity of the signals, only a limited quantity of identified peaks have been assigned so far, also because of the difficulty of a direct comparison with theoretical calculations. In the present work, we critically reconsider the comparison between differential vibrational spectroscopy and theoretical calculations performed on the structural models of the photosystem II active site and an inorganic structural mimic. Several factors are currently limiting the reliability of a quantitative comparison, such as intrinsic errors associated to theoretical methods, and most of all, the uncertainty attributed to the lack of knowledge about the localization of the underlying structural changes. Critical points in this comparison are extensively discussed. Comparing several computational data of differential S₂ /S₁ infrared spectroscopy, we have identified weak and strong points in their interpretation when compared with experimental spectra.

Proton Translocation via Tautomerization of Asn298 During the S2-S3 State Transition in the Oxygen-Evolving Complex of Photosystem II

In biological water oxidation, a redox-active tyrosine residue (D1-Tyr161 or Y_Z) mediates electron transfer between the Mn₄CaO₅ cluster of the oxygen-evolving complex and the charge-separation site of photosystem II (PSII), driving the cluster through progressively higher oxidation states S_i (i = 0–4). In contrast to lower S-states (S₀, S₁), in higher S-states (S₂, S₃) of the Mn₄CaO₅ cluster, Y_Z cannot be oxidized at cryogenic temperatures due to the accumulation of positive charge in the S₁ → S₂ transition. However, oxidation of Y_Z by illumination of S₂ at 77–190 K followed by rapid freezing and charge recombination between Y_Z^• and the plastoquinone radical Q_A^•– allows trapping of an S₂ variant, the so-called S₂^trapped state (S₂^t), that is capable of forming Y_Z^• at cryogenic temperature. To identify the differences between the S₂ and S₂^t states, we used the S₂^tY_Z^• intermediate as a probe for the S₂^t state and followed the S₂^tY_Z^•/Q_A^•– recombination kinetics at 10 K using time-resolved electron paramagnetic resonance spectroscopy in H₂O and D₂O. The results show that while S₂^tY_Z^•/Q_A^•– recombination can be described as pure electron transfer occurring in the Marcus inverted region, the S₂^t → S₂ reversion depends on proton rearrangement and exhibits a strong kinetic isotope effect. This suggests that Y_Z oxidation in the S₂^t state is facilitated by favorable proton redistribution in the vicinity of Y_Z, most likely within the hydrogen-bonded Y_Z–His190–Asn298 triad. Computational models show that tautomerization of Asn298 to its imidic acid form enables proton translocation to an adjacent asparagine-rich cavity of water molecules that functions as a proton reservoir and can further participate in proton egress to the lumen.

Spin, Valence, and Structural Isomerism in the S3 State of the Oxygen-Evolving Complex of Photosystem II as a Manifestation of Multimetallic Cooperativity

Photosynthetic water oxidation is catalyzed by a Mn₄CaO₅-cluster in photosystem II through an S-state cycle. Understanding the roles of heterogeneity in each S-state, as identified recently by the EPR spectroscopy, is very important to gain a complete description of the catalytic mechanism. We performed herein hybrid DFT calculations within the broken-symmetry formalism and associated analyses of Heisenberg spin models to study the electronic and spin structures of various isomeric structural motifs (hydroxo-oxo, oxyl-oxo, peroxo, and superoxo species) in the S₃ state. Our extensive study reveals several factors that affect the spin ground state: (1) (formal) Mn oxidation state; (2) metal-ligand covalency; (3) coordination geometry; and (4) structural change of the Mn cluster induced by alternations in Mn···Mn distances. Some combination of these effects could selectively stabilize/destabilize some spin states. We found that the high spin state ( S_total = 6) of the oxyl-oxo species can be causative for catalytic function, which manifests through mixing of the metal-ligand character in magnetic orbitals at relatively short O5···O6 distances (<2.0 Å) and long Mn_A···O5 distances (>2.0 Å). These results will serve as a basis to conceptually identify and rationalize the physicochemical synergisms that can be evoked by the unique "distorted chair" topology of the cluster through cooperative Jahn-Teller effects on multimetallic centers.

Tracing the Pathways of Waters and Protons in Photosystem II and Cytochrome c Oxidase

Photosystem II (PSII) uses water as the terminal electron donor, producing oxygen in the Mn4CaO5 oxygen evolving complex (OEC), while cytochrome c oxidase (CcO) reduces O2 to water in its heme–Cu binuclear center (BNC). Each protein is oriented in the membrane to add to the proton gradient. The OEC, which releases protons, is located near the P-side (positive, at low-pH) of the membrane. In contrast, the BNC is in the middle of CcO, so the protons needed for O2 reduction must be transferred from the N-side (negative, at high pH). In addition, CcO pumps protons from N- to P-side, coupled to the O2 reduction chemistry, to store additional energy. Thus, proton transfers are directly coupled to the OEC and BNC redox chemistry, as well as needed for CcO proton pumping. The simulations that study the changes in proton affinity of the redox active sites and the surrounding protein at different states of the reaction cycle, as well as the changes in hydration that modulate proton transfer paths, are described.

Thermodynamics of the S2-to-S3 state transition of the oxygen-evolving complex of photosystem II

The S₂ to S₃ transition in the OEC of PSII changes the structure of the Mn cluster. Monte Carlo sampling finds a Ca terminal water moves to form a bridge to Mn4 and the Mn1 ligand E189 can be replaced with a hydroxyl as a proton is lost.

Manganese-Hydroxido Complexes Supported by a Urea/Phosphinic Amide Tripodal Ligand

Hydrogen bonds (H-bonds) within the secondary coordination sphere are often invoked as essential noncovalent interactions that lead to productive chemistry in metalloproteins. Incorporating these types of effects within synthetic systems has proven a challenge in molecular design that often requires the use of rigid organic scaffolds to support H-bond donors or acceptors. We describe the preparation and characterization of a new hybrid tripodal ligand ([H₂pout]^3-) that contains two monodeprotonated urea groups and one phosphinic amide. The urea groups serve as H-bond donors, while the phosphinic amide group serves as a single H-bond acceptor. The [H₂pout]^3- ligand was utilized to stabilize a series of Mn-hydroxido complexes in which the oxidation state of the metal center ranges from 2+ to 4+. The molecular structure of the Mn^III-OH complex demonstrates that three intramolecular H-bonds involving the hydroxido ligand are formed. Additional evidence for the formation of intramolecular H-bonds was provided by vibrational spectroscopy in which the energy of the O-H vibration supports its assignment as an H-bond donor. The stepwise oxidation of [Mn^IIH₂pout(OH)]^2- to its higher oxidized analogs was further substantiated by electrochemical measurements and results from electronic absorbance and electron paramagnetic resonance spectroscopies. Our findings illustrate the utility of controlling both the primary and secondary coordination spheres to achieve structurally similar Mn-OH complexes with varying oxidation states.

Explaining the Different Geometries of the Water Oxidising Complex in the Nominal S3 State Crystal Structures of Photosystem II at 2.25 Å and 2.35 Å

Recently two atomic resolution crystal structures of Photosystem II, in the double flashed, nominal S3 intermediate state of its Mn4 Ca water oxidising complex (WOC), have been presented (Young et al., Nature 2016, 540, 453; Suga et al., Nature 2017, 543, 131). These structures are at 2.25 Å and 2.35 Å resolution, respectively. Although highly similar in most respects, the structures differ in a key region within the WOC catalytic site. In the 2.25 Å structure, one oxy species (O5) is observed within the WOC cavity, weakly associated with the Mn centres, similar to that seen earlier in the 1.95 Å XRD structure of the S1 intermediate (Suga et al., Nature, 2015, 517, 99). In the 2.35 Å structure, two such species are seen (O5, O6), with the Mn centres and O5 positioned as in the 2.25 Å structure and an O5-O6 separation of ∼1.5 Å, consistent with peroxo formation. This suggests O5 and O6 are substrate water derived species in this double flashed form. Recently we have presented (Petrie, et al., Chem. Phys. Chem., 2017) a large scale (220 atom) quantum chemical model of the Young et al. 2.25 Å structure, which quantitatively explains all significant features within the WOC region of that structure, particularly the positions of the metal centres and O5 group. Critical to this was our assumption of a 'low' Mn oxidation paradigm (mean S1 Mn oxidation level of +3.0, Petrie et al., Angew. Chem. Int. Ed., 2015), rather than a 'high' oxidation model (mean S1 oxidation level of +3.5), widely assumed in the literature. Here we show that our same oxidation state model predicts two classes of energetically close S3 structural forms, analogous to the S1 state, one with the metal centres and O5 positioned as in the 2.25 Å structure, and the other with the metals similarly placed, but with O5 located in the O6 position of the 2.35 Å structure. We show that the Suga et al. 2.35 Å structure is likely a superposition of two such forms, one from each class, which is consistent with reported atomic occupancies for that structure and the relative total energies we calculate for the two structural forms.

Thermodynamics of Proton and Electron Transfer in Tetranuclear Clusters with Mn-OH2/OH Motifs Relevant to H2O Activation by the Oxygen Evolving Complex in Photosystem II

We report the synthesis of site-differentiated heterometallic clusters with three Fe centers and a single Mn site that binds water and hydroxide in multiple cluster oxidation states. Deprotonation of Fe^III/II₃Mn^II-OH₂ clusters leads to internal reorganization resulting in formal oxidation at Mn to generate Fe^III/II₃Mn^III-OH. ⁵⁷Fe Mössbauer spectroscopy reveals that oxidation state changes (three for Fe^III/II₃Mn-OH₂ and four for Fe^III/II₃Mn-OH clusters) occur exclusively at the Fe centers; the Mn center is formally Mn^II when water is bound and Mn^III when hydroxide is bound. Experimentally determined p K_a (17.4) of the [Fe^III₂Fe^IIMn^II-OH₂] cluster and the reduction potentials of the [Fe₃Mn-OH₂] and [Fe₃Mn-OH] clusters were used to analyze the O-H bond dissociation enthalpies (BDE_O-H) for multiple cluster oxidation states. BDE_O-H increases from 69 to 78 and 85 kcal/mol for the [Fe^IIIFe^II₂Mn^II-OH₂], [Fe^III₂Fe^IIMn^II-OH₂], and [Fe^III₃Mn^II-OH₂] clusters, respectively. Further insight of the proton and electron transfer thermodynamics of the [Fe₃Mn-OH _x] system was obtained by constructing a potential-p K_a diagram; the shift in reduction potentials of the [Fe₃Mn-OH _x] clusters in the presence of different bases supports the BDE_O-H values reported for the [Fe₃Mn-OH₂] clusters. A lower limit of the p K_a for the hydroxide ligand of the [Fe₃Mn-OH] clusters was estimated for two oxidation states. These data suggest BDE_O-H values for the [Fe^III₂Fe^IIMn^III-OH] and [Fe^III₃Mn^III-OH] clusters are greater than 93 and 103 kcal/mol, which hints to the high reactivity expected of the resulting [Fe₃Mn═O] in this and related multinuclear systems.

Fourier Transform Infrared Analysis of the S-State Cycle of Water Oxidation in the Microcrystals of Photosystem II

Photosynthetic water oxidation is performed in photosystem II (PSII) through a light-driven cycle of intermediates called S states (S₀-S₄) at the water oxidizing center. Time-resolved serial femtosecond crystallography (SFX) has recently been applied to the microcrystals of PSII to obtain the structural information on these intermediates. However, it remains unanswered whether the reactions efficiently proceed throughout the S-state cycle retaining the native structures of the intermediates in PSII crystals. We investigated the water oxidation reactions in the PSII microcrystals using flash-induced Fourier transform infrared (FTIR) difference spectroscopy. In comparison with the FTIR spectra in solution, it was shown that all of the metastable intermediates in the microcrystals retained their native structures, and the efficiencies of the S-state transitions remained relatively high, although those of the S₂ → S₃ and S₃ → S₀ transitions were slightly lowered possibly due to some restriction of water movement in the crystals.

Computational studies of the Mn4/Ca cluster in photosystem II

Computational chemistry with the data from more detailed X-ray diffraction (XRD) oxygen evolving complex (OEC) structure has been used extensively of late in exploring the mechanisms of water oxidation in the OEC. The study reported in this paper involves density functional theory (DFT) calculations to investigate whether the data are in agreement with the four manganese ions in the OEC, being organized as a ‘3[Formula: see text]1’ (trimer plus one) model [Gatt et al. Angewandte Chemie International Edition, 51, 12025–12028, 2012; Petrie et al. Chemistry - A European Journal, 21, 6780–6792, 2015; Terrett et al. Chemical Communications, 50, 3187–3190, 2014] or ‘dimer of dimers’ model. [Terrett et al. Journal of Inorganic Biochemistry, 162, 178–189, 2016]. The data analysis method used involves quantum chemical DFT calculations on relevant models of the OEC cluster. DFT calculations were performed using both the so-called ‘open’ and ‘closed’ forms [Terrett et al. Journal of Inorganic Biochemistry, 162, 178–189, 2016] of the S2 OEC structure models with total spin ([Formula: see text]) 1/2, 7/2, 9/2 and 15/2 within the MnIII MnIV MnIII MnIII ‘low’ oxidation paradigm to examine exchange coupling within the OEC cluster. The results show that the [Formula: see text]-coupling in the ‘closed’ form: [Formula: see text][Formula: see text]cm[Formula: see text], [Formula: see text][Formula: see text]cm[Formula: see text], [Formula: see text][Formula: see text]cm[Formula: see text] and [Formula: see text]–[Formula: see text][Formula: see text]cm[Formula: see text]. In the ‘closed’ form, [Formula: see text] and [Formula: see text] represent the two largest exchange interactions within the manganese cluster, whereas [Formula: see text] and [Formula: see text] are small and almost net cancel. The magnetic coupling between the four Mn ions is close to ‘dimer of dimers’, with both dimers anti-ferromagnetically coupled internally and with weak inter-dimer net coupling.

The S2 to S3 transition for water oxidation in PSII (photosystem II), revisited

The formation of O₂ from water requires four transitions, each one after the absorption of one light flash.

D1-Asn-298 in photosystem II is involved in a hydrogen-bond network near the redox-active tyrosine YZ for proton exit during water oxidation

In photosynthetic water oxidation, two water molecules are converted into one oxygen molecule and four protons at the Mn₄CaO₅ cluster in photosystem II (PSII) via the S-state cycle. Efficient proton exit from the catalytic site to the lumen is essential for this process. However, the exit pathways of individual protons through the PSII proteins remain to be identified. In this study, we examined the involvement of a hydrogen-bond network near the redox-active tyrosine Y_Z in proton transfer during the S-state cycle. We focused on spectroscopic analyses of a site-directed variant of D1-Asn-298, a residue involved in a hydrogen-bond network near Y_Z We found that the D1-N298A mutant of Synechocystis sp. PCC 6803 exhibits an O₂ evolution activity of ∼10% of the wild-type. D1-N298A and the wild-type D1 had very similar features of thermoluminescence glow curves and of an FTIR difference spectrum upon Y_Z oxidation, suggesting that the hydrogen-bonded structure of Y_Z and electron transfer from the Mn₄CaO₅ cluster to Y_Z were little affected by substitution. In the D1-N298A mutant, however, the flash-number dependence of delayed luminescence showed a monotonic increase without oscillation, and FTIR difference spectra of the S-state cycle indicated partial and significant inhibition of the S₂ → S₃ and S₃ → S₀ transitions, respectively. These results suggest that the D1-N298A substitution inhibits the proton transfer processes in the S₂ → S₃ and S₃ → S₀ transitions. This in turn indicates that the hydrogen-bond network near Y_Z can be functional as a proton transfer pathway during photosynthetic water oxidation.

Tracking Reactive Water and Hydrogen-Bonding Networks in Photosynthetic Oxygen Evolution

In oxygenic photosynthesis, photosystem II (PSII) converts water to molecular oxygen through four photodriven oxidation events at a Mn₄CaO₅ cluster. A tyrosine, YZ (Y161 in the D1 polypeptide), transfers oxidizing equivalents from an oxidized, primary chlorophyll donor to the metal center. Calcium or its analogue, strontium, is required for activity. The Mn₄CaO₅ cluster and YZ are predicted to be hydrogen bonded in a water-containing network, which involves amide carbonyl groups, amino acid side chains, and water. This hydrogen-bonded network includes amino acid residues in intrinsic and extrinsic subunits. One of the extrinsic subunits, PsbO, is intrinsically disordered. This extensive (35 Å) network may be essential in facilitating proton release from substrate water. While it is known that some proteins employ internal water molecules to catalyze reactions, there are relatively few methods that can be used to study the role of water. In this Account, we review spectroscopic evidence from our group supporting the conclusion that the PSII hydrogen-bonding network is dynamic and that water in the network plays a direct role in catalysis. Two approaches, transient electron paramagnetic resonance (EPR) and reaction-induced FT-IR (RIFT-IR) spectroscopies, were used. The EPR experiments focused on the decay kinetics of YZ• via recombination at 190 K and the solvent isotope, pH, and calcium dependence of these kinetics. The RIFT-IR experiments focused on shifts in amide carbonyl frequencies, induced by photo-oxidation of the metal cluster, and on the isotope-based assignment of bands to internal, small protonated water clusters at 190, 263, and 283 K. To conduct these experiments, PSII was prepared in selected steps along the catalytic pathway, the S_n state cycle (n = 0-4). This cycle ultimately generates oxygen. In the EPR studies, S-state dependent changes were observed in the YZ• lifetime and in its solvent isotope effect. The YZ• lifetime depended on the presence of calcium at pH 7.5, but not at pH 6.0, suggesting a two-donor model for PCET. At pH 6.0 or 7.5, barium and ammonia both slowed the rate of YZ• recombination, consistent with disruption of the hydrogen-bonding network. In the RIFT-IR studies of the S state transitions, infrared bands associated with the transient protonation and deprotonation of internal waters were identified by D₂O and H₂¹⁸O labeling. The infrared bands of these protonated water clusters, W_n⁺ (or nH₂O(H₃O)⁺, n = 5-6), exhibited flash dependence and were produced during the S₁ to S₂ and S₃ to S₀ transitions. Calcium dependence was observed at pH 7.5, but not at pH 6.0. S-state induced shifts were observed in amide C═O frequencies during the S₁ to S₂ transition and attributed to alterations in hydrogen bonding, based on ammonia sensitivity. In addition, isotope editing of the extrinsic subunit, PsbO, established that amide vibrational bands of this lumenal subunit respond to the S state transitions and that PsbO is a structural template for the reaction center. Taken together, these spectroscopic results support the hypothesis that proton transfer networks, extending from YZ to PsbO, play a functional and dynamic role in photosynthetic oxygen evolution.