Computational insights into the O2-evolving complex of photosystem II

Covalent Metal-Metal-Bonded Mn4 Tetrahedron Inscribed within a Four-Coordinate Manganese Cubane Cluster, As Evidenced by Unexpected Temperature-Independent Diamagnetism

The electronic structures of the manganese(IV) cubane cluster Mn(μ₃-N^tBu)₄(N^tBu)₄ (1) and its one-electron-oxidized analogue, the 3:1 Mn^IV/Mn^V cluster [Mn(μ₃-N^tBu)₄(N^tBu)₄]⁺[PF₆]^- (1⁺[PF₆]), are described. The S = 0 spin quantum number of 1 is explained by a diamagnetic electronic structure where all metal-based d electrons are paired in Mn-Mn bonding orbitals. Temperature- and power-dependent studies of the S = ¹/₂ electron paramagnetic resonance signal of 1⁺ are consistent with an electronic structure described as a delocalized one-electron radical.

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.

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.

The role of nano-sized manganese oxides in the oxygen-evolution reactions by manganese complexes: towards a complete picture

Eighteen Mn complexes with N-donor and carboxylate ligands have been synthesized and characterized. Three Mn complexes among them are new and are reported for the first time. The reactions of oxygen evolution in the presence of oxone (2KHSO5·KHSO4·K2SO4) and cerium(iv) ammonium nitrate catalyzed by these complexes are studied and characterized by UV-visible spectroscopy, X-ray diffraction spectrometry, dynamic light scattering, Fourier transform infrared spectroscopy, electron paramagnetic resonance spectroscopy, transmission electron microscopy, scanning electron microscopy, membrane-inlet mass spectrometry and electrochemistry. Some of these complexes evolve oxygen in the presence of oxone as a primary oxidant. CO2 and MnO4(-) are other products of these reactions. Based on spectroscopic studies, the true catalysts for oxygen evolution in these reactions are different. We proposed that for the oxygen evolution reactions in the presence of oxone, the true catalysts are both high valent Mn complexes and Mn oxides, but for the reactions in the presence of cerium(iv) ammonium nitrate, the active catalyst is most probably a Mn oxide.

Mechanistic elucidation of the stepwise formation of a tetranuclear manganese pinned butterfly cluster via N-N bond cleavage, hydrogen atom transfer, and cluster rearrangement

A mechanistic pathway for the formation of the structurally characterized manganese-amide-hydrazide pinned butterfly complex, Mn4(μ3-PhN-NPh-κ(3)N,N')2(μ-PhN-NPh-κ(2)-N,N')(μ-NHPh)2L4 (L = THF, py), is proposed and supported by the use of labeling studies, kinetic measurements, kinetic competition experiments, kinetic isotope effects, and hydrogen atom transfer reagent substitution, and via the isolation and characterization of intermediates using X-ray diffraction and electron paramagnetic resonance spectroscopy. The data support a formation mechanism whereby bis[bis(trimethylsilyl)amido]manganese(II) (Mn(NR2)2, where R = SiMe3) reacts with N,N'-diphenylhydrazine (PhNHNHPh) via initial proton transfer, followed by reductive N-N bond cleavage to form a long-lived Mn(IV) imido multinuclear complex. Coordinating solvents activate this cluster for abstraction of hydrogen atoms from an additional equivalent of PhNHNHPh resulting in a Mn(II)phenylamido dimer, Mn2(μ-NHPh)2(NR2)2L2. This dimeric complex further assembles in fast steps with two additional equivalents of PhNHNHPh replacing the terminal silylamido ligands with η(1)-hydrazine ligands to give a dimeric Mn2(μ-NHPh)2(PhN-NHPh)2L4 intermediate, and finally, the addition of two additional equivalents of Mn(NR2)2 and PhNHNHPh gives the pinned butterfly cluster.

Preparation of a "twisted basket" Mn(4)N(8) cluster: a two-hydrogen-atom reduced analogue of the Mn(4)N(8) pinned butterfly

