All-Atom Molecular Dynamics Simulation of Photosystem II Embedded in Thylakoid Membrane

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

Theoretical elucidation of the structure, bonding, and reactivity of the CaMn4Ox clusters in the whole Kok cycle for water oxidation embedded in the oxygen evolving center of photosystem II. New molecular and quantum insights into the mechanism of the O-O bond formation

This paper reviews our historical developments of broken-symmetry (BS) and beyond BS methods that are applicable for theoretical investigations of metalloenzymes such as OEC in PSII. The BS hybrid DFT (HDFT) calculations starting from high-resolution (HR) XRD structure in the most stable S1 state have been performed to elucidate structure and bonding of whole possible intermediates of the CaMn4Ox cluster (1) in the Si (i = 0 ~ 4) states of the Kok cycle. The large-scale HDFT/MM computations starting from HR XRD have been performed to elucidate biomolecular system structures which are crucial for examination of possible water inlet and proton release pathways for water oxidation in OEC of PSII. DLPNO CCSD(T0) computations have been performed for elucidation of scope and reliability of relative energies among the intermediates by HDFT. These computations combined with EXAFS, XRD, XFEL, and EPR experimental results have elucidated the structure, bonding, and reactivity of the key intermediates, which are indispensable for understanding and explanation of the mechanism of water oxidation in OEC of PSII. Interplay between theory and experiments have elucidated important roles of four degrees of freedom, spin, charge, orbital, and nuclear motion for understanding and explanation of the chemical reactivity of 1 embedded in protein matrix, indicating the participations of the Ca(H2O)n ion and tyrosine(Yz)-O radical as a one-electron acceptor for the O-O bond formation. The Ca-assisted Yz-coupled O-O bond formation mechanisms for water oxidation are consistent with recent XES and very recent time-resolved SFX XFEL and FTIR results.

Evolutionary diversity of proton and water channels on the oxidizing side of photosystem II and their relevance to function

One of the reasons for the high efficiency and selectivity of biological catalysts arise from their ability to control the pathways of substrates and products using protein channels, and by modulating the transport in the channels using the interaction with the protein residues and the water/hydrogen-bonding network. This process is clearly demonstrated in Photosystem II (PS II), where its light-driven water oxidation reaction catalyzed by the Mn4CaO5 cluster occurs deep inside the protein complex and thus requires the transport of two water molecules to and four protons from the metal center to the bulk water. Based on the recent advances in structural studies of PS II from X-ray crystallography and cryo-electron microscopy, in this review we compare the channels that have been proposed to facilitate this mass transport in cyanobacteria, red and green algae, diatoms, and higher plants. The three major channels (O1, O4, and Cl1 channels) are present in all species investigated; however, some differences exist in the reported structures that arise from the different composition and arrangement of membrane extrinsic subunits between the species. Among the three channels, the Cl1 channel, including the proton gate, is the most conserved among all photosynthetic species. We also found at least one branch for the O1 channel in all organisms, extending all the way from Ca/O1 via the 'water wheel' to the lumen. However, the extending path after the water wheel varies between most species. The O4 channel is, like the Cl1 channel, highly conserved among all species while having different orientations at the end of the path near the bulk. The comparison suggests that the previously proposed functionality of the channels in T. vestitus (Ibrahim et al., Proc Natl Acad Sci USA 117:12624-12635, 2020; Hussein et al., Nat Commun 12:6531, 2021) is conserved through the species, i.e. the O1-like channel is used for substrate water intake, and the tighter Cl1 and O4 channels for proton release. The comparison does not eliminate the potential role of O4 channel as a water intake channel. However, the highly ordered hydrogen-bonded water wire connected to the Mn4CaO5 cluster via the O4 may strongly suggest that it functions in proton release, especially during the S0 → S1 transition (Saito et al., Nat Commun 6:8488, 2015; Kern et al., Nature 563:421-425, 2018; Ibrahim et al., Proc Natl Acad Sci USA 117:12624-12635, 2020; Sakashita et al., Phys Chem Chem Phys 22:15831-15841, 2020; Hussein et al., Nat Commun 12:6531, 2021).

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

Structural dynamics of water and its hydrogen-bonding networks play an important role in enzyme function via the transport of protons, ions, and substrates. To gain insights into these mechanisms in the water oxidation reaction in Photosystem II (PS II), we have performed crystalline molecular dynamics (MD) simulations of the dark-stable S1 state. Our MD model consists of a full unit cell with 8 PS II monomers in explicit solvent (861 894 atoms), enabling us to compute the simulated crystalline electron density and to compare it directly with the experimental density from serial femtosecond X-ray crystallography under physiological temperature collected at X-ray free electron lasers (XFELs). The MD density reproduced the experimental density and water positions with high fidelity. The detailed dynamics in the simulations provided insights into the mobility of water molecules in the channels beyond what can be interpreted from experimental B-factors and electron densities alone. In particular, the simulations revealed fast, coordinated exchange of waters at sites where the density is strong, and water transport across the bottleneck region of the channels where the density is weak. By computing MD hydrogen and oxygen maps separately, we developed a novel Map-based Acceptor-Donor Identification (MADI) technique that yields information which helps to infer hydrogen-bond directionality and strength. The MADI analysis revealed a series of hydrogen-bond wires emanating from the Mn cluster through the Cl1 and O4 channels; such wires might provide pathways for proton transfer during the reaction cycle of PS II. Our simulations provide an atomistic picture of the dynamics of water and hydrogen-bonding networks in PS II, with implications for the specific role of each channel in the water oxidation reaction.

