Water Networks in Photosystem II Using Crystalline Molecular Dynamics Simulations and Room-Temperature XFEL Serial Crystallography

Ammonia Binding to the Oxygen-Evolving Complex Probed by High-Energy Resolution Fluorescence Detected X-Ray Absorption Spectroscopy

The insertion pathways and binding sites of substrate water molecules at the catalytic Mn4CaO5 cluster of the oxygen-evolving complex (OEC) in photosystem II (PSII) remain a fundamentally unresolved question toward understanding biological water oxidation. To address this question, small molecules have been employed as "water analogues" to probe substrate binding to the OEC. In this context, the binding of ammonia has been extensively investigated and discussed using spectroscopic, structural, and quantum chemical methods, but a definitive answer regarding the ammonia binding site has not yet been achieved. Herein, we present high-energy resolution fluorescence detected (HERFD) Mn K-edge X-ray absorption spectroscopy (XAS) in ammonia-treated S2 state samples of the OEC. Pre-edge features were correlated with possible structural models with the aid of quantum chemical calculations. The comparison of calculated and experimental difference spectra between the native and ammonia-treated samples allows us to evaluate different modes of ammonia interaction with the OEC. The combined spectroscopic and theoretical investigation suggests the substitution of the terminal water ligand W2 on Mn4 as the most plausible ammonia binding mode, followed closely by the substitution of the second terminal water ligand (W1), and the coordination of ammonia on Mn1 as a sixth ligand. Our results are in line with the leading interpretations of other spectroscopic and kinetic studies, converging on the conclusion that the Mn4 ion is either the most accessible or the strongest binding site for substrate analogues.

Evaluating long-range orientational ordering of water around proteins: signature of a tug-of-war scenario

Long-range perturbations of water structure and dynamics by biomolecules are of great interest owing to their potential role in biomolecular recognition. In this article, we examined the local and long-range orientational structure of water molecules surrounding proteins with different total charges (+8, 0 and -8), both with and without the presence of a physiological salt environment. A prominent population of in-oriented water molecules was observed in the first hydration shell of the proteins, irrespective of their total charges. Starting from the third hydration layer, water molecules primarily reflected the total charge of the respective protein. This long-range ordering persisted up to the ninth hydration layer without a physiological salt environment and vanished beyond the fifth hydration shell in the presence of a physiological salt environment. Long-range orientational ordering around different types of surface atoms of a protein showed a particularly rich and heterogeneous behaviour. When the surface atom's charge and the protein's total charge were opposite, a clear signature of a tug-of-war was demonstrated in the long-range orientational ordering of water molecules. While water molecules reported the surface atom's charge at shorter distances, at longer distances, water molecules reported the total charge of the protein, with a crossover occurring around 10 Å. This phenomenon persisted even in the presence of a physiological salt environment. Evidence of destructive/constructive superposition of water-mediated orientation waves originating from two individual proteins with similar/opposite total charges was also demonstrated. These results are important for understanding long-range water-mediated recognition phenomena among biomolecules (e.g., protein-protein, protein-ligand, and protein-DNA interactions).

An Unexpected Water Channel in the Light-Harvesting Complex of a Diatom: Implications for the Switch between Light Harvesting and Photoprotection

Many important processes in cells depend on the transfer of protons through water wires embedded in transmembrane proteins. Herein, we have performed more than 55 μs all-atom simulations of the light-harvesting complex of a diatom, i.e., the fucoxanthin and chlorophyll a/c binding protein (FCP) from the marine diatom Phaeodactylum tricornutum. Diatoms are unique models to study natural photosynthesis as they exert an efficient light-harvesting machinery with a robust pH-dependent photoprotective mechanism. The present study reports on the dynamics of an FCP monomer, a dimer, and a tetramer at varying pH values. Surprisingly, we have identified at low pH a water channel across FCP that selectively hydrates and protonates the acrylate of a Chl-c2 pigment located in the middle of the membrane. These results are further supported by QM/MM calculations and steered MD simulations on the proton dynamics. It is shown that proton hopping events between the lumenal and stromal sides of the membrane through the observed water channel are highly disfavored. This hindrance is due to the presence of residues Arg31 and Lys82 close to the acrylate, along with an hydronium desolvation penalty that shows close similarities to the water conductance in aquaporins. Furthermore, we provide strong evidence that this identified water channel is governing the transition between light-harvesting and photoprotective states of the major FCP complex in the diatom P. tricornutum.

