Daniele Bovi

The S2 State of the Oxygen-Evolving Complex of Photosystem II Explored by QM/MM Dynamics: Spin Surfaces and Metastable States Suggest a Reaction Path Towards the S3 State†

One of the key steps in photosynthetic solar-energy conversion performed by plants, algae, and cyanobacteria is the splitting of water into molecular oxygen and hydrogen equivalents.[1] To achieve this challenging task photosynthetic organisms use a protein complex that remained almost unchanged during the evolution in the last two and a half billion years: the photosystem II (PSII). The reaction proceeds by the accumulation of four oxidizing equivalents on the {Mn4CaO5} cluster through five (S0–S4) oxidation states that are sequentially attained during water splitting (Kok cycle).[2] The deep understanding of the way nature has found to perform this difficult task efficiently has a great relevance not only for biology but also for inspiring the development of biomimetic artificial systems that can be used to store solar energy in an environmentally friendly way.[3] Atomic details of the structure of the oxygen-evolving complex (OEC) of PSII have been revealed by extended X-ray absorption fine structure (EXAFS) experiments and by X-ray crystallography at increasing resolution levels.[4] However, the accurate position of the {Mn4CaO5} cluster atoms and its ligands emerged only when a X-ray structure at 1.9 A resolution became accessible.[5] However, the effect of a possible X-ray photo-reduction, in particular on the characterization of the Kok’s state described by this structure and on the unrealistic bond lengths between the oxygen atom O5 and the two manganese ions Mn1 and Mn4, is matter of debate.[6] Additionally, important contributions to the structure refinement came from theoretical studies.[6b,7]

Pathway for Mn-cluster oxidation by tyrosine-Z in the S2 state of photosystem II

Significance A key step in natural photosynthesis is the water-splitting reaction into molecular oxygen and hydrogen equivalents. Understanding the molecular mechanisms behind this photoreaction will unravel the secrets of solar energy conversion in biochemistry and may inspire the design of artificial biomimetic materials for green energy production. Photosynthetic water oxidation occurs in the Mn4Ca core of the photosystem II complex and proceeds through five subsequent steps S0 – S4 of the Kok cycle. Four electrons are sequentially removed from the Mn4Ca core by a nearby tyrosine, which is in turn oxidized by the photoactivated chlorophyll special pair. Using first principles multiscale atomistic simulations we clarify the thermodynamics and the kinetics for such electron abstraction in the S2 state. Water oxidation in photosynthetic organisms occurs through the five intermediate steps S0–S4 of the Kok cycle in the oxygen evolving complex of photosystem II (PSII). Along the catalytic cycle, four electrons are subsequently removed from the Mn4CaO5 core by the nearby tyrosine Tyr-Z, which is in turn oxidized by the chlorophyll special pair P680, the photo-induced primary donor in PSII. Recently, two Mn4CaO5 conformations, consistent with the S2 state (namely, S2A and S2B models) were suggested to exist, perhaps playing a different role within the S2-to-S3 transition. Here we report multiscale ab initio density functional theory plus U simulations revealing that upon such oxidation the relative thermodynamic stability of the two previously proposed geometries is reversed, the S2B state becoming the leading conformation. In this latter state a proton coupled electron transfer is spontaneously observed at ∼100 fs at room temperature dynamics. Upon oxidation, the Mn cluster, which is tightly electronically coupled along dynamics to the Tyr-Z tyrosyl group, releases a proton from the nearby W1 water molecule to the close Asp-61 on the femtosecond timescale, thus undergoing a conformational transition increasing the available space for the subsequent coordination of an additional water molecule. The results can help to rationalize previous spectroscopic experiments and confirm, for the first time to our knowledge, that the water-splitting reaction has to proceed through the S2B conformation, providing the basis for a structural model of the S3 state.

Mechanism of Water Delivery to the Active Site of Photosystem II along the S2 to S3 Transition

The two water molecules serving as substrate for the oxygen evolution in Photosystem II are already bound in the S2 state of the Kok-Joliots cycle. Nevertheless, an additional water molecule is supposed to bind the cluster during the transition between the S2 and S3 states, which has been recently revealed to have the Mn4CaO5 catalytic cluster arranged in an open cubane fashion. In this Letter, by means of ab initio calculations, we investigated the possible pathways for the binding of the upcoming water molecule. Upon the four different possibilities checked in our calculations, the binding of the crystallographic water molecule, originally located nearby the Cl(-) binding site, showed the lowest activation energy barrier. Our findings therefore support the view in which the W2 hydroxyl group and the O5 oxygen act as substrates for the oxygen evolution. Within this framework the role of the open and closed Mn4CaO5 conformers is clarified as well as the exact mechanistic events occurring along the S2 to S3 transition.

