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

Intermediate Formation via Proton Release during the Photoassembly of the Water-Oxidizing Mn4CaO5 Cluster in Photosystem II

The early stages of the photoassembly of the water-oxidizing Mn4CaO5 cluster in spinach photosystem II (PSII) were monitored using rapid-scan time-resolved Fourier transform infrared (FTIR) spectroscopy. Carboxylate stretching and the amide I bands, which appeared upon the flash-induced oxidation of a Mn2+ ion, changed their features during the subsequent dark rearrangement process, indicating the relocation of the Mn3+ ion concomitant with protein conformational changes. Monitoring the isotope-edited FTIR signals of a Mes buffer estimated that nearly two protons are released upon the Mn2+ oxidation. Quantum chemical calculations for models of the Mn binding site suggested that the proton of a water ligand is transferred to D1-H332 through a hydrogen bond upon the Mn3+ formation and then released to the bulk as the Mn3+ shifts to bind to this histidine. Another Mn2+ ion may be inserted to form a binuclear Mn3+Mn2+ complex, whose structure was calculated to be stabilized by a μ-hydroxo bridge hydrogen-bonded with deprotonated D1-H337. Nearly one additional proton can thus be released from this histidine, assuming that it is mostly protonated before illumination. Alternatively, a proton could be released by further insertion of Ca2+, forming a Mn3+Mn2+Ca2+ complex with another hydroxo ligand connecting Ca2+ to the Mn3+Mn2+ complex.

Exploring the Deprotonation Process during Incorporation of a Ligand Water Molecule at the Dangling Mn Site in Photosystem II

The Mn4CaO5 cluster, featuring four ligand water molecules (W1 to W4), serves as the water-splitting site in photosystem II (PSII). X-ray free electron laser (XFEL) structures exhibit an additional oxygen site (O6) adjacent to the O5 site in the fourth lowest oxidation state, S3, forming Mn4CaO6. Here, we investigate the mechanism of the second water ligand molecule at the dangling Mn (W2) as a potential incorporating species, using a quantum mechanical/molecular mechanical (QM/MM) approach. Previous QM/MM calculations demonstrated that W1 releases two protons through a low-barrier H-bond toward D1-Asp61 and subsequently releases an electron during the S2 to S3 transition, resulting in O•- at W1 and protonated D1-Asp61. During the process of Mn4CaO6 formation, O•-, rather than H2O or OH-, best reproduced the O5···O6 distance. Although the catalytic cluster with O•- at O6 is more stable than that with O•- at W1 in S3, it does not occur spontaneously due to the significantly uphill deprotonation process. Assuming O•- at W2 incorporates into the O6 site, an exergonic conversion from Mn1(III)Mn2(IV)Mn3(IV)Mn4(IV) (equivalent to the open-cubane S2 valence state) to Mn1(IV)Mn2(IV)Mn3(IV)Mn4(III) (equivalent to the closed-cubane S2 valence state) occurs. These findings provide energetic insights into the deprotonation and structural conversion events required for the Mn4CaO6 formation.

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-).

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.

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 Mn₄CaO₅ cluster (S₁) leads to the S₂ state formation, Mn1(III)Mn2(IV)Mn3(IV)Mn4(IV) (open-cubane S₂) or Mn1(IV)Mn2(IV)Mn3(IV)Mn4(III) (closed-cubane S₂). 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 S₂ 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 [Mn₃(CaO₄)] 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 S₂ is less stable than the open-cubane S₂ in cyanobacterial PSII, which may explain why the g = 4.1 EPR signal is absent in cyanobacterial PSII.

Binding and functions of the two chloride ions in the oxygen-evolving center of photosystem II

Light-driven water oxidation in photosynthesis occurs at the oxygen-evolving center (OEC) of photosystem II (PSII). Chloride ions (Cl^-) are essential for oxygen evolution by PSII, and two Cl^- ions have been found to specifically bind near the Mn₄CaO₅ cluster in the OEC. The retention of these Cl^- ions within the OEC is critically supported by some of the membrane-extrinsic subunits of PSII. The functions of these two Cl^- ions and the mechanisms of their retention both remain to be fully elucidated. However, intensive studies performed recently have advanced our understanding of the functions of these Cl^- ions, and PSII structures from various species have been reported, aiding the interpretation of previous findings regarding Cl^- retention by extrinsic subunits. In this review, we summarize the findings to date on the roles of the two Cl^- ions bound within the OEC. Additionally, together with a short summary of the functions of PSII membrane-extrinsic subunits, we discuss the mechanisms of Cl^- retention by these extrinsic subunits.

