QM/MM Models of the O2-Evolving Complex of Photosystem II

An analysis of the structural changes of the oxygen evolving complex of Photosystem II in the S1 and S3 states revealed by serial femtosecond crystallography

Photosystem II (PSII) is a unique natural catalyst that converts solar energy into chemical energy using earth abundant elements in water at physiological pH. Understanding the reaction mechanism will aid the design of biomimetic artificial catalysts for efficient solar energy conversion. The Mn4O5Ca cluster cycles through five increasingly oxidized intermediates before oxidizing two water molecules into O2 and releasing protons to the lumen and electrons to drive PSII reactions. The Mn coordination and OEC electronic structure changes through these intermediates. Thus, obtaining a high-resolution structure of each catalytic intermediate would help reveal the reaction mechanism. While valuable structural information was obtained from conventional X-ray crystallography, time-resolution of conventional X-ray crystallography limits the analysis of shorted-lived reaction intermediates. Serial Femtosecond X-ray crystallography (SFX), which overcomes the radiation damage by using ultra short laser pulse for imaging, has been used extensively to study the water splitting intermediates in PSII. Here, we review the state of the art and our understanding of the water splitting reaction before and after the advent of SFX. Furthermore, we analyze the likely Mn coordination in multiple XFEL structures prepared in the dark-adapted S1 state and those following two-flashes which are poised in the penultimate S3 oxidation state based on Mn coordination chemistry. Finally, we summarize the major contributions of the SFX to our understanding of the structures of the S1 and S3 states.

Structural and quantum chemical basis for OCP-mediated quenching of phycobilisomes

Cyanobacteria use large antenna complexes called phycobilisomes (PBSs) for light harvesting. However, intense light triggers non-photochemical quenching, where the orange carotenoid protein (OCP) binds to PBS, dissipating excess energy as heat. The mechanism of efficiently transferring energy from phycocyanobilins in PBS to canthaxanthin in OCP remains insufficiently understood. Using cryo-electron microscopy, we unveiled the OCP-PBS complex structure at 1.6- to 2.1-angstrom resolution, showcasing its inherent flexibility. Using multiscale quantum chemistry, we disclosed the quenching mechanism. Identifying key protein residues, we clarified how canthaxanthin's transition dipole moment in its lowest-energy dark state becomes large enough for efficient energy transfer from phycocyanobilins. Our energy transfer model offers a detailed understanding of the atomic determinants of light harvesting regulation and antenna architecture in cyanobacteria.

Mapping the Oxygens in the Oxygen-Evolving Complex of Photosystem II by Their Nucleophilicity Using Quantum Descriptors

The oxygen-evolving complex (OEC) of Photosystem II catalyzes the water-splitting reaction using solar energy. Thus, understanding the reaction mechanism will inspire the design of biomimetic artificial catalysts that convert solar energy to chemical energy. Conceptual Density Functional Theory (CDFT) focuses on understanding the reactivity of molecules and the atomic contribution to the overall nucleophilicity and electrophilicity of the molecule using quantum descriptors. However, this method has not been applied to the OEC before. Here, we use Fukui functions and the dual descriptor to provide quantitative measures of the nucleophilicity and electrophilicity of oxygens in the OEC for different models in different S states. Our results show that the μ-oxo bridges connected to terminal Mn4 are nucleophilic, and those in the cube formed by Mn1, Mn2, and Mn3 are mostly electrophilic. The dual descriptors of the bridging oxygens in the OEC showed a similar reactivity to that of bridging oxygens in Mn model compounds. However, the terminal water W1, which is bound to Mn4, showed very strong reactivity in some of the S3 models. Thus, our calculations support the model that proposes the formation of the O2 molecule through nucleophilic attack by a terminal water.

Insights into PSII's S3YZ• State: An Electronic and Magnetic Analysis