Phase Transition and Phase Separation in Realistic Thylakoid Lipid Membrane of Marine Algae in All-Atom Simulations

Thylakoid membranes are specialized membranes predominantly composed of uncommon galacto- and sulfolipids, having distinct roles in photosynthesis. Large acyl chain variety and richness in polyunsaturated fatty acid (PUFA) content of thylakoid lipids further add to the compositional complexity. The function of these membrane systems is intimately dependent on the fluidity of its lipid matrix, which is strongly modulated by the lipid composition and temperature. The present work, employing extensive atomistic simulations, provides the first atomistic view of the phase transition and domain coexistence in a model membrane composed of thylakoid lipids of a commercially important red alga Gracilaria corticata between 10 and 40 °C. The growth and photosynthetic activity of marine algae are greatly influenced by the seawater temperature. So far, little is known about the molecular organization of lipids in thylakoid membranes, in particular their adaptive arrangements under temperature stress. Our simulations show that the algal thylakoid membrane undergoes a transition from a gel-like phase at a low temperature, 10-15 °C, to a homogeneous liquid-crystalline phase at a high temperature, 40 °C. Clear evidence of spontaneous phase separation into coexisting nanoscale domains is detected at intermediate temperatures nearing the optimal growth temperature range. Particularly, at 25-30 °C, we identified the formation of a stable ripple phase, where the gel-like domains rich in saturated and nearly hexagonally packed lipids were separated from fluid-like domains enriched in lipids containing PUFA chains. The phase separation is driven by the spontaneous and preferential segregation of lipids into differentially ordered domains, mainly depending on the acyl chain types. Cholesterol impairs the phase transition and the emergence of domains and induces a fairly uniform liquid-ordered phase in the membrane over the temperatures studied. This work improves the understanding of the properties and reorganization of lipids in the thylakoid membrane in response to temperature variation.

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.

Crystal structures of photosystem II from a cyanobacterium expressing psbA2 in comparison to psbA3 reveal differences in the D1 subunit

Three psbA genes (psbA₁, psbA₂, and psbA₃) encoding the D1 subunit of photosystem II (PSII) are present in the thermophilic cyanobacterium Thermosynechococcus elongatus and are expressed differently in response to changes in the growth environment. To clarify the functional differences of the D1 protein expressed from these psbA genes, PSII dimers from two strains, each expressing only one psbA gene (psbA₂ or psbA₃), were crystallized, and we analyzed their structures at resolutions comparable to previously studied PsbA1-PSII. Our results showed that the hydrogen bond between pheophytin/D1 (Pheo_D1) and D1-130 became stronger in PsbA2- and PsbA3-PSII due to change of Gln to Glu, which partially explains the increase in the redox potential of Pheo_D1 observed in PsbA3. In PsbA2, one hydrogen bond was lost in Pheo_D1 due to the change of D1-Y147F, which may explain the decrease in stability of Pheo_D1 in PsbA2. Two water molecules in the Cl-1 channel were lost in PsbA2 due to the change of D1-P173M, leading to the narrowing of the channel, which may explain the lower efficiency of the S-state transition beyond S₂ in PsbA2-PSII. In PsbA3-PSII, a hydrogen bond between D1-Ser270 and a sulfoquinovosyl-diacylglycerol molecule near Q_B disappeared due to the change of D1-Ser270 in PsbA1 and PsbA2 to D1-Ala270. This may result in an easier exchange of bound Q_B with free plastoquinone, hence an enhancement of oxygen evolution in PsbA3-PSII due to its high Q_B exchange efficiency. These results provide a structural basis for further functional examination of the three PsbA variants.

Role of D1-Glu65 in Proton Transfer during Photosynthetic Water Oxidation in Photosystem II

Photosynthetic water oxidation takes place at the Mn₄CaO₅ cluster in photosystem II (PSII) through a light-driven cycle of five intermediates called S states (S₀-S₄). Although the PSII structures have shown the presence of several channels around the Mn₄CaO₅ cluster leading to the lumen, the pathways for proton release in the individual S-state transitions remain unidentified. Here, we studied the involvement of the so-called Cl channel in proton transfer during water oxidation by examining the effect of the mutation of D1-Glu65, a key residue in this channel, to Ala using Fourier transform infrared difference and time-resolved infrared spectroscopies together with thermoluminescence and delayed luminescence measurements. It was shown that the structure and the redox property of the catalytic site were little affected by the D1-Glu65Ala mutation. In the S₂ → S₃ transition, the efficiency was still high and the transition rate was only moderately retarded in the D1-Glu65Ala mutant. In contrast, the S₃ → S₀ transition was significantly inhibited by this mutation. These results suggest that proton transfer in the S₂ → S₃ transition occurs through multiple pathways including the Cl channel, whereas this channel likely serves as a single pathway for proton exit in the S₃ → S₀ transition.

Theoretical Study on the Kinetics of the Rubisco Carboxylase Reaction by a Model Based on Quantum Chemistry and Absolute Reaction Rate Theory

