Distribution of the Cationic State over the Chlorophyll Pair of the Photosystem II Reaction Center

Temperature-Dependent Spin-Driven Dimerization Determines the Ultrafast Dynamics of a Copper(II)-Bound Tripyrrindione Radical

Radicals and other open-shell molecules play a central role in chemical transformations and redox chemistry. While radicals are often highly reactive, stable radical systems are desirable for a range of potential applications, ranging from materials chemistry and catalysis to spintronics and quantum information. Here we investigate the ultrafast properties of a stable radical system with temperature-dependent spin-tunable properties. This radical complex, Cu(II) hexaethyl tripyrrin-1,14-dione, accommodates unpaired electrons localized on both the copper metal center and the tripyrrolic ligand. The unusual combination of two unpaired electrons and high stability in this radical molecule enable switchable temperature-dependent spin coupling. Two-dimensional electronic spectroscopy measurements of Cu(II) hexaethyl tripyrrin-1,14-dione were collected at room temperature and at 77 K. At room temperature, the molecules are present as monomers and have short picosecond lifetimes. At 77 K, the molecules are present in a dimer form mediated by ferromagnetic and antiferromagnetic coupling. This reversible spin-driven dimerization changes the optical properties of the system, generating long-lived excitonic states.

Interplay of two low-barrier hydrogen bonds in long-distance proton-coupled electron transfer for water oxidation

D1-Tyr161 (TyrZ) forms a low-barrier H-bond with D1-His190 and functions as a redox-active group in photosystem II. When oxidized to the radical form (TyrZ-O•), it accepts an electron from the oxygen-evolving Mn4CaO5 cluster, facilitating an increase in the oxidation state (Sn; n = 0-3). In this study, we investigated the mechanism of how TyrZ-O• drives proton-coupled electron transfer during the S2 to S3 transition using a quantum mechanical/molecular mechanical approach. In response to TyrZ-O• formation and subsequent loss of the low-barrier H-bond, the ligand water molecule at the Ca2+ site (W4) reorients away from TyrZ and donates an H-bond to D1-Glu189 at Mn4 of Mn4CaO5 together with an adjacent water molecule. The H-bond donation to the Mn4CaO5 cluster triggers the release of the proton from the lowest pKa site (W1 at Mn4) along the W1…D1-Asp61 low-barrier H-bond, leading to protonation of D1-Asp61. The interplay of the two low-barrier H-bonds, involving the Ca2+ interface and forming the extended Grotthuss-like network [TyrZ…D1-His190]-[Mn4CaO5]-[W1…D1-Asp61], rather than the direct electrostatic interaction, is likely a basis of the apparent long-distance interaction (11.4 Å) between TyrZ-O• formation and D1-Asp61 protonation.

Triplet states in the reaction center of Photosystem II

In oxygenic photosynthesis sunlight is harvested and funneled as excitation energy into the reaction center (RC) of Photosystem II (PSII), the site of primary charge separation that initiates the photosynthetic electron transfer chain. The chlorophyll ChlD1 pigment of the RC is the primary electron donor, forming a charge-separated radical pair with the vicinal pheophytin PheoD1 (ChlD1+PheoD1-). To avert charge recombination, the electron is further transferred to plastoquinone QA, whereas the hole relaxes to a central pair of chlorophylls (PD1PD2), subsequently driving water oxidation. Spin-triplet states can form within the RC when forward electron transfer is inhibited or back reactions are favored. This can lead to formation of singlet dioxygen, with potential deleterious effects. Here we investigate the nature and properties of triplet states within the PSII RC using a multiscale quantum-mechanics/molecular-mechanics (QM/MM) approach. The low-energy spectrum of excited singlet and triplet states, of both local and charge-transfer nature, is compared using range-separated time-dependent density functional theory (TD-DFT). We further compute electron paramagnetic resonance properties (zero-field splitting parameters and hyperfine coupling constants) of relaxed triplet states and compare them with available experimental data. Moreover, the electrostatic modulation of excited state energetics and redox properties of RC pigments by the semiquinone QA- is described. The results provide a detailed electronic-level understanding of triplet states within the PSII RC and form a refined basis for discussing primary and secondary electron transfer, charge recombination pathways, and possible photoprotection mechanisms in PSII.

Redox Potentials of Iron-Sulfur Clusters in Type I Photosynthetic Reaction Centers

The electron transfer pathways in type I photosynthetic reaction centers, such as photosystem I (PSI) and reaction centers from green sulfur bacteria (GsbRC), are terminated by two Fe4S4 clusters, FA and FB. The protein structures are the basis of understanding how the protein electrostatic environment interacts with the Fe4S4 clusters and facilitates electron transfer. Using the protein structures, we calculated the redox potential (Em) values for FA and FB in PSI and GsbRC, solving the linear Poisson-Boltzmann equation. The FA-to-FB electron transfer is energetically downhill in the cyanobacterial PSI structure, while it is isoenergetic in the plant PSI structure. The discrepancy arises from differences in the electrostatic influences of conserved residues, including PsaC-Lys51 and PsaC-Arg52, located near FA. The FA-to-FB electron transfer is slightly downhill in the GsbRC structure. Em(FA) and Em(FB) exhibit similar levels upon isolation of the membrane-extrinsic PsaC and PscB subunits from the PSI and GsbRC reaction centers, respectively. The binding of the membrane-extrinsic subunit at the heterodimeric/homodimeric reaction center plays a key role in tuning Em(FA) and Em(FB).

Triplet Delocalization over the Reaction Center Chlorophylls in Photosystem II