Using BS-DFT (broken-symmetry density functional theory), the electronic and magnetic properties of the S3YZ• state of photosystem II were investigated and compared to those of the S3 state. While the O5 oxo-O6 hydroxo species presents little difference between the two states, a previously identified [O5O6]3- exhibits reduced stabilization of the O5-O6 shared spin. This species is shown to have some coupling with the YZ• center through Mn1 and O6. Similarly, a peroxo species is found to exhibit significant exchange couplings between the YZ• center and the Mn cluster through Mn1. Mechanistic changes in O-O bond formation in S3YZ• are highlighted by analysis of IBOs (intrinsic bonding orbitals) showing deviation for Mn1 and O6 centered IBOs. This change in coupling interactions throughout the complex as a result of S3YZ• formation presents implications for the determination of the mechanism spanning the end of the S3 and the start of the S4 states, affecting both electron movement and oxygen bond formation.

A dimeric chlorophyll electron acceptor differentiates type I from type II photosynthetic reaction centers

This research addresses one of the most compelling issues in the field of photosynthesis, namely, the role of the accessory chlorophyll molecules in primary charge separation. Using a combination of empirical and computational methods, we demonstrate that the primary acceptor of photosystem (PS) I is a dimer of accessory and secondary chlorophyll molecules, Chl_2A and Chl_3A, with an asymmetric electron charge density distribution. The incorporation of highly coupled donors and acceptors in PS I allows for extensive delocalization that prolongs the lifetime of the charge-separated state, providing for high quantum efficiency. The discovery of this motif has widespread implications ranging from the evolution of naturally occurring reaction centers to the development of a new generation of highly efficient artificial photosynthetic systems.

Video abstract

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

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

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.

The S3 State of the Oxygen-Evolving Complex: Overview of Spectroscopy and XFEL Crystallography with a Critical Evaluation of Early-Onset Models for O–O Bond Formation

The catalytic cycle of the oxygen-evolving complex (OEC) of photosystem II (PSII) comprises five intermediate states Si (i = 0–4), from the most reduced S0 state to the most oxidized S4, which spontaneously evolves dioxygen. The precise geometric and electronic structure of the Si states, and hence the mechanism of O–O bond formation in the OEC, remain under investigation, particularly for the final steps of the catalytic cycle. Recent advances in protein crystallography based on X-ray free-electron lasers (XFELs) have produced new structural models for the S3 state, which indicate that two of the oxygen atoms of the inorganic Mn4CaO6 core of the OEC are in very close proximity. This has been interpreted as possible evidence for “early-onset” O–O bond formation in the S3 state, as opposed to the more widely accepted view that the O–O bond is formed in the final state of the cycle, S4. Peroxo or superoxo formation in S3 has received partial support from computational studies. Here, a brief overview is provided of spectroscopic information, recent crystallographic results, and computational models for the S3 state. Emphasis is placed on computational S3 models that involve O–O formation, which are discussed with respect to their agreement with structural information, experimental evidence from various spectroscopic studies, and substrate exchange kinetics. Despite seemingly better agreement with some of the available crystallographic interpretations for the S3 state, models that implicate early-onset O–O bond formation are hard to reconcile with the complete line of experimental evidence, especially with X-ray absorption, X-ray emission, and magnetic resonance spectroscopic observations. Specifically with respect to quantum chemical studies, the inconclusive energetics for the possible isoforms of S3 is an acute problem that is probably beyond the capabilities of standard density functional theory.

Oxygenic photosynthesis: history, status and perspective

