Characterization of the Sr(2+)- and Cd(2+)-Substituted Oxygen-Evolving Complex of Photosystem II by Quantum Mechanics/Molecular Mechanics Calculations

Insights from Ca2+→Sr2+ substitution on the mechanism of O-O bond formation in photosystem II

In recent years, there has been a steady interest in unraveling the intricate mechanistic details of water oxidation mechanism in photosynthesis. Despite the substantial progress made over several decades, a comprehensive understanding of the precise kinetics underlying O-O bond formation and subsequent evolution remains elusive. However, it is well-established that the oxygen evolving complex (OEC), specifically the CaMn4O5 cluster, plays a crucial role in O-O bond formation, undergoing a series of four oxidative events as it progresses through the S-states of the Kok cycle. To gain further insights into the OEC, researchers have explored the substitution of the Ca2+ cofactor with strontium (Sr), the sole atomic replacement capable of retaining oxygen-evolving activity. Empirical investigations utilizing spectroscopic techniques such as XAS, XRD, EPR, FTIR, and XANES have been conducted to probe the structural consequences of Ca2+→Sr2+ substitution. In parallel, the development of DFT and QM/MM computational models has explored different oxidation and protonation states, as well as variations in ligand coordination at the catalytic center involving amino acid residues. In this review, we critically evaluate and integrate these computational and spectroscopic approaches, focusing on the structural and mechanistic implications of Ca2+→Sr2+ substitution in PS II. We contribute DFT modelling and simulate EXAFS Fourier transforms of Sr-substituted OEC, analyzing promising structures of the S3 state. Through the combination of computational modeling and spectroscopic investigations, valuable insights have been gained, developing a deeper understanding of the photosynthetic process.

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.

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.

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.

Substrate water exchange in the S2 state of photosystem II is dependent on the conformation of the Mn4Ca cluster

The structural flexibility of the Mn₄Ca cluster in photosystem II supports the exchange of the central O5 bridge.

Five-coordinate MnIV intermediate in the activation of nature's water splitting cofactor

Nature's water splitting cofactor passes through a series of catalytic intermediates (S₀-S₄) before O-O bond formation and O₂ release. In the second last transition (S₂ to S₃) cofactor oxidation is coupled to water molecule binding to Mn1. It is this activated, water-enriched all Mn^IV form of the cofactor that goes on to form the O-O bond, after the next light-induced oxidation to S₄ How cofactor activation proceeds remains an open question. Here, we report a so far not described intermediate (S₃') in which cofactor oxidation has occurred without water insertion. This intermediate can be trapped in a significant fraction of centers (>50%) in (i) chemical-modified cofactors in which Ca²⁺ is exchanged with Sr²⁺; the Mn₄O₅Sr cofactor remains active, but the S₂-S₃ and S₃-S₀ transitions are slower than for the Mn₄O₅Ca cofactor; and (ii) upon addition of 3% vol/vol methanol; methanol is thought to act as a substrate water analog. The S₃' electron paramagnetic resonance (EPR) signal is significantly broader than the untreated S₃ signal (2.5 T vs. 1.5 T), indicating the cofactor still contains a 5-coordinate Mn ion, as seen in the preceding S₂ state. Magnetic double resonance data extend these findings revealing the electronic connectivity of the S₃' cofactor is similar to the high spin form of the preceding S₂ state, which contains a cuboidal Mn₃O₄Ca unit tethered to an external, 5-coordinate Mn ion (Mn4). These results demonstrate that cofactor oxidation regulates water molecule insertion via binding to Mn4. The interaction of ammonia with the cofactor is also discussed.

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.

A Comparison of Experimental and Broken Symmetry Density Functional Theory (BS-DFT) Calculated Electron Paramagnetic Resonance (EPR) Parameters for Intermediates Involved in the S2 to S3 State Transition of Nature's Oxygen Evolving Complex

A broken symmetry density functional theory (BS-DFT) magnetic analysis of the S₂, S₂Y_Z^•, and S₃ states of Nature's oxygen evolving complex is performed for both the native Ca and Sr substituted forms. Good agreement with experiment is observed between the tyrosyl calculated g-tensor and ¹H hyperfine couplings for the native Ca form. Changes in the hydrogen bonding environment of the tyrosyl radical in S₂Y_Z^• caused by Sr substitution lead to notable changes in the calculated g-tensor of the tyrosyl radical. Comparison of calculated and experimental ⁵⁵Mn hyperfine couplings for the S₃ state presently favors an open cubane form of the complex with an additional OH ligand coordinating to Mn_D. In Ca models, this additional ligation can arise by closed-cubane form deprotonation of the Ca ligand W3 in the S₂Y_Z^• state accompanied by spontaneous movement to the vacant Mn coordination site or by addition of an external OH group. For the Sr form, no spontaneous movement of W3 to the vacant Mn coordination site is observed in contrast to the native Ca form, a difference which may lead to the reduced catalytic activity of the Sr substituted form. BS-DFT studies on peroxo models of S₃ as indicated by a recent X-ray free electron laser (XFEL) crystallography study give rise to a structural model compatible with experimental data and an S = 3 ground state compatible with EPR studies.

Quantum mechanics/molecular mechanics simulation of the ligand vibrations of the water-oxidizing Mn4CaO5 cluster in photosystem II

