Labile electronic and spin states of the CaMn4O5 cluster in the PSII system refined to the 1.9 Å X-ray resolution. UB3LYP computational results

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.

Electronic and spin structures of CaMn4Ox clusters in the S0 state of the oxygen evolving complex of photosystem II. Domain-based local pair natural orbital (DLPNO) coupled-cluster (CC) calculations using optimized geometries and natural orbitals (UNO) by hybrid density functional theory (HDFT) calculations

Domain-based local pair natural orbital (DLPNO) coupled cluster single and double (CCSD) with triple perturbation (T) correction methods were performed to elucidate the relative stabilities of ten different intermediate structures of the CaMn4Ox cluster in the S0 state of the oxygen evolving complex (OEC) of photosystem II (PSII). Full geometry optimizations of all the S0 intermediates were performed by the UB3LYP-D3/Def2-TZVP methods, providing the assumed geometrical structures and starting natural orbitals (UNO) for DLPNO-CCSD(T)/Def2TZVP calculations. The effective exchange integrals (J) for the spin Hamiltonian models for the ten intermediates were obtained by the UB3LYP/Def2-TZVP calculations followed by the general spin projections. DLPNO-CCSD(T) calculations followed by the CBS extrapolation procedure elucidated that the (II, III, IV, IV) and (III, III, III, IV) valence states in the CaMn4O5 cluster of the OEC of the PS II were nearly degenerated in energy in the S0 state, indicating an important role of dynamical electron correlation effects for the valence and spin fluctuations in strongly correlated electron systems (SCESs) consisting of 3d transition metals.

Theoretical and computational investigations of geometrical, electronic and spin structures of the CaMn4 OX (X = 5, 6) cluster in the Kok cycle Si (i = 0-3) of oxygen evolving complex of photosystem II

The optimized geometries of the CaMn₄ O_X (X = 5, 6) cluster in the oxygen evolving complex (OEC) of photosystem II (PSII) by large-scale quantum mechanics (QM) and molecular mechanics (MM) calculations are compared with recent serial femtosecond crystallography (SFX) results for the S_i (i = 0-3) states. The valence states of four Mn ions by the QM/MM calculations are also examined in relation to the experimental results by the X-ray emission spectroscopy (XES) for the S_i intermediates. Geometrical and valence structures of right-opened Mn-hydroxide, Mn-oxo and Mn-peroxide intermediates in the S₃ state are investigated in detail in relation to recent SFX and XES experiments for the S₃ state. Interplay between theory and experiment indicates that the Mn-oxo intermediate is a new possible candidate for the S₃ state. Implications of the computational results are discussed in relation to possible mechanisms of the oxygenoxygen bond formation for water oxidation in OEC of PSII.

Understanding Two Different Structures in the Dark Stable State of the Oxygen-Evolving Complex of Photosystem II: Applicability of the Jahn-Teller Deformation Formula

Tanaka et al. (J. Am. Chem. Soc., 2017, 139, 1718) recently reported the three-dimensional (3D) structure of the oxygen evolving complex (OEC) of photosystem II (PSII) by X-ray diffraction (XRD) using extremely low X-ray doses of 0.03 and 0.12 MGy. They observed two different 3D structures of the CaMn4O5 cluster with different hydrogen-bonding interactions in the S1 state of OEC keeping the surrounding polypeptide frameworks of PSII the same. Our Jahn-Teller (JT) deformation formula based on large-scale quantum mechanics/molecular mechanics (QM/MM) was applied for these low-dose XRD structures, elucidating important roles of JT effects of the MnIII ion for subtle geometric distortions of the CaMn4O5 cluster in OEC of PSII. The JT deformation formula revealed the similarity between the low-dose XRD and damage-free serial femtosecond X-ray diffraction (SFX) structures of the CaMn4O5 cluster in the dark stable state. The extremely low-dose XRD structures were not damaged by X-ray irradiation. Implications of the present results are discussed in relation to recent SFX results and a blue print for the design of artificial photocatalysts for water oxidation.

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.

Exchange Coupling Interactions from the Density Matrix Renormalization Group and N-Electron Valence Perturbation Theory: Application to a Biomimetic Mixed-Valence Manganese Complex