Cyanobacteria and plants carry out oxygenic photosynthesis. They use water to generate the atmospheric oxygen we breathe and carbon dioxide to produce the biomass serving as food, feed, fibre and fuel. This paper scans the emergence of structural and mechanistic understanding of oxygen evolution over the past 50 years. It reviews speculative concepts and the stepped insight provided by novel experimental and theoretical techniques. Driven by sunlight photosystem II oxidizes the catalyst of water oxidation, a hetero-metallic Mn4CaO5(H2O)4 cluster. Mn3Ca are arranged in cubanoid and one Mn dangles out. By accumulation of four oxidizing equivalents before initiating dioxygen formation it matches the four-electron chemistry from water to dioxygen to the one-electron chemistry of the photo-sensitizer. Potentially harmful intermediates are thereby occluded in space and time. Kinetic signatures of the catalytic cluster and its partners in the photo-reaction centre have been resolved, in the frequency domain ranging from acoustic waves via infra-red to X-ray radiation, and in the time domain from nano- to milli-seconds. X-ray structures to a resolution of 1.9 Å are available. Even time resolved X-ray structures have been obtained by clocking the reaction cycle by flashes of light and diffraction with femtosecond X-ray pulses. The terminal reaction cascade from two molecules of water to dioxygen involves the transfer of four electrons, two protons, one dioxygen and one water. A rigorous mechanistic analysis is challenging because of the kinetic enslaving at millisecond duration of six partial reactions (4e−, 1H+, 1O2). For the time being a peroxide-intermediate in the reaction cascade to dioxygen has been in focus, both experimentally and by quantum chemistry. Homo sapiens has relied on burning the products of oxygenic photosynthesis, recent and fossil. Mankind's total energy consumption amounts to almost one-fourth of the global photosynthetic productivity. If the average power consumption equalled one of those nations with the highest consumption per capita it was four times greater and matched the total productivity. It is obvious that biomass should be harvested for food, feed, fibre and platform chemicals rather than for fuel.

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.

Geometrical properties of the manganese(iv)/iron(iii) cofactor of Chlamydia trachomatis ribonucleotide reductase unveiled by simulations of XAS spectra

A combination of EXAFS simulations and DFT calculations, including a novel protocol to evaluate Debye–Waller factors, provide insights into the structure of the Mn(iv)/Fe(iii) cofactor ofCtR2.

Radiation Damage in XFEL: Case study from the oxygen-evolving complex of Photosystem II

Structural changes induced by radiation damage in X-ray crystallography hinder the ability to understand the structure/function relationship in chemical reactions. Serial femtosecond crystallography overcomes this problem by exposing the sample to very short and intense laser pulse leading to measurement before destruction. Here we use molecular modeling to map the radiation damage during the 10-50 fs to the intensity, the energy and the time duration of the laser pulse on the oxygen-evolving complex (OEC) of photosystem II. In the model, the nuclei move classically in a fully quantum potential created by electron density under the effect of strong laser pulse in the Ehrenfest dynamics regime. The results show that the Mn-Mn and Mn-Ca distances are less affected by radiation damage due to the their heavy masses, while one μ-oxo bridge (O5) moves significantly. The radiation damage may induce conformational changes of the water ligands but only bond elongation for the amino acids ligands. These effects are relatively intensity independent from 1016 to 1017 W/cm2, but changes increase dramatically if the beam intensity is increased to 1018 W/cm2. In addition, the self amplified spontaneous emission (SASE) nature of the laser beam does not affect the dynamics of the ions.

'Photosystem II: the water splitting enzyme of photosynthesis and the origin of oxygen in our atmosphere'

About 3 billion years ago an enzyme emerged which would dramatically change the chemical composition of our planet and set in motion an unprecedented explosion in biological activity. This enzyme used solar energy to power the thermodynamically and chemically demanding reaction of water splitting. In so doing it provided biology with an unlimited supply of reducing equivalents needed to convert carbon dioxide into the organic molecules of life while at the same time produced oxygen to transform our planetary atmosphere from an anaerobic to an aerobic state. The enzyme which facilitates this reaction and therefore underpins virtually all life on our planet is known as Photosystem II (PSII). It is a pigment-binding, multisubunit protein complex embedded in the lipid environment of the thylakoid membranes of plants, algae and cyanobacteria. Today we have detailed understanding of the structure and functioning of this key and unique enzyme. The journey to this level of knowledge can be traced back to the discovery of oxygen itself in the 18th-century. Since then there has been a sequence of mile stone discoveries which makes a fascinating story, stretching over 200 years. But it is the last few years that have provided the level of detail necessary to reveal the chemistry of water oxidation and O-O bond formation. In particular, the crystal structure of the isolated PSII enzyme has been reported with ever increasing improvement in resolution. Thus the organisational and structural details of its many subunits and cofactors are now well understood. The water splitting site was revealed as a cluster of four Mn ions and a Ca ion surrounded by amino-acid side chains, of which seven provide direct ligands to the metals. The metal cluster is organised as a cubane structure composed of three Mn ions and a Ca2+ linked by oxo-bonds with the fourth Mn ion attached to the cubane. This structure has now been synthesised in a non-protein environment suggesting that it is a totally inorganic precursor for the evolution of the photosynthetic oxygen-evolving complex. In summary, the overall structure of the catalytic site has given a framework on which to build a mechanistic scheme for photosynthetic dioxygen generation and at the same time provide a blue-print and incentive to develop catalysts for artificial photo-electrochemical systems to split water and generate renewable solar fuels.