The triplet state of chlorophyll formed by charge recombination in photosystem II (PSII) is a precursor of harmful singlet oxygen. Although main localization of the triplet state on the monomeric chlorophyll, Chl_D1, at cryogenic temperatures has been suggested, how the triplet state is delocalized on other chlorophylls remains unclear. Here, we investigated the distribution of the triplet state of chlorophyll in PSII using light-induced Fourier transform infrared (FTIR) difference spectroscopy. Measurements of triplet-minus-singlet FTIR difference spectra with PSII core complexes from cyanobacterial mutants, D1-V157H, D2-V156H, D2-H197A, and D1-H198A, in which the interactions of the 13¹-keto C═O groups of the reaction center chlorophylls, P_D1, P_D2, Chl_D1, and Chl_D2, respectively, were perturbed, identified the 13¹-keto C═O bands of the individual chlorophylls and showed that the triplet state is delocalized over all of these chlorophylls. It is suggested that the triplet delocalization plays important roles in the photoprotection and photodamage mechanisms in PSII.

Substitution of Ca2+ and changes in the H-bond network near the oxygen-evolving complex of photosystem II

Ca²⁺, which provides binding sites for ligand water molecules W3 and W4 in the Mn₄CaO₅ cluster, is a prerequisite for O₂ evolution in photosystem II (PSII). We report structural changes in the H-bond network and the catalytic cluster itself upon the replacement of Ca²⁺ with other alkaline earth metals, using a quantum mechanical/molecular mechanical approach. The small radius of Mg²⁺ makes W3 donate an H-bond to D1-Glu189 in Mg²⁺-PSII. If an additional water molecule binds at the large surface of Ba²⁺, it donates H-bonds to D1-Glu189 and the ligand water molecule at the dangling Mn, altering the H-bond network. The potential energy profiles of the H-bond between D1-Tyr161 (TyrZ) and D1-His190 and the interconversion between the open- and closed-cubane S₂ conformations remain substantially unaltered upon the replacement of Ca²⁺. Remarkably, the O5⋯Ca²⁺ distance is shortest among all O5⋯metal distances irrespective of the radius being larger than that of Mg²⁺. Furthermore, Ca²⁺ is the only alkaline earth metal that equalizes the O5⋯metal and O2⋯metal distances and facilitates the formation of the symmetric cubane structure.

Energetic Diversity in the Electron-Transfer Pathways of Type I Photosynthetic Reaction Centers

Photosynthetic reaction centers from heliobacteria (HbRC) and green sulfur bacteria (GsbRC) are homodimeric proteins and share a common ancestor with photosystem I (PSI), classified as type I reaction centers. Using the HbRC crystal structure, we calculated the redox potential (Em) values in the electron-transfer branches, solving the linear Poisson-Boltzmann equation and considering the protonation states of all titratable sites in the entire protein-pigment complex. Em(A-1) for bacteriochlorophyll g at the secondary site in HbRC (-1157 mV) is as low as Em(A-1) for chlorophyll a in PSI (-1173 mV). Em(A0/HbRC) is at the same level as Em(A0/GsbRC) and is 200 mV higher than Em(A0/PSI) due to the replacement of PsaA-Trp697/PsaB-Trp677 in PSI with PshA-Arg554 in HbRC. In contrast, Em(FX) for the Fe4S4 cluster in HbRC (-420 mV) is significantly higher than Em(FX) in GsbRC (-719 mV) and PSI (-705 mV) due to the absence of acidic residues that correspond to PscA-Asp634 in GsbRC and PsaB-Asp575 in PSI. It seems likely that type I reaction centers have evolved, adopting (bacterio)chlorophylls suitable for their light environments while maintaining electron-transfer cascades.

Alternative Fast and Slow Primary Charge-Separation Pathways in Photosystem II

Photosystem-II (PSII) is a multi-subunit protein complex that harvests sunlight to perform oxygenic photosynthesis. Initial light-activated charge separation takes place at a reaction centre consisting of four chlorophylls and two pheophytins. Understanding the processes following light excitation remains elusive due to spectral congestion, the ultrafast nature, and multi-component behaviour of the charge-separation process. Here, using advanced computational multiscale approaches which take into account the large-scale configurational flexibility of the system, we identify two possible primary pathways to radical-pair formation that differ by three orders of magnitude in their kinetics. The fast (short-range) pathway is dominant, but the existence of an alternative slow (long-range) charge-separation pathway hints at the evolution of redundancy that may serve other purposes, adaptive or protective, related to formation of the unique oxidative species that drives water oxidation in PSII.

Electron-Transfer Route in the Early Oxidation States of the Mn4CaO5 Cluster in Photosystem II

The electron transfer from the oxygen-evolving Mn4CaO5 cluster to the electron acceptor D1-Tyr161 (TyrZ) is a prerequisite for water oxidation and O2 evolution. Here, we analyzed the electronic coupling in the rate-limiting electron-transfer transitions using a combined quantum mechanical/molecular mechanical/polarizable continuum model approach. In the S0 to S1 transition, the electronic coupling between the electron-donor Mn3(III) and TyrZ is small (2 meV). In contrast, the electronic coupling between the dangling Mn4(III) and TyrZ is significantly large (172 meV), which suggests that the electron transfer proceeds from Mn3(III) to TyrZ via Mn4(III). In the S1 to S2 transition, the electronic coupling between Mn4(III) and TyrZ is also larger (124 meV) than that between Mn1(III) and TyrZ (1 meV), which favors the formation of the open-cubane S2 conformation with Mn4(IV) over the formation of the closed-cubane S2 conformation with Mn1(IV). In the S0 to S1 and S1 to S2 transitions, the Mn4 d-orbital and the TyrZ π-orbital are hybridized via D1-Asp170, which suggests that D1-Asp170 commonly provides a dominant electron-transfer route.

Release of a Proton and Formation of a Low-Barrier Hydrogen Bond between Tyrosine D and D2-His189 in Photosystem II