Release of a Proton and Formation of a Low-Barrier Hydrogen Bond between Tyrosine D and D2-His189 in Photosystem II

In photosystem II (PSII), the second-lowest oxidation state (S₁) of the oxygen-evolving Mn₄CaO₅ cluster is the most stable, as the radical form of the redox-active D2-Tyr160 is considered to be a candidate that accepts an electron from the lowest oxidation state (S₀) in the dark. Using quantum mechanical/molecular mechanical calculations, we investigated the redox potential (E _m) of TyrD and its H-bond partner, D2-His189. The potential energy profile indicates that the release of a proton from the TyrD...D2-His189 pair leads to the formation of a low-barrier H-bond. The E _m depends on the H⁺ position along the low-barrier H-bond, e.g., 680 mV when the H⁺ is at the D2-His189 moiety and 800 mV when the H⁺ is at the TyrD moiety, which can explain why TyrD mediates both the S₀ to S₁ oxidation and the S₂ to S₁ reduction.

Glycerol binding at the narrow channel of photosystem II stabilizes the low-spin S2 state of the oxygen-evolving complex

The oxygen-evolving complex (OEC) of photosystem II (PSII) cycles through redox intermediate states S_i (i = 0-4) during the photochemical oxidation of water. The S₂ state involves an equilibrium of two isomers including the low-spin S₂ (LS-S₂) state with its characteristic electron paramagnetic resonance (EPR) multiline signal centered at g = 2.0, and a high-spin S₂ (HS-S₂) state with its g = 4.1 EPR signal. The relative intensities of the two EPR signals change under experimental conditions that shift the HS-S₂/LS-S₂ state equilibrium. Here, we analyze the effect of glycerol on the relative stability of the LS-S₂ and HS-S₂ states when bound at the narrow channel of PSII, as reported in an X-ray crystal structure of cyanobacterial PSII. Our quantum mechanics/molecular mechanics (QM/MM) hybrid models of cyanobacterial PSII show that the glycerol molecule perturbs the hydrogen-bond network in the narrow channel, increasing the pK_a of D1-Asp61 and stabilizing the LS-S₂ state relative to the HS-S₂ state. The reported results are consistent with the absence of the HS-S₂ state EPR signal in native cyanobacterial PSII EPR spectra and suggest that the narrow water channel hydrogen-bond network regulates the relative stability of OEC catalytic intermediates during water oxidation.

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²⁺.

Experimental and Computational Insight into the Mechanism of NO Binding to Ferric Microperoxidase. The Likely Role of Tautomerization to Account for the pH Dependence

According to the current paradigm, the metal-hydroxo bond in a six-coordinate porphyrin complex is believed to be significantly less reactive in ligand substitution than the analogous metal-aqua bond, due to a much higher strength of the former bond. Here, we report kinetic studies for nitric oxide (NO) binding to a heme-protein model, acetylated microperoxidase-11 (AcMP-11), that challenge this paradigm. In the studied pH range 7.4-12.6, ferric AcMP-11 exists in three acid-base forms, assigned in the literature as [(AcMP-11)FeIII(H2O)(HisH)] (1), [(AcMP-11)FeIII(OH)(HisH)] (2), and [(AcMP-11)FeIII(OH)(His-)] (3). From the pH dependence of the second-order rate constant for NO binding (kon), we determined individual rate constants characterizing forms 1-3, revealing only a ca. 10-fold decrease in the NO binding rate on going from 1 (kon(1) = 3.8 × 106 M-1 s-1) to 2 (kon(2) = 4.0 × 105 M-1 s-1) and the inertness of 3. These findings lead to the abandonment of the dissociatively activated mechanism, in which the reaction rate can be directly correlated with the Fe-OH bond energy, as the mechanistic explanation for the process with regard to 2. The reactivity of 2 is accounted for through the existence of a tautomeric equilibrium between the major [(AcMP-11)FeIII(OH)(HisH)] (2a) and minor [(AcMP-11)FeIII(H2O)(His-)] (2b) species, of which the second one is assigned as the NO binding target due to its labile Fe-OH2 bond. The proposed mechanism is further substantiated by quantum-chemical calculations, which confirmed both the significant labilization of the Fe-OH2 bond in the [(AcMP-11)FeIII(H2O)(His-)] tautomer and the feasibility of the tautomer formation, especially after introducing empirical corrections to the computed relative acidities of the H2O and HisH ligands based on the experimental pKa values. It is shown that the "effective lability" of the axial ligand (OH-/H2O) in 2 may be comparable to the lability of the H2O ligand in 1.

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.

