Correlating Structure with Spectroscopy in Ascorbate Peroxidase Compound II

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.

Exploring Electrophilic Hydrophosphination via Metal Phosphenium Intermediates

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(PR2H)(PR2)]+ (R=Ph, Tolp), 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(PR2H)2(PR2)]+, 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.

Comprehensive Evaluation of Models for Ammonia Binding to the Oxygen Evolving Complex of Photosystem II

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 H2O/OH-/O2- group or addition as an extra ligand to a metal ion of the Mn4CaO5 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 S2 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.

Structure and Flexibility of Copper-Modified DNA G-Quadruplexes Investigated by 19 F ENDOR Experiments at 34 GHz

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.

Conformational energies of reference organic molecules: benchmarking of common efficient computational methods against coupled cluster theory

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.

Nature of S-States in the Oxygen-Evolving Complex Resolved by High-Energy Resolution Fluorescence Detected X-ray Absorption Spectroscopy

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 Mn4CaO5 cluster that cycles through five oxidation states Si (i = 0-4). The S3 state is the last metastable state before the O2 evolution. Its electronic structure and nature of the S2 → S3 transition are key topics of persisting controversy. Most spectroscopic studies suggest that the S3 state consists of four Mn(IV) ions, compared to the Mn(III)Mn(IV)3 of the S2 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 S3 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 S3 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 S3 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.

Reaction Center Excitation in Photosystem II: From Multiscale Modeling to Functional Principles

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 (ChlD1) 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 (PheoD1): a fast pathway with ChlD1 as the primary electron donor (short-range charge-separation) and a slow pathway with PD1PD2 as the initial donor (long-range charge separation). The low-energy spectrum is dominated by two states with significant charge-transfer character, ChlD1δ+PheoD1δ- and PD1δ+PheoD1δ-. 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.

Reactive high-spin iron(IV)-oxo sites through dioxygen activation in a metal-organic framework

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 O2 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 O2 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 O2, the framework is competent for catalytic oxygenation of cyclohexane and the stoichiometric conversion of ethane to ethanol.

Triplet states in the reaction center of Photosystem II

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 ChlD1 pigment of the RC is the primary electron donor, forming a charge-separated radical pair with the vicinal pheophytin PheoD1 (ChlD1+PheoD1-). To avert charge recombination, the electron is further transferred to plastoquinone QA, whereas the hole relaxes to a central pair of chlorophylls (PD1PD2), 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 QA- 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.

Does Serial Femtosecond Crystallography Depict State-Specific Catalytic Intermediates of the Oxygen-Evolving Complex?

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₁) to more oxidized intermediates (S₂ and S₃), eventually cycling back to the most reduced S₀. However, the interpretation of these models is controversial because geometric parameters within the Mn₄CaO₅ 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₁ → S₂, 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₂ state of the OEC. We show that the 1F/S₂ 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₂ state and with the nature of the S₁ → S₂ 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.

Alternative Fast and Slow Primary Charge-Separation Pathways in Photosystem II

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.

Synthesis and structural, magnetic and spectroscopic characterization of iron(III) complexes with in situ formed ligands from methyl-2-pyridyl ketone transformations

Reactions of methyl-2-pyridyl ketone, pyCOMe, with FeCl₃·6H₂O in various solvents gave complexes [Fe₄Cl₆(OMe)₂(L1)₂]·0.7MeCN·0.4MeOH (1·0.7MeCN·0.4MeOH) and [Fe₃Cl₄(bicine)(L2)]·Me₂CO·0.2H₂O (2·Me₂CO·0.2H₂O). The ligands (L1)^2- = pyCO(Me)CHCOpy (in 1) and (L2)^2- = pyCO(Me)CH₂CO(OMe)py (in 2) are formed in situ, through an aldol reaction-type mechanism between the carbanion pyC(O)CH₂^- (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₂- 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₂O₂] 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₂O₂] 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.

Copper Complexes with Diazoolefin Ligands and their Photochemical Conversion into Alkenylidene Complexes

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.

