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.

Introduction

Plants, algae, and cyanobacteria use the water-splitting enzyme photosystem II (PSII) for oxygen evolution. The oxygen evolution proceeds at the oxygen-evolving center, the Mn4CaO5 cluster. The cluster consists of a distorted cubane [Mn1, Mn2, Mn3, four oxygen atoms, and Ca2+] and “dangling” Mn4 (Fig. 1) (Umena et al. 2011). The Mn4CaO5 cluster has two ligand water molecules, W1 and W2, at the Mn4 site and another two ligand water molecules, W3 and W4, at the Ca2+ site (Fig. 1). The catalytic cycle moves through a series of oxidation states, denoted as S (n = 0, 1, 2, and 3). As electron transfer occurs, S increases. During the catalytic cycle, four electrons from two of the substrate water molecules are removed, and O2 evolves in the S3 to S0 transition (Shen 2015; Cardona and Rutherford 2019).

Fig. 1

Structure of the Mn4CaO5 cluster of PSII. Dotted lines indicate ligations to Ca2+

In the Mn4CaO5 cluster, a redox-inactive Ca2+ is essential for the oxygen evolution activity, as oxygen is not evolved when Ca2+ is removed (Ono and Inoue 1988) or replaced with Dy3+, Cu2+, Cd2+ (Lee et al. 2007), K+, Rb+, and Cs+ (Ono et al. 2001). The Mn4SrO5 cluster can evolve oxygen but the activity is lower than that of the native Mn4CaO5 cluster (Yachandra and Yano 2011). Koua et al. identified that the distance between Sr2+ and W3 (2.6 Å) was longer than that between Ca2+ and W3 (2.4 Å) (Koua et al. 2013) and proposed that the long Sr2+⋯W3 distance contributed to the decrease in the activity upon replacement of Ca2+ with Sr2+.

It was speculated that Ca2+ might be responsible for the distorted cubane structure of the Mn4CaO5 cluster (Kawakami et al. 2011). However, the removal of Ca2+ does not alter the Mn3CaO4 cubane structure (Saito and Ishikita 2014; Siegbahn 2014, 2017), as suggested by the extended X-ray absorption fine structure (EXAFS) and the electron paramagnetic resonance (EPR) measurements (Latimer et al. 1998; Yachandra and Yano 2011; Lohmiller et al. 2012). Note that the Jahn–Teller distortion for Mn(III) ions can be affected by Ca2+ (Yamaguchi et al. 2013). When Ca2+ is removed, the rearrangement of water molecules in the hydrogen-bond (H-bond) network of the redox-active tyrosine (TyrZ) is observed (Saito and Ishikita 2014; Saito et al. 2020a). TyrZ is involved in electron transfer from the Mn4CaO5 cluster to the reaction center chlorophyll PD1. The rearrangement of the H-bond network increases its redox potential (Em(TyrZ)) by ~ 300 mV and inhibits the formation of the downhill electron transfer pathway from the Mn4CaO5 cluster via TyrZ to PD1 (Saito et al. 2020a). Thus, Ca2+ is essential in both maintaining the H-bond network and optimizing electron transfer. The role of Ca2+ as the water binding site can be substituted with H3O+: recent theoretical studies showed that the H-bond network, including the low-barrier H-bond between TyrZ and D1-His190, remains unaltered upon the replacement of Ca2+ with H3O+ (Saito et al. 2020a).

Ca2+ is a prerequisite for the low-barrier H-bond between W1 and D1-Asp61: they form a low-barrier H-bond in native PSII (Kawashima et al. 2018b; Saito et al. 2020a), whereas they cannot form in the absence of Ca2+ (Saito et al. 2020a). That is, Ca2+ decreases pKa(W1) electrostatically to a level of pKa(D1-Asp61) in native PSII, thus forming the low-barrier H-bond and facilitating proton transfer from W1 to D1-Asp61.

So far, the role of the Ca2+ can be summarized as follows: (i) maintaining the TyrZ H-bond network (Saito et al. 2011, 2020a; Kawashima et al. 2018a), including the low-barrier H-bond between TyrZ and D1-His190 (Kawashima et al. 2018b; Saito et al. 2020a); (ii) optimizing Em(TyrZ) in the electron transfer cascade (Saito et al. 2020a); and (iii) electrostatically decreasing pKa(W1) and facilitating proton transfer via the low-barrier H-bond with D1-Asp61 (Saito et al. 2020a).

Althogh it was proposed that Ca2+ might electrostatically affect the properties of the cluster (e.g., pKa and Em) (McEvoy and Brudvig 2006), the replacements of Ca2+ with H2O and H3O+ lead to different Em(MnIII/IV) values due to different H-bond patterns in PSII (Saito et al. 2020a). Artificial Mn clusters with redox-inactive metals (Zhang et al. 2015; Mukherjee et al. 2012; Tsui et al. 2013; Kanady et al. 2013; Tsui and Agapie 2013; Lin et al. 2015) may serve as reference model systems since the corresponding H-bond network is absent. Tsui et al. synthesized artificial Mn3[M]O2 and clusters with redox-inactive metals [M] ([M] = Mg+, Ca2+, Zn2+, Sr2+, and Y3+) (Fig. 2), showing that Em(MnIII/IV) depends on the Lewis acidity of [M] (i.e., the pKa of aqua complexes of [M]) (Tsui et al. 2013). A similar correlation between the ligand-to-metal charge transfer energy (related to Em) and the Lewis acidity has also been reported in the Fe and Mn complexes (Bang et al. 2014; Krewald et al. 2016).