Mechanism of proton release during water oxidation in Photosystem II

Photosystem II (PSII) catalyzes light-driven water oxidation that releases dioxygen into our atmosphere and provides the electrons needed for the synthesis of biomass. The catalysis occurs in the oxygen-evolving oxo-manganese-calcium (Mn4O5Ca) cluster that drives the oxidation and deprotonation of substrate water molecules leading to the O2 formation. However, despite recent advances, the mechanism of these reactions remains unclear and much debated. Here, we show that the light-driven Tyr161D1 (Yz) oxidation adjacent to the Mn4O5Ca cluster, decreases the barrier for proton transfer from the putative substrate water molecule (W3/Wx) to Glu310D2, accessible to the luminal bulk. By combining hybrid quantum/classical (QM/MM) free energy calculations with atomistic molecular dynamics simulations, we probe the energetics of the proton transfer along the Cl1 pathway. We demonstrate that the proton transfer occurs via water molecules and a cluster of conserved carboxylates, driven by redox-triggered electric fields directed along the pathway. Glu65D1 establishes a local molecular gate that controls the proton transfer to the luminal bulk, while Glu312D2 acts as a local proton storage site. The identified gating region could be important in preventing backflow of protons to the Mn4O5Ca cluster. The structural changes, derived here based on the dark-state PSII structure, strongly support recent time-resolved X-ray free electron laser data of the S3 → S4 transition (Bhowmick et al. Nature 617, 2023) and reveal the mechanistic basis underlying deprotonation of the substrate water molecules. Our findings provide insight into the water oxidation mechanism of PSII and show how the interplay between redox-triggered electric fields, ion-pairs, and hydration effects control proton transport reactions.

Structure Function Studies of Photosystem II Using X-Ray Free Electron Lasers

The structure and mechanism of the water-oxidation chemistry that occurs in photosystem II have been subjects of great interest. The advent of X-ray free electron lasers allowed the determination of structures of the stable intermediate states and of steps in the transitions between these intermediate states, bringing a new perspective to this field. The room-temperature structures collected as the photosynthetic water oxidation reaction proceeds in real time have provided important novel insights into the structural changes and the mechanism of the water oxidation reaction. The time-resolved measurements have also given us a view of how this reaction-which involves multielectron, multiproton processes-is facilitated by the interaction of the ligands and the protein residues in the oxygen-evolving complex. These structures have also provided a picture of the dynamics occurring in the channels within photosystem II that are involved in the transport of the substrate water to the catalytic center and protons to the bulk.

Cryo-electron microscopy reveals hydrogen positions and water networks in photosystem II

Photosystem II starts the photosynthetic electron transport chain that converts solar energy into chemical energy and thus sustains life on Earth. It catalyzes two chemical reactions: water oxidation to molecular oxygen and plastoquinone reduction. Coupling of electron and proton transfer is crucial for efficiency; however, the molecular basis of these processes remains speculative owing to uncertain water binding sites and the lack of experimentally determined hydrogen positions. We thus collected high-resolution cryo-electron microscopy data of fully hydrated photosystem II from the thermophilic cyanobacterium Thermosynechococcus vestitus to a final resolution of 1.71 angstroms. The structure reveals several previously undetected partially occupied water binding sites and more than half of the hydrogen and proton positions. This clarifies the pathways of substrate water binding and plastoquinone B protonation.

Structural biology in the age of X-ray free-electron lasers and exascale computing