The accurate description of magnetic level energetics in oligonuclear exchange-coupled transition-metal complexes remains a formidable challenge for quantum chemistry. The density matrix renormalization group (DMRG) brings such systems for the first time easily within reach of multireference wave function methods by enabling the use of unprecedentedly large active spaces. But does this guarantee systematic improvement in predictive ability and, if so, under which conditions? We identify operational parameters in the use of DMRG using as a test system an experimentally characterized mixed-valence bis-μ-oxo/μ-acetato Mn(III,IV) dimer, a model for the oxygen-evolving complex of photosystem II. A complete active space of all metal 3d and bridge 2p orbitals proved to be the smallest meaningful starting point; this is readily accessible with DMRG and greatly improves on the unrealistic metal-only configuration interaction or complete active space self-consistent field (CASSCF) values. Orbital optimization is critical for stabilizing the antiferromagnetic state, while a state-averaged approach over all spin states involved is required to avoid artificial deviations from isotropic behavior that are associated with state-specific calculations. Selective inclusion of localized orbital subspaces enables probing the relative contributions of different ligands and distinct superexchange pathways. Overall, however, full-valence DMRG-CASSCF calculations fall short of providing a quantitative description of the exchange coupling owing to insufficient recovery of dynamic correlation. Quantitatively accurate results can be achieved through a DMRG implementation of second order N-electron valence perturbation theory (NEVPT2) in conjunction with a full-valence metal and ligand active space. Perspectives for future applications of DMRG-CASSCF/NEVPT2 to exchange coupling in oligonuclear clusters are discussed.

Large-scale QM/MM calculations of the CaMn4O5 cluster in the S3 state of the oxygen evolving complex of photosystem II. Comparison between water-inserted and no water-inserted structures

Large-scale QM/MM calculations were performed to elucidate an optimized geometrical structure of a CaMn₄O₅ cluster with and without water insertion in the S₃ state of the oxygen evolving complex (OEC) of photosystem II (PSII). The left (L)-opened structure was found to be stable under the assumption of no hydroxide anion insertion in the S₃ state, whereas the right (R)-opened structure became more stable if one water molecule is inserted to the Mn₄Ca cluster. The optimized Mn_a(4)–Mn_d(1) distance determined by QM/MM was about 5.0 Å for the S₃ structure without an inserted hydroxide anion, but this is elongated by 0.2–0.3 Å after insertion. These computational results are discussed in relation to the possible mechanisms of O–O bond formation in water oxidation by the OEC of PSII.

Generalized approximate spin projection calculations of effective exchange integrals of the CaMn4O5 cluster in the S1 and S3 states of the oxygen evolving complex of photosystem II

Full geometry optimizations followed by the vibrational analysis were performed for eight spin configurations of the CaMn4O4X(H2O)3Y (X = O, OH; Y = H2O, OH) cluster in the S1 and S3 states of the oxygen evolution complex (OEC) of photosystem II (PSII). The energy gaps among these configurations obtained by vertical, adiabatic and adiabatic plus zero-point-energy (ZPE) correction procedures have been used for computation of the effective exchange integrals (J) in the spin Hamiltonian model. The J values are calculated by the (1) analytical method and the (2) generalized approximate spin projection (AP) method that eliminates the spin contamination errors of UB3LYP solutions. Using J values derived from these methods, exact diagonalization of the spin Hamiltonian matrix was carried out, yielding excitation energies and spin densities of the ground and lower-excited states of the cluster. The obtained results for the right (R)- and left (L)-opened structures in the S1 and S3 states are found to be consistent with available optical and magnetic experimental results. Implications of the computational results are discussed in relation to (a) the necessity of the exact diagonalization for computations of reliable energy levels, (b) magneto-structural correlations in the CaMn4O5 cluster of the OEC of PSII, (c) structural symmetry breaking in the S1 and S3 states, and (d) the right- and left-handed scenarios for the O-O bond formation for water oxidation.

Metal oxidation states in biological water splitting

A central question in biological water splitting concerns the oxidation states of the manganese ions that comprise the oxygen-evolving complex of photosystem II.