In photosystem II (PSII), the second-lowest oxidation state (S₁) of the oxygen-evolving Mn₄CaO₅ cluster is the most stable, as the radical form of the redox-active D2-Tyr160 is considered to be a candidate that accepts an electron from the lowest oxidation state (S₀) in the dark. Using quantum mechanical/molecular mechanical calculations, we investigated the redox potential (E _m) of TyrD and its H-bond partner, D2-His189. The potential energy profile indicates that the release of a proton from the TyrD...D2-His189 pair leads to the formation of a low-barrier H-bond. The E _m depends on the H⁺ position along the low-barrier H-bond, e.g., 680 mV when the H⁺ is at the D2-His189 moiety and 800 mV when the H⁺ is at the TyrD moiety, which can explain why TyrD mediates both the S₀ to S₁ oxidation and the S₂ to S₁ reduction.

The initial charge separation step in oxygenic photosynthesis

Photosystem II is crucial for life on Earth as it provides oxygen as a result of photoinduced electron transfer and water splitting reactions. The excited state dynamics of the photosystem II-reaction center (PSII-RC) has been a matter of vivid debate because the absorption spectra of the embedded chromophores significantly overlap and hence it is extremely difficult to distinguish transients. Here, we report the two-dimensional electronic-vibrational spectroscopic study of the PSII-RC. The simultaneous resolution along both the visible excitation and infrared detection axis is crucial in allowing for the character of the excitonic states and interplay between them to be clearly distinguished. In particular, this work demonstrates that the mixed exciton-charge transfer state, previously proposed to be responsible for the far-red light operation of photosynthesis, is characterized by the Chl_D1⁺Phe radical pair and can be directly prepared upon photoexcitation. Further, we find that the initial electron acceptor in the PSII-RC is Phe, rather than P_D1, regardless of excitation wavelength.

The photosystem II reaction center (PSII-RC) is a model system to understand the initial steps of photosynthesis, but its excited state dynamics is difficult to disentangle with most spectroscopic methods. Here the authors perform a two-dimensional electronic-vibrational spectroscopic study of PSII-RC, providing detailed insight into such dynamics and into the mechanism of charge separation.

The Electronic Origin of Far-Red-Light-Driven Oxygenic Photosynthesis

Photosystem‐II uses sunlight to trigger charge separation and catalyze water oxidation. Intrinsic properties of chlorophyll a pigments define a natural “red limit” of photosynthesis at ≈680 nm. Nevertheless, charge separation can be triggered with far‐red photons up to 800 nm, without altering the nature of light‐harvesting pigments. Here we identify the electronic origin of this remarkable phenomenon using quantum chemical and multiscale simulations on a native Photosystem‐II model. We find that the reaction center is preorganized for charge separation in the far‐red region by specific chlorophyll–pheophytin pairs, potentially bypassing the light‐harvesting apparatus. Charge transfer can occur along two distinct pathways with one and the same pheophytin acceptor (Pheo_D1). The identity of the donor chlorophyll (Chl_D1 or P_D1) is wavelength‐dependent and conformational dynamics broaden the sampling of the far‐red region by the two charge‐transfer states. The two pathways rationalize spectroscopic observations and underpin designed extensions of the photosynthetically active radiation limit.

Effect of Hypoxia on Photosystem II of Tropical Seagrass Enhalus acoroides in the Dark

Hypoxia induced by eutrophication has become an important factor threatening the survival of coastal life such as Enhalus acoroides. The purpose of the current study was to explore the effect of hypoxia on photosystem II (PSII) of E. acoroides in the dark. The results showed that long-term dark hypoxia damages PSII activity of E. acoroides. The lower the oxygen content and the longer the hypoxic duration, the more seriously PSII was damaged and the less light-independent recovery parts of the damaged PSII. The damage to PSII caused by hypoxia was unrelated to ROS but related to respiration, because the respiration rate decreased with the decrease of oxygen content and PSII activity decreased significantly even at a normal oxygen content after the inhibition of aerobic respiratory pathway. Hypoxia reduced energy fluxes between the antennas and the RCs, and generated many inactive RCs, which significantly reduced the electron transfer efficiency of PSII. Severe hypoxia (2.65 mg L^-1 oxygen content) caused chlorophyll degradation. The study demonstrated that hypoxia damages PSII of E. acoroides and inhibits PSII recovery in the dark. We suggested that hypoxia together with other environment stressors would be the key reason for the decline of E. acoroides meadows.

Electronic Structure of Tyrosyl D Radical of Photosystem II, as Revealed by 2D-Hyperfine Sublevel Correlation Spectroscopy

The biological water oxidation takes place in Photosystem II (PSII), a multi-subunit protein located in thylakoid membranes of higher plant chloroplasts and cyanobacteria. The catalytic site of PSII is a Mn4Ca cluster and is known as the oxygen evolving complex (OEC) of PSII. Two tyrosine residues D1-Tyr161 (YZ) and D2-Tyr160 (YD) are symmetrically placed in the two core subunits D1 and D2 and participate in proton coupled electron transfer reactions. YZ of PSII is near the OEC and mediates electron coupled proton transfer from Mn4Ca to the photooxidizable chlorophyll species P680+. YD does not directly interact with OEC, but is crucial for modulating the various S oxidation states of the OEC. In PSII from higher plants the environment of YD• radical has been extensively characterized only in spinach (Spinacia oleracea) Mn-depleted non functional PSII membranes. Here, we present a 2D-HYSCORE investigation in functional PSII of spinach to determine the electronic structure of YD• radical. The hyperfine couplings of the protons that interact with the YD• radical are determined and the relevant assignment is provided. A discussion on the similarities and differences between the present results and the results from studies performed in non functional PSII membranes from higher plants and PSII preparations from other organisms is given.

The origin of unidirectional charge separation in photosynthetic reaction centers: nonadiabatic quantum dynamics of exciton and charge in pigment-protein complexes

