Theory of chemical bonds in metalloenzymes VI: Manganese–oxo bonds in the photosynthesis II system

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.

Controlling Chirogenic Effects in Porphyrin Based Supramolecular Systems: Theoretical Analysis Versus Experimental Observations

Electronic circular dichroism (ECD) spectroscopy is a widely employed method for studying chiral analysis, requiring the presence of a chromophore close to a chiral centre. Porphyrinoids are found to be one of the best chromophoric systems serving for this purpose and enabling the application of ECD spectroscopy for chirality determination across diverse classes of organic compounds. Consequently, it is crucial to understand the induction mechanisms of ECD in the porphyrin-based complexes. The present study explores systematically the influence of secondary chromophores, bonded to an achiral zinc porphyrin or to chiral guest molecules, on the B-region of ECD spectra using the time-dependent density functional theory (TD-DFT) calculations. The study analyses the impact of change in both the conformation of achiral porphyrin (host) and change in position and conformation of chiral organic molecule (guest) on the B-band of ECD spectra (energy, intensity, sign of Cotton effect). Finally, conclusions made on model complexes are applied to published experimental data, contributing to a deeper understanding of various factors influencing ECD spectra in chiral systems. In addition, a computer program aimed to help rationalise ECD spectra by visualizing corresponding orbital energies, rotatory strengths, electric and magnetic transition moments, and angles between them, is presented.

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.

Spin, Valence, and Structural Isomerism in the S3 State of the Oxygen-Evolving Complex of Photosystem II as a Manifestation of Multimetallic Cooperativity

Photosynthetic water oxidation is catalyzed by a Mn₄CaO₅-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 S₃ 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 ( S_total = 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 Mn_A···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.

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.

A light harvesting mononuclear manganese(ii) complex: synthesis, characterization, DFT and TDDFT calculations and photophysical profile

A new manganese(ii) [Mn^II(DEMP)(NCS)(H₂O)] (DEMP = Schiff base derived from salicylaldehyde and 2-diethylaminoethylamine) complex has been synthesized and characterized. The complex absorbs light ranging from 200–850 nm. Thus, the molecule is capable of harvesting the entire range of sunlight falling on earth.

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.

Theoretical evaluation of structural models of the S2 state in the oxygen evolving complex of Photosystem II: protonation states and magnetic interactions

Protonation states of water ligands and oxo bridges are intimately involved in tuning the electronic structures and oxidation potentials of the oxygen evolving complex (OEC) in Photosystem II, steering the mechanistic pathway, which involves at least five redox state intermediates S(n) (n = 0-4) resulting in the oxidation of water to molecular oxygen. Although protons are practically invisible in protein crystallography, their effects on the electronic structure and magnetic properties of metal active sites can be probed using spectroscopy. With the twin purpose of aiding the interpretation of the complex electron paramagnetic resonance (EPR) spectroscopic data of the OEC and of improving the view of the cluster at the atomic level, a complete set of protonation configurations for the S(2) state of the OEC were investigated, and their distinctive effects on magnetic properties of the cluster were evaluated. The most recent X-ray structure of Photosystem II at 1.9 Å resolution was used and refined to obtain the optimum structure for the Mn(4)O(5)Ca core within the protein pocket. Employing this model, a set of 26 structures was constructed that tested various protonation scenarios of the water ligands and oxo bridges. Our results suggest that one of the two water molecules that are proposed to coordinate the outer Mn ion (Mn(A)) of the cluster is deprotonated in the S(2) state, as this leads to optimal experimental agreement, reproducing the correct ground state spin multiplicity (S = 1/2), spin expectation values, and EXAFS-derived metal-metal distances. Deprotonation of Ca(2+)-bound water molecules is strongly disfavored in the S(2) state, but dissociation of one of the two water ligands appears to be facile. The computed isotropic hyperfine couplings presented here allow distinctions between models to be made and call into question the assumption that the largest coupling is always attributable to Mn(III). The present results impose limits for the total charge and the proton configuration of the OEC in the S(2) state, with implications for the cascade of events in the Kok cycle and for the water splitting mechanism.

Density functional theory for transition metals and transition metal chemistry

We introduce density functional theory and review recent progress in its application to transition metal chemistry. Topics covered include local, meta, hybrid, hybrid meta, and range-separated functionals, band theory, software, validation tests, and applications to spin states, magnetic exchange coupling, spectra, structure, reactivity, and catalysis, including molecules, clusters, nanoparticles, surfaces, and solids.