The MOD-QM/MM Method: Applications to Studies of Photosystem II and DNA G-Quadruplexes

Quantum mechanics/molecular mechanics (QM/MM) hybrid methods are currently the most powerful computational tools for studies of structure/function relations and catalytic sites embedded in macrobiomolecules (eg, proteins and nucleic acids). QM/MM methodologies are highly efficient since they implement quantum chemistry methods for modeling only the portion of the system involving bond-breaking/forming processes (QM layer), as influenced by the surrounding molecular environment described in terms of molecular mechanics force fields (MM layer). Some of the limitations of QM/MM methods when polarization effects are not explicitly considered include the approximate treatment of electrostatic interactions between QM and MM layers. Here, we review recent advances in the development of computational protocols that allow for rigorous modeling of electrostatic interactions in biomacromolecules and structural refinement, beyond the common limitations of QM/MM hybrid methods. We focus on photosystem II (PSII) with emphasis on the description of the oxygen-evolving complex (OEC) and its high-resolution extended X-ray absorption fine structure spectra (EXAFS) in conjunction with Monte Carlo structural refinement. Furthermore, we review QM/MM structural refinement studies of DNA G4 quadruplexes with embedded monovalent cations and direct comparisons to NMR data.

What computational chemistry and magnetic resonance reveal concerning the oxygen evolving centre in Photosystem II

Density Functional Theory (DFT) computational studies of the Mn₄/Ca Oxygen Evolving Complex (OEC) region of Photosystem II in the paramagnetic S₂ and S₃ states of the water oxdizing catalytic cycle are described. These build upon recent advances in computationally understanding the detailed S₁ state OEC geometries, revealed by the recent high resolution Photosystem II crystal structures of Shen et al., at 1.90Å and 1.95Å (Petrie et al., 2015, Angew. Chem. Int. Ed., 54, 7120). The models feature a 'Low Oxidation Paradigm' assumption for the mean Mn oxidation states in the functional enzyme, with the mean oxidation levels being 3.0, 3.25 and 3.5 in S₁, S₂ and S₃, respectively. These calculations are used to infer magnetic exchange interactions within the coupled OEC cluster, particularly in the Electron Paramagnetic Resonance (EPR)-visible S₂ and S₃ states. Detailed computational estimates of the intrinsic magnitudes and molecular orientations of the ⁵⁵Mn hyperfine tensors in the S₂ state are presented. These parameters, together with the resultant spin projected hyperfine values are compared with recent appropriate experimental EPR data (Continuous Wave (CW), Electron-Nuclear Double Resonance (ENDOR) and ELDOR (Electron-Electron Double Resonance)-Detected Nuclear Magnetic Resonance (EDNMR)) from the OEC. It is found that an effective Coupled Dimer magnetic organization of the four Mn in the OEC cluster in the S₂ and S₃ states is able to quantitatively rationalize the observed ⁵⁵Mn hyperfine data. This is consistent with structures we propose to represent the likely state of the OEC in the catalytically active form of the enzyme.

Redox-coupled substrate water reorganization in the active site of Photosystem II-The role of calcium in substrate water delivery