Exciton charge separation in photosynthetic reaction centers from purple bacteria (PbRC) and photosystem II (PSII) occurs exclusively along one of the two pseudo-symmetric branches (active branch) of pigment-protein complexes. The microscopic origin of unidirectional charge separation in photosynthesis remains controversial. Here we elucidate the essential factors leading to unidirectional charge separation in PbRC and PSII, using nonadiabatic quantum dynamics calculations in conjunction with time-dependent density functional theory (TDDFT) with the quantum mechanics/molecular mechanics/polarizable continuum model (QM/MM/PCM) method. This approach accounts for energetics, electronic coupling, and vibronic coupling of the pigment excited states under electrostatic interactions and polarization of whole protein environments. The calculated time constants of charge separation along the active branches of PbRC and PSII are similar to those observed in time-resolved spectroscopic experiments. In PbRC, Tyr-M210 near the accessary bacteriochlorophyll reduces the energy of the intermediate state and drastically accelerates charge separation overcoming the electron-hole interaction. Remarkably, even though both the active and inactive branches in PSII can accept excitons from light-harvesting complexes, charge separation in the inactive branch is prevented by a weak electronic coupling due to symmetry-breaking of the chlorophyll configurations. The exciton in the inactive branch in PSII can be transferred to the active branch via direct and indirect pathways. Subsequently, the ultrafast electron transfer to pheophytin in the active branch prevents exciton back transfer to the inactive branch, thereby achieving unidirectional charge separation.

Rigidly hydrogen-bonded water molecules facilitate proton transfer in photosystem II

In the water-splitting enzyme photosystem II (PSII), the proton is released from the catalytic site and transferred to the protein bulk surface via the proton-relay mechanism. Proton transfer occurs in a proton-conducting channel consisting of a series of water molecules connected by hydrogen-bonded (H-bonded) chains. The water-transport protein aquaporin (AQP) also contains a water chain with structure similar to that of the PSII proton channel, although the water chain does not transport protons. We compared the PSII proton channel with the AQP water channel from the following standpoints: (1) the energetics of proton transfer based on crystal structures obtained from quantum mechanical/molecular mechanical calculations, and (2) fluctuations in water molecules obtained from molecular dynamics simulations. The results showed that residues facing the channel and acting as H-bonded partners of water molecules predominantly determined the proton-transfer ability. In PSII, the water chain is surrounded by H-bond acceptors (e.g., carbonyl groups), and the water chain transports protons where the water molecules are rigidly fixed. In AQP, the water chain is surrounded by hydrophobic sidechains or H-bond donors (e.g., NH2 groups), and it does not transport protons where the water molecules are flexible and fluctuating.

Mechanism of absorption wavelength shifts in anion channelrhodopsin-1 mutants

Using a quantum mechanical/molecular mechanical approach, we show the mechanisms of how the protein environment of Guillardia theta anion channelrhodopsin-1 (GtACR1) can shift the absorption wavelength. The calculated absorption wavelengths for GtACR1 mutants, M105A, C133A, and C237A are in agreement with experimentally measured wavelengths. Among 192 mutant structures investigated, mutations at Thr101, Cys133, Pro208, and Cys237 are likely to increase the absorption wavelength. In particular, T101A GtACR1 was expressed in HEK293T cells. The measured absorption wavelength is 10 nm higher than that of wild type, consistent with the calculated wavelength. (i) Removal of a polar residue from the Schiff base moiety, (ii) addition of a polar or acidic residue to the β-ionone ring moiety, and (iii) addition of a bulky residue to increase the planarity of the β-ionone and Schiff base moieties are the basis of increasing the absorption wavelength.

Protein Matrix Control of Reaction Center Excitation in Photosystem II

Photosystem II (PSII) is a multisubunit pigment-protein complex that uses light-induced charge separation to power oxygenic photosynthesis. Its reaction center chromophores, where the charge transfer cascade is initiated, are arranged symmetrically along the D1 and D2 core polypeptides and comprise four chlorophyll (PD1, PD2, ChlD1, ChlD2) and two pheophytin molecules (PheoD1 and PheoD2). Evolution favored productive electron transfer only via the D1 branch, with the precise nature of primary excitation and the factors that control asymmetric charge transfer remaining under investigation. Here we present a detailed atomistic description for both. We combine large-scale simulations of membrane-embedded PSII with high-level quantum-mechanics/molecular-mechanics (QM/MM) calculations of individual and coupled reaction center chromophores to describe reaction center excited states. We employ both range-separated time-dependent density functional theory and the recently developed domain based local pair natural orbital (DLPNO) implementation of the similarity transformed equation of motion coupled cluster theory with single and double excitations (STEOM-CCSD), the first coupled cluster QM/MM calculations of the reaction center. We find that the protein matrix is exclusively responsible for both transverse (chlorophylls versus pheophytins) and lateral (D1 versus D2 branch) excitation asymmetry, making ChlD1 the chromophore with the lowest site energy. Multipigment calculations show that the protein matrix renders the ChlD1 → PheoD1 charge-transfer the lowest energy excitation globally within the reaction center, lower than any pigment-centered local excitation. Remarkably, no low-energy charge transfer states are located within the "special pair" PD1-PD2, which is therefore excluded as the site of initial charge separation in PSII. Finally, molecular dynamics simulations suggest that modulation of the electrostatic environment due to protein conformational flexibility enables direct excitation of low-lying charge transfer states by far-red light.

A TDDFT investigation of the Photosystem II reaction center: Insights into the precursors to charge separation

Examining the excited states of the Photosystem II reaction center furthers our understanding of available charge separation pathways that lead to successful photosynthesis. Our results comprise the largest complete model of the Photosystem II reaction center to be described using time-dependent density functional theory reported in the literature to date. We reveal the molecular orbitals contributing to the excited states that are precursors to charge separation. We demonstrate that our model can successfully predict the action of specific mutations, a valuable tool for the agricultural industry. These models may also be beneficial in informing the design of artificial photosynthetic complexes as well as enhanced bioengineered photosystems.

