Marius Retegan

Electronic structure of the oxygen-evolving complex in photosystem II prior to O-O bond formation

Setting the stage for release of oxygen Plants transform water into the oxygen we breathe using a protein-bound cluster of four manganese (Mn) ions and a calcium ion. Cox et al. now establish the precise electronic structure in that cluster immediately before formation of the O-O bond (see the Perspective by Britt and Oyala). Using the technique of electron paramagnetic resonance spectroscopy, they confirm a hypothesis that all four Mn ions are octahedrally coordinated and in the 4+ oxidation state. Such clues to the efficiency of the photosynthetic process, so essential to life on Earth, may also facilitate the development of artificial waters-plitting catalysts. Science, this issue p. 804; see also p. 736 Electron paramagnetic resonance spectroscopy reveals that photosynthetic production of O2 proceeds from four MnIV ions. [Also see Perspective by Britt and Oyala] The photosynthetic protein complex photosystem II oxidizes water to molecular oxygen at an embedded tetramanganese-calcium cluster. Resolving the geometric and electronic structure of this cluster in its highest metastable catalytic state (designated S3) is a prerequisite for understanding the mechanism of O-O bond formation. Here, multifrequency, multidimensional magnetic resonance spectroscopy reveals that all four manganese ions of the catalyst are structurally and electronically similar immediately before the final oxygen evolution step; they all exhibit a 4+ formal oxidation state and octahedral local geometry. Only one structural model derived from quantum chemical modeling is consistent with all magnetic resonance data; its formation requires the binding of an additional water molecule. O-O bond formation would then proceed by the coupling of two proximal manganese-bound oxygens in the transition state of the cofactor.

Toward Atomistic Resolution Structure of Phosphatidylcholine Headgroup and Glycerol Backbone at Different Ambient Conditions