A central question in biological water splitting concerns the oxidation states of the manganese ions that comprise the oxygen-evolving complex of photosystem II. Understanding the nature and order of oxidation events that occur during the catalytic cycle of five S_i states (i = 0–4) is of fundamental importance both for the natural system and for artificial water oxidation catalysts. Despite the widespread adoption of the so-called “high-valent scheme”—where, for example, the Mn oxidation states in the S₂ state are assigned as III, IV, IV, IV—the competing “low-valent scheme” that differs by a total of two metal unpaired electrons (i.e. III, III, III, IV in the S₂ state) is favored by several recent studies for the biological catalyst. The question of the correct oxidation state assignment is addressed here by a detailed computational comparison of the two schemes using a common structural platform and theoretical approach. Models based on crystallographic constraints were constructed for all conceivable oxidation state assignments in the four (semi)stable S states of the oxygen evolving complex, sampling various protonation levels and patterns to ensure comprehensive coverage. The models are evaluated with respect to their geometric, energetic, electronic, and spectroscopic properties against available experimental EXAFS, XFEL-XRD, EPR, ENDOR and Mn K pre-edge XANES data. New 2.5 K ⁵⁵Mn ENDOR data of the S₂ state are also reported. Our results conclusively show that the entire S state phenomenology can only be accommodated within the high-valent scheme by adopting a single motif and protonation pattern that progresses smoothly from S₀ (III, III, III, IV) to S₃ (IV, IV, IV, IV), satisfying all experimental constraints and reproducing all observables. By contrast, it was impossible to construct a consistent cycle based on the low-valent scheme for all S states. Instead, the low-valent models developed here may provide new insight into the over-reduced S states and the states involved in the assembly of the catalytically active water oxidizing cluster.

The Structure of Photosystem II and the Mechanism of Water Oxidation in Photosynthesis

Oxygenic photosynthesis forms the basis of aerobic life on earth by converting light energy into biologically useful chemical energy and by splitting water to generate molecular oxygen. The water-splitting and oxygen-evolving reaction is catalyzed by photosystem II (PSII), a huge, multisubunit membrane-protein complex located in the thylakoid membranes of organisms ranging from cyanobacteria to higher plants. The structure of PSII has been analyzed at 1.9-Å resolution by X-ray crystallography, revealing a clear picture of the Mn4CaO5 cluster, the catalytic center for water oxidation. This article provides an overview of the overall structure of PSII followed by detailed descriptions of the specific structure of the Mn4CaO5 cluster and its surrounding protein environment. Based on the geometric organization of the Mn4CaO5 cluster revealed by the crystallographic analysis, in combination with the results of a vast number of experimental studies involving spectroscopic and other techniques as well as various theoretical studies, the article also discusses possible mechanisms for water splitting that are currently under consideration.

Convergence of QM/MM and Cluster Models for the Spectroscopic Properties of the Oxygen-Evolving Complex in Photosystem II

The latest crystal structure of photosystem II at 1.9 Å resolution, which resolves the topology of the Mn4CaO5 oxygen evolving complex (OEC) at atomistic detail, enables a better correlation between structural features and spectroscopic properties than ever before. Building on the refined crystallographic model of the OEC and the protein, we present combined quantum mechanical/molecular mechanical (QM/MM) studies of the spectroscopic properties of the natural catalyst embedded in the protein matrix. Focusing on the S2 state of the catalytic cycle, we examine the convergence of not only structural parameters but also of the intracluster magnetic interactions in terms of exchange coupling constants and of experimentally relevant (55)Mn, (17)O, and (14)N hyperfine coupling constants with respect to QM/MM partitioning using five QM regions of increasing size. This enables us to assess the performance of the method and to probe second sphere effects by identifying amino acid residues that principally affect the spectroscopic properties of the OEC. Comparison between QM-only and QM/MM treatments reveals that whereas QM/MM models converge quickly to stable values, the QM cluster models need to incorporate significantly larger parts of the second coordination sphere and surrounding water molecules to achieve convergence for certain properties. This is mainly due to the sensitivity of the QM-only models to fluctuations in the hydrogen bonding network and ligand acidity. Additionally, a hydrogen bond that is typically omitted in QM-only treatments is shown to determine the hyperfine coupling tensor of the unique Mn(III) ion by regulating the rotation plane of the ligated D1-His332 imidazole ring, the only N-donor ligand of the OEC.