Photosystem II (PS II) captures solar energy and directs charge separation (CS) across the thylakoid membrane during photosynthesis. The highly oxidizing, charge-separated state generated within its reaction center (RC) drives water oxidation. Spectroscopic studies on PS II RCs are difficult to interpret due to large spectral congestion, necessitating modeling to elucidate key spectral features. Herein, we present results from time-dependent density functional theory (TDDFT) calculations on the largest PS II RC model reported to date. This model explicitly includes six RC chromophores and both the chlorin phytol chains and the amino acid residues <6 Å from the pigments’ porphyrin ring centers. Comparing our wild-type model results with calculations on mutant D1-His-198-Ala and D2-His-197-Ala RCs, our simulated absorption-difference spectra reproduce experimentally observed shifts in known chlorophyll absorption bands, demonstrating the predictive capabilities of this model. We find that inclusion of both nearby residues and phytol chains is necessary to reproduce this behavior. Our calculations provide a unique opportunity to observe the molecular orbitals that contribute to the excited states that are precursors to CS. Strikingly, we observe two high oscillator strength, low-lying states, in which molecular orbitals are delocalized over Chl_D1 and Phe_D1 as well as one weaker oscillator strength state with molecular orbitals delocalized over the P chlorophylls. Both these configurations are a match for previously identified exciton–charge transfer states (Chl_D1⁺Phe_D1⁻)* and (P_D2⁺P_D1⁻)*. Our results demonstrate the power of TDDFT as a tool, for studies of natural photosynthesis, or indeed future studies of artificial photosynthetic complexes.

Acquirement of water-splitting ability and alteration of the charge-separation mechanism in photosynthetic reaction centers

In photosynthetic reaction centers from purple bacteria (PbRC) and the water-oxidizing enzyme, photosystem II (PSII), charge separation occurs along one of the two symmetrical electron-transfer branches. Here we report the microscopic origin of the unidirectional charge separation, fully considering electron-hole interaction, electronic coupling of the pigments, and electrostatic interaction with the polarizable entire protein environments. The electronic coupling between the pair of bacteriochlorophylls is large in PbRC, forming a delocalized excited state with the lowest excitation energy (i.e., the special pair). The charge-separated state in the active branch is stabilized by uncharged polar residues in the transmembrane region and charged residues on the cytochrome c ₂ binding surface. In contrast, the accessory chlorophyll in the D1 protein (Chl_D1) has the lowest excitation energy in PSII. The charge-separated state involves Chl_D1 ^•+ and is stabilized predominantly by charged residues near the Mn₄CaO₅ cluster and the proceeding proton-transfer pathway. It seems likely that the acquirement of water-splitting ability makes Chl_D1 the initial electron donor in PSII.

Redox Potential of the Oxygen-Evolving Complex in the Electron Transfer Cascade of Photosystem II

In photosystem II (PSII), water oxidation occurs in the Mn₄CaO₅ cluster with the release of electrons via the redox-active tyrosine (TyrZ) to the reaction-center chlorophylls (P_D1/P_D2). Using a quantum mechanical/molecular mechanical approach, we report the redox potentials (E_m) of these cofactors in the PSII protein environment. The E_m values suggest that the Mn₄CaO₅ cluster, TyrZ, and P_D1/P_D2 form a downhill electron transfer pathway. E_m for the first oxidation step, E_m(S₀/S₁), is uniquely low (730 mV) and is ∼100 mV lower than that for the second oxidation step, E_m(S₁/S₂) (830 mV) only when the O4 site of the Mn₄CaO₅ cluster is protonated in S₀. The O4-water chain, which directly forms a low-barrier H-bond with the Mn₄CaO₅ cluster and mediates proton-coupled electron transfer in the S₀ to S₁ transition, explains why the second lowest oxidation state, S₁, is the most stable and S₀ is converted to S₁ even in the dark.

Redox potentials along the redox-active low-barrier H-bonds in electron transfer pathways

Local proton transfer along redox-active low-barrier H-bonds can alter the driving force or electronic coupling for electron transfer, as the redox potential values depend on the H⁺ position in low-barrier H-bonds.

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.

Reversible Aggregation of Chlorophyll Derivative Induced by Phase Transition of Lipid

Controlling the supramolecular organization of pigment molecules will provide innovative materials that exhibit variable optical properties. In nature, photosynthetic systems employ chlorophyllous supramolecules in which each pigment molecule is suitably organized in proteins, and their properties are adequately optimized by changing the structures of the surrounding amino acid residues. Here, we report a strategy for varying the aggregation behavior of a chlorophyll derivative by using a phase-transition phenomenon of lipid bilayers. Methyl pyropheophorbide a (MPP) was employed as a chlorophyllous pigment in our artificial system, and synthetic phosphatidylcholines with saturated acyl chain(s) were also used. The MPP molecules successfully accumulated within the lipid bilayer of liposomes without changing the vesicular structure. When the lipid bilayer was in a gel form (under the phase-transition temperature, T_m), the embedded MPP aggregated to yield a dimeric form showing red-shifted absorption bands and circular dichroism signals. When the solutions of MPP-containing liposomes were heated to higher temperatures than their T_m, MPP disaggregated to monomeric form as the absorption spectrum changed into its original fashion in dichloromethane. The reversible thermochromic (dis)aggregation of the MPP molecules had good cyclability. Additional careful examination of the phase transition in the MPP-lipid co-assemblies clarified that the critical temperatures of the MPP (dis)aggregation were in good agreement with the phase-transition temperatures of the pigment-containing bilayers. The reversible MPP aggregation in the lipid bilayers occurred in a wide range of temperatures (around 10-55 °C) by changing the length of the diacyl side chains of phospholipids. The reversible thermochromism of the chlorophyllous system was established by varying the nature of the surrounding lipid bilayer. This study can provide a useful strategy for making variable tetrapyrrolic aggregate systems induced by mild extrinsic stimuli.