Photosystem II (PSII) catalyzes light-driven water splitting in nature and is the key enzyme for energy input into the biosphere. Important details of its mechanism are not well understood. In order to understand the mechanism of water splitting, we perform here large-scale density functional theory (DFT) calculations on the active site of PSII in different oxidation, spin and ligand states. Prior to formation of the O-O bond, we find that all manganese atoms are oxidized to Mn(IV) in the S3 state, consistent with earlier studies. We find here, however, that the formation of the S3 state is coupled to the movement of a calcium-bound hydroxide (W3) from the Ca to a Mn (Mn1 or Mn4) in a process that is triggered by the formation of a tyrosyl radical (Tyr-161) and its protonated base, His-190. We find that subsequent oxidation and deprotonation of this hydroxide on Mn1 result in formation of an oxyl-radical that can exergonically couple with one of the oxo-bridges (O5), forming an O-O bond. When O(2) leaves the active site, a second Ca-bound water molecule reorients to bridge the gap between the manganese ions Mn1 and Mn4, forming a new oxo-bridge for the next reaction cycle. Our findings are consistent with experimental data, and suggest that the calcium ion may control substrate water access to the water oxidation sites.

Manganese Compounds as Water-Oxidizing Catalysts: From the Natural Water-Oxidizing Complex to Nanosized Manganese Oxide Structures

All cyanobacteria, algae, and plants use a similar water-oxidizing catalyst for water oxidation. This catalyst is housed in Photosystem II, a membrane-protein complex that functions as a light-driven water oxidase in oxygenic photosynthesis. Water oxidation is also an important reaction in artificial photosynthesis because it has the potential to provide cheap electrons from water for hydrogen production or for the reduction of carbon dioxide on an industrial scale. The water-oxidizing complex of Photosystem II is a Mn-Ca cluster that oxidizes water with a low overpotential and high turnover frequency number of up to 25-90 molecules of O2 released per second. In this Review, we discuss the atomic structure of the Mn-Ca cluster of the Photosystem II water-oxidizing complex from the viewpoint that the underlying mechanism can be informative when designing artificial water-oxidizing catalysts. This is followed by consideration of functional Mn-based model complexes for water oxidation and the issue of Mn complexes decomposing to Mn oxide. We then provide a detailed assessment of the chemistry of Mn oxides by considering how their bulk and nanoscale properties contribute to their effectiveness as water-oxidizing catalysts.

How does ammonia bind to the oxygen-evolving complex in the S2state of photosynthetic water oxidation? Theoretical support and implications for the W1 substitution mechanism

The mechanistic study shows that NH₃substitutes W1 rather than O5 of the OEC in the S₂state and leaves in the S₄′ state.

Metal oxidation states in biological water splitting

A central question in biological water splitting concerns the oxidation states of the manganese ions that comprise the oxygen-evolving complex of photosystem II.

A central question in biological water splitting concerns the oxidation states of the manganese ions that comprise the oxygen-evolving complex of photosystem II. Understanding the nature and order of oxidation events that occur during the catalytic cycle of five S_i states (i = 0–4) is of fundamental importance both for the natural system and for artificial water oxidation catalysts. Despite the widespread adoption of the so-called “high-valent scheme”—where, for example, the Mn oxidation states in the S₂ state are assigned as III, IV, IV, IV—the competing “low-valent scheme” that differs by a total of two metal unpaired electrons (i.e. III, III, III, IV in the S₂ state) is favored by several recent studies for the biological catalyst. The question of the correct oxidation state assignment is addressed here by a detailed computational comparison of the two schemes using a common structural platform and theoretical approach. Models based on crystallographic constraints were constructed for all conceivable oxidation state assignments in the four (semi)stable S states of the oxygen evolving complex, sampling various protonation levels and patterns to ensure comprehensive coverage. The models are evaluated with respect to their geometric, energetic, electronic, and spectroscopic properties against available experimental EXAFS, XFEL-XRD, EPR, ENDOR and Mn K pre-edge XANES data. New 2.5 K ⁵⁵Mn ENDOR data of the S₂ state are also reported. Our results conclusively show that the entire S state phenomenology can only be accommodated within the high-valent scheme by adopting a single motif and protonation pattern that progresses smoothly from S₀ (III, III, III, IV) to S₃ (IV, IV, IV, IV), satisfying all experimental constraints and reproducing all observables. By contrast, it was impossible to construct a consistent cycle based on the low-valent scheme for all S states. Instead, the low-valent models developed here may provide new insight into the over-reduced S states and the states involved in the assembly of the catalytically active water oxidizing cluster.