Fig. 2

Structure of the Mn3[M]O2 cluster for [M] = Ca2+, Sr2+, and Y3+ (Tsui et al. 2013) a Chemical structure. b Three-dimensional structure. For the detail of ligands, see Table S1

Em can be calculated as the free energy difference between the oxidized and reduced states, including the entropic effect of the solvent (Marenich et al. 2014; Pitari et al. 2015; Amin et al. 2013; Krewald et al. 2016). Em also correlates with the ionization potential as shown for various complexes, including Mn complexes (Marenich et al. 2014; Krewald et al. 2016). The ionization potential can be regarded as the free energy difference between the oxidized and the reduced states when the reorganization effect upon the redox reaction (including the electronic relaxation, the solvent reorganization, and the structural change of the molecule) is neglective. As the ionization potential (or the electronic affinity) is correlated with the energy levels of the lowest unoccupied molecular orbital (LUMO) and the highest occupied molecular orbital (HOMO) in density functional theory (DFT) (Kohn–Sham orbital) (Zhang and Musgrave 2007), Em should be calculated based on the HOMO or LUMO energy calculated using DFT. An electron releases from the HOMO upon oxidation, whereas an electron enters the LUMO upon reduction. Thus, the HOMO energy corresponds to the potential for one-electron oxidation, and the LUMO energy corresponds to the potential for one-electron reduction. When the redox reaction is reversible, the midpoint potential Em is located at the midpoint between the oxidation and reduction potentials, i.e., Em and the two potentials have the same tendency. Indeed, the Em of quinones can be determined based on the LUMO energy (Ishikita and Saito 2020) as accurately as the free energy difference (Kishi et al. 2017). For complexes that include transition metals, high correlations between the HOMO and/or LUMO energy and the experimentally measured Em value were observed [e.g., organic compounds (Mendez-Hernandez et al. 2013), β-diketones complexes, including Mn and Fe (Conradie 2015), and FeCo proteins (Dance 2006)]. For the natural Mn4CaO5 cluster in the PSII protein environment, Em can be determined based on the HOMO energy (Mandal et al. 2020; Saito et al. 2020a, b). Here we calculate the Em values of the artificial clusters (Tsui et al. 2013; Tsui and Agapie 2013) based on the HOMO energy and explain how the redox-inactive metal affects Em(MnIII/IV).

Computational details

The crystal structures of synthetic Mn3[M]O2 clusters ([M] = Na+, Sr2+, Ca2+, Zn2+, and Y3+) (Tsui et al. 2013) and Mn3[M]O4 clusters ([M] = Sr2+, Ca2+, Zn2+, Mn3+, Sc3+ and Y3+) (Tsui and Agapie 2013) were used as the basis for geometry optimization using unrestricted DFT (UDFT), with the B3LYP functional and LACVP* basis set (for optimized structures, see Supporting Information). For efficiency, the cluster was considered to comprise ferromagnetically coupled Mn atoms, (Tsui et al. 2013) where the total spin, S, was 12/2 for the Mn3[M]O2 cluster and 10/2 for the Mn3[M]O4 cluster. We note that, in the calculation of the native Mn4CaO5 of PSII, the difference in S (e.g., S = 1/2 in S2 (Zimmermann and Rutherford 1986), high, low, ferromagnetic, and antiferromagnetic) did not affect the (i) resulting geometry (Ames et al. 2011; Isobe et al. 2012), (ii) potential energy profile of proton transfer (Kawashima et al. 2018b), (iii) redox potential of each Mn site (Mandal et al. 2020), or (iv) pKa values of ligand water molecules W1–W4 (Saito et al. 2020c). The resulting oxidation states for three Mn atoms were Mn(III)3 and Mn(IV)Mn(III)2 for the Mn3[M]O2 and the Mn3[M]O4 clusters, respectively (for atomic spin density, see Table S2). Em(Mn1III/IV) was calculated from the HOMO energies, since the value of Em for one-electron oxidation is correlated with the energy of the highest occupied molecular orbital (HOMO) (Mendez-Hernandez et al. 2013; Igarashi and Seefeldt 2003; Mandal et al. 2020). Using the optimized geometries in vacuum, the HOMO energy (EHOMO) was calculated in dichloromethane (CH2Cl2, dielectric constant 8.93) using the polarizable continuum model (PCM). All calculations were performed with Jaguar program [Schrödinger, LLC, 2012, New York]. The initial-guess wavefunctions were obtained using the ligand field theory (Vacek et al. 1999) implemented in the Jaguar program.