Serial femtosecond X-ray crystallography has emerged as a powerful method for investigating biomolecular structure and dynamics. With the new generation of X-ray free-electron lasers, which generate ultrabright X-ray pulses at megahertz repetition rates, we can now rapidly probe ultrafast conformational changes and charge movement in biomolecules. Over the last year, another innovation has been the deployment of Frontier, the world's first exascale supercomputer. Synergizing extremely high repetition rate X-ray light sources and exascale computing has the potential to accelerate discovery in biomolecular sciences. Here we outline our perspective on each of these remarkable innovations individually, and the opportunities and challenges in yoking them within an integrated research infrastructure.

The biophysics of water in cell biology: perspectives on a keystone for both marine sciences and cancer research

The biophysics of water, has been debated over more than a century. Although its importance is still underestimated, significant breakthroughs occurred in recent years. The influence of protein condensation on water availability control was documented, new findings on water-transport proteins emerged, and the way water molecules rearrange to minimize free energy at interfaces was deciphered, influencing membrane thermodynamics. The state of knowledge continued to progress in the field of deep-sea marine biology, highlighting unknown effects of high hydrostatic pressure and/or temperature on interactions between proteins and ligands in extreme environments, and membrane structure adaptations. The role of osmolytes in protein stability control under stress is also discussed here in relation to fish egg hydration/buoyancy. The complexity of water movements within the cell is updated, all these findings leading to a better view of their impact on many cellular processes. The way water flow and osmotic gradients generated by ion transport work together to produce the driving force behind cell migration is also relevant to both marine biology and cancer research. Additional common points concern water dynamic changes during the neoplastic transformation of cells and tissues, or embryo development. This could improve imaging techniques, early cancer diagnosis, and understanding of the molecular and physiological basis of buoyancy for many marine species.

Independent Mutation of Two Bridging Carboxylate Ligands Stabilizes Alternate Conformers of the Photosynthetic O2-Evolving Mn4CaO5 Cluster in Photosystem II

The O2-evolving Mn4CaO5 cluster in photosystem II is ligated by six carboxylate residues. One of these is D170 of the D1 subunit. This carboxylate bridges between one Mn ion (Mn4) and the Ca ion. A second carboxylate ligand is D342 of the D1 subunit. This carboxylate bridges between two Mn ions (Mn1 and Mn2). D170 and D342 are located on opposite sides of the Mn4CaO5 cluster. Recently, it was shown that the D170E mutation perturbs both the intricate networks of H-bonds that surround the Mn4CaO5 cluster and the equilibrium between different conformers of the cluster in two of its lower oxidation states, S1 and S2, while still supporting O2 evolution at approximately 50% the rate of the wild type. In this study, we show that the D342E mutation produces much the same alterations to the cluster's FTIR and EPR spectra as D170E, while still supporting O2 evolution at approximately 20% the rate of the wild type. Furthermore, the double mutation, D170E + D342E, behaves similarly to the two single mutations. We conclude that D342E alters the equilibrium between different conformers of the cluster in its S1 and S2 states in the same manner as D170E and perturbs the H-bond networks in a similar fashion. This is the second identification of a Mn4CaO5 metal ligand whose mutation influences the equilibrium between the different conformers of the S1 and S2 states without eliminating O2 evolution. This finding has implications for our understanding of the mechanism of O2 formation in terms of catalytically active/inactive conformations of the Mn4CaO5 cluster in its lower oxidation states.

Stabilization of Charge-Transfer Excited States in Biological Systems: A Computational Focus on the Special Pair in Photosystem II Reaction Centers

Charge-transfer (CT) excited states play an important role in many biological processes. However, many computational approaches often inadequately address the equilibration effects of nuclear and environmental degrees of freedom on these states. One prominent example of systems in which CT states are of utmost importance is reaction centers (RC) in photosystems. Here we use a multiscale approach combined with time-dependent density functional theory to explore the lowest CT excited state of the special pair PD1-PD2 in the Photosystem II-RC of a cyanobacterium. We find that the nonequilibrium CT excited state resides near the Soret band, making an exciton the lowest-energy excited state. However, accounting for nuclear and state-specific dielectric equilibration along the CT potential energy surface (PES), the CT state PD1--PD2+ stabilizes energetically below the excitonic state. This underscores the crucial role of state-specific solvation in mapping the PES of CT states, as demonstrated in a simplified dimer model.