Similarities of artificial photosystems by ruthenium oxo complexes and native water splitting systems

The nature of chemical bonds of ruthenium(Ru)-quinine(Q) complexes, mononuclear [Ru(trpy)(3,5-t-Bu(2)Q)(OH(2))](ClO(4))(2) (trpy = 2,2':6',2''-terpyridine, 3,5-di-tert-butyl-1,2-benzoquinone) (1), and binuclear [Ru(2)(btpyan)(3,6-di-Bu(2)Q)(2)(OH(2))](2+) (btpyan = 1,8-bis(2,2':6',2''-terpyrid-4'-yl)anthracene, 3,6-t-Bu(2)Q = 3,6-di-tert-butyl-1,2-benzoquinone) (2), has been investigated by broken-symmetry (BS) hybrid density functional (DFT) methods. BS DFT computations for the Ru complexes have elucidated that the closed-shell structure (2b) Ru(II)-Q complex is less stable than the open-shell structure (2bb) consisting of Ru(III) and semiquinone (SQ) radical fragments. These computations have also elucidated eight different electronic and spin structures of tetraradical intermediates that may be generated in the course of water splitting reaction. The Heisenberg spin Hamiltonian model for these species has been derived to elucidate six different effective exchange interactions (J) for four spin systems. Six J values have been determined using total energies of the eight (or seven) BS solutions for different spin configurations. The natural orbital analyses of these BS DFT solutions have also been performed in order to obtain natural orbitals and their occupation numbers, which are useful for the lucid understanding of the nature of chemical bonds of the Ru complexes. Implications of the computational results are discussed in relation to the proposed reaction mechanisms of water splitting reaction in artificial photosynthesis systems and the similarity between artificial and native water splitting systems.

Influence of the axial ligand on the cationic properties of the chlorophyll pair in photosystem II from Thermosynechococcus vulcanus

Influence of the axial ligand of PD1 chlorophyll (D1-His-198) on the Em of monomer chlorophylls PD1 and PD2, and the PD1•+/PD2•+ charge ratio was investigated by theoretical calculations using the PSII crystal structure of Thermosynechococcus vulcanus analyzed at 1.9-Å resolution. It was found that the Em(PD1)/Em(PD2) values and PD1•+/PD2•+ ratio remained unchanged upon D1-H198Q mutation. However, Em(PD1) was increased in the D1-H198A mutant, resulting in a more even distribution of the positive charge over PD1/PD2. Introduction of a water molecule as an axial ligand resulted in equal Em values and PD1•+/PD2•+ ratios between the mutant and wild-type, thus confirming the presence of the water ligand in the mutant.

Cationic state distribution over the chlorophyll d-containing P(D1)/P(D2) pair in photosystem II

Most of the chlorophyll (Chl) cofactors in photosystem II (PSII) from Acaryochloris marina are Chld, although a few Chla molecules are also present. To evaluate the possibility that Chla may participate in the P(D1)/P(D2) Chl pair in PSII from A. marina, the P(D1)(•+)/P(D2)(•+) charge ratio was investigated using the PSII crystal structure analyzed at 1.9-Å resolution, while considering all possibilities for the Chld-containing P(D1)/P(D2) pair, i.e., Chld/Chld, Chla/Chld, and Chld/Chla pairs. Chld/Chld and Chla/Chld pairs resulted in a large P(D1)(•+) population relative to P(D2)(•+), as identified in Chla/Chla homodimer pairs in PSII from other species, e.g., Thermosynechococcus elongatus PSII. However, the Chld/Chla pair possessed a P(D1)(•+)/P(D2)(•+) ratio of approximately 50/50, which is in contrast to previous spectroscopic studies on A. marina PSII. The present results strongly exclude the possibility that the Chld/Chla pair serves as P(D1)/P(D2) in A. marina PSII. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.