During photosynthesis, the light-driven oxidation of water performed by photosystem II (PSII) provides electrons necessary to fix CO2, in turn supporting life on Earth by liberating molecular oxygen. Recent high-resolution X-ray images of PSII show that the water-oxidizing center (WOC) is composed of an Mn4CaO5 cluster with six carboxylate, one imidazole, and four water ligands. FTIR difference spectroscopy has shown significant structural changes of the WOC during the S-state cycle of water oxidation, especially within carboxylate groups. However, the roles that these carboxylate groups play in water oxidation as well as how they should be properly assigned in spectra are unresolved. In this study, we performed a normal mode analysis of the WOC using the quantum mechanics/molecular mechanics (QM/MM) method to simulate FTIR difference spectra on the S1 to S2 transition in the carboxylate stretching region. By evaluating WOC models with different oxidation and protonation states, we determined that models of high-oxidation states, Mn(III)2Mn(IV)2, satisfactorily reproduced experimental spectra from intact and Ca-depleted PSII compared with low-oxidation models. It is further suggested that the carboxylate groups bridging Ca and Mn ions within this center tune the reactivity of water ligands bound to Ca by shifting charge via their π conjugation.

Mechanism of Water Delivery to the Active Site of Photosystem II along the S(2) to S(3) Transition

The two water molecules serving as substrate for the oxygen evolution in Photosystem II are already bound in the S2 state of the Kok-Joliot's cycle. Nevertheless, an additional water molecule is supposed to bind the cluster during the transition between the S2 and S3 states, which has been recently revealed to have the Mn4CaO5 catalytic cluster arranged in an open cubane fashion. In this Letter, by means of ab initio calculations, we investigated the possible pathways for the binding of the upcoming water molecule. Upon the four different possibilities checked in our calculations, the binding of the crystallographic water molecule, originally located nearby the Cl(-) binding site, showed the lowest activation energy barrier. Our findings therefore support the view in which the W2 hydroxyl group and the O5 oxygen act as substrates for the oxygen evolution. Within this framework the role of the open and closed Mn4CaO5 conformers is clarified as well as the exact mechanistic events occurring along the S2 to S3 transition.

Redox-coupled substrate water reorganization in the active site of Photosystem II-The role of calcium in substrate water delivery

Photosystem II (PSII) catalyzes light-driven water splitting in nature and is the key enzyme for energy input into the biosphere. Important details of its mechanism are not well understood. In order to understand the mechanism of water splitting, we perform here large-scale density functional theory (DFT) calculations on the active site of PSII in different oxidation, spin and ligand states. Prior to formation of the O-O bond, we find that all manganese atoms are oxidized to Mn(IV) in the S3 state, consistent with earlier studies. We find here, however, that the formation of the S3 state is coupled to the movement of a calcium-bound hydroxide (W3) from the Ca to a Mn (Mn1 or Mn4) in a process that is triggered by the formation of a tyrosyl radical (Tyr-161) and its protonated base, His-190. We find that subsequent oxidation and deprotonation of this hydroxide on Mn1 result in formation of an oxyl-radical that can exergonically couple with one of the oxo-bridges (O5), forming an O-O bond. When O(2) leaves the active site, a second Ca-bound water molecule reorients to bridge the gap between the manganese ions Mn1 and Mn4, forming a new oxo-bridge for the next reaction cycle. Our findings are consistent with experimental data, and suggest that the calcium ion may control substrate water access to the water oxidation sites.

Redox potential tuning by redox-inactive cations in nature's water oxidizing catalyst and synthetic analogues

Fundamental differences between synthetic manganese clusters and the biological water oxidizing catalyst are demonstrated in the modulation of their redox potential by redox-inactive cations.

Magnetic Behavior of Heterometallic Wheels Having a [Mn(IV)6M2O9](10+) Core with M = Ca(2+) and Sr(2+)

Two new heterometallic Mn(IV)-M(2+) compounds with formula [Mn6M2O9(4-(t)BuC6H4COO)10(4-(t)BuC6H4COOH)5] (M = Ca(2+) (1), Sr(2+) (2)) have been crystallized. The core of both compounds consists of a planar Mn6 ring, where the Mn(IV) ions are alternatively bridged by (μ3-O)2(μ-RCOO) and (μ4-O)(μ-RCOO)2 ligands, and the two alkaline earth ions are located to both sides of the wheel, linked to the oxo bridges, generating three fused [Mn2M2O4](4+) cuboids. These compounds show a net antiferromagnetic behavior, more important for 2 (Sr(2+)) than for 1 (Ca(2+)). The fitting of the experimental data was performed with the support of DFT calculations, considering four different exchange pathways: two between adjacent Mn(IV) ions (J1 and J2) and two between nonadjacent Mn(IV) ions (J3 and J4). The results of the analysis show that J1 and J2 are of the opposite sign, the ferromagnetic contribution corresponding to the [Mn2(μ4-O)(μ-RCOO)2](4+) unit (J2). The influence of the M(2+) ions in the magnetic behavior is analyzed for 1 and 2 and for three hypothetical models with the structural parameters of 1 containing Mg(2+), Sr(2+) or without the M(2+) ions. In spite of the diamagnetic character of the alkaline earth ions, their influence on the magnetic behavior has been evidenced and correlated with their polarizing effect. Moreover, the magnetic interactions between nonadjacent ions are non-negligible.