Changes in an enzyme ensemble during catalysis observed by high-resolution XFEL crystallography

Enzymes populate ensembles of structures necessary for catalysis that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography at an x-ray free electron laser to observe catalysis in a designed mutant isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations, and formation of the thioimidate intermediate selects for catalytically competent substates. The influence of cysteine ionization on the ICH ensemble is validated by determining structures of the enzyme at multiple pH values. Large molecular dynamics simulations in crystallo and time-resolved electron density maps show that Asp17 ionizes during catalysis and causes conformational changes that propagate across the dimer, permitting water to enter the active site for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics.

Occupancy Analysis of Water Molecules inside Channels within 25 Å Radius of the Oxygen-Evolving Center of Photosystem II in Molecular Dynamics Simulations

At room temperature and neutral pH, the oxygen-evolving center (OEC) of photosystem II (PSII) catalyzes water oxidation. During this process, oxygen is released from the OEC, while substrate waters are delivered to the OEC and protons are passed from the OEC to the lumen through water channels known as the narrow or the O4 channel, broad or the Cl1 channel, and large or the O1 channel. Protein residues lining the surfaces of these channels play a critical role in stabilizing the hydrogen-bonding networks that assist in the process. We carried out an occupancy analysis to better understand the structural and possible substrate water dynamics in full PSII monomer molecular dynamics (MD) trajectories in both the S1 and S2 states. We find that the equilibrated positions of water molecules derived from MD-derived electron density maps largely match the experimentally observed positions in crystallography. Furthermore, the occupancy reduction in MD simulations of some water molecules inside the single-filed narrow channel also correlates well with the crystallographic data during a structural transition when the S1 state of the OEC advances to the S2 state. The overall reduced occupancies of water molecules are the source of their "vacancy-hopping" dynamic nature inside these channels, unlike water molecules inside an ice lattice where all water molecules have a fixed unit occupancy. We propose on the basis of findings in our structural and molecular dynamics analysis that the water molecule occupying a pocket formed by D1-D61, D1-S169, and O4 of the OEC could be the last steppingstone to enter into the OEC and that the broad channel may be favored for proton transfer.

Evolution in the Design of Water Oxidation Catalysts with Transition-Metals: A Perspective on Biological, Molecular, Supramolecular, and Hybrid Approaches

Increased demand for a carbon-neutral sustainable energy scheme augmented by climatic threats motivates the design and exploration of novel approaches that reserve intermittent solar energy in the form of chemical bonds in molecules and materials. In this context, inspired by biological processes, artificial photosynthesis has garnered significant attention as a promising solution to convert solar power into chemical fuels from abundantly found H2O. Among the two redox half-reactions in artificial photosynthesis, the four-electron oxidation of water according to 2H2O → O2 + 4H+ + 4e- comprises the major bottleneck and is a severe impediment toward sustainable energy production. As such, devising new catalytic platforms, with traditional concepts of molecular, materials and biological catalysis and capable of integrating the functional architectures of the natural oxygen-evolving complex in photosystem II would certainly be a value-addition toward this objective. In this review, we discuss the progress in construction of ideal water oxidation catalysts (WOCs), starting with the ingenuity of the biological design with earth-abundant transition metal ions, which then diverges into molecular, supramolecular and hybrid approaches, blurring any existing chemical or conceptual boundaries. We focus on the geometric, electronic, and mechanistic understanding of state-of-the-art homogeneous transition-metal containing molecular WOCs and summarize the limiting factors such as choice of ligands and predominance of environmentally unrewarding and expensive noble-metals, necessity of high-valency on metal, thermodynamic instability of intermediates, and reversibility of reactions that create challenges in construction of robust and efficient water oxidation catalyst. We highlight how judicious heterogenization of atom-efficient molecular WOCs in supramolecular and hybrid approaches put forth promising avenues to alleviate the existing problems in molecular catalysis, albeit retaining their fascinating intrinsic reactivities. Taken together, our overview is expected to provide guiding principles on opportunities, challenges, and crucial factors for designing novel water oxidation catalysts based on a synergy between conventional and contemporary methodologies that will incite the expansion of the domain of artificial photosynthesis.