Selective incorporation of 5-hydroxytryptophan blocks long range electron transfer in oxalate decarboxylase

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.

Water oxidation in oxygenic photosynthesis studied by magnetic resonance techniques

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₄ 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.

Functional Water Networks in Fully Hydrated Photosystem II

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.

Mechanism of the Aryl-F Bond-Forming Step from Bi(V) Fluorides

In this article, we describe a combined experimental and theoretical mechanistic investigation of the C(sp²)–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²) center, from a BF₄ 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.

Reconciling Local Coupled Cluster with Multireference Approaches for Transition Metal Spin-State Energetics

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₁) 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.

Decoding the Ambiguous Electron Paramagnetic Resonance Signals in the Lytic Polysaccharide Monooxygenase from Photorhabdus luminescens

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 PlAA10 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 PlAA10 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 PlAA10. 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 PlAA10 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.

The Electronic Origin of Far-Red-Light-Driven Oxygenic Photosynthesis

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.

Ionization Energies and Redox Potentials of Hydrated Transition Metal Ions: Evaluation of Domain-Based Local Pair Natural Orbital Coupled Cluster Approaches

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.

Characterization of a Triplet Vinylidene

Singlet vinylidenes (R₂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 ¹³C hyperfine coupling tensor was obtained by ¹³C-ENDOR measurements. Most strikingly, the isotropic ¹³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.

Electrostatic profiling of photosynthetic pigments: implications for directed spectral tuning

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.

Magnetic exchange coupling in Cu dimers studied with modern multireference methods and broken-symmetry coupled cluster theory

Spin-state energetics of exchange-coupled copper complexes pose a persistent challenge for applied quantum chemistry. Here, we provide a comprehensive comparison of all available theoretical approaches to the problem of exchange coupling in two antiferromagnetically coupled bis-μ-hydroxo Cu(II) dimers. The evaluated methods include multireference methods based on the density matrix renormalization group (DMRG), multireference methods that incorporate dynamic electron correlation either perturbatively, such as the N-electron valence state perturbation theory, or variationally, such as the difference-dedicated configuration interaction. In addition, we contrast the multireference results with those obtained using broken-symmetry approaches that utilize either density functional theory or, as demonstrated here for the first time in such systems, a local implementation of coupled cluster theory. The results show that the spin-state energetics of these copper dimers are dominated by dynamic electron correlation and represent an impossible challenge for multireference methods that rely on brute-force expansion of the active space to recover correlation energy. Therefore, DMRG-based methods even at the limit of their applicability cannot describe quantitatively the antiferromagnetic exchange coupling in these dimers, in contrast to dinuclear complexes of earlier transition metal ions. The convergence of the broken-symmetry coupled cluster approach is studied and shown to be a limiting factor for the practical application of the method. The advantages and disadvantages of all approaches are discussed, and recommendations are made for future developments.

Redox Isomerism in the S3 State of the Oxygen-Evolving Complex Resolved by Coupled Cluster Theory

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 Mn4 CaOx cofactor in the final metastable S3 state. These range from the widely accepted MnIV 4 oxo-hydroxo model, which presumes that O-O bond formation occurs in the ultimate transient intermediate (S4 ) of the catalytic cycle, to a MnIII 2 MnIV 2 peroxo model representative of the contrasting "early-onset" O-O bond formation hypothesis. Density functional theory energetics of suggested S3 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 S3 state. Our results enabled us to evaluate conflicting models for the S3 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.

Converged Structural and Spectroscopic Properties for Refined QM/MM Models of Azurin

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 33S 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.

Orientational Jahn-Teller Isomerism in the Dark-Stable State of Nature's Water Oxidase

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₂. 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₄CaO₅ cluster in the resting (dark‐stable) S₁ 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₁ state. This unique isomerism in the resting state is shown to be the electronic origin of valence isomerism in the S₂ state, establishing a functional role of orientational Jahn–Teller isomerism unprecedented in biological or artificial catalysis.

Structure-Spectroscopy Correlations for Intermediate Q of Soluble Methane Monooxygenase: Insights from QM/MM Calculations

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)₂O₂ closed-core conformation Fe(IV)₂O₂ 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 ¹⁶O–¹⁸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.