The rate of the Rubisco carboxylase reaction is evaluated by statistical mechanics and hybrid density functional theory (DFT). The Rubisco molecular model given by Kannappan et al. was modified and used in the present calculation. The activation energies of CO₂ addition reaction, H₂O addition reaction, C2-C3 bond scission, and C2 protonation are estimated. We calculated the turnover number (TON) for each of the four reaction steps based on a revised absolute reaction rate theory, which became applicable to soft matter reactions. The molecular parameters used in TON calculations were obtained by DFT calculations. The TON of the total Rubisco reaction was finally evaluated using rate equations. The calculation in a vacuum gave the total TON to be around 5 × 10^-5, which was much lower than the experimental value. The DFT calculation in water solvent gave the total TON to be around 0.1, which agreed reasonably well with experimentally reported values (∼2.71). The rate-limiting process was the scission reaction. The present calculation showed that both the phosphate groups in the substrate accelerate each reaction step. The present calculation showed that a more comprehensive molecular model including enolization and quantum chemical methods is necessary to make a more precise reaction model including the irreversibility of some reactions.

Binding of the substrate analog methanol in the oxygen-evolving complex of photosystem II in the D1-N87A genetic variant of cyanobacteria

The solar water-splitting protein complex, photosystem II (PSII), catalyzes one of the most energetically demanding reactions in nature by using light energy to drive a catalyst capable of oxidizing water. The water oxidation reaction is catalyzed at the Mn₄Ca-oxo cluster in the oxygen-evolving complex (OEC), which cycles through five light-driven S-state intermediates (S₀-S₄). A detailed mechanism of the reaction remains elusive as it requires knowledge of the delivery and binding of substrate water in the higher S-state intermediates. In this study, we use two-dimensional (2D) hyperfine sublevel correlation spectroscopy, in conjunction with quantum mechanics/molecular mechanics (QM/MM) and density functional theory (DFT), to probe the binding of the substrate analog, methanol, in the S₂ state of the D1-N87A variant of PSII from Synechocystis sp. PCC 6803. The results indicate that the size and specificity of the "narrow" channel is altered in D1-N87A PSII, allowing for the binding of deprotonated ¹³C-labeled methanol at the Mn4(IV) ion of the catalytic cluster in the S₂ state. This has important implications on the mechanistic models for water oxidation in PSII.

Structural and dynamic insights into Mn4Ca cluster-depleted Photosystem II

In the first steps of natural oxygenic photosynthesis, sunlight is used to oxidize water molecules to protons, electrons and molecular oxygen. This reaction takes place on the Mn₄Ca cluster located in the reaction centre of Photosystem II (PSII), where the cluster is assembled and continuously repaired through a process known as photoactivation. Understanding the molecular details of such a process has important implications in different fields, in particular inspiring synthesis and repair strategies for artificial photosynthesis devices. In this regard, a detailed structural and dynamic characterization of Photosystem II lacking a Mn₄Ca cluster, namely apo PSII, is a prerequisite for the full comprehension of the photoactivation. Recently, the structure of the apo PSII was resolved at 2.55 Å resolution [Zhang et al., eLife, 2017, 6, e26933], suggesting a pre-organized structure of the protein cavity hosting the cluster. Anyway, the question of whether these findings are a feature of the method used remains open. Here, by means of classical Molecular Dynamics simulations, we characterized the structural and dynamic features of the apo PSII for different protonation states of the cluster cavity. Albeit an overall conformational stability common to all investigated systems, we found significant deviations in the conformation of the side chains of the active site with respect to the X-ray positions. Our findings suggest that not all residues acting as Mn ligands are pre-organized prior to the Mn₄Ca formation and previous local conformational changes are required in order to bind the first Mn ion in the high-affinity binding site.

Impact of Increased Membrane Realism on Conformational Sampling of Proteins

The realism and accuracy of lipid bilayer simulations through molecular dynamics (MD) are heavily dependent on the lipid composition. While the field is pushing toward implementing more heterogeneous and realistic membrane compositions, a lack of high-resolution lipidomic data prevents some membrane protein systems from being modeled with the highest level of realism. Given the additional diversity of real-world cellular membranes and protein-lipid interactions, it is still not fully understood how altering membrane complexity affects modeled membrane protein functions or if it matters over long-timescale simulations. This is especially true for organisms whose membrane environments have little to no computational study, such as the plant plasma membrane. Tackling these issues in tandem, a generalized, realistic, and asymmetric plant plasma membrane with more than 10 different lipid species is constructed herein. Classical MD simulations of pure membrane constructs were performed to evaluate how altering the compositional complexity of the membrane impacted the plant membrane properties. The apo form of a plant sugar transporter, OsSWEET2b, was inserted into membrane models where lipid diversity was calculated in either a size-dependent or size-independent manner. An adaptive sampling simulation regime validated by Markov-state models was performed to capture the gating dynamics of OsSWEET2b in each of these membrane constructs. In comparison to previous OsSWEET2b simulations performed in a pure POPC bilayer, we confirm that simulations performed within a native-like membrane composition alter the stabilization of apo OsSWEET2b conformational states by ∼1 kcal/mol. The free-energy barriers of intermediate conformational states decrease when realistic membrane complexity is simplified, albeit roughly within sampling error, suggesting that protein-specific responses to membranes differ due to altered packing caused by compositional fluctuations. This work serves as a case study where a more realistic bilayer composition makes unbiased conformational sampling easier to achieve than with simplified bilayers.

Proton and Water Transfer Pathways in the S2 → S3 Transition of the Water-Oxidizing Complex in Photosystem II: Time-Resolved Infrared Analysis of the Effects of D1-N298A Mutation and NO3- Substitution

Photosynthetic water oxidation is performed through a light-driven cycle of five intermediates (S₀-S₄ states) in photosystem II (PSII). The S₂ → S₃ transition, which involves concerted water and proton transfer, is a key process for understanding the water oxidation mechanism. Here, to identify the water and proton transfer pathways during the S₂ → S₃ transition, we examined the effects of D1-N298A mutation and NO₃^- substitution for Cl^-, which perturbed the O1 and Cl channels, respectively, on the S₂ → S₃ kinetics using time-resolved infrared spectroscopy. The S₂ → S₃ transition was retarded both upon NO₃^- substitution and upon D1-N298A mutation, whereas it was unaffected by further NO₃^- substitution in N298A PSII. The H/D kinetic isotope effect in N298A PSII was relatively small, revealing that water transfer is a rate-limiting step in this mutant. From these results, it was suggested that during the S₂ → S₃ transition, water delivery and proton release occur through the O1 and Cl channels, respectively.

Interactions Determining the Structural Integrity of the Trimer of Plant Light Harvesting Complex in Lipid Membranes

The structural basis for the stability of the trimeric form of the light harvesting complex (LHCII), a pigmented protein from green plants pivotal for photosynthesis, remains elusive till date. The protein embedded in a dipalmitoylphosphatidylcholine (DPPC) lipid membrane is investigated using all-atom molecular dynamics simulations to find out the interactions responsible for the structural integrity of the trimer and its relation to antenna function. Central association of chlorophyll a (CLA) molecules near the LHCII chains is attributed to a conserved coordination between the Mg of CLA and the oxygen of a specific residue of the first helix of a chain. The residue forms a salt-bridge with the fourth helix of the same chain of the trimer, not of the monomer. In an earlier experiment, three residues (WYR) at each chain of the trimer have been found indispensable for the trimerization and referred to as trimerization motif. We find that the residues of the trimerization motif are connected to the lipids or pigments by a chain of interactions rather than a direct contact. Synergistic effects of sequentially located hydrogen bonds and salt-bridges within monomers of the trimer keep the trimer conformation stable in association with the pigments or the lipids. These interactions are exclusively present in the pigmented trimer and not present in the monomer or in the unpigmented trimer. Thus, our results provide a molecular basis for the inherent stability of the LHCII trimer in a lipid membrane and explain many pre-existing experimental data.

Rigidly hydrogen-bonded water molecules facilitate proton transfer in photosystem II

In the water-splitting enzyme photosystem II (PSII), the proton is released from the catalytic site and transferred to the protein bulk surface via the proton-relay mechanism. Proton transfer occurs in a proton-conducting channel consisting of a series of water molecules connected by hydrogen-bonded (H-bonded) chains. The water-transport protein aquaporin (AQP) also contains a water chain with structure similar to that of the PSII proton channel, although the water chain does not transport protons. We compared the PSII proton channel with the AQP water channel from the following standpoints: (1) the energetics of proton transfer based on crystal structures obtained from quantum mechanical/molecular mechanical calculations, and (2) fluctuations in water molecules obtained from molecular dynamics simulations. The results showed that residues facing the channel and acting as H-bonded partners of water molecules predominantly determined the proton-transfer ability. In PSII, the water chain is surrounded by H-bond acceptors (e.g., carbonyl groups), and the water chain transports protons where the water molecules are rigidly fixed. In AQP, the water chain is surrounded by hydrophobic sidechains or H-bond donors (e.g., NH2 groups), and it does not transport protons where the water molecules are flexible and fluctuating.

Proton transfer pathway from the oxygen-evolving complex in photosystem II substantiated by extensive mutagenesis

We report a structure-based biological approach to identify the proton-transfer pathway in photosystem II. First, molecular dynamics (MD) simulations were conducted to analyze the H-bond network that may serve as a Grotthuss-like proton conduit. MD simulations show that D1-Asp61, the H-bond acceptor of H₂O at the Mn₄CaO₅ cluster (W1), forms an H-bond via one water molecule with D1-Glu65 but not with D2-Glu312. Then, D1-Asp61, D1-Glu65, D2-Glu312, and the adjacent residues, D1-Arg334, D2-Glu302, and D2-Glu323, were thoroughly mutated to the other 19 residues, i.e., 114 Chlamydomonas chloroplast mutant cells were generated. Mutation of D1-Asp61 was most crucial. Only the D61E and D61C cells grew photoautotrophically and exhibit O₂-evolving activity. Mutations of D2-Glu312 were less crucial to photosynthetic growth than mutations of D1-Glu65. Quantum mechanical/molecular mechanical calculations indicated that in the PSII crystal structure, the proton is predominantly localized at D1-Glu65 along the H-bond with D2-Glu312, i.e., pK_a(D1-Glu65) > pK_a(D2-Glu312). The potential-energy profile shows that the release of the proton from D1-Glu65 leads to the formation of the two short H-bonds between D1-Asp61 and D1-Glu65, which facilitates downhill proton transfer along the Grotthuss-like proton conduit in the S₂ to S₃ transition. It seems possible that D1-Glu65 is involved in the dominant pathway that proceeds from W1 via D1-Asp61 toward the thylakoid lumen, whereas D2-Glu312 and D1-Arg334 may be involved in alternative pathways in some mutants.

Molecular dynamics simulations in photosynthesis

Photosynthesis is regulated by a dynamic interplay between proteins, enzymes, pigments, lipids, and cofactors that takes place on a large spatio-temporal scale. Molecular dynamics (MD) simulations provide a powerful toolkit to investigate dynamical processes in (bio)molecular ensembles from the (sub)picosecond to the (sub)millisecond regime and from the Å to hundreds of nm length scale. Therefore, MD is well suited to address a variety of questions arising in the field of photosynthesis research. In this review, we provide an introduction to the basic concepts of MD simulations, at atomistic and coarse-grained level of resolution. Furthermore, we discuss applications of MD simulations to model photosynthetic systems of different sizes and complexity and their connection to experimental observables. Finally, we provide a brief glance on which methods provide opportunities to capture phenomena beyond the applicability of classical MD.

Role of the Photosystem II as an Environment in the Oxidation Free Energy of the Mn Cluster from S1 to S2

The manganese cluster (CaMn₄O₅) in the photosystem II (PSII) is the reaction center of the light-driven oxidation reaction, which generates the molecular oxygen. In this paper, we address the issue of the effect of the environment on the free energy associated with the oxidation of the Mn cluster in S₁ state by conducting the large-scale quantum mechanical/molecular mechanical simulations, which involve the whole of the PSII monomer. It was found by the simulations at the level of the B3LYP functional that the environment surrounding the Mn cluster reduces the vertical oxidation free energy Δμ_vrt by 64.8 kcal/mol. A decomposition analysis of the free energy Δμ_vrt revealed that the system composed of peptide chains, ligands, lipids, and potassium ions contributes to lowering of Δμ_vrt by -98.0 kcal/mol, whereas the solvent water makes an opposite contribution of 38.9 kcal/mol. Reduction of the vertical oxidation free energy directly leads to the lowering of the activation free energy ΔG_ac for the electron transfer reaction from the Mn cluster in S₁ state to the neighboring Tyr_z⁺. Consequently, the electron transfer rate was found to be enhanced by a factor of 10¹² by virtue of the influence of the environment.

Water Network Dynamics Next to the Oxygen-Evolving Complex of Photosystem II

The influence of the environment on the functionality of the oxygen-evolving complex (OEC) of photosystem II has long been a subject of great interest. In particular, various water channels, which could serve as pathways for substrate water diffusion, or proton translocation, are thought to be critical to catalytic performance of the OEC. Here, we address the dynamical nature of hydrogen bonding along the water channels by performing molecular dynamics (MD) simulations of the OEC and its surrounding protein environment in the S1 and S2 states. Through the eigenvector centrality (EC) analysis, we are able to determine the characteristics of the water network and assign potential functions to the major channels, namely that the narrow and broad channels are likely candidates for proton/water transport, while the large channel may serve as a path for larger ions such as chloride and manganese thought to be essential during PSII assembly.

A theoretical study on the dynamics of light harvesting in the dimeric photosystem II core complex: regulation and robustness of energy transfer pathways

Here we present our theoretical investigations into the light reaction in the dimeric photosystem II (PSII) core complex. An effective model for excitation energy transfer (EET) and primary charge separation (CS) in the PSII core complex was developed, with model parameters constructed based on molecular dynamics (MD) simulation data. Compared to experimental results, we demonstrated that this model faithfully reproduces the absorption spectra of the RC and core light-harvesting complexes (CP43 and CP47) as well as the full EET dynamics among the chromophores in the PSII core complex. We then applied master equation simulations and network analysis to investigate detailed EET plus CS dynamics in the system, allowing us to identify key EET pathways and produce a coarse-grained cluster model for the light reaction in the dimeric PSII core complex. We show that non-equilibrium energy transfer channels play important roles in the efficient light harvesting process and that multiple EET pathways exist between subunits of PSII to ensure the robustness of light harvesting in the system. Furthermore, we revealed that inter-monomer energy transfer dominated by the coupling between the two CLA625 molecules enables efficient energy exchange between two CP47s in the dimeric PSII core complex, which leads to significant energy pooling in the CP47 domain during the light reaction. Our study provides a blueprint for the design of light harvesting in the PSII core and show that a structure-based approach using molecular dynamics simulations and quantum chemistry calculations can be effectively utilized to elucidate the dynamics of light harvesting in complex photosynthetic systems.

First principles calculations of surface dependent electronic structures: a study on β-FeOOH and γ-FeOOH

The surface dependent electronic structures of β-FeOOH and γ-FeOOH were studied using density functional calculations.

Fourier Transform Infrared Analysis of the S-State Cycle of Water Oxidation in the Microcrystals of Photosystem II

Photosynthetic water oxidation is performed in photosystem II (PSII) through a light-driven cycle of intermediates called S states (S₀-S₄) at the water oxidizing center. Time-resolved serial femtosecond crystallography (SFX) has recently been applied to the microcrystals of PSII to obtain the structural information on these intermediates. However, it remains unanswered whether the reactions efficiently proceed throughout the S-state cycle retaining the native structures of the intermediates in PSII crystals. We investigated the water oxidation reactions in the PSII microcrystals using flash-induced Fourier transform infrared (FTIR) difference spectroscopy. In comparison with the FTIR spectra in solution, it was shown that all of the metastable intermediates in the microcrystals retained their native structures, and the efficiencies of the S-state transitions remained relatively high, although those of the S₂ → S₃ and S₃ → S₀ transitions were slightly lowered possibly due to some restriction of water movement in the crystals.

Element strategy of oxygen evolution electrocatalysis based on in situ spectroelectrochemistry

Oxygen evolution electrocatalysis has received extensive attention due to its significance in biology, chemistry, and technology. However, it is still unclear how the abundant 3d-elements can be used to drive the four-electron oxidation of water as efficiently as in Nature. In this Feature Article, we will propose a design strategy concerning the optimization of the charge accumulation process based on our ongoing spectroelectrochemical study on Mn, Fe, and Ir oxygen evolution catalysts. Spectroscopic identification of the reaction intermediates showed that the activity of MnO₂ and Fe₂O₃ was dictated by the generation of Mn³⁺ and Fe⁴⁺, whereas in the case of IrO_x, the activity did not correlate with the valence change of Ir. The efficiency of charge accumulation through valence change is closely linked with the spin configuration of the metal center, because charge disproportionation, which was found to inhibit charge accumulation in the high-spin 3d metals, requires an electron in the e_g orbital. In addition to directly increasing the overpotential through the generation of an unstable intermediate, charge disproportionation inhibits charge accumulation by dissipating the total oxidative energy of the system. A favorable charge accumulation process may also be beneficial for electrode kinetics due to the enhanced coupling between reaction rates and electrochemical driving force. The model proposed in this study may help explain why low-spin 4d/5d rare metals are often more active than the abundant high-spin 3d materials for multi-electron transfer reactions in general, and provides new insight into how active 3d-metal catalysts can be synthesized by optimizing the energetics of both bond formation and charge accumulation.

In vivo NMR as a tool for probing molecular structure and dynamics in intact Chlamydomonas reinhardtii cells

We report the application of NMR dynamic spectral editing for probing the structure and dynamics of molecular constituents in fresh, intact cells and in freshly prepared thylakoid membranes of Chlamydomonas reinhardtii (Cr.) green algae. For isotope labeling, wild-type Cr. cells were grown on 13C acetate-enriched minimal medium. 1D 13C J-coupling based and dipolar-based MAS NMR spectra were applied to distinguish 13C resonances of different molecular components. 1D spectra were recorded over a physiological temperature range, and whole-cell spectra were compared to those taken from thylakoid membranes, evaluating their composition and dynamics. A theoretical model for NMR polarization transfer was used to simulate the relative intensities of direct, J-coupling, and dipolar-based polarization from which the degree of lipid segmental order and rotational dynamics of the lipid acyl chains were estimated. We observe that thylakoid lipid signals dominate the lipid spectral profile of whole algae cells, demonstrating that with our novel method, thylakoid membrane characteristics can be detected with atomistic precision inside intact photosynthetic cells. The experimental procedure is rapid and applicable to fresh cell cultures, and could be used as an original approach for detecting chemical profiles, and molecular structure and dynamics of photosynthetic membranes in vivo in functional states.

Prediction of Thylakoid Lipid Binding Sites on Photosystem II

The thylakoid membrane has a unique lipid composition, consisting mostly of galactolipids. These thylakoid lipids have important roles in photosynthesis. Here, we investigate to what extent these lipids bind specifically to the Photosystem II complex. To this end, we performed coarse-grain MD simulations of the Photosystem II complex embedded in a thylakoid membrane with realistic composition. Based on >85 μs simulation time, we find that monogalactosyldiacylglycerol and sulfoquinovosyldiacylglycerol lipids are enriched in the annular shell around the protein, and form distinct binding sites. From the analysis of residue contacts, we conclude that electrostatic interactions play an important role in stabilizing these binding sites. Furthermore, we find that chlorophyll a has a prevalent role in the coordination of the lipids. In addition, we observe lipids to diffuse in and out of the plastoquinone exchange cavities, allowing exchange of cocrystallized lipids with the bulk membrane and suggesting a more open nature of the plastoquinone exchange cavity. Together, our data provide a wealth of information on protein-lipid interactions for a key protein in photosynthesis.

Fine control of chlorophyll-carotenoid interactions defines the functionality of light-harvesting proteins in plants

Photosynthetic antenna proteins can be thought of as "programmed solvents", which bind pigments at specific mutual orientations, thus tuning the overall energetic landscape and ensuring highly efficient light-harvesting. While positioning of chlorophyll cofactors is well understood and rationalized by the principle of an "energy funnel", the carotenoids still pose many open questions. Particularly, their short excited state lifetime (<25 ps) renders them potential energy sinks able to compete with the reaction centers and drastically undermine light-harvesting efficiency. Exploration of the orientational phase-space revealed that the placement of central carotenoids minimizes their interaction with the nearest chlorophylls in the plant antenna complexes LHCII, CP26, CP29 and LHCI. At the same time we show that this interaction is highly sensitive to structural perturbations, which has a profound effect on the overall lifetime of the complex. This links the protein dynamics to the light-harvesting regulation in plants by the carotenoids.

D1-Asn-298 in photosystem II is involved in a hydrogen-bond network near the redox-active tyrosine YZ for proton exit during water oxidation

In photosynthetic water oxidation, two water molecules are converted into one oxygen molecule and four protons at the Mn₄CaO₅ cluster in photosystem II (PSII) via the S-state cycle. Efficient proton exit from the catalytic site to the lumen is essential for this process. However, the exit pathways of individual protons through the PSII proteins remain to be identified. In this study, we examined the involvement of a hydrogen-bond network near the redox-active tyrosine Y_Z in proton transfer during the S-state cycle. We focused on spectroscopic analyses of a site-directed variant of D1-Asn-298, a residue involved in a hydrogen-bond network near Y_Z We found that the D1-N298A mutant of Synechocystis sp. PCC 6803 exhibits an O₂ evolution activity of ∼10% of the wild-type. D1-N298A and the wild-type D1 had very similar features of thermoluminescence glow curves and of an FTIR difference spectrum upon Y_Z oxidation, suggesting that the hydrogen-bonded structure of Y_Z and electron transfer from the Mn₄CaO₅ cluster to Y_Z were little affected by substitution. In the D1-N298A mutant, however, the flash-number dependence of delayed luminescence showed a monotonic increase without oscillation, and FTIR difference spectra of the S-state cycle indicated partial and significant inhibition of the S₂ → S₃ and S₃ → S₀ transitions, respectively. These results suggest that the D1-N298A substitution inhibits the proton transfer processes in the S₂ → S₃ and S₃ → S₀ transitions. This in turn indicates that the hydrogen-bond network near Y_Z can be functional as a proton transfer pathway during photosynthetic water oxidation.

Differences in the Active Site of Water Oxidation among Photosynthetic Organisms

The site of biological water oxidation is highly conserved across photosynthetic organisms, but differences of unidentified structural and electronic origin exist between taxonomically discrete clades, revealed by distinct spectroscopic signatures of the oxygen-evolving Mn₄CaO₅ cluster and variations in active-site accessibility. Comparison of atomistic models of a native cyanobacterial form (Thermosynechococcus vulcanus) and a chimeric spinach-like form of photosystem II allows us to identify the precise atomic-level differences between organisms in the vicinity of the manganese cluster. Substitution of cyanobacterial D1-Asn87 by higher-plant D1-Ala87 is the principal discriminating feature: it drastically rearranges a network of proximal hydrogen bonds, modifying the local architecture of a water channel and the interaction of second coordination shell residues with the manganese cluster. The two variants explain species-dependent differences in spectroscopic properties and in the interaction of substrate analogues with the oxygen-evolving complex, enabling assignment of a substrate delivery channel to the active site.

Structurally conserved channels in cyanobacterial and plant photosystem II

In the cyanobacterial photosystem II (PSII), the O4-water chain in the D1 and CP43 proteins, a chain of water molecules that are directly H-bonded to O4 of the Mn₄Ca cluster, is linked with a channel that connects the protein bulk surface along with a membrane-extrinsic protein subunit, PsbU (O4-PsbU channel). The cyanobacterial PSII structure also shows that the O1 site of the Mn₄Ca cluster has a chain of H-bonded water molecules, which is linked with the channel that proceeds toward the bulk surface via PsbU and PsbV (O1-PsbU/V channel). Membrane-extrinsic protein subunits PsbU and PsbV in cyanobacterial PSII are replaced with PsbP and PsbQ in plant PSII. However, these four proteins have no structural similarity. It remains unknown whether the corresponding channels also exist in plant PSII, because water molecules are not identified in the plant PSII cryo-electron microscopy (cryo-EM) structure. Using the cyanobacterial and plant PSII structures, we analyzed the channels that proceed from the Mn₄Ca cluster. The cyanobacterial O4-PsbU and O1-PsbU/V channels were structurally conserved as the channel that proceeds along PsbP toward the protein bulk surface in the plant PSII (O4-PsbP and O1-PsbP channels, respectively). Calculated protonation states indicated that in contrast to the original geometry of the plant cryo-EM structure, protonated PsbP-Lys166 may form a salt-bridge with ionized D1-Glu329 and protonated PsbP-Lys173 may form a salt-bridge with ionized PsbQ-Asp28 near the O1-PsbP channel. The existence of these channels might explain the molecular mechanism of how PsbP can interact with the Mn₄Ca cluster.

Water permeation through the internal water pathway in activated GPCR rhodopsin

Rhodopsin is a light-driven G-protein-coupled receptor that mediates signal transduction in eyes. Internal water molecules mediate activation of the receptor in a rhodopsin cascade reaction and contribute to conformational stability of the receptor. However, it remains unclear how internal water molecules exchange between the bulk and protein inside, in particular through a putative solvent pore on the cytoplasmic. Using all-atom molecular dynamics simulations, we identified the solvent pore on cytoplasmic side in both the Meta II state and the Opsin. On the other hand, the solvent pore does not exist in the dark-adapted rhodopsin. We revealed two characteristic narrow regions located within the solvent pore in the Meta II state. The narrow regions distinguish bulk and the internal hydration sites, one of which is adjacent to the conserved structural motif "NPxxY". Water molecules in the solvent pore diffuse by pushing or sometimes jumping a preceding water molecule due to the geometry of the solvent pore. These findings revealed a total water flux between the bulk and the protein inside in the Meta II state, and suggested that these pathways provide water molecules to the crucial sites of the activated rhodopsin.

An Ab Initio Description of the Excitonic Properties of LH2 and Their Temperature Dependence

