Search

Iron-sulfur clusters play a crucial role in electron transfer for many essential enzymes, including [FeFe]-hydrogenases. This study focuses on the [4Fe4S] cluster ([4Fe]H) of the minimal [FeFe]-hydrogenase from Chlamydomonas reinhardtii (CrHydA1) and employs advanced spectroscopy, site-directed mutagenesis, molecular dynamics simulations, and QM/MM calculations. We provide insights into the complex electronic structure of [4Fe]H and its role in the catalytic reaction of CrHydA1, serving as paradigm for understanding [FeFe]-hydrogenases. We identified at least two distinct species within the apo-form of CrHydA1, designated 4Fe-R and 4Fe-A, with unique redox potentials and pH sensitivities. Our findings revealed that these species arise from a complex interplay of structural heterogeneity and valence isomer rearrangements, influenced by second-sphere residues. We propose that the interconversion between 4Fe-R and 4Fe-A could provide control over electron transfer in the absence of accessory FeS clusters typically found in other [FeFe]-hydrogenases. The insights gained from this study not only enhance our understanding of [FeFe]-hydrogenases but also provide a crucial foundation for future investigations into analysis of other FeS clusters across diverse biological systems., (© 2025 Wiley‐VCH GmbH.)

This study reports that the well-known Co(III) precursor Co(η 5 -Cp*)I 2 (CO) ( 1 ) is an effective precatalyst for dehydrocoupling reactions of phosphines. Reaction monitoring by NMR, complemented by ESI-MS and EPR, suggests that Co(III) complexes containing secondary and primary phosphine ligands play an important role in this catalysis but that redox chemistry is also facile for these complexes, resulting in paramagnetic species. Representative examples of the mono(phosphine) complexes Co(η 5 -Cp*)I 2 (PR 2 H) ( 2 ) and Co(η 5 -Cp*)I 2 (PRH 2 ) ( 4 ) and bis(phosphine) complexes [Co(η 5 -Cp*)I(PR 2 H) 2 ] I ( 3 ) have been prepared. The characterization of these complexes through structural, computational, spectroscopic, and electrochemical analysis is reported and serves as a foundation for considering their mechanistic significance in the observed catalytic dehydrocoupling chemistry.