MoD-QM/MM Structural Refinement Method: Characterization of Hydrogen Bonding in the Oxytricha nova G-Quadruplex

A generalization of the Moving-Domain Quantum Mechanics/Molecular Mechanics (MoD-QM/MM) hybrid method [Gascon, J. A.; Leung, S. S. F.; Batista, E. R.; Batista, V. S. J. Chem. Theory Comput. 2006, 2, 175-186] is introduced to provide a self-consistent computational protocol for structural refinement of extended systems. The method partitions the system into molecular domains that are iteratively optimized as quantum mechanical (QM) layers embedded in their surrounding molecular environment to obtain an ab initio quality description of the geometry and the molecular electrostatic potential of the extended system composed of those constituent fragments. The resulting methodology is benchmarked as applied to model systems that allow for full QM optimization as well as through refinement of the hydrogen bonding geometry in Oxytricha nova guanine quadruplex for which several studies have been reported, including the X-ray structure and NMR data. Calculations of (1)H NMR chemical shifts based on the gauge independent atomic orbital (GIAO) method and direct comparisons with experiments show that solvated MoD-QM/MM structures, sampled from explicit solvent molecular dynamics simulations, allow for NMR simulations in much improved agreement with experimental data than models based on the X-ray structure or those optimized using classical molecular mechanics force fields.

FTIR studies of metal ligands, networks of hydrogen bonds, and water molecules near the active site Mn₄CaO₅ cluster in Photosystem II

The photosynthetic conversion of water to molecular oxygen is catalyzed by the Mn₄CaO₅ cluster in Photosystem II and provides nearly our entire supply of atmospheric oxygen. The Mn₄CaO₅ cluster accumulates oxidizing equivalents in response to light-driven photochemical events within Photosystem II and then oxidizes two molecules of water to oxygen. The Mn₄CaO₅ cluster converts water to oxygen much more efficiently than any synthetic catalyst because its protein environment carefully controls the cluster's reactivity at each step in its catalytic cycle. This control is achieved by precise choreography of the proton and electron transfer reactions associated with water oxidation and by careful management of substrate (water) access and proton egress. This review describes the FTIR studies undertaken over the past two decades to identify the amino acid residues that are responsible for this control and to determine the role of each. In particular, this review describes the FTIR studies undertaken to characterize the influence of the cluster's metal ligands on its activity, to delineate the proton egress pathways that link the Mn₄CaO₅ cluster with the thylakoid lumen, and to characterize the influence of specific residues on the water molecules that serve as substrate or as participants in the networks of hydrogen bonds that make up the water access and proton egress pathways. This information will improve our understanding of water oxidation by the Mn₄CaO₅ catalyst in Photosystem II and will provide insight into the design of new generations of synthetic catalysts that convert sunlight into useful forms of storable energy. This article is part of a Special Issue entitled: Vibrational spectroscopies and bioenergetic systems.

Molecular mechanisms for generating transmembrane proton gradients

Membrane proteins use the energy of light or high energy substrates to build a transmembrane proton gradient through a series of reactions leading to proton release into the lower pH compartment (P-side) and proton uptake from the higher pH compartment (N-side). This review considers how the proton affinity of the substrates, cofactors and amino acids are modified in four proteins to drive proton transfers. Bacterial reaction centers (RCs) and photosystem II (PSII) carry out redox chemistry with the species to be oxidized on the P-side while reduction occurs on the N-side of the membrane. Terminal redox cofactors are used which have pKas that are strongly dependent on their redox state, so that protons are lost on oxidation and gained on reduction. Bacteriorhodopsin is a true proton pump. Light activation triggers trans to cis isomerization of a bound retinal. Strong electrostatic interactions within clusters of amino acids are modified by the conformational changes initiated by retinal motion leading to changes in proton affinity, driving transmembrane proton transfer. Cytochrome c oxidase (CcO) catalyzes the reduction of O2 to water. The protons needed for chemistry are bound from the N-side. The reduction chemistry also drives proton pumping from N- to P-side. Overall, in CcO the uptake of 4 electrons to reduce O2 transports 8 charges across the membrane, with each reduction fully coupled to removal of two protons from the N-side, the delivery of one for chemistry and transport of the other to the P-side.

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.

