Elucidating substrate binding in the light-dependent protochlorophyllide oxidoreductase.

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

Excitation landscape of the CP43 photosynthetic antenna complex from multiscale simulations.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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 (S1) to more oxidized intermediates (S2 and S3), eventually cycling back to the most reduced S0. However, the interpretation of these models is controversial because geometric parameters within the Mn4CaO5 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, S1 → S2, 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 S2 state of the OEC. We show that the 1F/S2 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 S2 state and with the nature of the S1 → S2 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 FeCl3·6H2O in various solvents gave complexes [Fe4Cl6(OMe)2(L1)2]·0.7MeCN·0.4MeOH (1·0.7MeCN·0.4MeOH) and [Fe3Cl4(bicine)(L2)]·Me2CO·0.2H2O (2·Me2CO·0.2H2O). The ligands (L1)2- = pyCO(Me)CHCOpy (in 1) and (L2)2- = pyCO(Me)CH2CO(OMe)py (in 2) are formed in situ, through an aldol reaction-type mechanism between the carbanion pyC(O)CH2- (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 -CH2- 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 [Fe2O2] rhombic units in cis-arrangement. The two terminal FeIII ions display distorted square pyramidal geometry and the two central FeIII ions are distorted octahedral. The molecular structure of 2 consists of two corner-sharing [Fe2O2] rhombic units, with the two terminal FeIII ions in distorted square pyramidal geometry and the central FeIII in distorted octahedral. The differentiation in the coordination environment of the FeIII 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 CuI 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 Mn4 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 QB. 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.