Here we report the reaction of aryl thianthrenium salts that allows selective functionalization of the meta position of arenes. The combination of a site-selective thianthrenation with a Catellani reaction provides access to 3,5-dimethylated arenes. The developed reaction is complementary to the previously discovered reductive ipso-alkylation of aryl thianthrenium salts and extends the possibilities for late-stage methylation of arenes with a single aryl thianthrenium salt., (© 2024 The Author(s). Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Mössbauer spectroscopy is a powerful technique for probing the local electronic structure of iron compounds, because it reports in an element-selective manner on both the oxidation state and coordination environment of the Fe ion. Computational prediction of the two main Mössbauer parameters, isomer shift ( δ ) and quadrupole splitting (Δ E Q ), has long been targeted by quantum chemical studies, and useful protocols based on density functional theory have been proposed. Here we present an extensive curated reference set of Fe compounds that is considerably larger and more diverse than literature precedents. We make a distinction between low-temperature and high-temperature experimental subgroups. This set is employed for optimizing a refined computational protocol utilizing the scalar version of the exact 2-component (X2C) Hamiltonian with the finite nucleus approximation. Attention is devoted to having an accurate and flexible all-electron basis set for Fe. We assess the performance of several DFT methods that cover all representative families and rungs of functionals and find that hybrid functionals with ca. 25-30% exact exchange offer the best accuracy for isomer shifts. The work establishes a refined general protocol of wide applicability that achieves good performance for the prediction of isomer shifts in a wider variety of systems than before, but the limitations of DFT for quadrupole splittings are also highlighted. Finally, comparison of calculated values with high-temperature experimental results shows that the use of an empirical correction factor is required to account for the second-order Doppler shift and to achieve the same quality of correlation as with the low-temperature data.

A first-principles description of the primary photochemical processes that drive photosynthesis and sustain life on our planet remains one of the grand challenges of modern science. Recent research established that explicit incorporation of protein electrostatics in excited-state calculations of photosynthetic pigments, achieved for example with quantum-mechanics/molecular-mechanics (QM/MM) approaches, is essential for a meaningful description of the properties and function of pigment-protein complexes. Although time-dependent density functional theory has been used productively so far in QM/MM approaches for the study of such systems, this methodology has limitations. Here we pursue for the first time a QM/MM description of the reaction center in the principal enzyme of oxygenic photosynthesis, Photosystem II, using multireference wave function theory for the high-level QM region. We identify best practices and establish guidelines regarding the rational choice of active space and appropriate state-averaging for the efficient and reliable use of complete active space self-consistent field (CASSCF) and the N-electron valence state perturbation theory (NEVPT2) in the prediction of low-lying excited states of chlorophyll and pheophytin pigments. Given that the Gouterman orbitals are inadequate as a minimal active space, we define specific minimal and extended active spaces for the NEVPT2 description of electronic states that fall within the Q and B bands. Subsequently, we apply our multireference-QM/MM protocol to the description of all pigments in the reaction center of Photosystem II. The calculations reproduce the electrochromic shifts induced by the protein matrix and the ordering of site energies consistent with the identity of the primary donor (Chl D1 ) and the experimentally known asymmetric and directional electron transfer. The optimized protocol sets the stage for future multireference treatments of multiple pigments, and hence for multireference studies of charge separation, while it is transferable to the study of any photoactive embedded tetrapyrrole system.

The Mn 4 CaO 5(6) cluster in photosystem II catalyzes water splitting through the S i state cycle (i = 0-4). Molecular O 2 is formed and the natural catalyst is reset during the final S 3  → (S 4 ) → S 0 transition. Only recently experimental breakthroughs have emerged for this transition but without explicit information on the S 0 -state reconstitution, thus the progression after O 2 release remains elusive. In this report, our molecular dynamics simulations combined with density functional calculations suggest a likely missing link for closing the cycle, i.e., restoring the first catalytic state. Specifically, the formation of closed-cubane intermediates with all hexa-coordinate Mn is observed, which would undergo proton release, water dissociation, and ligand transfer to produce the open-cubane structure of the S 0 state. Thereby, we theoretically identify the previously unknown structural isomerism in the S 0 state that acts as the origin of the proposed structural flexibility prevailing in the cycle, which may be functionally important for nature's water oxidation catalysis., (© 2024. The Author(s).)

The Light-Dependent Protochlorophyllide Oxidoreductase (LPOR) catalyzes a crucial step in chlorophyll biosynthesis: the rare biological photocatalytic reduction of the double C[double bond, length as m-dash]C bond in the precursor, protochlorophyllide (Pchlide). Despite its fundamental significance, limited structural insights into the active complex have hindered understanding of its reaction mechanism. Recently, a high-resolution cryo-EM structure of LPOR in its active conformation challenged our view of pigment binding, residue interactions, and the catalytic process. Surprisingly, this structure contrasts markedly with previous assumptions, particularly regarding the orientation of the bound Pchlide. To gain insights into the substrate binding puzzle, we conducted molecular dynamics simulations, quantum-mechanics/molecular-mechanics (QM/MM) calculations, and site-directed mutagenesis. Two Pchlide binding modes were considered, one aligning with historical proposals (mode A) and another consistent with the recent experimental data (mode B). Binding energy calculations revealed that in contrast to the non-specific interactions found for mode A, mode B exhibits distinct stabilizing interactions that support more thermodynamically favorable binding. A comprehensive analysis incorporating QM/MM-based local energy decomposition unraveled a complex interaction network involving Y177, H319, and the C13 1 carboxy group, influencing the pigment's excited state energy and potentially contributing to substrate specificity. Importantly, our results uniformly favor mode B, challenging established interpretations and emphasizing the need for a comprehensive re-evaluation of the LPOR reaction mechanism in a way that incorporates accurate structural information on pigment interactions and substrate-cofactor positioning in the binding pocket. The results shed light on the intricacies of LPOR's catalytic mechanism and provide a solid foundation for further elucidating the secrets of chlorophyll biosynthesis., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Structural and spectroscopic investigations of compound II in ascorbate peroxidase (APX) have yielded conflicting conclusions regarding the protonation state of the crucial Fe(IV) intermediate. Neutron diffraction and crystallographic data support an iron(IV)-hydroxo formulation, whereas Mössbauer, X-ray absorption (XAS), and nuclear resonance vibrational spectroscopy (NRVS) studies appear consistent with an iron(IV)-oxo species. Here we examine APX with spectroscopy-oriented QM/MM calculations and extensive exploration of the conformational space for both possible formulations of compound II. We establish that irrespective of variations in the orientation of a vicinal arginine residue and potential reorganization of proximal water molecules and hydrogen bonding, the Fe-O distances for the oxo and hydroxo forms consistently fall within distinct, narrow, and nonoverlapping ranges. The accuracy of geometric parameters is validated by coupled-cluster calculations with the domain-based local pair natural orbital approach, DLPNO-CCSD(T). QM/MM calculations of spectroscopic properties are conducted for all structural variants, encompassing Mössbauer, optical, X-ray absorption, and X-ray emission spectroscopies and NRVS. All spectroscopic observations can be assigned uniquely to an Fe(IV)═O form. A terminal hydroxy group cannot be reconciled with the spectroscopic data. Under no conditions can the Fe(IV)═O distance be sufficiently elongated to approach the crystallographically reported Fe-O distance. The latter is consistent only with a hydroxo species, either Fe(IV) or Fe(III). Our findings strongly support the Fe(IV)═O formulation of APX-II and highlight unresolved discrepancies in the nature of samples used across different experimental studies.

Photosystem II (PSII), the principal enzyme of oxygenic photosynthesis, contains two integral light harvesting proteins (CP43 and CP47) that bind chlorophylls and carotenoids. The two intrinsic antennae play crucial roles in excitation energy transfer and photoprotection. CP43 interacts most closely with the reaction center of PSII, specifically with the branch of the reaction center (D1) that is responsible for primary charge separation and electron transfer. Deciphering the function of CP43 requires detailed atomic-level insights into the properties of the embedded pigments. To advance this goal, we employ a range of multiscale computational approaches to determine the site energies and excitonic profile of CP43 chlorophylls, using large all-atom models of a membrane-bound PSII monomer. In addition to time-dependent density functional theory (TD-DFT) used in the context of a quantum-mechanics/molecular-mechanics setup (QM/MM), we present a thorough analysis using the perturbed matrix method (PMM), which enables us to utilize information from long-timescale molecular dynamics simulations of native PSII-complexed CP43. The excited state energetics and excitonic couplings have both similarities and differences compared with previous experimental fits and theoretical calculations. Both static TD-DFT and dynamic PMM results indicate a layered distribution of site energies and reveal specific groups of chlorophylls that have shared contributions to low-energy excitations. Importantly, the contribution to the lowest energy exciton does not arise from the same chlorophylls at each system configuration, but rather changes as a function of conformational dynamics. An unexpected finding is the identification of a low-energy charge-transfer excited state within CP43 that involves a lumenal (C2) and the central (C10) chlorophyll of the complex. The results provide a refined basis for structure-based interpretation of spectroscopic observations and for further deciphering excitation energy transfer in oxygenic photosynthesis., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Two Mo(0) phosphenium complexes containing ancillary secondary phosphine ligands have been investigated with respect to their ability to participate in electrophilic addition at unsaturated substrates and subsequent P-H hydride transfer to "quench" the resulting carbocations. These studies provide stoichiometric "proof of concept" for a proposed new metal-catalyzed electrophilic hydrophosphination mechanism. The more strongly Lewis acidic phosphenium complex, [Mo(CO) 4 (PR 2 H)(PR 2 )] + (R=Ph, Tol p ), cleanly hydrophosphinates 1,1-diphenylethylene, benzophenone, and ethylene, while other substrates react rapidly to give products resulting from competing electrophilic processes. A less Lewis acidic complex, [Mo(CO) 3 (PR 2 H) 2 (PR 2 )] + , generally reacts more slowly but participates in clean hydrophosphination of a wider range of unsaturated substrates, including styrene, indene, 1-hexene, and cyclohexanone, in addition to 1,1-diphenylethylene, benzophenone, and ethylene. Mechanistic studies are described, including stoichiometric control reactions and computational and kinetic analyses, which probe whether the observed P-H addition actually does occur by the proposed electrophilic mechanism, and whether hydridic P-H transfer in this system is intra- or intermolecular. Preliminary reactivity studies indicate challenges that must be addressed to exploit these promising results in catalysis., (© 2024 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

The identity and insertion pathway of the substrate oxygen atoms that are coupled to dioxygen by the oxygen-evolving complex (OEC) remains a central question toward understanding Nature's water oxidation mechanism. In several studies, ammonia has been used as a small "water analogue" to elucidate the pathway of substrate access to the OEC and to aid in determining which of the oxygen ligands of the tetramanganese cluster are substrates for O-O bond formation. On the basis of structural and spectroscopic investigations, five first-sphere binding modes of ammonia have been suggested, involving either substitution of an existing H 2 O/OH - /O 2- group or addition as an extra ligand to a metal ion of the Mn 4 CaO 5 cluster. Some of these modes, specifically the ones involving substitution, have already been subject to spectroscopy-oriented quantum chemical investigations, whereas more recent suggestions that postulate the addition of ammonia have not been examined so far with quantum chemistry for their agreement with spectroscopic data. Herein, we use a common structural framework and theoretical methodology to evaluate structural models of the OEC that represent all proposed modes of first-sphere ammonia interaction with the OEC in its S 2 state. Criteria include energetic, magnetic, kinetic, and spectroscopic properties compared against available experimental EPR, ENDOR, ESEEM, and EDNMR data. Our results show that models featuring ammonia replacing one of the two terminal water ligands on Mn4 align best with experimental data, while they definitively exclude substitution of a bridging μ-oxo ligand as well as incorporation of ammonia as a sixth ligand on Mn1 or Mn4.

DNA G-quadruplexes (GQs) are of great interest due to their involvement in crucial biological processes such as telomerase maintenance and gene expression. Furthermore, they are reported as catalytically active DNAzymes and building blocks in bio-nanotechnology. GQs exhibit remarkable structural diversity and conformational heterogeneity, necessitating precise and reliable tools to unravel their structure-function relationships. Here, we present insights into the structure and conformational flexibility of a unimolecular GQ with high spatial resolution via electron-nuclear double resonance (ENDOR) experiments combined with Cu(II) and fluorine labeling. These findings showcase the successful application of the 19 F-ENDOR methodology at 34 GHz, overcoming the limitations posed by the complexity and scarcity of higher-frequency spectrometers. Importantly, our approach retains both sensitivity and orientational resolution. This integrated study not only enhances our understanding of GQs but also expands the methodological toolbox for studying other macromolecules., (© 2023 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

We selected 145 reference organic molecules that include model fragments used in computer-aided drug design. We calculated 158 conformational energies and barriers using force fields, with wide applicability in commercial and free softwares and extensive application on the calculation of conformational energies of organic molecules, e.g. the UFF and DREIDING force fields, the Allinger's force fields MM3-96, MM3-00, MM4-8, the MM2-91 clones MMX and MM+, the MMFF94 force field, MM4, ab initio Hartree-Fock (HF) theory with different basis sets, the standard density functional theory B3LYP, the second-order post-HF MP2 theory and the Domain-based Local Pair Natural Orbital Coupled Cluster DLPNO-CCSD(T) theory, with the latter used for accurate reference values. The data set of the organic molecules includes hydrocarbons, haloalkanes, conjugated compounds, and oxygen-, nitrogen-, phosphorus- and sulphur-containing compounds. We reviewed in detail the conformational aspects of these model organic molecules providing the current understanding of the steric and electronic factors that determine the stability of low energy conformers and the literature including previous experimental observations and calculated findings. While progress on the computer hardware allows the calculations of thousands of conformations for later use in drug design projects, this study is an update from previous classical studies that used, as reference values, experimental ones using a variety of methods and different environments. The lowest mean error against the DLPNO-CCSD(T) reference was calculated for MP2 (0.35 kcal mol -1 ), followed by B3LYP (0.69 kcal mol -1 ) and the HF theories (0.81-1.0 kcal mol -1 ). As regards the force fields, the lowest errors were observed for the Allinger's force fields MM3-00 (1.28 kcal mol -1 ), ΜΜ3-96 (1.40 kcal mol -1 ) and the Halgren's MMFF94 force field (1.30 kcal mol -1 ) and then for the MM2-91 clones MMX (1.77 kcal mol -1 ) and MM+ (2.01 kcal mol -1 ) and MM4 (2.05 kcal mol -1 ). The DREIDING (3.63 kcal mol -1 ) and UFF (3.77 kcal mol -1 ) force fields have the lowest performance. These model organic molecules we used are often present as fragments in drug-like molecules. The values calculated using DLPNO-CCSD(T) make up a valuable data set for further comparisons and for improved force field parameterization., (© 2023. The Author(s).)

Photosystem II, the water splitting enzyme of photosynthesis, utilizes the energy of sunlight to drive the four-electron oxidation of water to dioxygen at the oxygen-evolving complex (OEC). The OEC harbors a Mn 4 CaO 5 cluster that cycles through five oxidation states S i ( i = 0-4). The S 3 state is the last metastable state before the O 2 evolution. Its electronic structure and nature of the S 2 → S 3 transition are key topics of persisting controversy. Most spectroscopic studies suggest that the S 3 state consists of four Mn(IV) ions, compared to the Mn(III)Mn(IV) 3 of the S 2 state. However, recent crystallographic data have received conflicting interpretations, suggesting either metal- or ligand-based oxidation, the latter leading to an oxyl radical or a peroxo moiety in the S 3 state. Herein, we utilize high-energy resolution fluorescence detected (HERFD) X-ray absorption spectroscopy to obtain a highly resolved description of the Mn K pre-edge region for all S-states, paying special attention to use chemically unperturbed S 3 state samples. In combination with quantum chemical calculations, we achieve assignment of specific spectroscopic features to geometric and electronic structures for all S-states. These data are used to confidently discriminate between the various suggestions concerning the electronic structure and the nature of oxidation events in all observable catalytic intermediates of the OEC. Our results do not support the presence of either peroxo or oxyl in the active configuration of the S 3 state. This establishes Mn-centered storage of oxidative equivalents in all observable catalytic transitions and constrains the onset of the O-O bond formation until after the final light-driven oxidation event.

Oxygenic photosynthesis is the fundamental energy-converting process that utilizes sunlight to generate molecular oxygen and the organic compounds that sustain life. Protein-pigment complexes harvest light and transfer excitation energy to specialized pigment assemblies, reaction centers (RC), where electron transfer cascades are initiated. A molecular-level understanding of the primary events is indispensable for elucidating the principles of natural photosynthesis and enabling development of bioinspired technologies. The primary enzyme in oxygenic photosynthesis is Photosystem II (PSII), a membrane-embedded multisubunit complex, that catalyzes the light-driven oxidation of water. The RC of PSII consists of four chlorophyll a and two pheophytin a pigments symmetrically arranged along two core polypeptides; only one branch participates in electron transfer. Despite decades of research, fundamental questions remain, including the origin of this functional asymmetry, the nature of primary charge-transfer states and the identity of the initial electron donor, the origin of the capability of PSII to enact charge separation with far-red photons, i.e., beyond the "red limit" where individual chlorophylls absorb, and the role of protein conformational dynamics in modulating charge-separation pathways.In this Account, we highlight developments in quantum-chemistry based excited-state computations for multipigment assemblies and the refinement of protocols for computing protein-induced electrochromic shifts and charge-transfer excitations calibrated with modern local correlation coupled cluster methods. We emphasize the importance of multiscale atomistic quantum-mechanics/molecular-mechanics and large-scale molecular dynamics simulations, which enabled direct and accurate modeling of primary processes in RC excitation at the quantum mechanical level.Our findings show how differential protein electrostatics enable spectral tuning of RC pigments and generate functional asymmetry in PSII. A chlorophyll pigment on the active branch (Chl D1 ) has the lowest site energy in PSII and is the primary electron donor. The complete absence of low-lying charge-transfer states within the central pair of chlorophylls excludes a long-held assumption about the initial charge separation. Instead, we identify two primary charge separation pathways, both with the same pheophytin acceptor (Pheo D1 ): a fast pathway with Chl D1 as the primary electron donor (short-range charge-separation) and a slow pathway with P D1 P D2 as the initial donor (long-range charge separation). The low-energy spectrum is dominated by two states with significant charge-transfer character, Chl D1 δ+ Pheo D1 δ- and P D1 δ+ Pheo D1 δ- . The conformational dynamics of PSII allows these charge-transfer states to span wide energy ranges, pushing oxygenic photosynthesis beyond the "red limit". These results provide a quantum mechanical picture of the primary events in the RC of oxygenic photosynthesis, forming a solid basis for interpreting experimental observations and for extending photosynthesis research in new directions.

In nature, nonheme iron enzymes use dioxygen to generate high-spin iron(IV)=O species for a variety of oxygenation reactions. Although synthetic chemists have long sought to mimic this reactivity, the enzyme-like activation of O 2 to form high-spin iron(IV) = O species remains an unrealized goal. Here, we report a metal-organic framework featuring iron(II) sites with a local structure similar to that in α-ketoglutarate-dependent dioxygenases. The framework reacts with O 2 at low temperatures to form high-spin iron(IV) = O species that are characterized using in situ diffuse reflectance infrared Fourier transform, in situ and variable-field Mössbauer, Fe Kβ x-ray emission, and nuclear resonance vibrational spectroscopies. In the presence of O 2 , the framework is competent for catalytic oxygenation of cyclohexane and the stoichiometric conversion of ethane to ethanol.

In oxygenic photosynthesis sunlight is harvested and funneled as excitation energy into the reaction center (RC) of Photosystem II (PSII), the site of primary charge separation that initiates the photosynthetic electron transfer chain. The chlorophyll Chl D1 pigment of the RC is the primary electron donor, forming a charge-separated radical pair with the vicinal pheophytin Pheo D1 (Chl D1 + Pheo D1 - ). To avert charge recombination, the electron is further transferred to plastoquinone Q A , whereas the hole relaxes to a central pair of chlorophylls (P D1 P D2 ), subsequently driving water oxidation. Spin-triplet states can form within the RC when forward electron transfer is inhibited or back reactions are favored. This can lead to formation of singlet dioxygen, with potential deleterious effects. Here we investigate the nature and properties of triplet states within the PSII RC using a multiscale quantum-mechanics/molecular-mechanics (QM/MM) approach. The low-energy spectrum of excited singlet and triplet states, of both local and charge-transfer nature, is compared using range-separated time-dependent density functional theory (TD-DFT). We further compute electron paramagnetic resonance properties (zero-field splitting parameters and hyperfine coupling constants) of relaxed triplet states and compare them with available experimental data. Moreover, the electrostatic modulation of excited state energetics and redox properties of RC pigments by the semiquinone Q A - is described. The results provide a detailed electronic-level understanding of triplet states within the PSII RC and form a refined basis for discussing primary and secondary electron transfer, charge recombination pathways, and possible photoprotection mechanisms in PSII., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Recent advances in serial femtosecond crystallography (SFX) of photosystem II (PSII), enabled by X-ray free electron lasers (XFEL), provided the first geometric models of distinct intermediates in the catalytic S-state cycle of the oxygen-evolving complex (OEC). These models are obtained by flash-advancing the OEC from the dark-stable state (S 1 ) to more oxidized intermediates (S 2 and S 3 ), eventually cycling back to the most reduced S 0 . However, the interpretation of these models is controversial because geometric parameters within the Mn 4 CaO 5 cluster of the OEC do not exactly match those expected from coordination chemistry for the spectroscopically verified manganese oxidation states of the distinct S-state intermediates. Here we focus on the first catalytic transition, S 1 → S 2 , which represents a one-electron oxidation of the OEC. Combining geometric and electronic structure criteria, including a novel effective oxidation state approach, we analyze existing 1-flash (1F) SFX-XFEL crystallographic models that should depict the S 2 state of the OEC. We show that the 1F/S 2 equivalence is not obvious, because the Mn oxidation states and total unpaired electron counts encoded in these models are not fully consistent with those of a pure S 2 state and with the nature of the S 1 → S 2 transition. Furthermore, the oxidation state definition in two-flashed (2F) structural models is practically impossible to elucidate. Our results advise caution in the extraction of electronic structure information solely from the literal interpretation of crystallographic models and call for re-evaluation of structural and mechanistic interpretations that presume exact correspondence of such models to specific catalytic intermediates of the OEC.

Photosystem-II (PSII) is a multi-subunit protein complex that harvests sunlight to perform oxygenic photosynthesis. Initial light-activated charge separation takes place at a reaction centre consisting of four chlorophylls and two pheophytins. Understanding the processes following light excitation remains elusive due to spectral congestion, the ultrafast nature, and multi-component behaviour of the charge-separation process. Here, using advanced computational multiscale approaches which take into account the large-scale configurational flexibility of the system, we identify two possible primary pathways to radical-pair formation that differ by three orders of magnitude in their kinetics. The fast (short-range) pathway is dominant, but the existence of an alternative slow (long-range) charge-separation pathway hints at the evolution of redundancy that may serve other purposes, adaptive or protective, related to formation of the unique oxidative species that drives water oxidation in PSII., (© 2023 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Reactions of methyl-2-pyridyl ketone, pyCOMe, with FeCl 3 ·6H 2 O in various solvents gave complexes [Fe 4 Cl 6 (OMe) 2 (L1) 2 ]·0.7MeCN·0.4MeOH (1·0.7MeCN·0.4MeOH) and [Fe 3 Cl 4 (bicine)(L2)]·Me 2 CO·0.2H 2 O (2·Me 2 CO·0.2H 2 O). The ligands (L1) 2- = pyCO(Me)CHCOpy (in 1) and (L2) 2- = pyCO(Me)CH 2 CO(OMe)py (in 2) are formed in situ , through an aldol reaction-type mechanism between the carbanion pyC(O)CH 2 - (formed by the nucleophilic attack of the MeO - in pyCOMe) and pyCOMe which results in the formation of a new C-C bond. The intermediate compound undergoes attack in the -CH 2 - or -CO- group by a MeO - group, and the new ligands (L1) 2- and (L2) 2- , respectively, are formed. The molecular structure of 1 consists of three corner-sharing [Fe 2 O 2 ] rhombic units in cis -arrangement. The two terminal Fe III ions display distorted square pyramidal geometry and the two central Fe III ions are distorted octahedral. The molecular structure of 2 consists of two corner-sharing [Fe 2 O 2 ] rhombic units, with the two terminal Fe III ions in distorted square pyramidal geometry and the central Fe III in distorted octahedral. The differentiation in the coordination environment of the Fe III ions in 1-2 is reflected in the values of the Mössbauer hyperfine parameters. In agreement with theoretical calculations, the square pyramidal sites exhibit a smaller isomer shift value in comparison to the octahedral sites. Magnetic studies indicate antiferromagnetic interactions leading to an S = 0 ground state in 1 and to an S = 5/2 ground state in 2, consistent with Electron Paramagnetic Resonance spectroscopy. Mössbauer spectra of 2 indicate the onset of relaxation effects below 80 K. At 1.5 K the spectrum of 2 consists of magnetic sextets. The determined hyperfine magnetic fields are consistent with the exchange coupling scheme imposed by the crystal structure of 2. Theoretical calculations shed light on the differences in the electronic structure between the square pyramidal and the octahedral sites.

Homometallic copper complexes with alkenylidene ligands are discussed as intermediates in catalysis but the isolation of such complexes has remained elusive. Herein, we report the structural characterization of copper complexes with bridging and terminal alkenylidene ligands. The compounds were obtained by irradiation of Cu I complexes with N-heterocyclic diazoolefin ligands. The complex with a terminal alkenylidene ligand required isolation in a crystalline matrix, and its structural characterization was enabled by in crystallo photolysis at low temperature., (© 2022 Wiley-VCH GmbH.)

Oxalate decarboxylase from Bacillus subtilis is a binuclear Mn-dependent acid stress response enzyme that converts the mono-anion of oxalic acid into formate and carbon dioxide in a redox neutral unimolecular disproportionation reaction. A π-stacked tryptophan dimer, W96 and W274, at the interface between two monomer subunits facilitates long-range electron transfer between the two Mn ions and plays an important role in the catalytic mechanism. Substitution of W96 with the unnatural amino acid 5-hydroxytryptophan leads to a persistent EPR signal which can be traced back to the neutral radical of 5-hydroxytryptophan with its hydroxyl proton removed. 5-Hydroxytryptophan acts as a hole sink preventing the formation of Mn(III) at the N-terminal active site and strongly suppresses enzymatic activity. The lower boundary of the standard reduction potential for the active site Mn(II)/Mn(III) couple can therefore be estimated as 740 mV against the normal hydrogen electrode at pH 4, the pH of maximum catalytic efficiency. Our results support the catalytic importance of long-range electron transfer in oxalate decarboxylase while at the same time highlighting the utility of unnatural amino acid incorporation and specifically the use of 5-hydroxytryptophan as an energetic sink for hole hopping to probe electron transfer in redox proteins., (© 2022 The Protein Society.)

The understanding of light-induced biological water oxidation in oxygenic photosynthesis is of great importance both for biology and (bio)technological applications. The chemically difficult multistep reaction takes place at a unique protein-bound tetra-manganese/calcium cluster in photosystem II whose structure has been elucidated by X-ray crystallography (Umena et al. Nature 2011, 473, 55). The cluster moves through several intermediate states in the catalytic cycle. A detailed understanding of these intermediates requires information about the spatial and electronic structure of the Mn 4 Ca complex; the latter is only available from spectroscopic techniques. Here, the important role of Electron Paramagnetic Resonance (EPR) and related double resonance techniques (ENDOR, EDNMR), complemented by quantum chemical calculations, is described. This has led to the elucidation of the cluster's redox and protonation states, the valence and spin states of the manganese ions and the interactions between them, and contributed substantially to the understanding of the role of the protein surrounding, as well as the binding and processing of the substrate water molecules, the O-O bond formation and dioxygen release. Based on these data, models for the water oxidation cycle are developed., (© 2022 The Authors. FEBS Letters published by John Wiley & Sons Ltd on behalf of Federation of European Biochemical Societies.)

Water channels and networks within photosystem II (PSII) of oxygenic photosynthesis are critical for enzyme structure and function. They control substrate delivery to the oxygen-evolving center and mediate proton transfer at both the oxidative and reductive endpoints. Current views on PSII hydration are derived from protein crystallography, but structural information may be compromised by sample dehydration and technical limitations. Here, we simulate the physiological hydration structure of a cyanobacterial PSII model following a thorough hydration procedure and large-scale unconstrained all-atom molecular dynamics enabled by massively parallel simulations. We show that crystallographic models of PSII are moderately to severely dehydrated and that this problem is particularly acute for models derived from X-ray free electron laser (XFEL) serial femtosecond crystallography. We present a fully hydrated representation of cyanobacterial PSII and map all water channels, both static and dynamic, associated with the electron donor and acceptor sides. Among them, we describe a series of transient channels and the attendant conformational gating role of protein components. On the acceptor side, we characterize a channel system that is absent from existing crystallographic models but is likely functionally important for the reduction of the terminal electron acceptor plastoquinone Q B . The results of the present work build a foundation for properly (re)evaluating crystallographic models and for eliciting new insights into PSII structure and function.

In this article, we describe a combined experimental and theoretical mechanistic investigation of the C(sp 2 )-F bond formation from neutral and cationic high-valent organobismuth(V) fluorides, featuring a dianionic bis-aryl sulfoximine ligand. An exhaustive assessment of the substitution pattern in the ligand, the sulfoximine, and the reactive aryl on neutral triarylbismuth(V) difluorides revealed that formation of dimeric structures in solution promotes facile Ar-F bond formation. Noteworthy, theoretical modeling of reductive elimination from neutral bismuth(V) difluorides agrees with the experimentally determined kinetic and thermodynamic parameters. Moreover, the addition of external fluoride sources leads to inactive octahedral anionic Bi(V) trifluoride salts, which decelerate reductive elimination. On the other hand, a parallel analysis for cationic bismuthonium fluorides revealed the crucial role of tetrafluoroborate anion as fluoride source. Both experimental and theoretical analyses conclude that C-F bond formation occurs through a low-energy five-membered transition-state pathway, where the F anion is delivered to a C(sp 2 ) center, from a BF 4 anion, reminiscent of the Balz-Schiemann reaction. The knowledge gathered throughout the investigation permitted a rational assessment of the key parameters of several ligands, identifying the simple sulfone-based ligand family as an improved system for the stoichiometric and catalytic fluorination of arylboronic acid derivatives.

Spin-state energetics of transition metal complexes remain one of the most challenging targets for electronic structure methods. Among single-reference wave function approaches, local correlation approximations to coupled cluster theory, most notably the domain-based local pair natural orbital (DLPNO) approach, hold the promise of bringing the accuracy of coupled cluster theory with single, double, and perturbative triple excitations, CCSD(T), to molecular systems of realistic size with acceptable computational cost. However, recent studies on spin-state energetics of iron-containing systems raised doubts about the ability of the DLPNO approach to adequately and systematically approximate energetics obtained by the reference-quality complete active space second-order perturbation theory with coupled-cluster semicore correlation, CASPT2/CC. Here, we revisit this problem using a diverse set of iron complexes and examine several aspects of the application of the DLPNO approach. We show that DLPNO-CCSD(T) can accurately reproduce both CASPT2/CC and canonical CCSD(T) results if two basic principles are followed. These include the consistent use of the improved iterative (T 1 ) versus the semicanonical perturbative triple corrections and, most importantly, a simple two-point extrapolation to the PNO space limit. The latter practically eliminates errors arising from the default truncation of electron-pair correlation spaces and should be viewed as standard practice in applications of the method to transition metal spin-state energetics. Our results show that reference-quality results can be readily achieved with DLPNO-CCSD(T) if these principles are followed. This is important also in view of the applicability of the method to larger single-reference systems and multinuclear clusters, whose treatment of dynamic correlation would be challenging for multireference-based approaches.

Understanding the structure and function of lytic polysaccharide monooxygenases (LPMOs), copper enzymes that degrade recalcitrant polysaccharides, requires the reliable atomistic interpretation of electron paramagnetic resonance (EPR) data on the Cu(II) active site. Among various LPMO families, the chitin-active Pl AA10 shows an intriguing phenomenology with distinct EPR signals, a major rhombic and a minor axial signal. Here, we combine experimental and computational investigations to uncover the structural identity of these signals. X-band EPR spectra recorded at different pH values demonstrate pH-dependent population inversion: the major rhombic signal at pH 6.5 becomes minor at pH 8.5, where the axial signal dominates. This suggests that a protonation change is involved in the interconversion. Precise structural interpretations are pursued with quantum chemical calculations. Given that accurate calculations of Cu g -tensors remain challenging for quantum chemistry, we first address this problem via a thorough calibration study. This enables us to define a density functional that achieves accurate and reliable prediction of g -tensors, giving confidence in our evaluation of Pl AA10 LPMO models. Large models were considered that include all parts of the protein matrix surrounding the Cu site, along with the characteristic second-sphere features of Pl AA10. The results uniquely identify the rhombic signal with a five-coordinate Cu ion bearing two water molecules in addition to three N-donor ligands. The axial signal is attributed to a four-coordinate Cu ion where only one of the waters remains bound, as hydroxy. Alternatives that involve decoordination of the histidine brace amino group are unlikely based on energetics and spectroscopy. These results provide a reliable spectroscopy-consistent view on the plasticity of the resting state in Pl AA10 LPMO as a foundation for further elucidating structure-property relationships and the formation of catalytically competent species. Our strategy is generally applicable to the study of EPR parameters of mononuclear copper-containing metalloenzymes.

An isostructural series of heavy Group 14 E(I) radical anions (Ge, Sn, Pb), stabilized by a bulky xanthene-based diamido ligand are reported. The radical anions were synthesised by the one-electron reduction of their corresponding E(II) precursor complexes with sodium naphthalenide in THF, yielding the radical anions as charge-separated sodium salts. The series of main group radicals have been comprehensively characterized by EPR spectroscopy, X-ray crystallography and DFT analysis, which reveal that in all cases, the spin density of the unpaired electron almost exclusively resides in a p-orbital of π symmetry located on the Group 14 center., (© 2022 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Photosystem-II uses sunlight to trigger charge separation and catalyze water oxidation. Intrinsic properties of chlorophyll a pigments define a natural "red limit" of photosynthesis at ≈680 nm. Nevertheless, charge separation can be triggered with far-red photons up to 800 nm, without altering the nature of light-harvesting pigments. Here we identify the electronic origin of this remarkable phenomenon using quantum chemical and multiscale simulations on a native Photosystem-II model. We find that the reaction center is preorganized for charge separation in the far-red region by specific chlorophyll-pheophytin pairs, potentially bypassing the light-harvesting apparatus. Charge transfer can occur along two distinct pathways with one and the same pheophytin acceptor (Pheo D1 ). The identity of the donor chlorophyll (Chl D1 or P D1 ) is wavelength-dependent and conformational dynamics broaden the sampling of the far-red region by the two charge-transfer states. The two pathways rationalize spectroscopic observations and underpin designed extensions of the photosynthetically active radiation limit., (© 2022 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Hydrated transition metal ions are prototypical systems that can be used to model properties of transition metals in complex chemical environments. These seemingly simple systems present challenges for computational chemistry and are thus crucial in evaluations of quantum chemical methods for spin-state and redox energetics. In this work, we explore the applicability of the domain-based pair natural orbital implementation of coupled cluster (DLPNO-CC) theory to the calculation of ionization energies and redox potentials for hydrated ions of all first transition row (3d) metals in the 2+/3+ oxidation states, in connection with various solvation approaches. In terms of model definition, we investigate the construction of a minimally explicitly hydrated quantum cluster with a first and second hydration layer. We report on the convergence with respect to the coupled cluster expansion and the PNO space, as well as on the role of perturbative triple excitations. A recent implementation of the conductor-like polarizable continuum model (CPCM) for the DLPNO-CC approach is employed to determine self-consistent redox potentials at the coupled cluster level. Our results establish conditions for the convergence of DLPNO-CCSD(T) energetics and stress the absolute necessity to explicitly consider the second solvation sphere even when CPCM is used. The achievable accuracy for redox potentials of a practical DLPNO-based approach is, on average, 0.13 V. Furthermore, multilayer approaches that combine a higher-level DLPNO-CCSD(T) description of the first solvation sphere with a lower-level description of the second solvation layer are investigated. The present work establishes optimal and transferable methodological choices for employing DLPNO-based coupled cluster theory, the associated CPCM implementation, and cost-efficient multilayer derivatives of the approach for open-shell transition metal systems in complex environments.

Singlet vinylidenes (R 2 C═C:) are proposed as intermediates in a series of organic reactions, and very few have been studied by matrix isolation or gas-phase spectroscopy. Triplet vinylidenes, however, featuring two unpaired electrons at a monosubstituted carbon atom are thus far only predicted as electronically excited-state species and represent an unexplored class of carbon-centered diradicals. We report the photochemical generation and low-temperature EPR/ENDOR characterization of the first ground-state high-spin (triplet) vinylidene. The zero-field splitting parameters ( D = 0.377 cm -1 and | E| / D = 0.028) were determined, and the 13 C hyperfine coupling tensor was obtained by 13 C-ENDOR measurements. Most strikingly, the isotropic 13 C hyperfine coupling constant (50 MHz) is far smaller than the characteristic values of triplet carbenes, demonstrating a unique electronic structure which is supported by quantum chemical calculations.

Photosynthetic pigment-protein complexes harvest solar energy with a high quantum efficiency. Protein scaffolds are known to tune the spectral properties of embedded pigments principally through structured electrostatic environments. Although the physical nature of electrostatic tuning is straightforward, the precise spatial principles of electrostatic preorganization remain poorly explored for different protein matrices and incompletely characterized with respect to the intrinsic properties of different photosynthetic pigments. In this work, we study the electronic structure features associated with the lowest excited state of a series of eight naturally occurring (bacterio)chlorophylls and pheophytins to describe the precise topological differences in electrostatic potentials and hence determine intrinsic differences in the expected mode and impact of electrostatic tuning. The difference electrostatic potentials between the ground and first excited states are used as fingerprints. Both the spatial profile and the propensity for spectral tuning are found to be unique for each pigment, indicating spatially and directionally distinct modes of electrostatic tuning. The results define a specific partitioning of the protein matrix around each pigment as an aid to identify regions with a maximal impact on spectral tuning and have direct implications for dimensionality reduction in protein design and engineering. Thus, a quantum mechanical basis is provided for understanding, predicting, and ultimately designing sequence-modified or pigment-exchanged biological systems, as suggested for selected examples of pigment-reconstituted proteins.

The electronic and geometric structures of the water-oxidizing complex of photosystem II in the steps of the catalytic cycle that precede dioxygen evolution remain hotly debated. Recent structural and spectroscopic investigations support contradictory redox formulations for the active-site Mn 4 CaO x cofactor in the final metastable S 3 state. These range from the widely accepted Mn IV 4 oxo-hydroxo model, which presumes that O-O bond formation occurs in the ultimate transient intermediate (S 4 ) of the catalytic cycle, to a Mn III 2 Mn IV 2 peroxo model representative of the contrasting "early-onset" O-O bond formation hypothesis. Density functional theory energetics of suggested S 3 redox isomers are inconclusive because of extreme functional dependence. Here, we use the power of the domain-based local pair natural orbital approach to coupled cluster theory, DLPNO-CCSD(T), to present the first correlated wave function theory calculations of relative stabilities for distinct redox-isomeric forms of the S 3 state. Our results enabled us to evaluate conflicting models for the S 3 state of the oxygen-evolving complex (OEC) and to quantify the accuracy of lower-level theoretical approaches. Our assessment of the relevance of distinct redox-isomeric forms for the mechanism of biological water oxidation strongly disfavors the scenario of early-onset O-O formation advanced by literal interpretations of certain crystallographic models. This work serves as a case study in the application of modern coupled cluster implementations to redox isomerism problems in oligonuclear transition metal systems., (© 2021 The Authors. Chemistry - A European Journal published by Wiley-VCH GmbH.)

The tetramanganese-calcium cluster of the oxygen-evolving complex of photosystem II adopts electronically and magnetically distinct but interconvertible valence isomeric forms in its first light-driven oxidized catalytic state, S 2 . This bistability is implicated in gating the final catalytic states preceding O-O bond formation, but it is unknown how the biological system enables its emergence and controls its effect. Here we show that the Mn 4 CaO 5 cluster in the resting (dark-stable) S 1 state adopts orientational Jahn-Teller isomeric forms arising from a directional change in electronic configuration of the "dangler" Mn III ion. The isomers are consistent with available structural data and explain previously unresolved electron paramagnetic resonance spectroscopic observations on the S 1 state. This unique isomerism in the resting state is shown to be the electronic origin of valence isomerism in the S 2 state, establishing a functional role of orientational Jahn-Teller isomerism unprecedented in biological or artificial catalysis., (© 2021 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

Most functional groups, especially those consisting of the abundant elements of organic matter-carbon, nitrogen and oxygen-have been extensively studied and only very few remain speculative due to their high intrinsic reactivity. In contrast to the well-explored chemistry of diazoalkanes (R 2 C=N 2 ), diazoalkenes (R 2 C=C=N 2 ) have been postulated in several organic transformations, but remain elusive long-sought intermediates. Here, we present a room-temperature stable diazoalkene, utilizing a dinitrogen transfer from nitrous oxide. This functional group shows dual-site nucleophilicity (C and N atoms) and features a bent C-C-N entity (124°) and a long N-N bond together with a remarkable low infrared absorption (1,944 cm -1 ). Substitution of N 2 by an isocyanide leads to a vinylidene ketenimine. Furthermore, photochemically triggered loss of dinitrogen might proceed through a transient triplet vinylidene. We anticipate the existence of a stable diazoalkene functional group to pave an exciting avenue into the chemistry of low-valent carbon and unsaturated carbenes.

Blue copper proteins continue to challenge experiment and theory with their electronic structure and spectroscopic properties that respond sensitively to the coordination environment of the copper ion. In this work, we report state-of-the art electronic structure studies for geometric and spectroscopic properties of the archetypal "Type I" copper protein azurin in its Cu(II) state. A hybrid quantum mechanics/molecular mechanics (QM/MM) approach is used, employing both density functional theory (DFT) and coupled cluster with singles, doubles, and perturbative triples (CCSD(T)) methods for the QM region, the latter method making use of the domain-based local pair natural orbital (DLPNO) approach. Models of increasing QM size are employed to investigate the convergence of critical geometric parameters. It is shown that convergence is slow and that a large QM region is critical for reproducing the short experimental Cu-SCys112 distance. The study of structural convergence is followed by investigation of spectroscopic parameters using both DFT and DLPNO-CC methods and comparing these to the experimental spectrum using simulations. The results allow us to examine for the first time the distribution of spin densities and hyperfine coupling constants at the coupled cluster level, leading us to revisit the experimental assignment of the 33 S hyperfine splitting. The wavefunction-based approach to obtain spin-dependent properties of open-shell systems demonstrated here for the case of azurin is transferable and applicable to a large array of bioinorganic systems.

The determination of the diiron core intermediate structures involved in the catalytic cycle of soluble methane monooxygenase (sMMO), the enzyme that selectively catalyzes the conversion of methane to methanol, has been a subject of intense interest within the bioinorganic scientific community. Particularly, the specific geometry and electronic structure of the intermediate that precedes methane binding, known as intermediate Q (or MMOH Q ), has been debated for over 30 years. Some reported studies support a bis-μ-oxo-bridged Fe(IV) 2 O 2 closed-core conformation Fe(IV) 2 O 2 core, whereas others favor an open-core geometry, with a longer Fe-Fe distance. The lack of consensus calls for a thorough re-examination and reinterpretation of the spectroscopic data available on the MMOH Q intermediate. Herein, we report extensive simulations based on a hybrid quantum mechanics/molecular mechanics approach (QM/MM) approach that takes into account the complete enzyme to explore possible conformations for intermediates MMOH ox and MMOH Q of the sMMOH catalytic cycle. High-level quantum chemical approaches are used to correlate specific structural motifs with geometric parameters for comparison with crystallographic and EXAFS data, as well as with spectroscopic data from Mössbauer spectroscopy, Fe K-edge high-energy resolution X-ray absorption spectroscopy (HERFD XAS), and resonance Raman 16 O- 18 O difference spectroscopy. The results provide strong support for an open-core-type configuration in MMOH Q , with the most likely topology involving mono-oxo-bridged Fe ions and alternate terminal Fe-oxo and Fe-hydroxo groups that interact via intramolecular hydrogen bonding. The implications of an open-core intermediate Q on the reaction mechanism of sMMO are discussed.

Quantum chemical approaches today are a powerful tool to study the properties and reactivity of metalloenzymes. In the field of solar fuels research these involve predominantly photosystem II and hydrogenases, which catalyze water oxidation and hydrogen evolution, as well as related biomimetic and bio-inspired models. Theoretical methods are extensively used to better comprehend the nature of catalytic intermediates, establish important structure-function and structure-property correlations, elucidate functional principles, and uncover the catalytic activity of these complex systems by unravelling the key steps of their reaction mechanism. Computations in the field of water oxidation and hydrogen evolution are used as predictive tools to elucidate structures, explain and synthesize complex experimental observations from advanced spectroscopic techniques, rationalize reactivity on the basis of atomistic models and electronic structure, and guide the design of new synthetic targets. This feature article covers recent advances in the application of quantum chemical methods for understanding the nature of catalytic intermediates and the mechanism by which photosystem II and hydrogenases achieve their function, and points at essential questions that remain only partly answered and at challenges that will have to be met by future advances and applications of quantum and computational chemistry.

Protein-embedded chromophores are responsible for light harvesting, excitation energy transfer, and charge separation in photosynthesis. A critical part of the photosynthetic apparatus are reaction centers (RCs), which comprise groups of (bacterio)chlorophyll and (bacterio)pheophytin molecules that transform the excitation energy derived from light absorption into charge separation. The lowest excitation energies of individual pigments (site energies) are key for understanding photosynthetic systems, and form a prime target for quantum chemistry. A major theoretical challenge is to accurately describe the electrochromic (Stark) shifts in site energies produced by the inhomogeneous electric field of the protein matrix. Here, we present large-scale quantum mechanics/molecular mechanics calculations of electrochromic shifts for the RC chromophores of photosystem II (PSII) using various quantum chemical methods evaluated against the domain-based local pair natural orbital (DLPNO) implementation of the similarity-transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD). We show that certain range-separated density functionals (ωΒ97, ωΒ97X-V, ωΒ2PLYP, and LC-BLYP) correctly reproduce RC site energy shifts with time-dependent density functional theory (TD-DFT). The popular CAM-B3LYP functional underestimates the shifts and is not recommended. Global hybrid functionals are too insensitive to the environment and should be avoided, while nonhybrid functionals are strictly nonapplicable. Among the applicable approximate coupled cluster methods, the canonical versions of CC2 and ADC(2) were found to deviate significantly from the reference results both for the description of the lowest excited state and for the electrochromic shifts. By contrast, their spin-component-scaled (SCS) and particularly the scale-opposite-spin (SOS) variants compare well with the reference DLPNO-STEOM-CCSD and the best range-separated DFT methods. The emergence of RC excitation asymmetry is discussed in terms of intrinsic and protein electrostatic potentials. In addition, we evaluate a minimal structural scaffold of PSII, the D1-D2-Cyt B559 RC complex often employed in experimental studies, and show that it would have the same site energy distribution of RC chromophores as the full PSII supercomplex, but only under the unlikely conditions that the core protein organization and cofactor arrangement remain identical to those of the intact enzyme.

Natural photosynthesis relies on light harvesting and excitation energy transfer by specialized pigment-protein complexes. Their structure and the electronic properties of the embedded chromophores define the mechanisms of energy transfer. An important example of a pigment-protein complex is CP47, one of the integral antennae of the oxygen-evolving photosystem II (PSII) that is responsible for efficient excitation energy transfer to the PSII reaction center. The charge-transfer excitation induced among coupled reaction center chromophores resolves into charge separation that initiates the electron transfer cascade driving oxygenic photosynthesis. Mapping the distribution of site energies among the 16 chlorophyll molecules of CP47 is essential for understanding excitation energy transfer and overall antenna function. In this work, we demonstrate a multiscale quantum mechanics/molecular mechanics (QM/MM) approach utilizing full time-dependent density functional theory with modern range-separated functionals to compute for the first time the excitation energies of all CP47 chlorophylls in a complete membrane-embedded cyanobacterial PSII dimer. The results quantify the electrostatic effect of the protein on the site energies of CP47 chlorophylls, providing a high-level quantum chemical excitation profile of CP47 within a complete computational model of "near-native" cyanobacterial PSII. The ranking of site energies and the identity of the most red-shifted chlorophylls (B3, followed by B1) differ from previous hypotheses in the literature and provide an alternative basis for evaluating past approaches and semiempirically fitted sets. Given that a lot of experimental studies on CP47 and other light-harvesting complexes utilize extracted samples, we employ molecular dynamics simulations of isolated CP47 to identify which parts of the polypeptide are most destabilized and which pigments are most perturbed when the antenna complex is extracted from PSII. We demonstrate that large parts of the isolated complex rapidly refold to non-native conformations and that certain pigments (such as chlorophyll B1 and β-carotene h1) are so destabilized that they are probably lost upon extraction of CP47 from PSII. The results suggest that the properties of isolated CP47 are not representative of the native complexed antenna. The insights obtained from CP47 are generalizable, with important implications for the information content of experimental studies on biological light-harvesting antenna systems., Competing Interests: There are no conflicts to declare., (This journal is © The Royal Society of Chemistry.)

Among the intermediate catalytic steps of the water-oxidizing Mn 4 CaO 5 cluster of photosystem II (PSII), the final metastable S 3 state is critically important because it binds one substrate and precedes O 2 evolution. Herein, we combine X- and Q-band EPR experiments on native and methanol-treated PSII of Spinacia oleracea and show that methanol-treated PSII preparations of the S 3 state correspond to a previously uncharacterized high-spin (S=6) species. This is confirmed as a major component also in intact photosynthetic membranes, coexisting with the previously known intermediate-spin conformation (S=3). The high-spin intermediate is assigned to a water-unbound form, with a Mn IV 3 subunit interacting ferromagnetically via anisotropic exchange with a coordinatively unsaturated Mn IV ion. These results resolve and define the structural heterogeneity of the S 3 state, providing constraints on the S 3 to S 4 transition, on substrate identity and delivery pathways, and on the mechanism of O-O bond formation., (© 2020 The Authors. Angewandte Chemie International Edition published by Wiley-VCH GmbH.)

The Mo=PR 2 π* orbital in a Mo phosphenium complex acts as acceptor in a new P III -based Lewis superacid. This Lewis acid (LA) participates in electrophilic Si-H abstraction from E 3 SiH to give a Mo-bound secondary phosphine ligand, Mo-PR 2 H. The resulting Et 3 Si + ion remains associated with the Mo complex, stabilized by η 1 -P-H donation, yet undergoes rapid exchange with an η 1 -Si-H adduct of free silane in solution. The equilibrium between these two adducts presents an opportunity to assess the role of this new LA in catalytic reactions of silanes: is the LA acting as a catalyst or as an initiator? Preliminary results suggest that a cycle including the Mo-bound phosphine-silylium adduct dominates in the catalytic hydrosilylation of acetophenone, relative to a putative cycle involving the silane-silylium adduct or "free" silylium., (© 2020 Wiley-VCH GmbH.)

The reliable prediction of Cu(II) hyperfine coupling constants remains a challenge for quantum chemistry. Until recently only density functional theory (DFT) could target this property for systems of realistic size. However, wave function based methods become increasingly applicable. In the present work, we define a large set of Cu(II) complexes with experimentally known hyperfine coupling constants and use it to investigate the performance of modern quantum chemical methods for the prediction of this challenging spectroscopic parameter. DFT methods are evaluated against orbital-optimized second-order Møller-Plesset (OO-MP2) theory and coupled cluster calculations including singles and doubles excitations, driven by the domain-based local pair natural orbital approach (DLPNO-CCSD). Special attention is paid to the definition of a basis set that converges adequately toward the basis set limit for the given property for all methods considered in this study, and a specifically optimized basis set is proposed for this purpose. The results suggest that wave function based methods can supplant but do not outcompete DFT for the calculation of Cu(II) hyperfine coupling constants. Mainstream hybrid functionals such as B3PW91 remain on average the best choice., (© 2020 The Authors. ChemPhysChem published by Wiley-VCH GmbH.)

Photosystem II (PSII) is a multisubunit pigment-protein complex that uses light-induced charge separation to power oxygenic photosynthesis. Its reaction center chromophores, where the charge transfer cascade is initiated, are arranged symmetrically along the D1 and D2 core polypeptides and comprise four chlorophyll (P D1 , P D2 , Chl D1 , Chl D2 ) and two pheophytin molecules (Pheo D1 and Pheo D2 ). Evolution favored productive electron transfer only via the D1 branch, with the precise nature of primary excitation and the factors that control asymmetric charge transfer remaining under investigation. Here we present a detailed atomistic description for both. We combine large-scale simulations of membrane-embedded PSII with high-level quantum-mechanics/molecular-mechanics (QM/MM) calculations of individual and coupled reaction center chromophores to describe reaction center excited states. We employ both range-separated time-dependent density functional theory and the recently developed domain based local pair natural orbital (DLPNO) implementation of the similarity transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD), the first coupled cluster QM/MM calculations of the reaction center. We find that the protein matrix is exclusively responsible for both transverse (chlorophylls versus pheophytins) and lateral (D1 versus D2 branch) excitation asymmetry, making Chl D1 the chromophore with the lowest site energy. Multipigment calculations show that the protein matrix renders the Chl D1 → Pheo D1 charge-transfer the lowest energy excitation globally within the reaction center, lower than any pigment-centered local excitation. Remarkably, no low-energy charge transfer states are located within the "special pair" P D1 -P D2 , which is therefore excluded as the site of initial charge separation in PSII. Finally, molecular dynamics simulations suggest that modulation of the electrostatic environment due to protein conformational flexibility enables direct excitation of low-lying charge transfer states by far-red light.

The ability to accurately compute low-energy excited states of chlorophylls is critically important for understanding the vital roles they play in light harvesting, energy transfer, and photosynthetic charge separation. The challenge for quantum chemical methods arises both from the intrinsic complexity of the electronic structure problem and, in the case of biological models, from the need to account for protein-pigment interactions. In this work, we report electronic structure calculations of unprecedented accuracy for the low-energy excited states in the Q and B bands of chlorophyll a . This is achieved by using the newly developed domain-based local pair natural orbital (DLPNO) implementation of the similarity transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD) in combination with sufficiently large and flexible basis sets. The results of our DLPNO-STEOM-CCSD calculations are compared with more approximate approaches. The results demonstrate that, in contrast to time-dependent density functional theory, the DLPNO-STEOM-CCSD method provides a balanced performance for both absorption bands. In addition to vertical excitation energies, we have calculated the vibronic spectrum for the Q and B bands through a combination of DLPNO-STEOM-CCSD and ground-state density functional theory frequency calculations. These results serve as a basis for comparison with gas-phase experiments.

Segmented all-electron relativistically contracted (SARC) basis sets are presented for the elements 37 Rb- 54 Xe, for use with the second-order Douglas-Kroll-Hess approach and the zeroth-order regular approximation. The basis sets have a common set of exponents produced with established heuristic procedures, but have contractions optimized individually for each scalar relativistic Hamiltonian. Their compact size and loose segmented contraction, which is in line with the construction of SARC basis sets for heavier elements, makes them suitable for routine calculations on large systems and when core spectroscopic properties are of interest. The basis sets are of triple-zeta quality and come in singly or doubly polarized versions, which are appropriate for both density functional theory and correlated wave function theory calculations. The quality of the basis sets is assessed against large decontracted reference basis sets for a number of atomic and ionic properties, while their general applicability is demonstrated with selected molecular examples., (© 2020 The Authors. Journal of Computational Chemistry published by Wiley Periodicals, Inc.)

The investigation of water oxidation in photosynthesis has remained a central topic in biochemical research for the last few decades due to the importance of this catalytic process for technological applications. Significant progress has been made following the 2011 report of a high-resolution X-ray crystallographic structure resolving the site of catalysis, a protein-bound Mn 4 CaO x complex, which passes through ≥5 intermediate states in the water-splitting cycle. Spectroscopic techniques complemented by quantum chemical calculations aided in understanding the electronic structure of the cofactor in all (detectable) states of the enzymatic process. Together with isotope labeling, these techniques also revealed the binding of the two substrate water molecules to the cluster. These results are described in the context of recent progress using X-ray crystallography with free-electron lasers on these intermediates. The data are instrumental for developing a model for the biological water oxidation cycle.

Co(iii) complexes are increasingly prevalent in homogeneous catalysis. Catalytic cycles involve multiple intermediates, many of which will feature unsaturated metal centres. This raises the possibility of multi-state character along reaction pathways and so requires an accurate approach to calculating spin-state energetics. Here we report an assessment of the performance of DLPNO-CCSD(T) (domain-based local pair natural orbital approximation to coupled cluster theory) against experimental 1Co to 3Co spin splitting energies for a series of pseudo-octahedral Co(iii) complexes. The alternative NEVPT2 (strongly-contracted n-electron valence perturbation theory) and a range of density functionals are also assessed. DLPNO-CCSD(T) is identified as a highly promising method, with mean absolute deviations (MADs) as small as 1.3 kcal mol-1 when Kohn-Sham reference orbitals are used. DLPNO-CCSD(T) out-performs NEVPT2 for which a MAD of 3.5 kcal mol-1 can be achieved when a (10,12) active space is employed. Of the nine DFT methods investigated TPSS is the leading functional, with a MAD of 1.9 kcal mol-1. Our results show how DLPNO-CCSD(T) can provide accurate spin state energetics for Co(iii) species in particular and first row transition metal systems in general. DLPNO-CCSD(T) is therefore a promising method for applications in the burgeoning field of homogeneous catalysis based on Co(iii) species.