Reorganization of substrate waters between the closed and open cubane conformers during the S2 to S3 transition in the oxygen evolving complex.

A crucial step in the mechanism for oxygen evolution in the Photosystem II complex resides in the transition from the S2 state to the S3 state of Kok–Joliot’s cycle, in which an additional water molecule binds to the cluster. On the basis of computational chemistry calculations on Photosystem II models, we propose a reorganization mechanism involving a hydroxyl (W2) and a μ2-oxo bridge (O5) that is able to link the closed cubane S2B intermediate conformer to the S3 open cubane structure. This mechanism can reconcile the apparent conflict between recently reported water exchange and electron paramagnetic resonance experiments, and theoretical studies.

Fermi resonance as a tool for probing peridinin environment.

In the present paper, we provide an extended study of the vibrational signature of a butenolide carotenoid, peridinin, in various solvents by combining resonance Raman spectroscopy (RRS) with theoretical calculations. The presence of a Fermi resonance due to coupling between the lactonic C═O stretching and the overtone of the wagging of the C-H in the lactonic ring provides a spectroscopic way of differentiating between peridinins lying in different environments. This is a significant achievement, given that simultaneous presence of several peridinins (each with a peculiar photophysical role) in different environments occurs in the most important peridinin containing proteins, the peridinin-chlorophyll proteins (PCPs) and the Chl a-c2-peridinin binding proteins. In RRS, small modifications of solvent polarity can give rise to large differences in the intensity and splitting between the two bands, resulting from the Fermi resonance. By changing the polarity, we can tune the frequency of stretching of the C═O and, while the C-H wagging frequency is almost always constant in different solvents, move the system from a perfect resonance condition to off-resonance ones. We have corroborated our spectroscopic findings with a quasi-classical dynamical model of two coupled oscillators, and DFT calculations on peridinin in different solvents; we have also used calculations to complete the peridinin vibrational mode assignments in the 800-1600 cm(-1) region of RRS spectra, corresponding to polyene chain motion. Finally, the presence of Fermi resonance has been used to reinterpret previous vibrational spectroscopic experiments in PCPs.

Magnetic interactions in the catalyst used by nature to split water: a DFT+U multiscale study on the Mn 4 CaO 5 core in photosystem II

An important approach in the design of new environmentally friendly materials is represented by the study of analogous systems already existing in nature. In the search for new water splitting catalysts, the corresponding natural analogue is represented by the oxygen-evolving complex of photosystem II, which is a large membrane protein complex present in photosynthetic organisms. The understanding of the catalytic strategy of its active Mn4CaO5 core is important to unravel the mechanisms of water oxidation in photosynthesis and can serve as an inspiring model for the design of biomimetic catalysts based on largely non-toxic, earth abundant elements. The magnetic interactions between Mn ions are studied in the present work by means of DFT+U broken symmetry ab initio molecular dynamics within a quantum mechanics/molecular mechanics framework. The room temperature dynamics of two different structural models (i.e. with total high-spin and total low-spin ground states) was stable during the simulated time. We observed large fluctuations of the magnetic coupling

Characterization of the Sr(2+)- and Cd(2+)-Substituted Oxygen-Evolving Complex of Photosystem II by Quantum Mechanics/Molecular Mechanics Calculations.

The Mn4CaO5 cluster in the oxygen-evolving complex is the catalytic core of the Photosystem II (PSII) enzyme, responsible for the water splitting reaction in oxygenic photosynthesis. The role of the redox-inactive ion in the cluster has not yet been fully clarified, although several experimental data are available on Ca2+-depleted and Ca2+-substituted PSII complexes, indicating Sr2+-substituted PSII as the only modification that preserves oxygen evolution. In this work, we investigated the structural and electronic properties of the PSII catalytic core with Ca2+ replaced with Sr2+ and Cd2+ in the S2 state of the Kok−Joliot cycle by means of density functional theory and ab initio molecular dynamics based on a quantum mechanics/ molecular mechanics approach. Our calculations do not reveal significant differences between the substituted and wild-type systems in terms of geometries, thermodynamics, and kinetics of two previously identified intermediate states along the S2 to S3 transition, namely, the open cubane S2 A and closed cubane S2 B conformers. Conversely, our calculations show different pKa values for the water molecule bound to the three investigated heterocations. Specifically, for Cd-substituted PSII, the pKa value is 5.3 units smaller than the respective value in wild type Ca-PSII. On the basis of our results, we conclude that, assuming all the cations sharing the same binding site, the induced difference in the acidity of the binding pocket might influence the hydrogen bonding network and the redox levels to prevent the further evolution of the cycle toward the S3 state.