A Computational Study of the S2 State in the Oxygen-Evolving Complex of Photosystem II by Electron Paramagnetic Resonance Spectroscopy

The S2 state produces two basic electron paramagnetic resonance signal types due to the manganese cluster in oxygen-evolving complex, which are influenced by the solvents, and cryoprotectant added to the photosystem II samples. It is presumed that a single manganese center oxidation occurs on S1 → S2 state transition. The S2 state has readily visible multiline and g4.1 electron paramagnetic resonance signals and hence it has been the most studied of all the Kok cycle intermediates due to the ease of experimental preparation and stability. The S2 state was studied using electron paramagnetic resonance spectroscopy at X-band frequencies. The aim of this study was to determine the spin states of the g4.1 signal. The multiline signal was observed to arise from a ground state spin ½ centre while the g4.1 signal generated at ≈140 K NIR illumination was proposed to arise from a spin 52 center with rhombic distortion. The 'ground' state g4.1 signal was generated solely or by conversion from the multiline. The data analysis methods used involved numerical simulations of the experimental spectra on relevant models of the oxygen-evolving complex cluster. A strong focus in this paper was on the 'ground' state g4.1 signal, whether it is a rhombic 52 spin state signal or an axial 32 spin state signal. The data supported an X-band CW-EPR-generated g4.1 signal as originating from a near rhombic spin 5/2 of the S2 state of the PSII manganese cluster.

The Nature of the Short Oxygen-Oxygen Distance in the Mn4CaO6 Complex of Photosystem II Crystals

The O···O distance for a typical H-bond is ∼2.8 Å, whereas the radiation-damage-free structures of photosystem II (PSII), obtained using the X-ray free electron laser (XFEL), shows remarkably short O···O distances of ∼2 Å in the oxygen-evolving Mn₄CaO_5/6 complex. Herein, we report the protonation/oxidation states of the short O···O atoms in the XFEL structures using a quantum mechanical/molecular mechanical approach. The O5···O6 distance of 1.9 Å is reproduced only when O6 is an unprotonated O radical (O^•-) with Mn(IV)₃Mn(III), i.e., the S₃ state. The potential energy profile shows a barrier-less energy minimum region when O5···O6 = 1.90-2.05 Å (O^•- ↓) or 2.05-2.20 Å (O^•- ↑). Formation of such a short O5···O6 distance is not possible when O6 is OH^- with Mn(IV)₄. In the case in which the O5···O6 distance is 1.9 Å, it seems likely that the O radical species exists in the oxygen-evolving complex of the XFEL-S₃ crystals.

Proton transfer pathway from the oxygen-evolving complex in photosystem II substantiated by extensive mutagenesis

We report a structure-based biological approach to identify the proton-transfer pathway in photosystem II. First, molecular dynamics (MD) simulations were conducted to analyze the H-bond network that may serve as a Grotthuss-like proton conduit. MD simulations show that D1-Asp61, the H-bond acceptor of H₂O at the Mn₄CaO₅ cluster (W1), forms an H-bond via one water molecule with D1-Glu65 but not with D2-Glu312. Then, D1-Asp61, D1-Glu65, D2-Glu312, and the adjacent residues, D1-Arg334, D2-Glu302, and D2-Glu323, were thoroughly mutated to the other 19 residues, i.e., 114 Chlamydomonas chloroplast mutant cells were generated. Mutation of D1-Asp61 was most crucial. Only the D61E and D61C cells grew photoautotrophically and exhibit O₂-evolving activity. Mutations of D2-Glu312 were less crucial to photosynthetic growth than mutations of D1-Glu65. Quantum mechanical/molecular mechanical calculations indicated that in the PSII crystal structure, the proton is predominantly localized at D1-Glu65 along the H-bond with D2-Glu312, i.e., pK_a(D1-Glu65) > pK_a(D2-Glu312). The potential-energy profile shows that the release of the proton from D1-Glu65 leads to the formation of the two short H-bonds between D1-Asp61 and D1-Glu65, which facilitates downhill proton transfer along the Grotthuss-like proton conduit in the S₂ to S₃ transition. It seems possible that D1-Glu65 is involved in the dominant pathway that proceeds from W1 via D1-Asp61 toward the thylakoid lumen, whereas D2-Glu312 and D1-Arg334 may be involved in alternative pathways in some mutants.