Phospholipids are essential building blocks of biological membranes. Despite a vast amount of very accurate experimental data, the atomistic resolution structures sampled by the glycerol backbone and choline headgroup in phoshatidylcholine bilayers are not known. Atomistic resolution molecular dynamics simulations have the potential to resolve the structures, and to give an arrestingly intuitive interpretation of the experimental data, but only if the simulations reproduce the data within experimental accuracy. In the present work, we simulated phosphatidylcholine (PC) lipid bilayers with 13 different atomistic models, and compared simulations with NMR experiments in terms of the highly structurally sensitive C–H bond vector order parameters. Focusing on the glycerol backbone and choline headgroups, we showed that the order parameter comparison can be used to judge the atomistic resolution structural accuracy of the models. Accurate models, in turn, allow molecular dynamics simulations to be used as an interpretation tool that translates these NMR data into a dynamic three-dimensional representation of biomolecules in biologically relevant conditions. In addition to lipid bilayers in fully hydrated conditions, we reviewed previous experimental data for dehydrated bilayers and cholesterol-containing bilayers, and interpreted them with simulations. Although none of the existing models reached experimental accuracy, by critically comparing them we were able to distill relevant chemical information: (1) increase of choline order parameters indicates the P–N vector tilting more parallel to the membrane, and (2) cholesterol induces only minor changes to the PC (glycerol backbone) structure. This work has been done as a fully open collaboration, using nmrlipids.blogspot.fi as a communication platform; all the scientific contributions were made publicly on this blog. During the open research process, the repository holding our simulation trajectories and files (https://zenodo.org/collection/user-nmrlipids) has become the most extensive publicly available collection of molecular dynamics simulation trajectories of lipid bilayers.

Structure, ligands and substrate coordination of the oxygen-evolving complex of photosystem II in the S2 state: a combined EPR and DFT study

The S2 state of the oxygen-evolving complex of photosystem II, which consists of a Mn4O5Ca cofactor, is EPR-active, typically displaying a multiline signal, which arises from a ground spin state of total spin ST = 1/2. The precise appearance of the signal varies amongst different photosynthetic species, preparation and solvent conditions/compositions. Over the past five years, using the model species Thermosynechococcus elongatus, we have examined modifications that induce changes in the multiline signal, i.e. Ca(2+)/Sr(2+)-substitution and the binding of ammonia, to ascertain how structural perturbations of the cluster are reflected in its magnetic/electronic properties. This refined analysis, which now includes high-field (W-band) data, demonstrates that the electronic structure of the S2 state is essentially invariant to these modifications. This assessment is based on spectroscopies that examine the metal centres themselves (EPR, (55)Mn-ENDOR) and their first coordination sphere ligands ((14)N/(15)N- and (17)O-ESEEM, -HYSCORE and -EDNMR). In addition, extended quantum mechanical models from broken-symmetry DFT now reproduce all EPR, (55)Mn and (14)N experimental magnetic observables, with the inclusion of second coordination sphere ligands being crucial for accurately describing the interaction of NH3 with the Mn tetramer. These results support a mechanism of multiline heterogeneity reported for species differences and the effect of methanol [Biochim. Biophys. Acta, Bioenerg., 2011, 1807, 829], involving small changes in the magnetic connectivity of the solvent accessible outer MnA4 to the cuboidal unit Mn3O3Ca, resulting in predictable changes of the measured effective (55)Mn hyperfine tensors. Sr(2+) and NH3 replacement both affect the observed (17)O-EDNMR signal envelope supporting the assignment of O5 as the exchangeable μ-oxo bridge and it acting as the first site of substrate inclusion.

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.

Principles of Natural Photosynthesis

Nature relies on a unique and intricate biochemical setup to achieve sunlight-driven water splitting. Combined experimental and computational efforts have produced significant insights into the structural and functional principles governing the operation of the water-oxidizing enzyme Photosystem II in general, and of the oxygen-evolving manganese-calcium cluster at its active site in particular. Here we review the most important aspects of biological water oxidation, emphasizing current knowledge on the organization of the enzyme, the geometric and electronic structure of the catalyst, and the role of calcium and chloride cofactors. The combination of recent experimental work on the identification of possible substrate sites with computational modeling have considerably limited the possible mechanistic pathways for the critical O-O bond formation step. Taken together, the key features and principles of natural photosynthesis may serve as inspiration for the design, development, and implementation of artificial systems.

Azurin as a Protein Scaffold for a Low-coordinate Nonheme Iron Site with a Small-molecule Binding Pocket.

The apoprotein of Pseudomonas aeruginosa azurin binds iron(II) to give a 1:1 complex, which has been characterized by electronic absorption, Mössbauer, and NMR spectroscopies, as well as X-ray crystallography and quantum-chemical computations. Despite potential competition by water and other coordinating residues, iron(II) binds tightly to the low-coordinate site. The iron(II) complex does not react with chemical redox agents to undergo oxidation or reduction. Spectroscopically calibrated quantum-chemical computations show that the complex has high-spin iron(II) in a pseudotetrahedral coordination environment, which features interactions with side chains of two histidines and a cysteine as well as the C═O of Gly45. In the (5)A(1) ground state, the d(z(2)) orbital is doubly occupied. Mutation of Met121 to Ala leaves the metal site in a similar environment but creates a pocket for reversible binding of small anions to the iron(II) center. Specifically, azide forms a high-spin iron(II) complex and cyanide forms a low-spin iron(II) complex.

A First-Principles Approach to the Calculation of the on-Site Zero-Field Splitting in Polynuclear Transition Metal Complexes

The interpretation of electron paramagnetic resonance spectra of polynuclear transition metal complexes in terms of individual contributions from each paramagnetic center can be greatly facilitated by the availability of theoretical methods that enable the reliable prediction of local spectroscopic parameters. In this work we report an approach that enables the application of multireference ab initio methods for the calculation of local zero field splitting tensors, one of the leading terms in the spin Hamiltonian for exchange-coupled systems of high nuclearity. The method referred to as local complete active space configuration interaction (L-CASCI) represents a multireference calculation with an active space composed of local orbitals of the center of interest. By successive permutation of the active space to include the localized orbitals corresponding to a particular center of the complex, all on-site parameters can be easily obtained at a high-level of theory with a corresponding low computational cost. Benchmark calculations on synthetic complexes confirm the validity of the approach. As an example of the applicability of the L-CASCI method to large systems, we determine the local anisotropy of the Mn(III) ion of the tetranuclear manganese cluster of photosystem II in both structural forms of its S2 state.

Differences in the Active Site of Water Oxidation among Photosynthetic Organisms

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

The iron–sulfur core in Rieske proteins is not symmetric

At variance with ferredoxins, Rieske-type proteins contain a chemically asymmetric iron–sulfur cluster. Nevertheless, X-ray crystallography apparently finds their [2Fe–2S] cores to be structurally symmetric or very close to symmetric (i.e. the four iron–sulfur bonds in the [2Fe–2S] core are equidistant). Using a combination of advanced density-based approaches, including finite-temperature molecular dynamics to access thermal fluctuations and free-energy profiles, in conjunction with correlated wavefunction-based methods we clearly predict an asymmetric core structure. This reveals a fundamental and intrinsic difference within the iron–sulfur clusters hosted by Rieske proteins and ferredoxins and thus opens up a new dimension for the ongoing efforts in understanding the role of Rieske-type [2Fe–2S] cluster in electron transfer processes that occur in almost all biological systems.

Time-Resolved Electron Paramagnetic Resonance and Theoretical Investigations of Metal-Free Room-Temperature Triplet Emitters

Utilization of triplets is important for preparing organic light-emitting diodes with high efficiency. Very recently, both electrophosphorescence and electrofluorescence could be observed at room temperature for thienyl-substituted phenazines without any heavy metals ( Ratzke et al. J. Phys. Chem. Lett. , 2016 , 7 , 4802 ). It was found that the phosphorescence efficiency depends on the orientation of fused thiophenes. In this work, the thienyl-substituted phenazines are investigated in more detail by time-resolved electron paramagnetic resonance (EPR) and quantum chemical calculations. Spin dynamics, zero-field splitting constants, and electron-spin structures of the excited triplet states for the metal-free room-temperature triplet emitters are correlated with phosphorescence efficiency. Complete active space self-consistent field (CASSCF) calculations clearly show that the electron spin density distributions of the first excited triplet states are strongly affected by the molecular geometry. For the phosphorescent molecules, the electron spins are localized on the phenazine unit, in which the sulfur atom of the fused thiophene points upward. The electron spins are delocalized onto the thiophene unit just by changing the orientation of the fused thiophenes from upward to downward, resulting in the suppression of phosphorescence. Time-resolved EPR measurements and time-dependent density functional theory (TD-DFT) calculations demonstrate that the electron spins delocalized onto the thiophene unit lead to the acceleration of nonradiative decays, in conjunction with the narrowing of the singlet-triplet energy gap.