Microsolvation of the Redox-Active Tyrosine-D in Photosystem II: Correlation of Energetics with EPR Spectroscopy and Oxidation-Induced Proton Transfer

Photosystem II (PSII) of oxygenic photosynthesis captures sunlight to drive the catalytic oxidation of water and the reduction of plastoquinone. Among the several redox-active cofactors that participate in intricate electron transfer pathways there are two tyrosine residues, Y_Z and Y_D. They are situated in symmetry-related electron transfer branches but have different environments and play distinct roles. Y_Z is the immediate oxidant of the oxygen-evolving Mn₄CaO₅ cluster, whereas Y_D serves regulatory and protective functions. The protonation states and hydrogen-bond network in the environment of Y_D remain debated, while the role of microsolvation in stabilizing different redox states of Y_D and facilitating oxidation or mediating deprotonation, as well the fate of the phenolic proton, is unclear. Here we present detailed structural models of Y_D and its environment using large-scale quantum mechanical models and all-atom molecular dynamics of a complete PSII monomer. The energetics of water distribution within a hydrophobic cavity adjacent to Y_D are shown to correlate directly with electron paramagnetic resonance (EPR) parameters such as the tyrosyl g-tensor, allowing us to map the correspondence between specific structural models and available experimental observations. EPR spectra obtained under different conditions are explained with respect to the mode of interaction of the proximal water with the tyrosyl radical and the position of the phenolic proton within the cavity. Our results revise previous models of the energetics and build a detailed view of the role of confined water in the oxidation and deprotonation of Y_D. Finally, the model of microsolvation developed in the present work rationalizes in a straightforward way the biphasic oxidation kinetics of Y_D, offering new structural insights regarding the function of the radical in biological photosynthesis.

Energetic insights into two electron transfer pathways in light-driven energy-converting enzymes

We report E_m values of (bacterio-)chlorophylls for one-electron reduction in both electron-transfer branches of PbRC, PSI, and PSII.

We report redox potentials (E_m) for one-electron reduction for all chlorophylls in the two electron-transfer branches of water-oxidizing enzyme photosystem II (PSII), photosystem I (PSI), and purple bacterial photosynthetic reaction centers (PbRC). In PSI, E_m values for the accessory chlorophylls were similar in both electron-transfer branches. In PbRC, the corresponding E_m value was 170 mV less negative in the active L-branch (B_L) than in the inactive M-branch (B_M), favoring B_L˙^– formation. This contrasted with the corresponding chlorophylls, Chl_D1 and Chl_D2, in PSII, where E_m(Chl_D1) was 120 mV more negative than E_m(Chl_D2), implying that to rationalize electron transfer in the D1-branch, Chl_D1 would need to serve as the primary electron donor. Residues that contributed to E_m(Chl_D1) < E_m(Chl_D2) simultaneously played a key role in (i) releasing protons from the substrate water molecules and (ii) contributing to the larger cationic population on the chlorophyll closest to the Mn₄CaO₅ cluster (P_D1), favoring electron transfer from water molecules. These features seem to be the nature of PSII, which needs to possess the proton-exit pathway to use a protonated electron source—water molecules.

Genetically introduced hydrogen bond interactions reveal an asymmetric charge distribution on the radical cation of the special-pair chlorophyll P680

The special-pair chlorophyll (Chl) P680 in photosystem II has an extremely high redox potential (E_m ) to enable water oxidation in photosynthesis. Significant positive-charge localization on one of the Chl constituents, P_D1 or P_D2, in P680⁺ has been proposed to contribute to this high E_m To identify the Chl molecule on which the charge is mainly localized, we genetically introduced a hydrogen bond to the 13¹-keto C=O group of P_D1 and P_D2 by changing the nearby D1-Val-157 and D2-Val-156 residues to His, respectively. Successful hydrogen bond formation at P_D1 and P_D2 in the obtained D1-V157H and D2-V156H mutants, respectively, was monitored by detecting 13¹-keto C=O vibrations in Fourier transfer infrared (FTIR) difference spectra upon oxidation of P680 and the symmetrically located redox-active tyrosines Y_Z and Y_D, and they were simulated by quantum-chemical calculations. Analysis of the P680⁺/P680 FTIR difference spectra of D1-V157H and D2-V156H showed that upon P680⁺ formation, the 13¹-keto C=O frequency upshifts by a much larger extent in P_D1 (23 cm^-1) than in P_D2 (<9 cm^-1). In addition, thermoluminescence measurements revealed that the D1-V157H mutation increased the E_m of P680 to a larger extent than did the D2-V156H mutation. These results, together with the previous results for the mutants of the His ligands of P_D1 and P_D2, lead to a definite conclusion that a charge is mainly localized to P_D1 in P680<sup/>.

Regulated Electron Tunneling of Photoinduced Primary Charge-Separated State in the Photosystem II Reaction Center

In initial events of the photosynthesis by higher plants, the photosystem II (PSII) generates photoinduced primary charge-separated (CS) state composed of reduced pheophytin (Pheo_D1^-•) and oxidized special pair (P^+•) in chlorophyll a (Chla) P_D1/P_D2 in the D1/D2 heterodimer, ultimately leading to the water oxidation at the oxygen-evolving Mn₄CaO₅ cluster by P^+•. To understand the mechanism of the efficient generation of initially localized CS state (P_D1^+• Pheo_D1^-•), we have characterized cofactor geometries and electronic coupling of the photoinduced primary CS state in quinone prereduced membrane of PSII from spinach using the time-resolved electron paramagnetic resonance method. It has been revealed that the electronic coupling between the charges is significantly weak in the CS state separated by 1.5 nm, showing the importance of regulated cofactor-cofactor electronic interaction between a vinyl substituent in Pheo_D1 and an accessory chlorophyll to inhibit the energy-wasting charge recombination after the primary electron-transfer processes.