Mechanism of Proton Transfer through the D1-E65/D2-E312 Gate during Photosynthetic Water Oxidation

In photosystem II, the D1-E65/D2-E312 dyad in the Cl-1 channel has been proposed to play a pivotal role in proton transfer during water oxidation. However, the precise mechanism remains elusive. Here, the proton transfer mechanism within the Cl-1 channel was investigated using quantum mechanics/molecular mechanics calculations. The molecular vibration of the E65/E312 dyad and its deuteration effect revealed that the recently suggested stepwise proton transfer, i.e., initial proton release from the dyad followed by slow reprotonation, does not occur in the Cl-1 channel. Instead, proton transfer is proposed to take place via a conformational change at the E65/E312 dyad, acting as a gate. In its closed form, a proton is trapped within the dyad, preventing forward proton transfer. This closed form converts into the open form, where protonated D1-E65 provides a hydrogen bond to the water network, thereby facilitating fast Grotthuss-type proton transfer.

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.

A breath of sunshine: oxygenic photosynthesis by functional molecular architectures

The conversion of light into chemical energy is the game-changer enabling technology for the energetic transition to renewable and clean solar fuels. The photochemistry of interest includes the overall reductive/oxidative splitting of water into hydrogen and oxygen and alternatives based on the reductive conversion of carbon dioxide or nitrogen, as primary sources of energy-rich products. Devices capable of performing such transformations are based on the integration of three sequential core functions: light absorption, photo-induced charge separation, and the photo-activated breaking/making of molecular bonds via specific catalytic routes. The key to success does not rely simply on the individual components' performance, but on their optimized integration in terms of type, number, geometry, spacing, and linkers dictating the photosynthetic architecture. Natural photosynthesis has evolved along this concept, by integrating each functional component in one specialized "body" (from the Greek word "soma") to enable the conversion of light quanta with high efficiency. Therefore, the natural "quantasome" represents the key paradigm to inspire man-made constructs for artificial photosynthesis. The case study presented in this perspective article deals with the design of artificial photosynthetic systems for water oxidation and oxygen production, engineered as molecular architectures then rendered on electrodic surfaces. Water oxidation to oxygen is indeed the pervasive oxidative reaction used by photosynthetic organisms, as the source of reducing equivalents (electrons and protons) to be delivered for the processing of high-energy products. Considering the vast and abundant supply of water (including seawater) as a renewable source on our planet, this is also a very appealing option for photosynthetic energy devices. We will showcase the progress in the last 15 years (2009-2023) in the strategies for integrating functional building blocks as molecular photosensitizers, multi-redox water oxidation catalysts and semiconductor materials, highlighting how additional components such as redox mediators, hydrophilic/hydrophobic pendants, and protective layers can impact on the overall photosynthetic performance. Emerging directions consider the modular tuning of the multi-component device, in order to target a diversity of photocatalytic oxidations, expanding the scope of the primary electron and proton sources while enhancing the added-value of the oxidation product beyond oxygen: the selective photooxidation of organics combines the green chemistry vision with renewable energy schemes and is expected to explode in coming years.

Changes in an Enzyme Ensemble During Catalysis Observed by High Resolution XFEL Crystallography

Enzymes populate ensembles of structures with intrinsically different catalytic proficiencies that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography (MISC) at an X-ray free electron laser (XFEL) to observe catalysis in a designed mutant (G150T) isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations and formation of the thioimidate catalytic intermediate selects for catalytically competent substates. A prior proposal for active site cysteine charge-coupled conformational changes in ICH is validated by determining structures of the enzyme over a range of pH values. A combination of large molecular dynamics simulations of the enzyme in crystallo and time-resolved electron density maps shows that ionization of the general acid Asp17 during catalysis causes additional conformational changes that propagate across the dimer interface, connecting the two active sites. These ionization-linked changes in the ICH conformational ensemble permit water to enter the active site in a location that is poised for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics.