QM/MM computational studies of substrate water binding to the oxygen-evolving centre of photosystem II

Glycerol binding at the narrow channel of photosystem II stabilizes the low-spin S2 state of the oxygen-evolving complex

The oxygen-evolving complex (OEC) of photosystem II (PSII) cycles through redox intermediate states S_i (i = 0-4) during the photochemical oxidation of water. The S₂ state involves an equilibrium of two isomers including the low-spin S₂ (LS-S₂) state with its characteristic electron paramagnetic resonance (EPR) multiline signal centered at g = 2.0, and a high-spin S₂ (HS-S₂) state with its g = 4.1 EPR signal. The relative intensities of the two EPR signals change under experimental conditions that shift the HS-S₂/LS-S₂ state equilibrium. Here, we analyze the effect of glycerol on the relative stability of the LS-S₂ and HS-S₂ states when bound at the narrow channel of PSII, as reported in an X-ray crystal structure of cyanobacterial PSII. Our quantum mechanics/molecular mechanics (QM/MM) hybrid models of cyanobacterial PSII show that the glycerol molecule perturbs the hydrogen-bond network in the narrow channel, increasing the pK_a of D1-Asp61 and stabilizing the LS-S₂ state relative to the HS-S₂ state. The reported results are consistent with the absence of the HS-S₂ state EPR signal in native cyanobacterial PSII EPR spectra and suggest that the narrow water channel hydrogen-bond network regulates the relative stability of OEC catalytic intermediates during water oxidation.

High-resolution cryo-electron microscopy structure of photosystem II from the mesophilic cyanobacterium, Synechocystis sp. PCC 6803

Photosystem II (PSII) enables global-scale, light-driven water oxidation. Genetic manipulation of PSII from the mesophilic cyanobacterium Synechocystis sp. PCC 6803 has provided insights into the mechanism of water oxidation; however, the lack of a high-resolution structure of oxygen-evolving PSII from this organism has limited the interpretation of biophysical data to models based on structures of thermophilic cyanobacterial PSII. Here, we report the cryo-electron microscopy structure of PSII from Synechocystis sp. PCC 6803 at 1.93-Å resolution. A number of differences are observed relative to thermophilic PSII structures, including the following: the extrinsic subunit PsbQ is maintained, the C terminus of the D1 subunit is flexible, some waters near the active site are partially occupied, and differences in the PsbV subunit block the Large (O1) water channel. These features strongly influence the structural picture of PSII, especially as it pertains to the mechanism of water oxidation.

The open-cubane oxo–oxyl coupling mechanism dominates photosynthetic oxygen evolution: a comprehensive DFT investigation on O–O bond formation in the S4state

How is O₂created in nature? Comprehensive DFT investigations determine the dominance of the open-cubane oxo–oxyl coupling mechanism over alternative possibilities.

Photosynthetic water oxidation: binding and activation of substrate waters for O-O bond formation

Photosynthetic water oxidation occurs at the oxygen-evolving complex (OEC) of Photosystem II (PSII). The OEC, which contains a Mn4CaO5 inorganic cluster ligated by oxides, waters and amino-acid residues, cycles through five redox intermediates known as S(i) states (i = 0-4). The electronic and structural properties of the transient S4 intermediate that forms the O-O bond are not well understood. In order to gain insight into how water is activated for O-O bond formation in the S4 intermediate, we have performed a detailed analysis of S-state dependent substrate water binding kinetics taking into consideration data from Mn coordination complexes. This analysis supports a model in which the substrate waters are both bound as terminal ligands and react via a water-nucleophile attack mechanism.

Energetics of proton release on the first oxidation step in the water-oxidizing enzyme

In photosystem II (PSII), the Mn₄CaO₅ cluster catalyses the water splitting reaction. The crystal structure of PSII shows the presence of a hydrogen-bonded water molecule directly linked to O4. Here we show the detailed properties of the H-bonds associated with the Mn₄CaO₅ cluster using a quantum mechanical/molecular mechanical approach. When O4 is taken as a μ-hydroxo bridge acting as a hydrogen-bond donor to water539 (W539), the S₀ redox state best describes the unusually short O4–O_W539 distance (2.5 Å) seen in the crystal structure. We find that in S₁, O4 easily releases the proton into a chain of eight strongly hydrogen-bonded water molecules. The corresponding hydrogen-bond network is absent for O5 in S₁. The present study suggests that the O4-water chain could facilitate the initial deprotonation event in PSII. This unexpected insight is likely to be of real relevance to mechanistic models for water oxidation.

The availability of crystal structures of photosystem II opens up the possibility of gaining insights into its mechanism. Here, the authors use a computational approach and propose a deprotonation event at O4 followed by long-range proton-transfer along a chain of strongly bonded water molecules.

Photosynthetic water oxidation: insights from manganese model chemistry

Catalysts for light-driven water oxidation are a critical component for development of solar fuels technology. The multielectron redox chemistry required for this process has been successfully deployed on a global scale in natural photosynthesis by green plants and cyanobacteria using photosystem II (PSII). PSII employs a conserved, cuboidal Mn4CaOX cluster called the O2-evolving complex (OEC) that offers inspiration for artificial O2-evolution catalysts. In this Account, we describe our work on manganese model chemistry relevant to PSII, particularly the functional model [Mn(III/IV)2(terpy)2(μ-O)2(OH2)2](NO3)3 complex (terpy = 2,2';6',2″-terpyridine), a mixed-valent di-μ-oxo Mn dimer with two terminal aqua ligands. In the presence of oxo-donor oxidants such as HSO5(-), this complex evolves O2 by two pathways, one of which incorporates solvent water in an O-O bond-forming reaction. Deactivation pathways of this catalyst include comproportionation to form an inactive Mn(IV)Mn(IV) dimer and also degradation to MnO2, a consequence of ligand loss when the oxidation state of the complex is reduced to labile Mn(II) upon release of O2. The catalyst's versatility has been shown by its continued catalytic activity after direct binding to the semiconductor titanium dioxide. In addition, after binding to the surface of TiO2 via a chromophoric linker, the catalyst can be oxidized by a photoinduced electron-transfer mechanism, mimicking the natural PSII process. Model oxomanganese complexes have also aided in interpreting biophysical and computational studies on PSII. In particular, the μ-oxo exchange rates of the Mn-terpy dimer have been instrumental in establishing that the time scale for μ-oxo exchange of high-valent oxomanganese complexes with terminal water ligands is slower than O2 evolution in the natural photosynthetic system. Furthermore, computational studies on the Mn-terpy dimer and the OEC point to similar Mn(IV)-oxyl intermediates in the O-O bond-forming mechanism. Comparison between the OEC and the Mn-terpy dimer indicates that challenges remain in the development of synthetic Mn water-oxidation catalysts. These include redox leveling to couple multielectron reactions with one-electron steps, avoiding labile Mn(II) species during the catalytic cycle, and protecting the catalyst active site from undesired side reactions. As the first example of a functional manganese O2-evolution catalyst, the Mn-terpy dimer exemplifies the interrelatedness of biomimetic chemistry with biophysical studies. The design of functional model complexes enriches the study of the natural photosynthetic system, while biology continues to provide inspiration for artificial photosynthetic technologies to meet global energy demand.

Alternative mechanisms for O2release and O–O bond formation in the oxygen evolving complex of photosystem II

A new transition state for O₂release has been found. An alternative, nearly competitive, mechanism for O–O bond formation is described, which is very similar to the one previously suggested.

Substrate-water exchange in photosystem II is arrested before dioxygen formation

Light-driven oxidation of water into dioxygen, catalysed by the oxygen-evolving complex (OEC) in photosystem II, is essential for life on Earth and provides the blueprint for devices for producing fuel from sunlight. Although the structure of the OEC is known at atomic level for its dark-stable state, the mechanism by which water is oxidized remains unsettled. Important mechanistic information was gained in the past two decades by mass spectrometric studies of the H₂¹⁸O/H₂¹⁶O substrate–water exchange in the four (semi) stable redox states of the OEC. However, until now such data were not attainable in the transient states formed immediately before the O–O bond formation. Using modified photosystem II complexes displaying up to 40-fold slower O₂ production rates, we show here that in the transient state the substrate–water exchange is dramatically slowed as compared with the earlier S states. This further constrains the possible sites for substrate–water binding in photosystem II.

The oxygen-evolving complex of photosystem II converts water into oxygen during photosynthesis, but how this process occurs is not yet fully understood. Here, the authors use modified complexes with reduced reaction rates to study the process of oxygen evolution in more detail.

The structure and activation of substrate water molecules in Sr2+-substituted photosystem II

An EPR spectroscopy study with direct evidence that the Ca²⁺ ion plays a structural role in maintaining the hydrogen-bond network in photosystem II.

Competitive oxygen-18 kinetic isotope effects expose O–O bond formation in water oxidation catalysis by monomeric and dimeric ruthenium complexes

Competitive ¹⁸O KIEs on water oxidation catalysis provide a probe of transition states for O–O bond formation.

Convergence of QM/MM and Cluster Models for the Spectroscopic Properties of the Oxygen-Evolving Complex in Photosystem II

The latest crystal structure of photosystem II at 1.9 Å resolution, which resolves the topology of the Mn4CaO5 oxygen evolving complex (OEC) at atomistic detail, enables a better correlation between structural features and spectroscopic properties than ever before. Building on the refined crystallographic model of the OEC and the protein, we present combined quantum mechanical/molecular mechanical (QM/MM) studies of the spectroscopic properties of the natural catalyst embedded in the protein matrix. Focusing on the S2 state of the catalytic cycle, we examine the convergence of not only structural parameters but also of the intracluster magnetic interactions in terms of exchange coupling constants and of experimentally relevant (55)Mn, (17)O, and (14)N hyperfine coupling constants with respect to QM/MM partitioning using five QM regions of increasing size. This enables us to assess the performance of the method and to probe second sphere effects by identifying amino acid residues that principally affect the spectroscopic properties of the OEC. Comparison between QM-only and QM/MM treatments reveals that whereas QM/MM models converge quickly to stable values, the QM cluster models need to incorporate significantly larger parts of the second coordination sphere and surrounding water molecules to achieve convergence for certain properties. This is mainly due to the sensitivity of the QM-only models to fluctuations in the hydrogen bonding network and ligand acidity. Additionally, a hydrogen bond that is typically omitted in QM-only treatments is shown to determine the hyperfine coupling tensor of the unique Mn(III) ion by regulating the rotation plane of the ligated D1-His332 imidazole ring, the only N-donor ligand of the OEC.

A reaction center-dependent photoprotection mechanism in a highly robust photosystem II from an extremophilic red alga, Cyanidioschyzon merolae

Members of the rhodophytan order Cyanidiales are unique among phototrophs in their ability to live in extremely low pH levels and moderately high temperatures. The photosynthetic apparatus of the red alga Cyanidioschyzon merolae represents an intermediate type between cyanobacteria and higher plants, suggesting that this alga may provide the evolutionary link between prokaryotic and eukaryotic phototrophs. Although we now have a detailed structural model of photosystem II (PSII) from cyanobacteria at an atomic resolution, no corresponding structure of the eukaryotic PSII complex has been published to date. Here we report the isolation and characterization of a highly active and robust dimeric PSII complex from C. merolae. We show that this complex is highly stable across a range of extreme light, temperature, and pH conditions. By measuring fluorescence quenching properties of the isolated C. merolae PSII complex, we provide the first direct evidence of pH-dependent non-photochemical quenching in the red algal PSII reaction center. This type of quenching, together with high zeaxanthin content, appears to underlie photoprotection mechanisms that are efficiently employed by this robust natural water-splitting complex under excess irradiance. In order to provide structural details of this eukaryotic form of PSII, we have employed electron microscopy and single particle analyses to obtain a 17 Å map of the C. merolae PSII dimer in which we locate the position of the protein mass corresponding to the additional extrinsic protein stabilizing the oxygen-evolving complex, PsbQ'. We conclude that this lumenal subunit is present in the vicinity of the CP43 protein, close to the membrane plane.

Reflections on substrate water and dioxygen formation

This brief article aims at presenting a concise summary of all experimental findings regarding substrate water-binding to the Mn4CaO5 cluster in photosystem II. Mass spectrometric and spectroscopic results are interpreted in light of recent structural information of the water oxidizing complex obtained by X-ray crystallography, spectroscopy and theoretical modeling. Within this framework current proposals for the mechanism of photosynthetic water-oxidation are evaluated. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

Similarities of artificial photosystems by ruthenium oxo complexes and native water splitting systems

The nature of chemical bonds of ruthenium(Ru)-quinine(Q) complexes, mononuclear [Ru(trpy)(3,5-t-Bu(2)Q)(OH(2))](ClO(4))(2) (trpy = 2,2':6',2''-terpyridine, 3,5-di-tert-butyl-1,2-benzoquinone) (1), and binuclear [Ru(2)(btpyan)(3,6-di-Bu(2)Q)(2)(OH(2))](2+) (btpyan = 1,8-bis(2,2':6',2''-terpyrid-4'-yl)anthracene, 3,6-t-Bu(2)Q = 3,6-di-tert-butyl-1,2-benzoquinone) (2), has been investigated by broken-symmetry (BS) hybrid density functional (DFT) methods. BS DFT computations for the Ru complexes have elucidated that the closed-shell structure (2b) Ru(II)-Q complex is less stable than the open-shell structure (2bb) consisting of Ru(III) and semiquinone (SQ) radical fragments. These computations have also elucidated eight different electronic and spin structures of tetraradical intermediates that may be generated in the course of water splitting reaction. The Heisenberg spin Hamiltonian model for these species has been derived to elucidate six different effective exchange interactions (J) for four spin systems. Six J values have been determined using total energies of the eight (or seven) BS solutions for different spin configurations. The natural orbital analyses of these BS DFT solutions have also been performed in order to obtain natural orbitals and their occupation numbers, which are useful for the lucid understanding of the nature of chemical bonds of the Ru complexes. Implications of the computational results are discussed in relation to the proposed reaction mechanisms of water splitting reaction in artificial photosynthesis systems and the similarity between artificial and native water splitting systems.

Oxomanganese complexes for natural and artificial photosynthesis

The oxygen-evolving complex (OEC) of Photosystem II (PSII) is an oxomanganese complex that catalyzes water-splitting into O2, protons and electrons. Recent breakthroughs in X-ray crystallography have resolved the cuboidal OEC structure at 1.9 Å resolution, stimulating significant interest in studies of structure/function relations. This article summarizes recent advances on studies of the OEC along with studies of synthetic oxomanganese complexes for artificial photosynthesis. Quantum mechanics/molecular mechanics hybrid methods have enabled modeling the S1 state of the OEC, including the ligation proposed by the most recent X-ray data where D170 is bridging Ca and the Mn center outside the CaMn3 core. Molecular dynamics and Monte Carlo simulations have explored the structural/functional roles of chloride, suggesting that it regulates the electrostatic interactions between D61 and K317 that might be critical for proton abstraction. Furthermore, structural studies of synthetic oxomanganese complexes, including the [H2O(terpy)MnIII(μ-O)2MnIV(terpy)OH2]3+ (1, terpy=2,2':6',2″-terpyridine) complex, provided valuable insights on the mechanistic influence of carboxylate moieties in close contact with the Mn catalyst during oxygen evolution. Covalent attachment of 1 to TiO2 has been achieved via direct deposition and by using organic chromophoric linkers. The (III,IV) oxidation state of 1 attached to TiO2 can be advanced to (IV,IV) by visible-light photoexcitation, leading to photoinduced interfacial electron transfer. These studies are particularly relevant to the development of artificial photosynthetic devices based on inexpensive materials.

Theoretical evaluation of structural models of the S2 state in the oxygen evolving complex of Photosystem II: protonation states and magnetic interactions

Protonation states of water ligands and oxo bridges are intimately involved in tuning the electronic structures and oxidation potentials of the oxygen evolving complex (OEC) in Photosystem II, steering the mechanistic pathway, which involves at least five redox state intermediates S(n) (n = 0-4) resulting in the oxidation of water to molecular oxygen. Although protons are practically invisible in protein crystallography, their effects on the electronic structure and magnetic properties of metal active sites can be probed using spectroscopy. With the twin purpose of aiding the interpretation of the complex electron paramagnetic resonance (EPR) spectroscopic data of the OEC and of improving the view of the cluster at the atomic level, a complete set of protonation configurations for the S(2) state of the OEC were investigated, and their distinctive effects on magnetic properties of the cluster were evaluated. The most recent X-ray structure of Photosystem II at 1.9 Å resolution was used and refined to obtain the optimum structure for the Mn(4)O(5)Ca core within the protein pocket. Employing this model, a set of 26 structures was constructed that tested various protonation scenarios of the water ligands and oxo bridges. Our results suggest that one of the two water molecules that are proposed to coordinate the outer Mn ion (Mn(A)) of the cluster is deprotonated in the S(2) state, as this leads to optimal experimental agreement, reproducing the correct ground state spin multiplicity (S = 1/2), spin expectation values, and EXAFS-derived metal-metal distances. Deprotonation of Ca(2+)-bound water molecules is strongly disfavored in the S(2) state, but dissociation of one of the two water ligands appears to be facile. The computed isotropic hyperfine couplings presented here allow distinctions between models to be made and call into question the assumption that the largest coupling is always attributable to Mn(III). The present results impose limits for the total charge and the proton configuration of the OEC in the S(2) state, with implications for the cascade of events in the Kok cycle and for the water splitting mechanism.

Structural-functional role of chloride in photosystem II

Chloride binding in photosystem II (PSII) is essential for photosynthetic water oxidation. However, the functional roles of chloride and possible binding sites, during oxygen evolution, remain controversial. This paper examines the functions of chloride based on its binding site revealed in the X-ray crystal structure of PSII at 1.9 Å resolution. We find that chloride depletion induces formation of a salt bridge between D2-K317 and D1-D61 that could suppress the transfer of protons to the lumen.

Structural and mechanistic investigations of photosystem II through computational methods

The advent of oxygenic photosynthesis through water oxidation by photosystem II (PSII) transformed the planet, ultimately allowing the evolution of aerobic respiration and an explosion of ecological diversity. The importance of this enzyme to life on Earth has ironically been paralleled by the elusiveness of a detailed understanding of its precise catalytic mechanism. Computational investigations have in recent years provided more and more insights into the structural and mechanistic details that underlie the workings of PSII. This review will present an overview of some of these studies, focusing on those that have aimed at elucidating the mechanism of water oxidation at the CaMn₄ cluster in PSII, and those exploring the features of the structure and dynamics of this enzyme that enable it to catalyse this energetically demanding reaction. This article is part of a Special Issue entitled: Photosystem II.

Water exchange on manganese(III) porphyrins. Mechanistic insights relevant for oxygen evolving complex and superoxide dismutation catalysis

In this work the rate constants (k(ex)) and the activation parameters (DeltaH(double dagger), DeltaS(double dagger), and DeltaV(double dagger)) for the water exchange process on Mn(III) centers have experimentally been determined using temperature and pressure dependent (17)O NMR techniques. For the investigations the Mn(III) porphyrin complexes [Mn(III)(TPPS)S(2)](n-) and [Mn(III)(TMpyP)S(2)](n+) (S = H(2)O and/or OH(-)) have been selected due to their high solution stability in a wide pH range, enabling the measurements of water exchange in the case of both diaqua and aqua-hydroxo complexes. We have experimentally demonstrated that the water exchange on Mn(III) porphyrins is a fast process (k(ex) approximately = 10(7) s(-1)) of an I(d) to I mechanism, strongly influenced by a Jahn-Teller effect and as such almost independent of a porphyrin charge and a trans ligand. This is also supported by our DFT calculations which show only a slight difference in an average Mn(III)-OH(2) bond found for a positively charged model porphyrin with protonated pyridine groups (2.446 A) and for a simple model without any substituents on the porphyrin ring (2.437 A). The calculated effective charge on the Mn center, which is significantly lower than its formal +3 charge (ca. +1.5 for diaqua; +1.4 for aqua-hydroxo), also contributes to its substitution lability. The herein presented results are discussed in connection to a possible fast exchanging substrate binding site in photosystem II and corresponding inorganic model complexes, as well as in the context of a possible inner-sphere catalytic pathway for superoxide dismutation on Mn centers.

Is Mn-Bound Substrate Water Protonated in the S(2) State of Photosystem II?

In spite of great progress in resolving the geometric structure of the water-splitting Mn(4)O(x)Ca cluster in photosystem II, the binding sites and modes of the two substrate water molecules are still insufficiently characterized. While time-resolved membrane-inlet mass spectrometry measurements indicate that both substrate water molecules are bound to the oxygen-evolving complex (OEC) in the S(2) and S(3) states (Hendry and Wydrzynski in Biochemistry 41:13328-13334, 2002), it is not known (1) if they are both Mn-bound, (2) if they are terminal or bridging ligands, and (3) in what protonation state they are bound in the different oxidation states S(i) (i = 0, 1, 2, 3, 4) of the OEC. By employing (17)O hyperfine sublevel correlation (HYSCORE) spectroscopy we recently demonstrated that in the S(2) state there is only one (type of) Mn-bound oxygen that is water exchangeable. We therefore tentatively identified this oxygen as one substrate 'water' molecule, and on the basis of the finding that it has a hyperfine interaction of about 10 MHz with the electron spin of the Mn(4)O(x)Ca cluster, we suggest that it is bound as a Mn-O-Mn bridge within a bis-mu(2) oxo-bridged unit (Su et al. in J Am Chem Soc 130:786-787, 2008). Employing pulse electron paramagnetic resonance, (1)H/(2)H Mims electron-nuclear double resonance and (2)H-HYSCORE spectroscopies together with (1)H/(2)H-exchange here, we test this hypothesis by probing the protonation state of this exchangeable oxygen. We conclude that this oxygen is fully deprotonated. This result is discussed in the light of earlier reports in the literature.

Structures and energetics for O2 formation in photosystem II

Water oxidation, forming O(2) from water and sunlight, is a fundamental process for life on earth. In nature, the enzyme photosystem II (PSII) catalyzes this reaction. The oxygen evolving complex (OEC), the complex within PSII that catalyzes the actual formation of the O-O bond, contains four manganese atoms and one calcium atom connected by oxo bonds. Seven amino acid side chains in the structure, mostly carboxylates, are ligated to the metal atoms. In the study of many enzyme mechanisms, theoretical modeling using density functional theory has served as an indispensable tool. This Account summarizes theoretical research to elucidate the mechanism for water oxidation in photosynthesis, including the most recent findings. The development of successively larger models, ranging from 50 atoms in the active site up to the present model size of 170 atoms, has revealed the mechanism of O(2) formation with increasing detail. The X-ray crystal structures of PSII have provided a framework for optimizing the theoretical models. By constraint of the backbone atoms to be at the same positions as those in the X-ray structures, the theoretical structures are in good agreement with both the measured electron density and extended X-ray absorption fine structure (EXAFS) interpretations. By following the structural and energetic changes in those structures through the different steps in the catalytic process, we have modeled the oxidation of the catalytic complex, the binding of the two substrate water molecules, and the subsequent deprotonations of those substrate molecules. In these models, the OEC forms a basin into which the water molecules naturally fit. These findings demonstrate that the binding of the second water molecule causes a reconstruction, results that are consistent with earlier EXAFS measurements. Most importantly, this Account describes a low-barrier mechanism for formation of the O-O bond, involving an oxygen radical that reacts with a mu-oxo ligand of the OEC. Further research revealed that the oxygen radical is bound in the Mn(3)Ca cube rather than to the outside manganese. This Account provides detailed diagrams of the energetics of the different S-transitions both without and with a membrane gradient. An interesting detail of these reactions concerns the role of the tyrosine (Tyr(Z)), which appears as an intermediate radical in the oxidation of the OEC. By simple electrostatic arguments, these results show that the initial oxidation of Tyr(Z) is downhill for the first two transitions but uphill for the final ones. In these later transitions, the oxidation of the OEC is coupled to deprotonations of water.

Density functional theory

Density functional theory (DFT) finds increasing use in applications related to biological systems. Advancements in methodology and implementations have reached a point where predicted properties of reasonable to high quality can be obtained. Thus, DFT studies can complement experimental investigations, or even venture with some confidence into experimentally unexplored territory. In the present contribution, we provide an overview of the properties that can be calculated with DFT, such as geometries, energies, reaction mechanisms, and spectroscopic properties. A wide range of spectroscopic parameters is nowadays accessible with DFT, including quantities related to infrared and optical spectra, X-ray absorption and Mössbauer, as well as all of the magnetic properties connected with electron paramagnetic resonance spectroscopy except relaxation times. We highlight each of these fields of application with selected examples from the recent literature and comment on the capabilities and limitations of current methods.

The MoD-QM/MM methodology for structural refinement of photosystem II and other biological macromolecules

Quantum mechanics/molecular mechanics (QM/MM) hybrid methods are currently the most powerful computational tools for studies of structure/function relations and structural refinement of macrobiomolecules (e.g., proteins and nucleic acids). These methods are highly efficient, since they implement quantum chemistry techniques for modeling only the small part of the system (QM layer) that undergoes chemical modifications, charge transfer, etc., under the influence of the surrounding environment. The rest of the system (MM layer) is described in terms of molecular mechanics force fields, assuming that its influence on the QM layer can be roughly decomposed in terms of electrostatic interactions and steric hindrance. Common limitations of QM/MM methods include inaccuracies in the MM force fields, when polarization effects are not explicitly considered, and the approximate treatment of electrostatic interactions at the boundaries between QM and MM layers. This article reviews recent advances in the development of computational protocols that allow for rigorous modeling of electrostatic interactions in extended systems beyond the common limitations of QM/MM hybrid methods. We focus on the moving-domain QM/MM (MoD-QM/MM) methodology that partitions the system into many molecular domains and obtains the electrostatic and structural properties of the whole system from an iterative self-consistent treatment of the constituent molecular fragments. We illustrate the MoD-QM/MM method as applied to the description of photosystem II as well as in conjunction with the application of spectroscopically constrained QM/MM optimization methods, based on high-resolution spectroscopic data (extended X-ray absorption fine structure spectra, and exchange coupling constants).

Evidence that D1-His332 in photosystem II from Thermosynechococcus elongatus interacts with the S3-state and not with the S2-state

Oxygen evolution by Photosystem II (PSII) is catalyzed by a Mn(4)Ca cluster. Thus far, from the crystallographic three-dimensional (3D) structures, seven amino acid residues have been identified as possible ligands of the Mn(4)Ca cluster. Among them, there is only one histidine, His332, which belongs to the D1 polypeptide. The relationships of the D1-His332 amino acid with kinetics and thermodynamic properties of the Mn(4)Ca cluster in the S(2)- and S(3)-states of the catalytic cycle were investigated in purified PSII from Thermosynechococcus elongatus. This was done by examining site-directed D1-His332Gln and D1-His332Ser mutants by a variety of spectroscopic techniques such as time-resolved UV-visible absorption change spectroscopy, cw- and pulse-EPR, thermoluminescence, and measurement of substrate water exchange. Both mutants grew photo-autotrophically and active PSII could be purified. On the basis of the parameters assessed in this work, the D1-His332(Gln, Ser) mutations had no effect in the S(2)-state. Electron spin-echo envelope modulation (ESEEM) spectroscopy also showed that possible interactions between the nuclear spin of the nitrogen(s) of D1-His332 with the electronic spin S = 1/2 of the Mn(4)Ca cluster in the S(2)-state were not detectable and that the D1-His332Ser mutation did not affect the detected hyperfine couplings. In contrast, the following changes were observed in the S(3)-state of the D1-His332 mutants: (1) The redox potential of the S(3)/S(2) couple was slightly increased by < or = 20 meV, (2) The S(3)-EPR spectrum was slightly modified, (3) The D1-His332Gln mutation resulted in a approximately 3 fold decrease of the slow (tightly bound) exchange rate and a approximately 2 fold increase of the fast exchange rate of the water substrate molecules. All these results suggest that the D1-His332 would be more involved in S(3) than in S(2). This could be one element of the conformational changes put forward in the S(2) to S(3) transition.

A new quantum chemical approach to the magnetic properties of oligonuclear transition-metal complexes: application to a model for the tetranuclear manganese cluster of photosystem II

The reliable correlation of structural features and magnetic or spectroscopic properties of oligonuclear transition-metal complexes is a critical requirement both for research into innovative magnetic materials and for elucidating the structure and function of many metalloenzymes. We have developed a novel method that for the first time enables the extraction of hyperfine coupling constants (HFCs) from broken-symmetry density functional theory (BS-DFT) calculations on clusters. Using the geometry-optimized tetranuclear manganese complex [Mn(4)O(6)(bpy)(6)](4+/3+) as a model, we first examine in detail the calculation of exchange coupling constants J through the BS-DFT approach. Complications arising from the indeterminacy of experimentally fitted J constants are identified and analyzed. It is found that only the energy levels derived from Hamiltonian diagonalization are a physically meaningful basis for comparing theory and experiment. Subsequently, the proposed theoretical scheme is applied to the calculation of (55)Mn HFCs of the Mn(III,IV,IV,IV) state of the complex, which is similar to the S(2) state of the oxygen-evolving complex (OEC) in photosystem II of oxygenic photosynthesis. The new approach performs reliably and accurately, and yields calculated HFCs that can be directly compared with experimental data. Finally, we carefully examine the dependence of HFC on the J value and draw attention to the sensitivity of the calculated values to the exchange coupling parameters. The proposed strategy extends naturally to hetero-oligonuclear clusters of arbitrary shape and nuclearity, and hence is of general validity and usefulness in the study of magnetic metal clusters. The successful application of the new approach presented here is a first step in the effort to establish correlations between the available spectroscopic information and the structural features of complex metalloenzymes like OEC.

A structure-consistent mechanism for dioxygen formation in photosystem II

In recent DFT studies a new mechanism for O-O bond formation at the oxygen evolving center (OEC) in photosystem II has been suggested. With the structure of the S(4) state required for that mechanism, the structures of the lower S states are investigated herein by adding protons and electrons. A model was used including the full amino acids for the ones ligating the OEC, and in which the backbone positions were held fixed from the X-ray structure. The only charged second-shell ligand Arg357 was also included. An optimized structure for the S(1) state was reached with a large similarity to one of those suggested by EXAFS. A full catalytic cycle was derived which can rationalize the structural relaxation in the S(2) to S(3) transition, and the fact that only an electron leaves in the transition before. Water is suggested to bind to the OEC in the S(2) to S(3), and S(4) to S(0) transitions. A new possibility for water exchange is suggested from the final energy diagram. The optimal O-O bond formation occurs between an oxygen radical and an oxo ligand. The alternative mechanism, where the oxygen radical reacts with an external water, has a barrier about 20 kcal mol(-1) higher.

Substrate water binding and oxidation in photosystem II

This mini review presents a general introduction to photosystem II with an emphasis on the oxygen evolving complex. An attempt is made to summarise what is currently known about substrate interaction in the oxygen evolving complex of photosystem II in terms of the nature of the substrate, the timing and the location of its binding. As the nature of substrate water binding has a direct bearing on the mechanism of O-O bond formation in PSII, a discussion of O-O bond formation follows the summary of current opinion in substrate interaction.

Uncovering channels in photosystem II by computer modelling: current progress, future prospects, and lessons from analogous systems

Even prior to the publication of the crystal structures for photosystem II (PSII), it had already been suggested that water, O(2) and H(+) channels exist in PSII to achieve directed transport of these molecules, and to avoid undesirable side reactions. Computational efforts to uncover these channels and investigate their properties are still at early stages, and have so far only been based on the static PSII structure. The rationale behind the proposals for such channels and the computer modelling studies thus far are reviewed here. The need to take the dynamic protein into account is then highlighted with reference to the specific issues and techniques applicable to the simulation of each of the three channels. In particular, lessons are drawn from simulation studies on other protein systems containing similar channels.

Computational insights into the O2-evolving complex of photosystem II

Mechanistic investigations of the water-splitting reaction of the oxygen-evolving complex (OEC) of photosystem II (PSII) are fundamentally informed by structural studies. Many physical techniques have provided important insights into the OEC structure and function, including X-ray diffraction (XRD) and extended X-ray absorption fine structure (EXAFS) spectroscopy as well as mass spectrometry (MS), electron paramagnetic resonance (EPR) spectroscopy, and Fourier transform infrared spectroscopy applied in conjunction with mutagenesis studies. However, experimental studies have yet to yield consensus as to the exact configuration of the catalytic metal cluster and its ligation scheme. Computational modeling studies, including density functional (DFT) theory combined with quantum mechanics/molecular mechanics (QM/MM) hybrid methods for explicitly including the influence of the surrounding protein, have proposed chemically satisfactory models of the fully ligated OEC within PSII that are maximally consistent with experimental results. The inorganic core of these models is similar to the crystallographic model upon which they were based, but comprises important modifications due to structural refinement, hydration, and proteinaceous ligation which improve agreement with a wide range of experimental data. The computational models are useful for rationalizing spectroscopic and crystallographic results and for building a complete structure-based mechanism of water-splitting in PSII as described by the intermediate oxidation states of the OEC. This review summarizes these recent advances in QM/MM modeling of PSII within the context of recent experimental studies.

A model of the oxygen-evolving center of photosystem II predicted by structural refinement based on EXAFS simulations

A refined computational structural model of the oxygen-evolving complex (OEC) of photosystem II (PSII) is introduced. The model shows that the cuboidal core Mn3CaO4 with a "dangler" Mn ligated to a corner mu4-oxide ion is maximally consistent with the positioning of the amino acids around the metal cluster as characterized by XRD models and high-resolution spectroscopic data, including polarized EXAFS of oriented single crystals and isotropic EXAFS. It is, therefore, natural to expect that the proposed structural model should be particularly useful to establish the structure of the OEC, consistently with high-resolution spectroscopic data, and for elucidating the mechanism of water-splitting in PSII as described by the intermediate oxidation states of the EC along the catalytic cycle.