Conformational control of cofactors in nature - the influence of protein-induced macrocycle distortion on the biological function of tetrapyrroles

Tetrapyrrole-containing proteins are one of the most fundamental classes of enzymes in nature and it remains an open question to give a chemical rationale for the multitude of biological reactions that can be catalyzed by these pigment-protein complexes. There are many fundamental processes where the same (i.e., chemically identical) porphyrin cofactor is involved in chemically quite distinct reactions. For example, heme is the active cofactor for oxygen transport and storage (hemoglobin, myoglobin) and for the incorporation of molecular oxygen in organic substrates (cytochrome P450). It is involved in the terminal oxidation (cytochrome c oxidase) and the metabolism of H2O2 (catalases and peroxidases) and catalyzes various electron transfer reactions in cytochromes. Likewise, in photosynthesis the same chlorophyll cofactor may function as a reaction center pigment (charge separation) or as an accessory pigment (exciton transfer) in light harvesting complexes (e.g., chlorophyll a). Whilst differences in the apoprotein sequences alone cannot explain the often drastic differences in physicochemical properties encountered for the same cofactor in diverse protein complexes, a critical factor for all biological functions must be the close structural interplay between bound cofactors and the respective apoprotein in addition to factors such as hydrogen bonding or electronic effects. Here, we explore how nature can use the same chemical molecule as a cofactor for chemically distinct reactions using the concept of conformational flexibility of tetrapyrroles. The multifaceted roles of tetrapyrroles are discussed in the context of the current knowledge on distorted porphyrins. Contemporary analytical methods now allow a more quantitative look at cofactors in protein complexes and the development of the field is illustrated by case studies on hemeproteins and photosynthetic complexes. Specific tetrapyrrole conformations are now used to prepare bioengineered designer proteins with specific catalytic or photochemical properties.

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.

Energetics of the Proton Transfer Pathway for Tyrosine D in Photosystem II

The proton transfer pathway for redox active tyrosine D (TyrD) in photosystem II is a hydrogen-bond network that involves D2-Arg180 and a series of water molecules. Using quantum mechanical/molecular mechanical calculations, the detailed properties of the energetics and structural geometries were investigated. The potential-energy profile of all hydrogen bonds along the proton transfer pathway indicates that the overall proton transfer from TyrD is energetically downhill. D2-Arg180 plays a key role in the proton transfer pathway, providing a driving force for proton transfer, maintaining the hydrogen-bond network structure, stabilising P680•+, and thus deprotonating TyrD-OH to TyrD-O•. A hydrophobic environment near TyrD enhances the electrostatic interactions between TyrD and redox active groups, e.g. P680 and the catalytic Mn4CaO5 cluster: the redox states of those groups are linked with the protonation state of TyrD, i.e. release of the proton from TyrD. Thus, the proton transfer pathway from TyrD may ultimately contribute to the conversion of S0 into S1 in the dark in order to stabilise the Mn4CaO5 cluster when the photocycle is interrupted in S0.

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.

Design of dinuclear manganese cofactors for bacterial reaction centers

A compelling target for the design of electron transfer proteins with novel cofactors is to create a model for the oxygen-evolving complex, a Mn4Ca cluster, of photosystem II. A mononuclear Mn cofactor can be added to the bacterial reaction center, but the addition of multiple metal centers is constrained by the native protein architecture. Alternatively, metal centers can be incorporated into artificial proteins. Designs for the addition of dinuclear metal centers to four-helix bundles resulted in three artificial proteins with ligands for one, two, or three dinuclear metal centers able to bind Mn. The three-dimensional structure determined by X-ray crystallography of one of the Mn-proteins confirmed the design features and revealed details concerning coordination of the Mn center. Electron transfer between these artificial Mn-proteins and bacterial reaction centers was investigated using optical spectroscopy. After formation of a light-induced, charge-separated state, the experiments showed that the Mn-proteins can donate an electron to the oxidized bacteriochlorophyll dimer of modified reaction centers, with the Mn-proteins having additional metal centers being more effective at this electron transfer reaction. Modeling of the structure of the Mn-protein docked to the reaction center showed that the artificial protein likely binds on the periplasmic surface similarly to cytochrome c2, the natural secondary donor. Combining reaction centers with exogenous artificial proteins provides the opportunity to create ligands and investigate the influence of inhomogeneous protein environments on multinuclear redox-active metal centers. This article is part of a Special Issue entitled Biodesign for Bioenergetics--the design and engineering of electronic transfer cofactors, proteins and protein networks, edited by Ronald L. Koder and J.L. Ross Anderson.

Does the singlet minus triplet spectrum with major photobleaching band near 680-682 nm represent an intact reaction center of Photosystem II?

We use both frequency- and time-domain low-temperature (5-20 K) spectroscopies to further elucidate the shape and spectral position of singlet minus triplet (triplet-bottleneck) spectra in the reaction centers (RCs) of Photosystem II (PSII) isolated from wild-type Chlamydomonas reinhardtii and spinach. It is shown that the shape of the nonresonant transient hole-burned spectrum in destabilized RCs from C. reinhardtii is very similar to that typically observed for spinach. This suggests that the previously observed difference in transient spectra between RCs from C. reinhardtii and spinach is not due to the sample origin but most likely due to a partial destabilization of the D1 and D2 polypeptides. This supports our previous assignments that destabilized RCs (referred to as RC680) (Acharya, K. et al. J. Phys. Chem. B 2012, 116, 4860-4870), with a major photobleaching band near 680-682 nm and the absence of a photobleaching band near 673 nm, do not represent the intact RC residing within the PSII core complex. Time-resolved absorption difference spectra obtained for partially destabilized RCs of C. reinhardtii and for typical spinach RCs support the above conclusions. The absence of clear photobleaching bands near 673 and 684 nm (where the PD1 chlorophyll and the active pheophytin (PheoD1) contribute, respectively) in picosecond transient absorption spectra in both RCs studied in this work indicates that the cation can move from the primary electron donor (ChlD1) to PD1 (i.e., PD1ChlD1(+)PheoD1(-) → PD1(+)ChlD1PheoD1(-)). Therefore, we suggest that ChlD1 is the major electron donor in usually studied destabilized RCs (with a major photobleaching near 680-682 nm), although the PD1 path (where PD1 serves as the primary electron donor) is likely present in intact RCs, as discussed in Acharya, K. et al. J. Phys. Chem. B 2012, 116, 4860-4870.