The spectroscopic properties of light-harvesting (LH) antennae in photosyntehtic organisms represent a fingerprint that is unique for each specific pigment-protein complex. Because of that, spectroscopic observations are generally combined with structural data from X-ray crystallography to obtain an indirect representation of the excitonic properties of the system. Here, an alternative strategy is presented which goes beyond this empirical approach and introduces an ab initio computational description of both structural and electronic properties and their dependence on the temperature. The strategy is applied to the peripheral light-harvesting antenna complex (LH2) present in purple bacteria. By comparing this model with the one based on the crystal structure, a detailed, molecular level explanation of the absorption and circular dichroism (CD) spectra and their temperature dependence is achieved. The agreement obtained with the experiments at both low and room temperature lays the groundwork for an atomistic understanding of the excitation dynamics in the LH2 system.

Molecular Dynamics of Photosystem II Embedded in the Thylakoid Membrane

Photosystem II (PSII) is one of the key protein complexes in photosynthesis. We introduce a coarse grained model of PSII and present the analysis of 60 μs molecular dynamics simulations of PSII in both monomeric and dimeric form, embedded in a thylakoid membrane model that reflects its native lipid composition. We describe in detail the setup of the protein complex and the many natural cofactors and characterize their mobility. Overall we find that the protein subunits and cofactors are more flexible toward the periphery of the complex as well as near the PLQ exchange cavity and at the dimer interface. Of all cofactors, β-carotenes show the highest mobility. Some of the β-carotenes diffuse in and out of the protein complex via the thylakoid membrane. In contrast with the PSII dimer, the monomeric form adopts a tilted conformation in the membrane, with strong interactions between the soluble PsbO subunit and the glycolipid headgroups. Interestingly, the tilted conformation causes buckling of the membrane. Together, our results provide an unprecedented view of PSII dynamics on a microsecond time scale. Our data may be used as basis for the interpretation of experimental data as well as for theoretical models describing exciton energy transfer.

Interaction of methanol with the oxygen-evolving complex: atomistic models, channel identification, species dependence, and mechanistic implications

Methanol has long being used as a substrate analogue to probe access pathways and investigate water delivery at the oxygen-evolving complex (OEC) of photosystem-II. In this contribution we study the interaction of methanol with the OEC by assembling available spectroscopic data into a quantum mechanical treatment that takes into account the local channel architecture of the active site. The effect on the magnetic energy levels of the Mn₄Ca cluster in the S₂ state of the catalytic cycle can be explained equally well by two models that involve either methanol binding to the calcium ion of the cluster, or a second-sphere interaction in the vicinity of the "dangler" Mn4 ion. However, consideration of the latest ¹³C hyperfine interaction data shows that only one model is fully consistent with experiment. In contrast to previous hypotheses, methanol is not a direct ligand to the OEC, but is situated at the end-point of a water channel associated with the O4 bridge. Its effect on magnetic properties of plant PS-II results from disruption of hydrogen bonding between O4 and proximal channel water molecules, thus enhancing superexchange (antiferromagnetic coupling) between the Mn3 and Mn4 ions. The same interaction mode applies to the dark-stable S₁ state and possibly to all other states of the complex. Comparison of protein sequences from cyanobacteria and plants reveals a channel-altering substitution (D1-Asn87 versus D1-Ala87) in the proximity of the methanol binding pocket, explaining the species-dependence of the methanol effect. The water channel established as the methanol access pathway is the same that delivers ammonia to the Mn4 ion, supporting the notion that this is the only directly solvent-accessible manganese site of the OEC. The results support the pivot mechanism for water binding at a component of the S₃ state and would be consistent with partial inhibition of water delivery by methanol. Mechanistic implications for enzymatic regulation and catalytic progression are discussed.

pK(a) of a Proton-Conducting Water Chain in Photosystem II

Recent high-resolution crystal structures of the water-oxidizing enzyme photosystem II (PSII) show that O4 of the catalytic Mn4CaO5 cluster forms an H-bond with a water molecule W539, which belongs to a chain of water molecules (O4-water chain). Oxidation of Mn4CaO5 to S1 resulted in elongation of the O-H bonds and decrease in pKa(O-H/O(-)) in the [O4-H···OW539-H···OW538-H···OW393] region along the O4-water chain. In S1, removal of all water molecules from the O4-water chain, except W539, resulted in a significant pKa upshift at O4; this suggests that the proton-conducting water chain serves as a conducting media for protons and significantly decreases the donor pKa, leading to a downhill proton transfer. The absence of a corresponding proton-conducting channel is disadvantageous for release of protons from the proton-releasing site, as in the case of O5 that has no H-bond partner.

Combining classical molecular dynamics and quantum mechanical methods for the description of electronic excitations: The case of carotenoids

Carotenoids are important actors both in light-harvesting (LH) and in photoprotection functions of photosynthetic pigment-protein complexes. A deep theoretical investigation of this multiple role is still missing owing to the difficulty of describing the delicate interplay between electronic and nuclear degrees of freedom. A possible strategy is to combine accurate quantum mechanical (QM) methods with classical molecular dynamics. To do this, however, accurate force-fields (FF) are necessary. This article presents a new FF for the different carotenoids present in LH complexes of plants. The results show that all the important structural properties described by the new FF are in very good agreement with QM reference values. This increased accuracy in the simulation of the structural fluctuations is also reflected in the description of excited states. Both the energy order and the different nature of the lowest singlet states are preserved during the dynamics when the new FF is used, whereas an unphysical mixing is found when a standard FF is used.