How Can We Predict Accurate Electrochromic Shifts for Biochromophores? A Case Study on the Photosynthetic Reaction Center

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-CytB559 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.

Chlorophyll excitation energies and structural stability of the CP47 antenna of photosystem II: a case study in the first-principles simulation of light-harvesting complexes

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.

Arrested Substrate Binding Resolves Catalytic Intermediates in Higher-Plant Water Oxidation

Among the intermediate catalytic steps of the water-oxidizing Mn4 CaO5 cluster of photosystem II (PSII), the final metastable S3 state is critically important because it binds one substrate and precedes O2 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 S3 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 MnIV3 subunit interacting ferromagnetically via anisotropic exchange with a coordinatively unsaturated MnIV ion. These results resolve and define the structural heterogeneity of the S3 state, providing constraints on the S3 to S4 transition, on substrate identity and delivery pathways, and on the mechanism of O-O bond formation.

Reversible Silylium Transfer between P-H and Si-H Donors

The Mo=PR₂ π* 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₃ SiH to give a Mo-bound secondary phosphine ligand, Mo-PR₂ H. The resulting Et₃ Si⁺ ion remains associated with the Mo complex, stabilized by η¹ -P-H donation, yet undergoes rapid exchange with an η¹ -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.

Successes, challenges, and opportunities for quantum chemistry in understanding metalloenzymes for solar fuels research

Overview of the rich and diverse contributions of quantum chemistry to understanding the structure and function of the biological archetypes for solar fuel research, photosystem II and hydrogenases.

Comparison of Density Functional and Correlated Wave Function Methods for the Prediction of Cu(II) Hyperfine Coupling Constants

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.

Protein Matrix Control of Reaction Center Excitation in Photosystem II

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 (PD1, PD2, ChlD1, ChlD2) and two pheophytin molecules (PheoD1 and PheoD2). 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 ChlD1 the chromophore with the lowest site energy. Multipigment calculations show that the protein matrix renders the ChlD1 → PheoD1 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" PD1-PD2, 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.

Accurate Computation of the Absorption Spectrum of Chlorophyll a with Pair Natural Orbital Coupled Cluster Methods

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.

Current Understanding of the Mechanism of Water Oxidation in Photosystem II and Its Relation to XFEL Data

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 Mn4CaOxcomplex, 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.

Unusual 31P Hyperfine Strain Effects in a Conformationally Flexible Cu(II) Complex Revealed by Two-Dimensional Pulse EPR Spectroscopy

Strain effects on g and metal hyperfine coupling tensors, A, are often manifested in Electron Paramagnetic Resonance (EPR) spectra of transition metal complexes, as a result of their intrinsic and/or solvent-mediated structural variations. Although distributions of these tensors are quite common and well understood in continuous-wave (cw) EPR spectroscopy, reported strain effects on ligand hyperfine coupling constants are rather scarce. Here we explore the case of a conformationally flexible Cu(II) complex, [Cu{Ph₂P(O)NP(O)Ph₂-κ²O,O'}₂], bearing P atoms in its second coordination sphere and exhibiting two structurally distinct CuO₄ coordination spheres, namely a square planar and a tetrahedrally distorted one, as revealed by X-ray crystallography. The Hyperfine Sublevel Correlation (HYSCORE) spectra of this complex exhibit ³¹P correlation ridges that have unusual inverse or so-called "boomerang" shapes and features that cannot be reproduced by standard simulation procedures assuming only one set of magnetic parameters. Our work shows that a distribution of isotropic hyperfine coupling constants (hfc) spanning a range between negative and positive values is necessary in order to describe in detail the unusual shapes of HYSCORE spectra. By employing DFT calculations we show that these hfc correspond to molecules showing variable distortions from square planar to tetrahedral geometry, and we demonstrate that line shape analysis of such HYSCORE spectra provides new insight into the conformation-dependent spectroscopic response of the spin system under investigation.