Dynamic protein conformations preferentially drive energy transfer along the active chain of the photosystem II reaction centre

One longstanding puzzle concerning photosystem II, a core component of photosynthesis, is that only one of the two symmetric branches in its reaction centre is active in electron transfer. To investigate the effect of the photosystem II environment on the preferential selection of the energy transfer pathway (a prerequisite for electron transfer), we have constructed an exciton model via extensive molecular dynamics simulations and quantum mechanics/molecular mechanics calculations based on a recent X-ray structure. Our results suggest that it is essential to take into account an ensemble of protein conformations to accurately compute the site energies. We identify the cofactor CLA606 of active chain as the most probable site for the energy excitation. We further pinpoint a number of charged protein residues that collectively lower the CLA606 site energy. Our work provides insights into the understanding of molecular mechanisms of the core machinery of the green-plant photosynthesis.

Cofactor-mediated energy and electron transfer in photosystem II occurs preferentially through one branch of the reaction centre, despite there being a symmetric path available. Here, the authors use computational methods to determine the influence of protein conformation on this selectivity.

Mechanism of tyrosine D oxidation in Photosystem II

Using quantum mechanics/molecular mechanics calculations and the 1.9-Å crystal structure of Photosystem II [Umena Y, Kawakami K, Shen J-R, Kamiya N (2011) Nature 473(7345):55-60], we investigated the H-bonding environment of the redox-active tyrosine D (TyrD) and obtained insights that help explain its slow redox kinetics and the stability of TyrD(•). The water molecule distal to TyrD, located ~4 Å away from the phenolic O of TyrD, corresponds to the presence of the tyrosyl radical state. The water molecule proximal to TyrD, in H-bonding distance to the phenolic O of TyrD, corresponds to the presence of the unoxidized tyrosine. The H(+) released on oxidation of TyrD is transferred to the proximal water, which shifts to the distal position, triggering a concerted proton transfer pathway involving D2-Arg180 and a series of waters, through which the proton reaches the aqueous phase at D2-His61. The water movement linked to the ejection of the proton from the hydrophobic environment near TyrD makes oxidation slow and quasiirreversible, explaining the great stability of the TyrD(•). A symmetry-related proton pathway associated with tyrosine Z is pointed out, and this is associated with one of the Cl(-) sites. This may represent a proton pathway functional in the water oxidation cycle.

Influence of protein interactions on oxidation/reduction midpoint potentials of cofactors in natural and de novo metalloproteins

As discussed throughout this special issue, oxidation and reduction reactions play critical roles in the function of many organisms. In photosynthetic organisms, the conversion of light energy drives oxidation and reduction reactions through the transfer of electrons and protons in order to create energy-rich compounds. These reactions occur in proteins such as cytochrome c, a heme-containing water-soluble protein, the bacteriochlorophyll-containing reaction center, and photosystem II where water is oxidized at the manganese cluster. A critical measure describing the ability of cofactors in proteins to participate in such reactions is the oxidation/reduction midpoint potential. In this review, the basic concepts of oxidation/reduction reactions are reviewed with a summary of the experimental approaches used to measure the midpoint potential of metal cofactors. For cofactors in proteins, the midpoint potential not only depends upon the specific chemical characteristics of cofactors but also upon interactions with the surrounding protein, such as the nature of the coordinating ligands and protein environment. These interactions can be tailored to optimize an oxidation/reduction reaction carried out by the protein. As examples, the midpoint potentials of hemes in cytochromes, bacteriochlorophylls in reaction centers, and the manganese cluster of photosystem II are discussed with an emphasis on the influence that protein interactions have on these potentials. This article is part of a Special Issue entitled: Metals in Bioenergetics and Biomimetics Systems.

Mechanism of proton-coupled quinone reduction in Photosystem II

Photosystem II uses light to drive water oxidation and plastoquinone (PQ) reduction. PQ reduction involves two PQ cofactors, Q(A) and Q(B), working in series. Q(A) is a one-electron carrier, whereas Q(B) undergoes sequential reduction and protonation to form Q(B)H(2). Q(B)H(2) exchanges with PQ from the pool in the membrane. Based on the atomic coordinates of the Photosystem II crystal structure, we analyzed the proton transfer (PT) energetics adopting a quantum mechanical/molecular mechanical approach. The potential-energy profile suggests that the initial PT to Q(B)(•-) occurs from the protonated, D1-His252 to Q(B)(•)(-) via D1-Ser264. The second PT is likely to occur from D1-His215 to Q(B)H(-) via an H-bond with an energy profile with a single well, resulting in the formation of Q(B)H(2) and the D1-His215 anion. The pathway for reprotonation of D1-His215(-) may involve bicarbonate, D1-Tyr246 and water in the Q(B) site. Formate ligation to Fe(2+) did not significantly affect the protonation of reduced Q(B), suggesting that formate inhibits Q(B)H(2) release rather than its formation. The presence of carbonate rather than bicarbonate seems unlikely because the calculations showed that this greatly perturbed the potential of the nonheme iron, stabilizing the Fe(3+) state in the presence of Q(B)(•-), a situation not encountered experimentally. H-bonding from D1-Tyr246 and D2-Tyr244 to the bicarbonate ligand of the nonheme iron contributes to the stability of the semiquinones. A detailed mechanistic model for Q(B) reduction is presented.