Mn4(μ-NHPh)4(μ-PhNNPh-κ(2)N,N')2(py)4 () is synthesized via self assembly from dimeric Mn2(μ-NHPh)2(NR2)2 and PhNHNHPh (R = SiMe3). This cluster represents the N-N cleaved version of the previously-reported Mn4(μ-NHPh)2(μ3-PhNNPh-κ(3)N,N')2(μ-PhNNPh-κ(2)N,N')(py)4 "pinned butterfly" cluster (), formally reduced by two hydrogen atoms. Cluster may be converted to by addition of N,N'-diphenylhydrazine as a two-electron reductant.

Taking snapshots of photosynthetic water oxidation using femtosecond X-ray diffraction and spectroscopy

The dioxygen we breathe is formed by light-induced oxidation of water in photosystem II. O2 formation takes place at a catalytic manganese cluster within milliseconds after the photosystem II reaction centre is excited by three single-turnover flashes. Here we present combined X-ray emission spectra and diffraction data of 2-flash (2F) and 3-flash (3F) photosystem II samples, and of a transient 3F' state (250 μs after the third flash), collected under functional conditions using an X-ray free electron laser. The spectra show that the initial O-O bond formation, coupled to Mn reduction, does not yet occur within 250 μs after the third flash. Diffraction data of all states studied exhibit an anomalous scattering signal from Mn but show no significant structural changes at the present resolution of 4.5 Å. This study represents the initial frames in a molecular movie of the structural changes during the catalytic reaction in photosystem II.

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.

Molecular dynamics simulations reveal highly permeable oxygen exit channels shared with water uptake channels in photosystem II

Photosystem II (PSII) catalyzes the oxidation of water in the conversion of light energy into chemical energy in photosynthesis. Water delivery and oxygen removal from the oxygen evolving complex (OEC), buried deep within PSII, are critical requirements to facilitate the reaction and minimize reactive oxygen damage. It has often been assumed that water and oxygen travel through separate channels within PSII, as demonstrated in cytochrome c oxidase. This study describes all-atom molecular dynamics simulations of PSII designed to investigate channels by fully characterizing the distribution and permeation of both water and oxygen. Interestingly, most channels found in PSII were permeable to both oxygen and water, however individual channels exhibited different energetic barriers for the two solutes. Several routes for oxygen diffusion within PSII with low energy permeation barriers were found, ensuring its fast removal from the OEC. In contrast, all routes for water showed significant energy barriers, corresponding to a much slower permeation rate for water through PSII. Two major factors were responsible for this selectivity: (1) hydrogen bonds between water and channel amino acids, and (2) steric restraints. Our results reveal the presence of a shared network of channels in PSII optimized to both facilitate the quick removal of oxygen and effectively restrict the water supply to the OEC to help stabilize and protect it from small water soluble inhibitors.

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.

Non-heme manganese catalase--the 'other' catalase

Non-heme manganese catalases are widely distributed over microbial life and represent an environmentally important alternative to heme-containing catalases in antioxidant defense. Manganese catalases contain a binuclear manganese complex as their catalytic active site rather than a heme, and cycle between Mn(2)(II,II) and Mn(2)(III,III) states during turnover. X-ray crystallography has revealed the key structural elements of the binuclear manganese active site complex that can serve as the starting point for computational studies on the protein. Four manganese catalase enzymes have been isolated and characterized, and the enzyme appears to have a broad phylogenetic distribution including both bacteria and archae. More than 100 manganese catalase genes have been annotated in genomic databases, although the assignment of many of these putative manganese catalases needs to be experimentally verified. Iron limitation, exposure to low levels of peroxide stress, thermostability and cyanide resistance may provide the biological and environmental context for the occurrence of manganese catalases.

Structural models of the manganese complex of photosystem II and mechanistic implications

Photosynthetic water oxidation and O₂ formation are catalyzed by a Mn₄Ca complex bound to the proteins of photosystem II (PSII). The catalytic site, including the inorganic Mn₄CaO(n)H(x) core and its protein environment, is denoted as oxygen-evolving complex (OEC). Earlier and recent progress in the endeavor to elucidate the structure of the OEC is reviewed, with focus on recent results obtained by (i) X−ray spectroscopy (specifically by EXAFS analyses), and (ii) X-ray diffraction (XRD, protein crystallography). Very recently, an impressive resolution of 1.9Å has been achieved by XRD. Most likely however, all XRD data on the Mn₄CaO(n)H(x) core of the OEC are affected by X-ray induced modifications (radiation damage). Therefore and to address (important) details of the geometric and electronic structure of the OEC, a combined analysis of XRD and XAS data has been approached by several research groups. These efforts are reviewed and extended using an especially comprehensive approach. Taking into account XRD results on the protein environment of the inorganic core of the Mn complex, 12 alternative OEC models are considered and evaluated by quantitative comparison to (i) extended-range EXAFS data, (ii) polarized EXAFS of partially oriented PSII membrane particles, and (iii) polarized EXAFS of PSII crystals. We conclude that there is a class of OEC models that is in good agreement with both the recent crystallographic models and the XAS data. On these grounds, mechanistic implications for the O−O bond formation chemistry are discussed. This article is part of a Special Issue entitled: Photosystem II.

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.

Participation of glutamate-354 of the CP43 polypeptide in the ligation of manganese and the binding of substrate water in photosystem II

In the current X-ray crystallographic structural models of photosystem II, Glu354 of the CP43 polypeptide is the only amino acid ligand of the oxygen-evolving Mn(4)Ca cluster that is not provided by the D1 polypeptide. To further explore the influence of this structurally unique residue on the properties of the Mn(4)Ca cluster, the CP43-E354Q mutant of the cyanobacterium Synechocystis sp. PCC 6803 was characterized with a variety of biophysical and spectroscopic methods, including polarography, EPR, X-ray absorption, FTIR, and mass spectrometry. The kinetics of oxygen release in the mutant were essentially unchanged from those in wild type. In addition, the oxygen flash yields exhibited normal period four oscillations having normal S state parameters, although the yields were lower, correlating with the mutant's lower steady-state rate (approximately 20% compared to wild type). Experiments conducted with H(2)(18)O showed that the fast and slow phases of substrate water exchange in CP43-E354Q thylakoid membranes were accelerated 8.5- and 1.8-fold, respectively, in the S(3) state compared to wild type. Purified oxygen-evolving CP43-E354Q PSII core complexes exhibited a slightly altered S(1) state Mn-EXAFS spectrum, a slightly altered S(2) state multiline EPR signal, a substantially altered S(2)-minus-S(1) FTIR difference spectrum, and an unusually long lifetime for the S(2) state (>10 h) in a substantial fraction of reaction centers. In contrast, the S(2) state Mn-EXAFS spectrum was nearly indistinguishable from that of wild type. The S(2)-minus-S(1) FTIR difference spectrum showed alterations throughout the amide and carboxylate stretching regions. Global labeling with (15)N and specific labeling with l-[1-(13)C]alanine revealed that the mutation perturbs both amide II and carboxylate stretching modes and shifts the symmetric carboxylate stretching modes of the α-COO(-) group of D1-Ala344 (the C-terminus of the D1 polypeptide) to higher frequencies by 3-4 cm(-1) in both the S(1) and S(2) states. The EPR and FTIR data implied that 76-82% of CP43-E354Q PSII centers can achieve the S(2) state and that most of these can achieve the S(3) state, but no evidence for advancement beyond the S(3) state was observed in the FTIR data, at least not in a majority of PSII centers. Although the X-ray absorption and EPR data showed that the CP43-E354Q mutation only subtly perturbs the structure and spin state of the Mn(4)Ca cluster in the S(2) state, the FTIR and H(2)(18)O exchange data show that the mutation strongly influences other properties of the Mn(4)Ca cluster, altering the response of numerous carboxylate and amide groups to the increased positive charge that develops on the cluster during the S(1) to S(2) transition and weakening the binding of both substrate water molecules (or water-derived ligands), especially the one that exchanges rapidly in the S(3) state. The FTIR data provide evidence that CP43-Glu354 coordinates to the Mn(4)Ca cluster in the S(1) state as a bridging ligand between two metal ions but provide no compelling evidence that this residue changes its coordination mode during the S(1) to S(2) transition. The H(2)(18)O exchange data provide evidence that CP43-Glu354 interacts with the Mn ion that ligates the substrate water molecule (or water-derived ligand) that is in rapid exchange in the S(3) state.

Solid-state (55)Mn NMR spectroscopy of bis(μ-oxo)dimanganese(IV) [Mn(2)O(2)(salpn)(2)], a model for the oxygen evolving complex in photosystem II

We have examined the antiferromagneticly coupled bis(μ-oxo)dimanganese(IV) complex [Mn(2)O(2)(salpn)(2)] (1) with (55)Mn solid-state NMR at cryogenic temperatures and first-principle theory. The extracted values of the (55)Mn quadrupole coupling constant, C(Q), and its asymmetry parameter, η(Q), for 1 are 24.7 MHz and 0.43, respectively. Further, there was a large anisotropic contribution to the shielding of each Mn(4+), i.e. a Δσ of 3375 ppm. Utilizing broken symmetry density functional theory, the predicted values of the electric field gradient (EFG) or equivalently the C(Q) and η(Q) at ZORA, PBE QZ4P all electron level of theory are 23.4 MHz and 0.68, respectively, in good agreement with experimental observations.

Evidence from FTIR difference spectroscopy of an extensive network of hydrogen bonds near the oxygen-evolving Mn(4)Ca cluster of photosystem II involving D1-Glu65, D2-Glu312, and D1-Glu329

Analyses of the refined X-ray crystallographic structures of photosystem II (PSII) at 2.9-3.5 A have revealed the presence of possible channels for the removal of protons from the catalytic Mn(4)Ca cluster during the water-splitting reaction. As an initial attempt to verify these channels experimentally, the presence of a network of hydrogen bonds near the Mn(4)Ca cluster was probed with FTIR difference spectroscopy in a spectral region sensitive to the protonation states of carboxylate residues and, in particular, with a negative band at 1747 cm(-1) that is often observed in the S(2)-minus-S(1) FTIR difference spectrum of PSII from the cyanobacterium Synechocystis sp. PCC 6803. On the basis of its 4 cm(-1) downshift in D(2)O, this band was assigned to the carbonyl stretching vibration (C horizontal lineO) of a protonated carboxylate group whose pK(a) decreases during the S(1) to S(2) transition. The positive charge that forms on the Mn(4)Ca cluster during the S(1) to S(2) transition presumably causes structural perturbations that are transmitted to this carboxylate group via electrostatic interactions and/or an extended network of hydrogen bonds. In an attempt to identify the carboxylate group that gives rise to this band, the FTIR difference spectra of PSII core complexes from the mutants D1-Asp61Ala, D1-Glu65Ala, D1-Glu329Gln, and D2-Glu312Ala were examined. In the X-ray crystallographic models, these are the closest carboxylate residues to the Mn(4)Ca cluster that do not ligate Mn or Ca and all are highly conserved. The 1747 cm(-1) band is present in the S(2)-minus-S(1) FTIR difference spectrum of D1-Asp61Ala but absent from the corresponding spectra of D1-Glu65Ala, D2-Glu312Ala, and D1-Glu329Gln. The band is also sharply diminished in magnitude in the wild type when samples are maintained at a relative humidity of </=85%. It is proposed that D1-Glu65, D2-Glu312, and D1-Glu329 participate in a common network of hydrogen bonds that includes water molecules and the carboxylate group that gives rise to the 1747 cm(-1) band. It is further proposed that the mutation of any of these three residues, or partial dehydration caused by maintaining samples at a relative humidity of <or=85%, disrupts the network sufficiently that the structural perturbations associated with the S(1) to S(2) transition are no longer transmitted to the carboxylate group that gives rise to the 1747 cm(-1) band. Because D1-Glu329 is located approximately 20 A from D1-Glu65 and D2-Glu312, the postulated network of hydrogen bonds must extend for at least 20 A across the lumenal face of the Mn(4)Ca cluster. The D1-Asp61Ala, D1-Glu65Ala, and D2-Glu312Ala mutations also appear to substantially decrease the fraction of PSII reaction centers that undergo the S(3) to S(0) transition in response to a saturating flash. This behavior is consistent with D1-Asp61, D1-Glu65, and D2-Glu312 participating in a dominant proton egress channel that links the Mn(4)Ca cluster with the thylakoid lumen.

Tracking the flow of water through photosystem II using molecular dynamics and streamline tracing

The CaMn(4) cluster of the oxygen-evolving complex (OEC) of photosynthesis catalyzes the light-driven splitting of water into molecular oxygen, protons, and electrons. The OEC is buried within photosystem II (PSII), a multisubunit integral membrane protein complex, and water must find its way to the CaMn(4) cluster by moving through protein. Channels for water entrance, and proton and oxygen exit, have previously been proposed following the analysis of cavities found within X-ray structures of PSII. However, these analyses do not account for the dynamic motion of proteins and cannot track the movement of water within PSII. To study water dynamics in PSII, we performed molecular dynamics simulations and developed a novel approach for the visualization of water diffusion within protein based on a streamline tracing algorithm used in fluid dynamics and diffusion tensor imaging. We identified a system of branching pathways of water diffusion in PSII leading to the OEC that connect to a number of distinct entrance points on the lumenal surface. We observed transient changes in the connections between channels and entrance points that served to moderate both the flow of water near the OEC and the exchange of water inside and outside of the protein. Water flow was significantly altered in simulations lacking the OEC which were characterized by a simpler and wider channel with only two openings, consistent with the creation of an ion channel that allows entry of Mn(2+), Ca(2+), and Cl(-) as required for construction of the CaMn(4) cluster.

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.

Computational analysis of photosynthetic systems

The use of various computational techniques for the study of photosynthetic systems is described ranging from genome analysis to density functional simulations of the oxygen evolving complex of PSII. The use of simulations for analyzing protein structures can aid in clarifying ambiguous and incomplete experimental results to identifying underlying rules to create efficient light-initiated charge separation at high efficiency.