Electrostatic effects on proton coupled electron transfer in oxomanganese complexes inspired by the oxygen-evolving complex of photosystem II

The influence of electrostatic interactions on the free energy of proton coupled electron transfer in biomimetic oxomanganese complexes inspired by the oxygen-evolving complex (OEC) of photosystem II (PSII) are investigated. The reported study introduces an enhanced multiconformer continuum electrostatics (MCCE) model, parametrized at the density functional theory (DFT) level with a classical valence model for the oxomanganese core. The calculated pKa's and oxidation midpoint potentials (E(m)'s) match experimental values for eight complexes, indicating that purely electrostatic contributions account for most of the observed couplings between deprotonation and oxidation state transitions. We focus on pKa's of terminal water ligands in [Mn(II/III)(H2O)6](2+/3+) (1), [Mn(III)(P)(H2O)2](3-) (2, P = 5,10,15,20-tetrakis(2,6-dichloro-3-sulfonatophenyl)porphyrinato), [Mn2(IV,IV)(μ-O)2(terpy)2(H2O)2](4+) (3, terpy = 2,2':6',2″-terpyridine), and [Mn3(IV,IV,IV)(μ-O)4(phen)4(H2O)2](4+) (4, phen = 1,10-phenanthroline) and the pKa's of μ-oxo bridges and Mn E(m)'s in [Mn2(μ-O)2(bpy)4] (5, bpy = 2,2'-bipyridyl), [Mn2(μ-O)2(salpn)2] (6, salpn = N,N'-bis(salicylidene)-1,3-propanediamine), [Mn2(μ-O)2(3,5-di(Cl)-salpn)2] (7), and [Mn2(μ-O)2(3,5-di(NO2)-salpn)2] (8). The analysis of complexes 6-8 highlights the strong coupling between electron and proton transfers, with any Mn oxidation lowering the pKa of an oxo bridge by 10.5 ± 0.9 pH units. The model also accounts for changes in the E(m)'s by ligand substituents, such as found in complexes 6-8, due to the electron withdrawing Cl (7) and NO2 (8). The reported study provides the foundation for analysis of electrostatic effects in other oxomanganese complexes and metalloenzymes, where proton coupled electron transfer plays a fundamental role in redox-leveling mechanisms.

On the magnetic and spectroscopic properties of high-valent Mn3CaO4 cubanes as structural units of natural and artificial water-oxidizing catalysts

The Mn(IV)3CaO4 cubane is a structural motif present in the oxygen-evolving complex (OEC) of photosystem II and in water-oxidizing Mn/Ca layered oxides. This work investigates the magnetic and spectroscopic properties of two recently synthesized complexes and a series of idealized models that incorporate this structural unit. Magnetic interactions, accessible spin states, and (55)Mn isotropic hyperfine couplings are computed with quantum chemical methods and form the basis for structure-property correlations. Additionally, the effects of oxo-bridge protonation and one-electron reduction are examined. The calculated properties are found to be in excellent agreement with available experimental data. It is established that all synthetic and model Mn(IV)3CaO4 cubane complexes have the same high-spin S = (9)/2 ground state. The magnetic coupling conditions under which different ground spin states can be accessed are determined. Substitution of Mn(IV) magnetic centers by diamagnetic ions [e.g., Ge(IV)] allows one to "switch off" specific spin sites in order to examine the magnetic orbitals along individual Mn-Mn exchange pathways, which confirms the predominance of ferromagnetic interactions within the cubane framework. The span of the Heisenberg spin ladder is found to correlate inversely with the number of protonated oxo bridges. Energetic comparisons for protonated models show that the tris-μ-oxo bridge connecting only Mn ions in the cubane has the lowest proton affinity and that the average relaxation energy per additional proton is on the order of 18 kcal·mol(-1), thus making access to ground states other than the high-spin S = (9)/2 state in these cubanes unlikely. The relevance of these cubanes for the OEC and synthetic oxides is discussed.

From natural to artificial photosynthesis

Demand for energy is projected to increase at least twofold by mid-century relative to the present global consumption because of predicted population and economic growth. This demand could be met, in principle, from fossil energy resources, particularly coal. However, the cumulative nature of carbon dioxide (CO(2)) emissions demands that stabilizing the atmospheric CO(2) levels to just twice their pre-anthropogenic values by mid-century will be extremely challenging, requiring invention, development and deployment of schemes for carbon-neutral energy production on a scale commensurate with, or larger than, the entire present-day energy supply from all sources combined. Among renewable and exploitable energy resources, nuclear fusion energy or solar energy are by far the largest. However, in both cases, technological breakthroughs are required with nuclear fusion being very difficult, if not impossible on the scale required. On the other hand, 1 h of sunlight falling on our planet is equivalent to all the energy consumed by humans in an entire year. If solar energy is to be a major primary energy source, then it must be stored and despatched on demand to the end user. An especially attractive approach is to store solar energy in the form of chemical bonds as occurs in natural photosynthesis. However, a technology is needed which has a year-round average conversion efficiency significantly higher than currently available by natural photosynthesis so as to reduce land-area requirements and to be independent of food production. Therefore, the scientific challenge is to construct an 'artificial leaf' able to efficiently capture and convert solar energy and then store it in the form of chemical bonds of a high-energy density fuel such as hydrogen while at the same time producing oxygen from water. Realistically, the efficiency target for such a technology must be 10 per cent or better. Here, we review the molecular details of the energy capturing reactions of natural photosynthesis, particularly the water-splitting reaction of photosystem II and the hydrogen-generating reaction of hydrogenases. We then follow on to describe how these two reactions are being mimicked in physico-chemical-based catalytic or electrocatalytic systems with the challenge of creating a large-scale robust and efficient artificial leaf technology.

Photosystem II and the unique role of bicarbonate: a historical perspective

In photosynthesis, cyanobacteria, algae and plants fix carbon dioxide (CO(2)) into carbohydrates; this is necessary to support life on Earth. Over 50 years ago, Otto Heinrich Warburg discovered a unique stimulatory role of CO(2) in the Hill reaction (i.e., O(2) evolution accompanied by reduction of an artificial electron acceptor), which, obviously, does not include any carbon fixation pathway; Warburg used this discovery to support his idea that O(2) in photosynthesis originates in CO(2). During the 1960s, a large number of researchers attempted to decipher this unique phenomenon, with limited success. In the 1970s, Alan Stemler, in Govindjee's lab, perfected methods to get highly reproducible results, and observed, among other things, that the turnover of Photosystem II (PSII) was stimulated by bicarbonate ions (hydrogen carbonate): the effect would be on the donor or the acceptor, or both sides of PSII. In 1975, Thomas Wydrzynski, also in Govindjee's lab, discovered that there was a definite bicarbonate effect on the electron acceptor (the plastoquinone) side of PSII. The most recent 1.9Å crystal structure of PSII, unequivocally shows HCO(3)(-) bound to the non-heme iron that sits in-between the bound primary quinone electron acceptor, Q(A), and the secondary quinone electron acceptor Q(B). In this review, we focus on the historical development of our understanding of this unique bicarbonate effect on the electron acceptor side of PSII, and its mechanism as obtained by biochemical, biophysical and molecular biological approaches in many laboratories around the World. We suggest an atomic level model in which HCO(3)(-)/CO(3)(2-) plays a key role in the protonation of the reduced Q(B). In addition, we make comments on the role of bicarbonate on the donor side of PSII, as has been extensively studied in the labs of Alan Stemler (USA) and Vyacheslav Klimov (Russia). We end this review by discussing the uniqueness of bicarbonate's role in oxygenic photosynthesis and its role in the evolutionary development of O(2)-evolving PSII. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.

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.