The Impact of Nitric Oxide Toxicity on the Evolution of the Glutathione Transferase Superfamily: A PROPOSAL FOR AN EVOLUTIONARY DRIVING FORCE*

Background: Why do ancestral GSTs utilize cysteine/serine as catalytic residues, whereas more recently evolved GSTs utilize tyrosine? Results: Only the more recently evolved GSTs display enough affinity to bind and make harmless the toxic DNDGIC (a natural NO carrier). Conclusion: GST evolution could be linked to the defense against NO. Significance: This represents a further piece in the puzzle of evolutive adaptation to NO toxicity. Glutathione transferases (GSTs) are protection enzymes capable of conjugating glutathione (GSH) to toxic compounds. During evolution an important catalytic cysteine residue involved in GSH activation was replaced by serine or, more recently, by tyrosine. The utility of these replacements represents an enigma because they yield no improvements in the affinity toward GSH or in its reactivity. Here we show that these changes better protect the cell from nitric oxide (NO) insults. In fact the dinitrosyl·diglutathionyl·iron complex (DNDGIC), which is formed spontaneously when NO enters the cell, is highly toxic when free in solution but completely harmless when bound to GSTs. By examining 42 different GSTs we discovered that only the more recently evolved Tyr-based GSTs display enough affinity for DNDGIC (KD < 10−9 m) to sequester the complex efficiently. Ser-based GSTs and Cys-based GSTs show affinities 102–104 times lower, not sufficient for this purpose. The NO sensitivity of bacteria that express only Cys-based GSTs could be related to the low or null affinity of their GSTs for DNDGIC. GSTs with the highest affinity (Tyr-based GSTs) are also over-represented in the perinuclear region of mammalian cells, possibly for nucleus protection. On the basis of these results we propose that GST evolution in higher organisms could be linked to the defense against NO.

Dynamics of the Special Pair of Chlorophylls of Photosystem II

Cholophylls are at the basis of the photosynthetic energy conversion mechanisms in algae, plants, and cyanobacteria. In photosystem II, the photoproduced electrons leave a special pair of chlorophylls (namely, P(D1) and P(D2)) that becomes cationic. This oxidizing pair [P(D1),P(D2)](+), in turn, triggers a cascade of oxidative events, eventually leading to water splitting and oxygen evolution. In the present work, using quantum mechanics/molecular mechanics calculations, we investigate the electronic structure and the dynamics of the P(D1)P(D2) special pair in both its oxidized and reduced states. In agreement with previously reported static calculations, the symmetry between the two chlorophylls was found to be broken, the positive charge being preferentially located on P(D1). Nevertheless, this study reveals for the first time that large charge fluctuations occur along dynamics, temporarily inverting the charge preference for the two branches. Finally, a vibrational analysis pinpointed that such charge fluctuations are strongly coupled to specific modes of the special pair.

Effects of Static Correlation between Spin Centers in Multicenter Transition Metal Complexes

Multicenter transition metal complexes are the key moieties of many processes in chemistry, biochemistry, and materials science such as in the active sites of enzymes, molecular catalysts, and biological electron carriers. Their electronic structure, often characterized by high-spin-polarized metal sites, is a challenge for theoretical chemists because of their high degree of dynamical and static correlation. Static correlation is necessary both for the appropriate description of the metal-ligand bonding and for a correct description of the multideterminant character arising from the magnetic interactions between spin centers. Density functional theory (DFT) is usually applied using a single-determinant broken-symmetry state that is lacking the correct spin symmetry when the ground state has total low-spin character. To alleviate this drawback, we use the extended broken-symmetry (EBS) approach to derive approximate ground-state energies and, for the first time, forces for the correctly symmetric ground state of an arbitrary number of spin centers within the framework of the Heisenberg-Dirac-van Vleck Hamiltonian. Remarkably, the proposed procedure supplies relaxed geometries that are fully consistent with the calculated J-coupling constants. We apply the method to investigate the relaxed geometrical structure of the low-spin ground state of iron-sulfur clusters with two, three, and four iron centers. We observed significant differences in both geometrical parameters and coupling constant J between the symmetrized ground state, the high-spin, and the broken-symmetry optimized structures. These changes are often comparable with the differences observed by using different functionals, and the use of EBS always improves the description of the studied systems. It will be therefore important to include it in any DFT attempt to quantitatively describe multicenter transition metal complexes in the future.