ABSTRACT
Photosystem II (PSII) catalyses the oxidation of water through a four-step cycle of Si states (i = 0-4) at the Mn4CaO5 cluster1-3, during which an extra oxygen (O6) is incorporated at the S3 state to form a possible dioxygen4-7. Structural changes of the metal cluster and its environment during the S-state transitions have been studied on the microsecond timescale. Here we use pump-probe serial femtosecond crystallography to reveal the structural dynamics of PSII from nanoseconds to milliseconds after illumination with one flash (1F) or two flashes (2F). YZ, a tyrosine residue that connects the reaction centre P680 and the Mn4CaO5 cluster, showed structural changes on a nanosecond timescale, as did its surrounding amino acid residues and water molecules, reflecting the fast transfer of electrons and protons after flash illumination. Notably, one water molecule emerged in the vicinity of Glu189 of the D1 subunit of PSII (D1-E189), and was bound to the Ca2+ ion on a sub-microsecond timescale after 2F illumination. This water molecule disappeared later with the concomitant increase of O6, suggesting that it is the origin of O6. We also observed concerted movements of water molecules in the O1, O4 and Cl-1 channels and their surrounding amino acid residues to complete the sequence of electron transfer, proton release and substrate water delivery. These results provide crucial insights into the structural dynamics of PSII during S-state transitions as well as O-O bond formation.
Subject(s)
Oxygen , Photosystem II Protein Complex , Biocatalysis/radiation effects , Calcium/metabolism , Crystallography , Electron Transport/radiation effects , Electrons , Manganese/metabolism , Oxidation-Reduction/radiation effects , Oxygen/chemistry , Oxygen/metabolism , Photosystem II Protein Complex/chemistry , Photosystem II Protein Complex/metabolism , Photosystem II Protein Complex/radiation effects , Protons , Time Factors , Tyrosine/metabolism , Water/chemistry , Water/metabolismABSTRACT
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.
ABSTRACT
This review article describes a historical perspective of elucidation of the nature of the chemical bonds of the high-valent transition metal oxo (M=O) and peroxo (M-O-O) compounds in chemistry and biology. The basic concepts and theoretical backgrounds of the broken-symmetry (BS) method are revisited to explain orbital symmetry conservation and orbital symmetry breaking for the theoretical characterization of four different mechanisms of chemical reactions. Beyond BS methods using the natural orbitals (UNO) of the BS solutions, such as UNO CI (CC), are also revisited for the elucidation of the scope and applicability of the BS methods. Several chemical indices have been derived as the conceptual bridges between the BS and beyond BS methods. The BS molecular orbital models have been employed to explain the metal oxyl-radical character of the M=O and M-O-O bonds, which respond to their radical reactivity. The isolobal and isospin analogy between carbonyl oxide R2C-O-O and metal peroxide LFe-O-O has been applied to understand and explain the chameleonic chemical reactivity of these compounds. The isolobal and isospin analogy among Fe=O, O=O, and O have also provided the triplet atomic oxygen (3O) model for non-heme Fe(IV)=O species with strong radical reactivity. The chameleonic reactivity of the compounds I (Cpd I) and II (Cpd II) is also explained by this analogy. The early proposals obtained by these theoretical models have been examined based on recent computational results by hybrid DFT (UHDFT), DLPNO CCSD(T0), CASPT2, and UNO CI (CC) methods and quantum computing (QC).
ABSTRACT
The primary coordination sphere of the multinuclear cofactor (Mn4CaOx) in the oxygen-evolving complex (OEC) of photosystem II is absolutely conserved to maintain its structure and function. Recent time-resolved serial femtosecond crystallography identified large reorganization of the primary coordination sphere in the S2 to S3 transition, which elicits a cascade of events involving Mn oxidation and water molecule binding to a putative catalytic Mn site. We examined how the crystallographic fields, created by transient conformational states of the OEC at various time points, affect the thermodynamics of various isomers of the Mn cluster using DFT calculations, with an aim of comprehending the functional roles of the flexible primary coordination sphere in the S2 to S3 transition and in the recovery of the S2 state. The results show that the relative movements of surrounding residues change the size and shape of the cavity of the cluster and thereby affect the thermodynamics of various catalytic intermediates as well as the ability to capture a new water molecule at a coordinatively unsaturated site. The implication of these findings is that the protein dynamics may serve to gate the catalytic reaction efficiently by controlling the sequence of Mn oxidation/reduction and water binding/release. This interpretation is consistent with EPR experiments; g â¼ 5 and g â¼ 3 signals obtained after near-infrared (NIR) excitation of the S3 state at 4 K and a g â¼ 5 only signal produced after prolonged incubation of the S3 state at 77 K can be best explained as originating from water-bound S2 clusters (Stotal = 7/2) under a S3 ligand field, i.e., the immediate one-electron reduction products of the oxyl-oxo (Stotal = 6) and hydroxo-oxo (Stotal = 3) species in the S3 state.
Subject(s)
Oxygen , Photosystem II Protein Complex , Electron Spin Resonance Spectroscopy , Ligands , Oxidation-Reduction , Oxygen/chemistry , Photosystem II Protein Complex/chemistry , Water/chemistryABSTRACT
Photosynthetic water oxidation is catalyzed by the Mn4CaO5 cluster of photosystem II (PSII) with linear progression through five S-state intermediates (S0 to S4). To reveal the mechanism of water oxidation, we analyzed structures of PSII in the S1, S2, and S3 states by x-ray free-electron laser serial crystallography. No insertion of water was found in S2, but flipping of D1 Glu189 upon transition to S3 leads to the opening of a water channel and provides a space for incorporation of an additional oxygen ligand, resulting in an open cubane Mn4CaO6 cluster with an oxyl/oxo bridge. Structural changes of PSII between the different S states reveal cooperative action of substrate water access, proton release, and dioxygen formation in photosynthetic water oxidation.
Subject(s)
Oxygen/chemistry , Photosystem II Protein Complex/chemistry , Photosystem II Protein Complex/metabolism , Water/chemistry , Calcium/chemistry , Crystallography, X-Ray , Fourier Analysis , Hydrogen/chemistry , Hydrogen Bonding , Lasers , Ligands , Manganese/chemistry , Models, Molecular , Oxidation-Reduction , Oxygen/metabolism , Protein Conformation , Water/metabolismABSTRACT
The optimized geometries of the CaMn4 OX (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 Si (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 Si intermediates. Geometrical and valence structures of right-opened Mn-hydroxide, Mn-oxo and Mn-peroxide intermediates in the S3 state are investigated in detail in relation to recent SFX and XES experiments for the S3 state. Interplay between theory and experiment indicates that the Mn-oxo intermediate is a new possible candidate for the S3 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.
Subject(s)
Oxygen/metabolism , Photosystem II Protein Complex/metabolism , Crystallography , Manganese , Oxidation-Reduction , Spectrometry, X-Ray EmissionABSTRACT
Photosynthetic water oxidation is catalyzed by a Mn4CaO5-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 S3 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 ( Stotal = 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 MnA···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.
Subject(s)
Bacterial Proteins/chemistry , Manganese/chemistry , Oxygen/chemistry , Photosystem II Protein Complex/chemistry , Synechococcus/chemistry , Electron Spin Resonance Spectroscopy , Isomerism , Models, Molecular , Oxidation-Reduction , Water/chemistryABSTRACT
We applied our analysis, based on a linear response function of density and spin density, to two typical transition metal complex systems-the reaction centers of P450, and oxygen evolving center in Photosystem II, both of which contain open-shell transition metal ions. We discuss the relationship between LRF of electron density and spin density and the types of units and interactions of the systems. The computational results are discussed in relation to quantum mechanics (QM) cluster and quantum mechanics/molecular mechanics (QM/MM) modeling that are employed to compute the reaction centers of enzymes.
Subject(s)
Coordination Complexes/chemistry , Oxygen/chemistry , Transition Elements/chemistry , Catalysis , Electron Transport , Molecular Dynamics Simulation , Molecular Structure , Oxidation-Reduction , Photochemical Processes , Photosystem II Protein Complex/chemistry , Quantum Theory , WaterABSTRACT
Both direct exchange and super-exchange interactions cooperate to realize inter-spin magnetic interaction in binuclear manganese complex Mn(IV)2 O2 (NHCHCO2 )4 with a di-µ-oxo path. We revisited this spin system using DMRG CAS methods and CAS selection procedures. Our results indicate that our previous "dynamically extended spin polarization" (DE-SP) procedure for organic polyradicals and so forth does not work well. Thus, we have examined another selection procedure, the "dynamically extended super-exchange" (DE-SE) procedure. DMRG CASCI [18,18] by UB3LYP(HS)-UNO(DE-SE) can realize antiferromagnetic J values similar to experimental ones (-87 cm-1 ). In addition, all J values between all spin states (HS[septet],IS[quintet],IS[triplet],LS[singlet])were also shown to be correct under sufficiently large M values. © 2018 Wiley Periodicals, Inc.
ABSTRACT
The O2 release of the oxygen-evolving complex of the photosystem II (PSII) is one of the essential processes responsible for the highly efficient O2 production. Despite its importance, the detailed molecular mechanism is still unsolved. In the present study, we show that the O2 release is directly coupled with water insertion into the Mn cluster based on the quantum mechanics/molecular mechanics (QM/MM) calculations. In this mechanism, the O2 molecule first dissociates from the Mn sites in order, that is, the O atom coordinating to the Mn3 (O5a) first dissociates, then the other O atom coordinating to the Mn1 (O5d) dissociates in the next step in the late S4 state (1 â 2). Next, the O2 migrates to a space surrounded by the Val185 and His332 side chains as one water molecule coordinating to the Ca2+ ion (W3) comes into the O2 bonded site (2 â 3). Finally, a pre-S0 state (4) is formed after a proton transfer from the inserted water to the other proton acceptor site (W2) (3 â 4). The highest activation barrier during these reactions was found at the O2 release step (2 â 3) that only requires E⧧ = 12.7 kcal mol-1 ( G⧧ = 10.4 kcal mol-1). A series of the reactions (2 â 3) look like a chain crash of billiard balls because the W3 is inserted into the catalytic center from the water-abundant side (Ca2+ ion side), and then the O2 moiety is pushed out to the opposite side (Val185 side). The hydrophobic residue of Val185 covers the active O5 site and forms an O2-specific permeation tunnel. The present sequential reactions clearly demonstrate the efficient removal of the toxic O2 from the catalytic center and implications of the essential roles of Val185, Ca2+ ions and water molecules, which are all present in the active site of PSII as the indispensable constituents.
Subject(s)
Oxygen/chemistry , Photosystem II Protein Complex/chemistry , Water/chemistry , Catalytic Domain , Manganese/chemistry , Molecular Dynamics Simulation , Oxygen/metabolism , Photosystem II Protein Complex/metabolism , Quantum Theory , Static Electricity , ThermodynamicsABSTRACT
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.
ABSTRACT
Large-scale QM/MM calculations were performed to elucidate an optimized geometrical structure of a CaMn4O5 cluster with and without water insertion in the S3 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 S3 state, whereas the right (R)-opened structure became more stable if one water molecule is inserted to the Mn4Ca cluster. The optimized Mna(4)-Mnd(1) distance determined by QM/MM was about 5.0 Å for the S3 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.
ABSTRACT
Water oxidation by photosystem II (PSII) converts light energy into chemical energy with the concomitant production of molecular oxygen, both of which are indispensable for sustaining life on Earth. This reaction is catalyzed by an oxygen-evolving complex (OEC) embedded in the huge PSII complex, and its mechanism remains elusive in spite of the extensive studies of the geometric and electronic structures. In order to elucidate the water-splitting mechanism, synthetic approaches have been extensively employed to mimic the native OEC. Very recently, a synthetic complex [Mn4CaO4(Bu(t)COO)8(py)(Bu(t)COOH)2] (1) closely mimicking the structure of the native OEC was obtained. In this study, we extensively examined the geometric, electronic and spin structures of 1 using the density functional theory method. Our results showed that the geometric structure of 1 can be accurately reproduced by theoretical calculations, and revealed many similarities in the ground valence and spin states between 1 and the native OEC. We also revealed two different valence states in the one-electron oxidized state of 1 (corresponding to the S2 state), which lie in the lower and higher ground spin states (S = 1/2 and S = 5/2), respectively. One remarkable difference between 1 and the native OEC is the presence of a non-negligible antiferromagnetic interaction between the Mn1 and Mn4 sites, which slightly influenced their ground spin structures (spin alignments). The major reason causing the difference can be attributed to the short Mn1-O5 and Mn1-Mn4 distances in 1. The introduction of the missing O4 atom and the reorientation of the Ca coordinating ligands improved the Mn1-O5 and Mn1-Mn4 distances comparable to the native OEC. These modifications will therefore be important for the synthesis of further advanced model complexes more closely mimicking the native OEC beyond 1.
Subject(s)
Calcium Compounds/chemistry , Manganese Compounds/chemistry , Models, Chemical , Oxides/chemistry , Oxygen/chemistry , Photosynthesis , Electron Spin Resonance Spectroscopy , Quantum TheoryABSTRACT
We have performed hybrid density functional theory (DFT) calculations to investigate how chemical equilibria can be described in the S3 state of the oxygen-evolving complex in photosystem II. For a chosen 340-atom model, 1 stable and 11 metastable intermediates have been identified within the range of 13 kcal mol(-1) that differ in protonation, charge, spin, and conformational states. The results imply that reversible interconversion of these intermediates gives rise to dynamic equilibria that involve processes with relocations of protons and electrons residing in the Mn4CaO5 cluster, as well as bound water ligands, with concomitant large changes in the cluster geometry. Such proton tautomerism and redox isomerism are responsible for reversible activation/deactivation processes of substrate oxygen species, through which Mn-O and O-O bonds are transiently ruptured and formed. These results may allow for a tentative interpretation of kinetic data on substrate water exchange on the order of seconds at room temperature, as measured by time-resolved mass spectrometry. The reliability of the hybrid DFT method for the multielectron redox reaction in such an intricate system is also addressed.
Subject(s)
Models, Chemical , Oxygen/chemistry , Photosystem II Protein Complex/chemistry , Crystallography, X-Ray , Isomerism , Oxidation-Reduction , Proteins/chemistry , ThermodynamicsABSTRACT
We have studied the early phase of the S2 â S3 transition in the oxygen-evolving complex (OEC) of photosystem II using the hybrid density functional theory with a quantum mechanical model composed of 338-341 atoms. Special attention is given to the vital role of water molecules in the vicinity of the Mn4CaO5 core. Our results demonstrate how important the dynamic behavior of surrounding water molecules is in mediating critical chemical transformations such as binding and deprotonation of substrates and hydration of the catalytic site and identify a strong coupling of water-chain relocation near the redox-active tyrosine residue Tyr161 (TyrZ) with oxidation of the Mn4CaO5 cluster by TyrZ(â¢+). The oxidation reaction is further promoted when the catalytic site is more solvated by water. These results indicate the importance of surrounding water molecules in biological catalysts as they ultimately lead to effective catalytic function and/or favorable electron-transfer dynamics.
Subject(s)
Oxygen/chemistry , Oxygen/metabolism , Photosystem II Protein Complex/chemistry , Photosystem II Protein Complex/metabolism , Electron Transport , Hydrogen Bonding , Models, Molecular , Oxidation-Reduction , Quantum TheoryABSTRACT
We examined nearsightedness of electronic matter (NEM) of finite systems on the basis of linear response function (LRF). From the computational results of a square-well model system, the behavior of responses obviously depends on the number of electrons (N): as N increases, LRF, δρ(r)/δv(r'), decays rapidly for the distance, |r-r'|. This exemplifies that the principle suggested by Kohn and Prodan holds even for finite systems: the cause of NEM is destructive interference among electron density amplitudes. In addition, we examined double-well model systems, which have low-lying degenerate levels. In this case, there are two types of LRF: the cases of the half-filled and of full-filled in low-lying degenerate levels. The response for the former is delocalized, while that of the later is localized. These behaviors of model systems are discussed in relation to the molecular systems' counterparts, H2, He22+, and He2 systems. We also see that NEM holds for the dissociated limit of H2, of which the mechanism is similar to that of the insulating state of solids as suggested by Kohn. We also examined LRF of alanine tripeptide system as well as butane and butadiene molecules, showing that NEM of the polypeptide system is caused by sp3 junctions at Cα atoms that prevent propagation of amplitudes of LRF, which is critically different from that of NEM for finite and infinite homogeneous systems.
Subject(s)
Computer Simulation , Algorithms , Linear Models , Models, Molecular , Oligopeptides/chemistry , Protein Structure, Secondary , Quantum TheoryABSTRACT
We investigated theoretically the catalytic mechanism of electrochemical water oxidation in aqueous solution by a dinuclear ruthenium complex containing redox-active quinone ligands, [Ru2(X)(Y)(3,6-tBu2Q)2(btpyan)](m+) [X, Y = H2O, OH, O, O2; 3,6-tBu2Q = 3,6-di-tert-butyl-1,2-benzoquinone; btpyan =1,8-bis(2,2':6',2â³-terpyrid-4'-yl)anthracene] (m = 2, 3, 4) (1). The reaction involves a series of electron and proton transfers to achieve redox leveling, with intervening chemical transformations in a mesh scheme, and the entire molecular structure and motion of the catalyst 1 work together to drive the catalytic cycle for water oxidation. Two substrate water molecules can bind to 1 with simultaneous loss of one or two proton(s), which allows pH-dependent variability in the proportion of substrate-bound structures and following pathways for oxidative activation of the aqua/hydroxo ligands at low thermodynamic and kinetic costs. The resulting bis-oxo intermediates then undergo endothermic O-O radical coupling between two Ru(III)-O(â¢) units in an anti-coplanar conformation leading to bridged µ-peroxo or µ-superoxo intermediates. The µ-superoxo species can liberate oxygen with the necessity for the preceding binding of a water molecule, which is possible only after four-electron oxidation is completed. The magnitude of catalytic current would be limited by the inherent sluggishness of the hinge-like bending motion of the bridged µ-superoxo complex that opens up the compact, hydrophobic active site of the catalyst and thereby allows water entry under dynamic conditions. On the basis of a newly proposed mechanism, we rationalize the experimentally observed behavior of electrode kinetics with respect to potential and discuss what causes a high overpotential for water oxidation by 1.
Subject(s)
Benzoquinones/chemistry , Organometallic Compounds/chemistry , Ruthenium/chemistry , Water/chemistry , Ligands , Molecular Structure , Oxidation-Reduction , Quantum Theory , ThermodynamicsABSTRACT
Electronic structures, magnetic properties and electron conductivities of linearly aligned one-dimensional (1-D) Ni(II)3, Ni(II)5 and Ni(II)7 complexes, i.e. [Ni3(dpa)4NCS2], [Ni5(tpda)4X2] (X = Cl, CN, N3, NCS) and [Ni7(teptra)4Cl2], are systematically investigated by the broken-symmetry B3LYP calculations and simulations based on an elastic scattering Green's function theory. Calculated spin densities appear only at terminal Ni ions, while inner Ni ions are the closed-shell. The calculated effective exchange integrals (J(ab)) values reproduce well the experimental results that indicate anti-ferromagnetic (AF) interactions between two terminal Ni ions. Natural orbitals and their occupation numbers show that a change in the weak AF couplings by axial ligands in penta-nickel complexes originates in σ-type orbitals. Simulated electron conductivities of [Ni3(dpa)4NCS2] and [Ni5(tpda)4NCS2] semi-quantitatively correspond to the experimental results. By the analyses, it is elucidated that electrons are mainly transmitted by σ-type orbitals, but the bonds between Au and axial ligands are also dominant factors for conductivity.
ABSTRACT
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.