Results and discussion

The calculated EHOMO values for the artificial Mn3[M]O2 clusters ([M] = Na+, Sr2+, Ca2+, Zn2+, and Y3+) show a correlation with the experimentally measured Em(MnIII/IV) values (Fig. 3) in CH2Cl2 and are best fitted to Eq. (1).where Fc/Fc+ denotes ferrocene electrode. A similar correlation is also observed in the Mn3[M]O4 cubane clusters ([M] = Sr2+, Ca2+, Zn2+, Mn3+, Sc3+, and Y3+) (Tsui and Agapie 2013) (Fig. S1). The coefficient of −0.302 in Eq. (1) is the conversion factor from MO energy to Em, which may be associated with the solvation effect (Schmidt am Busch and Knapp 2005), whereas the offset of −1.710 V is associated with a difference between the absolute electrode potential and the Fc/Fc+ electrode potential and liquid junction potential (Kishi et al. 2017). These factors depend on the size and the net charge of the QM system, the solvent, and the reference electrode (e.g., see the caption of Fig. S1). Thus, Eq. (1) is applicable only to similar molecular groups (e.g., artificial Mn3[M]O2 clusters).

Fig. 3

Experimentally measured Em(MnIII/IV) values in CH2Cl2 (Tsui and Agapie 2013) and calculated HOMO energy levels (EHOMO) of the Mn3[M]O2 clusters. The coefficient of determination (R2) is 0.96

Using Eq. (1), EHOMO can be converted to Em. The calculated Em values correlate with the experimentally measured Em values (Fig. 4, blue diamonds). In the Mn3[M]O2 clusters synthesized by Tsui, each [M] has different ligand groups (Table S1). Accordingly, the absolute Em values are affected by [M] and the ligand groups. To evaluate the direct influence of electrostatic and the van der Waals interactions with [M] on the Em of the Mn3[M]O2 cluster, the metal [M] was removed from the geometry-optimized Mn3[M]O2 cluster. The calculated Em values for the metal-removed clusters do not correlate with the experimentally measured Em values (Fig. 4, red circles)..

Fig. 4

Calculated Em(MnIII/IV) values of the Mn3[M]O2 clusters (blue diamonds) and calculated Em(MnIII/IV) values of the metal-removed Mn3O2 clusters (red circles) plotted with experimentally measured Em(MnIII/IV) values in CH2Cl2 (Tsui and Agapie 2013)

The removal of Y3+ resulted in an increase of 1.5 V in Em, whereas the removal of Na+ resulted in an increase of 0.5 V in Em (Fig. 5a). These results suggest that the valence of [M] is the main factor determining Em. In addition, the removal of a metal with a large radius (e.g., Sr2+) resulted in a large increase in Em, whereas the removal of a metal with a small radius (e.g., Zn2+) resulted in a small increase in Em (Fig. 5b). For metals with the same valence (e.g., Sr2+, Ca2+, and Zn2+), the difference in Em can be explained by the difference in the ionic radius of the redox-inactive metal [M]. As the radius of [M] increases, the distance between [M]2+ and Mn increases, leading to weak electrostatic interactions between Mn and [M]2+. Thus, of the [M]2+ metals, [M] with large radii have a smaller influence on Em. The effect that the ionic radius has on the difference in Em can be explained in terms of the Lewis acidity of [M] because the Lewis acidity decreases with an increase in the ionic radius (Lin et al. 2015).

Fig. 5

Shift in Em upon the removal of [M] ([M] = Na+ (red), Sr2+, Ca2+, Zn2+ (green), and Y3+ (blue)) a with respect to the valence (R2 = 0.93) and b with respect to the ionic radius (Shannon 1976) of [M] ([M] = Sr2+, Ca2+, and Zn2+) (R2 = 0.99)

In PSII, Em(Mn1III/IV) of the Mn4CaO5 cluster changes by only ~40 mV even upon the replacement ion of Ca2+ with H3O+ irrespective of the loss of a +1 charge (Saito et al. 2020a). In contrast, Em(MnIII/IV) of the Mn3[M]O2 cluster changes by 450 mV upon the loss of a +1 charge (Fig. 5a). These indicate that the protein environment including the H-bond network (e.g., D1-Asp61, TyrZ, D1-His190, and water molecules) plays a key role in determining the Em(Mn4CaO5) in PSII.

In summary, the quantum-chemically calculated HOMO energies of artificial Mn3[M]O2 clusters ([M] = Na+, Sr2+, Ca2+, Zn2+, and Y3+) correlate with the experimentally measured Em(MnIII/IV) values (Fig. 3). The Em calculation for the metal-deleted Mn3O2 clusters shows that the valence of [M] predominantly affects Em (Fig. 5a), whereas the ionic radius of [M] affects Em only slightly (Fig. 5b).

Below is the link to the electronic supplementary material.

Supplementary file1 (PDF 200 kb)

Supplementary file3 (PDB 55 kb)

Supplementary file4 (PDB 67 kb)