Electrostatic Asymmetry in the Reaction Center of Photosystem II

Femtosecond optical studies of the primary charge separation reactions in far-red photosystem II from Synechococcus sp. PCC 7335

Primary processes of light energy conversion by Photosystem II (PSII) were studied using femtosecond broadband pump-probe absorption difference spectroscopy. Transient absorption changes of core complexes isolated from the cyanobacterium Synechococcus sp. PCC 7335 grown under far-red light (FRL-PSII) were compared with the canonical Chl a containing spinach PSII core complexes upon excitation into the red edge of the Qy band. Absorption changes of FRL-PSII were monitored at 278 K in the 400-800 nm spectral range on a timescale of 0.1-500 ps upon selective excitation at 740 nm of four chlorophyll (Chl) f molecules in the light harvesting antenna, or of one Chl d molecule at the ChlD1 position in the reaction center (RC) upon pumping at 710 nm. Numerical analysis of absorption changes and assessment of the energy levels of the presumed ion-radical states made it possible to identify PD1+ChlD1- as the predominant primary charge-separated radical pair, the formation of which upon selective excitation of Chl d has an apparent time of ∼1.6 ps. Electron transfer to the secondary acceptor pheophytin PheoD1 has an apparent time of ∼7 ps with a variety of excitation wavelengths. The energy redistribution between Chl a and Chl f in the antenna occurs within 1 ps, whereas the energy migration from Chl f to the RC occurs mostly with lifetimes of 60 and 400 ps. Potentiometric analysis suggests that in canonical PSII, PD1+ChlD1- can be partially formed from the excited (PD1ChlD1)* state.

PigmentHunter: A point-and-click application for automated chlorophyll-protein simulations

Chlorophyll proteins (CPs) are the workhorses of biological photosynthesis, working together to absorb solar energy, transfer it to chemically active reaction centers, and control the charge-separation process that drives its storage as chemical energy. Yet predicting CP optical and electronic properties remains a serious challenge, driven by the computational difficulty of treating large, electronically coupled molecular pigments embedded in a dynamically structured protein environment. To address this challenge, we introduce here an analysis tool called PigmentHunter, which automates the process of preparing CP structures for molecular dynamics (MD), running short MD simulations on the nanoHUB.org science gateway, and then using electrostatic and steric analysis routines to predict optical absorption, fluorescence, and circular dichroism spectra within a Frenkel exciton model. Inter-pigment couplings are evaluated using point-dipole or transition-charge coupling models, while site energies can be estimated using both electrostatic and ring-deformation approaches. The package is built in a Jupyter Notebook environment, with a point-and-click interface that can be used either to manually prepare individual structures or to batch-process many structures at once. We illustrate PigmentHunter's capabilities with example simulations on spectral line shapes in the light harvesting 2 complex, site energies in the Fenna-Matthews-Olson protein, and ring deformation in photosystems I and II.

Tuning of the ChlD1 and ChlD2 properties in photosystem II by site-directed mutagenesis of neighbouring amino acids

Photosystem II is the water/plastoquinone photo-oxidoreductase of photosynthesis. The photochemistry and catalysis occur in a quasi-symmetrical heterodimer, D1D2, that evolved from a homodimeric ancestor. Here, we studied site-directed mutants in PSII from the thermophilic cyanobacterium Thermosynechoccocus elongatus, focusing on the primary electron donor chlorophyll a in D1, ChlD1, and on its symmetrical counterpart in D2, ChlD2, which does not play a direct photochemical role. The main conserved amino acid specific to ChlD1 is D1/T179, which H-bonds the water ligand to its Mg2+, while its counterpart near ChlD2 is the non-H-bonding D2/I178. The symmetrical-swapped mutants, D1/T179I and D2/I178T, and a second ChlD2 mutant, D2/I178H, were studied. The D1 mutations affected the 686 nm absorption attributed to ChlD1, while the D2 mutations affected a 663 nm feature, tentatively attributed to ChlD2. The mutations had little effect on enzyme activity and forward electron transfer, reflecting the robustness of the overall enzyme function. In contrast, the mutations significantly affected photodamage and protective mechanisms, reflecting the importance of redox tuning in these processes. In D1/T179I, the radical pair recombination triplet on ChlD1 was shared onto a pheophytin, presumably PheD1 and the detection of 3PheD1 supports the proposed mechanism for the anomalously short lifetime of 3ChlD1; e.g. electron transfer quenching by QA- of 3PheD1 after triplet transfer from 3ChlD1. In D2/I178T, a charge separation could occur between ChlD2 and PheD2, a reaction that is thought to occur in ancestral precursors of PSII. These mutants help understand the evolution of asymmetry in PSII.

Reinvestigation on primary processes of PSII-dimer from Thermosynechococcus vulcanus by femtosecond pump-probe spectroscopy

Cyanobacterial photosynthetic apparatus efficiently capture sunlight, and the energy is subsequently transferred to photosystem I (PSI) and II (PSII), to produce electrochemical potentials. PSII is a unique membrane protein complex that photo-catalyzes oxidation of water and majorly contains photosynthetic pigments of chlorophyll a and carotenoids. In the present study, the ultrafast energy transfer and charge separation dynamics of PSII from a thermophilic cyanobacterium Thermosynechococcus vulcanus were reinvestigated by femtosecond pump-probe spectroscopic measurements under low temperature and weak intensity excitation condition. The results imply the two possible models of the energy transfers and subsequent charge separation in PSII. One is the previously suggested "transfer-to-trapped limit" model. Another model suggests that the energy transfers from core CP43 and CP47 antennas to the primary electron donor ChlD1 with time-constants of 0.71 ps and 3.28 ps at 140 K (0.17 and 1.33 ps at 296 K), respectively and that the pheophytin anion (PheoD1-) is generated with the time-constant of 43.0 ps at 140 K (14.8 ps at 296 K) upon excitation into the Qy band of chlorophyll a at 670 nm. The secondary electron transfer to quinone QA: PheoD1-QA → PheoD1QA- is observed with the time-constant of 650 ps only at 296 K. On the other hand, an inefficient β-carotene → chlorophyll a energy transfer (33%) occurred after excitation to the S2 state of β-carotene at 500 nm. Instead, the carotenoid triplet state appeared in an ultrafast timescale after excitation at 500 nm.

Short-Range Effects in the Special Pair of Photosystem II Reaction Centers: The Nonconservative Nature of Circular Dichroism

Photosystem II reaction centers extract electrons from water, providing the basis of oxygenic life on earth. Among the light-sensitive pigments of the reaction center, a central chlorophyll a dimer, known as the special pair, so far has escaped a complete theoretical characterization of its excited state properties. The close proximity of the special pair pigments gives rise to short-range effects that comprise a coupling between local and charge transfer (CT) excited states as well as other intermolecular quantum effects. Using a multiscale simulation and a diabatization technique, we show that the coupling to CT states is responsible for 45% of the excitonic coupling in the special pair. The other short-range effects cause a nonconservative nature of the circular dichroism spectrum of the reaction center by effectively rotating the electric transition dipole moments of the special pair pigments inverting and strongly enhancing their intrinsic rotational strength.

Reaction Center Excitation in Photosystem II: From Multiscale Modeling to Functional Principles

Oxygenic photosynthesis is the fundamental energy-converting process that utilizes sunlight to generate molecular oxygen and the organic compounds that sustain life. Protein-pigment complexes harvest light and transfer excitation energy to specialized pigment assemblies, reaction centers (RC), where electron transfer cascades are initiated. A molecular-level understanding of the primary events is indispensable for elucidating the principles of natural photosynthesis and enabling development of bioinspired technologies. The primary enzyme in oxygenic photosynthesis is Photosystem II (PSII), a membrane-embedded multisubunit complex, that catalyzes the light-driven oxidation of water. The RC of PSII consists of four chlorophyll a and two pheophytin a pigments symmetrically arranged along two core polypeptides; only one branch participates in electron transfer. Despite decades of research, fundamental questions remain, including the origin of this functional asymmetry, the nature of primary charge-transfer states and the identity of the initial electron donor, the origin of the capability of PSII to enact charge separation with far-red photons, i.e., beyond the "red limit" where individual chlorophylls absorb, and the role of protein conformational dynamics in modulating charge-separation pathways.In this Account, we highlight developments in quantum-chemistry based excited-state computations for multipigment assemblies and the refinement of protocols for computing protein-induced electrochromic shifts and charge-transfer excitations calibrated with modern local correlation coupled cluster methods. We emphasize the importance of multiscale atomistic quantum-mechanics/molecular-mechanics and large-scale molecular dynamics simulations, which enabled direct and accurate modeling of primary processes in RC excitation at the quantum mechanical level.Our findings show how differential protein electrostatics enable spectral tuning of RC pigments and generate functional asymmetry in PSII. A chlorophyll pigment on the active branch (ChlD1) has the lowest site energy in PSII and is the primary electron donor. The complete absence of low-lying charge-transfer states within the central pair of chlorophylls excludes a long-held assumption about the initial charge separation. Instead, we identify two primary charge separation pathways, both with the same pheophytin acceptor (PheoD1): a fast pathway with ChlD1 as the primary electron donor (short-range charge-separation) and a slow pathway with PD1PD2 as the initial donor (long-range charge separation). The low-energy spectrum is dominated by two states with significant charge-transfer character, ChlD1δ+PheoD1δ- and PD1δ+PheoD1δ-. The conformational dynamics of PSII allows these charge-transfer states to span wide energy ranges, pushing oxygenic photosynthesis beyond the "red limit". These results provide a quantum mechanical picture of the primary events in the RC of oxygenic photosynthesis, forming a solid basis for interpreting experimental observations and for extending photosynthesis research in new directions.

New foundations for the physical mechanism of variable chlorophyll a fluorescence. Quantum efficiency versus the light-adapted state of photosystem II

Photosystem II (PSII) uses solar energy to oxidize water and delivers electrons to fix CO2. Although the structure at atomic resolution and the basic photophysical and photochemical functions of PSII are well understood, many important questions remain. The activity of PSII in vitro and in vivo is routinely monitored by recording the induction kinetics of chlorophyll a fluorescence (ChlF). According to the 'mainstream' model, the rise from the minimum level (Fo) to the maximum (Fm) of ChlF of dark-adapted PSII reflects the closure of all functionally active reaction centers, and the Fv/Fm ratio is equated with the maximum photochemical quantum yield of PSII (where Fv=Fm-Fo). However, this model has never been free of controversies. Recent experimental data from a number of studies have confirmed that the first single-turnover saturating flash (STSF), which generates the closed state (PSIIC), produces F1 Collapse

Key Words

• F v/Fm
• Chlorophyll a fluorescence induction
• QA-model
• conformational changes
• dielectric relaxation
• electric field effects
• light-adapted charge-separated state
• photochemical quantum efficiency
• photosystem II
• purple bacterial reaction center

Collapse

MESH Headings Collapse

Grants

• Hungarian Ministry of Innovation and Technology
• ANN-144012 National Research, Development and Innovation Fund
• GA ČR 23-07744S Czech Science Foundation
• ELKH KÖ-36/2021 Eötvös Loránd Research Network

Collapse

Affiliation(s)

Győző Garab Institute of Plant Biology, Biological Research Centre, Szeged, Hungary Department of Physics, Faculty of Science, University of Ostrava, Ostrava, Czech Republic
Melinda Magyar Institute of Plant Biology, Biological Research Centre, Szeged, Hungary
Gábor Sipka Institute of Plant Biology, Biological Research Centre, Szeged, Hungary
Petar H Lambrev Institute of Plant Biology, Biological Research Centre, Szeged, Hungary

Collapse

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.

Charge separation in the photosystem II reaction center resolved by multispectral two-dimensional electronic spectroscopy

The photosystem II reaction center (PSII RC) performs the primary energy conversion steps of oxygenic photosynthesis. While the PSII RC has been studied extensively, the similar time scales of energy transfer and charge separation and the severely overlapping pigment transitions in the Q_y region have led to multiple models of its charge separation mechanism and excitonic structure. Here, we combine two-dimensional electronic spectroscopy (2DES) with a continuum probe and two-dimensional electronic vibrational spectroscopy (2DEV) to study the cyt b559-D1D2 PSII RC at 77 K. This multispectral combination correlates the overlapping Q_y excitons with distinct anion and pigment-specific Q_x and mid-infrared transitions to resolve the charge separation mechanism and excitonic structure. Through extensive simultaneous analysis of the multispectral 2D data, we find that charge separation proceeds on multiple time scales from a delocalized excited state via a single pathway in which Pheo_D1 is the primary electron acceptor, while Chl_D1 and P_D1 act in concert as the primary electron donor.

Energetics and proton release in photosystem II from Thermosynechococcus elongatus with a D1 protein encoded by either the psbA2 or psbA3 gene

In the cyanobacterium Thermosynechococcus elongatus, there are three psbA genes coding for the Photosystem II (PSII) D1 subunit that interacts with most of the main cofactors involved in the electron transfers. Recently, the 3D crystal structures of both PsbA2-PSII and PsbA3-PSII have been solved [Nakajima et al., J. Biol. Chem. 298 (2022) 102668.]. It was proposed that the loss of one hydrogen bond of Phe_D1 due to the D1-Y147F exchange in PsbA2-PSII resulted in a more negative E_m of Phe_D1 in PsbA2-PSII when compared to PsbA3-PSII. In addition, the loss of two water molecules in the Cl-1 channel was attributed to the D1-P173M substitution in PsbA2-PSII. This exchange, by narrowing the Cl-1 proton channel, could be at the origin of a slowing down of the proton release. Here, we have continued the characterization of PsbA2-PSII by measuring the thermoluminescence from the S₂Q_A^-/DCMU charge recombination and by measuring proton release kinetics using time-resolved absorption changes of the dye bromocresol purple. It was found that i) the E_m of Phe_D1^-•/Phe_D1 was decreased by ~30 mV in PsbA2-PSII when compared to PsbA3-PSII and ii) the kinetics of the proton release into the bulk was significantly slowed down in PsbA2-PSII in the S₂Tyr_Z^• to S₃Tyr_Z and S₃Tyr_Z^• → (S₃Tyr_Z^•)' transitions. This slowing down was partially reversed by the PsbA2/M173P mutation and induced by the PsbA3/P173M mutation thus confirming a role of the D1-173 residue in the egress of protons trough the Cl-1 channel.

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.

Quantum mechanical analysis of excitation energy transfer couplings in photosystem II

We evaluated excitation energy transfer (EET) coupling (J) between all pairs of chlorophylls (Chls) and pheophytins (Pheos) in the protein environment of photosystem II based on the time-dependent density functional theory with a quantum mechanical/molecular mechanics approach. In the reaction center, the EET coupling between Chls PD1 and PD2 is weaker (|J(PD1/PD2)| = 79 cm-1), irrespective of a short edge-to-edge distance of 3.6 Å (Mg-to-Mg distance of 8.1 Å), than the couplings between PD1 and the accessory ChlD1 (|J(PD1/ChlD2)| = 104 cm-1) and between PD2 and ChlD2 (|J(PD2/ChlD1)| = 101 cm-1), suggesting that PD1 and PD2 are two monomeric Chls rather than a "special pair". There exist strongly coupled Chl pairs (|J| > ∼100 cm-1) in the CP47 and CP43 core antennas, which may be candidates for the red-shifted Chls observed in spectroscopic studies. In CP47 and CP43, Chls ligated to CP47-His26 and CP43-His56, which are located in the middle layer of the thylakoid membrane, play a role in the "hub" that mediates the EET from the lumenal to stromal layers. In the stromal layer, Chls ligated to CP47-His466, CP43-His441, and CP43-His444 mediate the EET from CP47 to ChlD2/PheoD2 and from CP43 to ChlD1/PheoD1 in the reaction center. Thus, the excitation energy from both CP47 and CP43 can always be utilized for the charge-separation reaction in the reaction center.

Functional Water Networks in Fully Hydrated Photosystem II

Water channels and networks within photosystem II (PSII) of oxygenic photosynthesis are critical for enzyme structure and function. They control substrate delivery to the oxygen-evolving center and mediate proton transfer at both the oxidative and reductive endpoints. Current views on PSII hydration are derived from protein crystallography, but structural information may be compromised by sample dehydration and technical limitations. Here, we simulate the physiological hydration structure of a cyanobacterial PSII model following a thorough hydration procedure and large-scale unconstrained all-atom molecular dynamics enabled by massively parallel simulations. We show that crystallographic models of PSII are moderately to severely dehydrated and that this problem is particularly acute for models derived from X-ray free electron laser (XFEL) serial femtosecond crystallography. We present a fully hydrated representation of cyanobacterial PSII and map all water channels, both static and dynamic, associated with the electron donor and acceptor sides. Among them, we describe a series of transient channels and the attendant conformational gating role of protein components. On the acceptor side, we characterize a channel system that is absent from existing crystallographic models but is likely functionally important for the reduction of the terminal electron acceptor plastoquinone Q_B. The results of the present work build a foundation for properly (re)evaluating crystallographic models and for eliciting new insights into PSII structure and function.

From antenna to reaction center: Pathways of ultrafast energy and charge transfer in photosystem II

The photosystem II core complex (PSII-CC) is a photosynthetic complex that contains antenna proteins, which collect energy from sunlight, and a reaction center, which converts the collected energy to redox potential. Understanding the interplay between the antenna proteins and the reaction center will facilitate the development of more efficient solar energy conversion technologies. Here, we study the sub-100-ps dynamics of PSII-CC with two-dimensional electronic-vibrational spectroscopy, which connects energy flows with physical space, allowing a direct mapping of energy transfer pathways. Our results reveal a complex dynamical scheme which includes a specific pathway that connects CP43 to the reaction center. Resolving this pathway experimentally provides insights into the energy conversion processes in natural photosynthesis.

The photosystem II core complex (PSII-CC) is the smallest subunit of the oxygenic photosynthetic apparatus that contains core antennas and a reaction center, which together allow for rapid energy transfer and charge separation, ultimately leading to efficient solar energy conversion. However, there is a lack of consensus on the interplay between the energy transfer and charge separation dynamics of the core complex. Here, we report the application of two-dimensional electronic-vibrational (2DEV) spectroscopy to the spinach PSII-CC at 77 K. The simultaneous temporal and spectral resolution afforded by 2DEV spectroscopy facilitates the separation and direct assignment of coexisting dynamical processes. Our results show that the dominant dynamics of the PSII-CC are distinct in different excitation energy regions. By separating the excitation regions, we are able to distinguish the intraprotein dynamics and interprotein energy transfer. Additionally, with the improved resolution, we are able to identify the key pigments involved in the pathways, allowing for a direct connection between dynamical and structural information. Specifically, we show that C505 in CP43 and the peripheral chlorophyll Chlz_D1 in the reaction center are most likely responsible for energy transfer from CP43 to the reaction center.

Structure-based model of fucoxanthin-chlorophyll protein complex: Calculations of chlorophyll electronic couplings

Diatoms are a group of marine algae that are responsible for a significant part of global oxygen production. Adapted to life in an aqueous environment dominated by the blue-green light, their major light-harvesting antennae-fucoxanthin-chlorophyll protein complexes (FCPs)-exhibit different pigment compositions than of plants. Despite extensive experimental studies, until recently the theoretical description of excitation energy dynamics in these complexes was limited by the lack of high-resolution structural data. In this work, we use the recently resolved crystallographic information of the FCP complex from Phaeodactylum tricornutum diatom [Wang et al., Science 363, 6427 (2019)] and quantum chemistry-based calculations to evaluate the chlorophyll transition dipole moments, atomic transition charges from electrostatic potential, and the inter-chlorophyll couplings in this complex. The obtained structure-based excitonic couplings form the foundation for any modeling of stationary or time-resolved spectroscopic data. We also calculate the inter-pigment Förster energy transfer rates and identify two quickly equilibrating chlorophyll clusters.

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.

Searching for a Unique Exciton Model of Photosynthetic Pigment–Protein Complexes: Photosystem II Reaction Center Study by Differential Evolution

Studying the optical properties of photosynthetic pigment–protein complexes (PPCs) in the visible light range, both experimentally and theoretically, is one of the ways of gaining knowledge about the function of the photosynthetic machinery of living species. To simulate the PPC optical response, it is necessary to use semiclassical theories describing the effect of external fields–matter interaction, energy migration in molecular crystals, and electron–phonon coupling. In this paper, we report the results of photosystem II reaction center (PSIIRC) linear optical response simulations. Applying the multimode Brownian oscillator model and the theory of molecular excitons, we have demonstrated that the absorption, circular and linear dichroism, and steady-state fluorescence of PSIIRC can be accurately fitted with the help of differential evolution (DE), the multiparametric evolutionary optimization algorithm. To explore the effectiveness of DE, we used the simulated experimental data as the target functions instead of those actually measured. Only 2 of 10 DE strategies have shown the best performance of the optimization algorithm. With the best tuning parameters of DE/rand-to-best/1/exp strategy determined from the strategy tests, we found the exact solution for the PSIIRC exciton model and fitted the spectra with a reasonable convergence rate.

Electrostatic profiling of photosynthetic pigments: implications for directed spectral tuning

Photosynthetic pigment-protein complexes harvest solar energy with a high quantum efficiency. Protein scaffolds are known to tune the spectral properties of embedded pigments principally through structured electrostatic environments. Although the physical nature of electrostatic tuning is straightforward, the precise spatial principles of electrostatic preorganization remain poorly explored for different protein matrices and incompletely characterized with respect to the intrinsic properties of different photosynthetic pigments. In this work, we study the electronic structure features associated with the lowest excited state of a series of eight naturally occurring (bacterio)chlorophylls and pheophytins to describe the precise topological differences in electrostatic potentials and hence determine intrinsic differences in the expected mode and impact of electrostatic tuning. The difference electrostatic potentials between the ground and first excited states are used as fingerprints. Both the spatial profile and the propensity for spectral tuning are found to be unique for each pigment, indicating spatially and directionally distinct modes of electrostatic tuning. The results define a specific partitioning of the protein matrix around each pigment as an aid to identify regions with a maximal impact on spectral tuning and have direct implications for dimensionality reduction in protein design and engineering. Thus, a quantum mechanical basis is provided for understanding, predicting, and ultimately designing sequence-modified or pigment-exchanged biological systems, as suggested for selected examples of pigment-reconstituted proteins.

Modeling and Simulating the Excited-State Dynamics of a System with Condensed Phases: A Machine Learning Approach

Simulating the irreversible quantum dynamics of exciton- and electron-transfer problems poses a nontrivial challenge. Because the irreversibility of the system dynamics is a result of quantum thermal activation and dissipation caused by the surrounding environment, it is necessary to include infinite environmental degrees of freedom in the simulation. Because the capabilities of full quantum dynamics simulations that include the surrounding molecular degrees of freedom are limited, employing a system-bath model is a practical approach. In such a model, the dynamics of excitons or electrons are described by a system Hamiltonian, while the other degrees of freedom that arise from the environmental molecules are described by a harmonic oscillator bath (HOB) and system-bath interaction parameters. By extending on a previous study of molecular liquids [ J. Chem. Theory Comput. 2020, 16, 2099], here, we construct a system-bath model for exciton- and electron-transfer problems by means of a machine learning approach. We determine both the system and system-bath interaction parameters, including the spectral distribution of the bath, using the electronic excitation energies obtained from a quantum mechanics/molecular mechanics (QM/MM) simulation that is conducted as a function of time. Using the analytical expressions of optical response functions, we calculate linear and two-dimensional electronic spectra (2DES) for indocarbocyanine dimers in methanol. From these results, we demonstrate the capability of our approach to elucidate the nonequilibrium exciton dynamics of a quantum system in a nonintuitive manner.

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.

How Can We Predict Accurate Electrochromic Shifts for Biochromophores? A Case Study on the Photosynthetic Reaction Center

Protein-embedded chromophores are responsible for light harvesting, excitation energy transfer, and charge separation in photosynthesis. A critical part of the photosynthetic apparatus are reaction centers (RCs), which comprise groups of (bacterio)chlorophyll and (bacterio)pheophytin molecules that transform the excitation energy derived from light absorption into charge separation. The lowest excitation energies of individual pigments (site energies) are key for understanding photosynthetic systems, and form a prime target for quantum chemistry. A major theoretical challenge is to accurately describe the electrochromic (Stark) shifts in site energies produced by the inhomogeneous electric field of the protein matrix. Here, we present large-scale quantum mechanics/molecular mechanics calculations of electrochromic shifts for the RC chromophores of photosystem II (PSII) using various quantum chemical methods evaluated against the 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). We show that certain range-separated density functionals (ωΒ97, ωΒ97X-V, ωΒ2PLYP, and LC-BLYP) correctly reproduce RC site energy shifts with time-dependent density functional theory (TD-DFT). The popular CAM-B3LYP functional underestimates the shifts and is not recommended. Global hybrid functionals are too insensitive to the environment and should be avoided, while nonhybrid functionals are strictly nonapplicable. Among the applicable approximate coupled cluster methods, the canonical versions of CC2 and ADC(2) were found to deviate significantly from the reference results both for the description of the lowest excited state and for the electrochromic shifts. By contrast, their spin-component-scaled (SCS) and particularly the scale-opposite-spin (SOS) variants compare well with the reference DLPNO-STEOM-CCSD and the best range-separated DFT methods. The emergence of RC excitation asymmetry is discussed in terms of intrinsic and protein electrostatic potentials. In addition, we evaluate a minimal structural scaffold of PSII, the D1-D2-CytB559 RC complex often employed in experimental studies, and show that it would have the same site energy distribution of RC chromophores as the full PSII supercomplex, but only under the unlikely conditions that the core protein organization and cofactor arrangement remain identical to those of the intact enzyme.

Probing the role of arginine 323 of the D1 protein in photosystem II function

The Mn₄ CaO₅ cluster of photosystem II (PSII) advances sequentially through five oxidation states (S₀ to S₄ ). Under the enzyme cycle, two water molecules are oxidized, O₂ is generated and four protons are released into the lumen. Umena et al. (2011) have proposed that, with other charged amino acids, the R323 residue of the D1 protein could contribute to regulate a proton egress pathway from the Mn₄ CaO₅ cluster and Tyr_Z via a proton channel identified from the 3D structure. To test this suggestion, a PsbA3/R323E site-directed mutant has been constructed and the properties of its PSII have been compared to those of the PsbA3-PSII by using EPR spectroscopy, polarography, thermoluminescence and time-resolved UV-visible absorption spectroscopy. Neither the oscillations with a period four nor the kinetics and S-state-dependent stoichiometry of the proton release were affected. However, several differences have been found: (1) the P₆₈₀ ⁺ decay in the hundreds of ns time domain was much slower in the mutant, (2) the S₂ Q_A ^- /DCMU and S₃ Q_A ^- /DCMU radiative charge recombination occurred at higher temperatures and (3) the S₀ Tyr_Z ^• , S₁ Tyr_Z ^• , S₂ Tyr_Z ^• split EPR signals induced at 4.2 K by visible light from the S₀ Tyr_Z , S₁ Tyr_Z , S₂ Tyr_Z , respectively, and the (S₂ Tyr_Z ^• )' induced by NIR illumination at 4.2 K of the S₃ Tyr_Z state differed. It is proposed that the R323 residue of the D1 protein interacts with Tyr_Z likely via the H-bond network previously proposed to be a proton channel. Therefore, rather than participating in the egress of protons to the lumen, this channel could be involved in the relaxations of the H-bonds around Tyr_Z by interacting with the bulk, thus tuning the driving force required for Tyr_Z oxidation.

Normal mode analysis of spectral density of FMO trimers: Intra- and intermonomer energy transfer

The intermolecular contribution to the spectral density of the exciton-vibrational coupling of the homotrimeric Fenna-Matthews-Olson (FMO) light-harvesting protein of green sulfur bacteria P. aestuarii is analyzed by combining a normal mode analysis of the protein with the charge density coupling method for the calculation of local transition energies of the pigments. Correlations in site energy fluctuations across the whole FMO trimer are found at low vibrational frequencies. Including, additionally, the high-frequency intrapigment part of the spectral density, extracted from line-narrowing spectra, we study intra- and intermonomer exciton transfer. Whereas the intrapigment part of the spectral density is important for fast intramonomer exciton relaxation, the intermolecular contributions (due to pigment-environment coupling) determine the intermonomer exciton transfer. Neither the variations of the local Huang-Rhys factors nor the correlations in site energy fluctuations have a critical influence on energy transfer. At room temperature, the intermonomer transfer in the FMO protein occurs on a 10 ps time scale, whereas intramonomer exciton equilibration is roughly two orders of magnitude faster. At cryogenic temperatures, intermonomer transfer limits the lifetimes of the lowest exciton band. The lifetimes are found to increase between 20 ps in the center of this band up to 100 ps toward lower energies, which is in very good agreement with the estimates from hole burning data. Interestingly, exciton delocalization in the FMO monomers is found to slow down intermonomer energy transfer, at both physiological and cryogenic temperatures.

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

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 ChlD1 and PheD1 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 (ChlD1+PheD1-)* and (PD2+PD1-)*. Our results demonstrate the power of TDDFT as a tool, for studies of natural photosynthesis, or indeed future studies of artificial photosynthetic complexes.

Photosynthetic inner antenna CP47 plays important roles in ephemeral plants in adapting to high light stress

Throughout 3.5 billion years of evolution, photosynthesis of land plants has developed a complicated antenna system to cope with the ever-changing environments. The antenna system of photosystem (PS) II includes the outer antennae and inner antennae. The inner antennae CP43 and CP47, located in the closest peripheral of PSII reaction center (RC), play important roles in facilitating excitation energy transport from the outer antennae to the PSII RC. Although PSII RC is the last station of energy transport, it is the inner antenna CP47, not the RC, which possesses the lowest energy level in PSII. Berteroa incana (B. incana), which is a vascular plant grown in the Gobi region, can sustain very high photosynthesis capacity under very high light conditions. It has been discovered that the thylakoid membrane of B. incana possesses a unique low fluorescence emission spectrum because the fluorescence emission of CP47 (695 nm) is the main fluorescence emission peak of PSII. In this paper, the thylakoid membrane, isolated from B. incana grown under a light condition of 100 μM photons m^-2 s^-1 and subjected to high light treatment (1600 μM photons m^-2 s^-1 for 1.5 h or 3 h) was employed for the research. It has been found that the high fluorescence emission of CP47 decreased remarkably upon exposure to the high light treatment. Further observation revealed that the drastic changes in the CP47 fluorescence emission were accompanied by a slight reduction in the amount of CP47 and an enhancement of the CP29-LHCII-CP24 assembly. It is proposed that CP47 enables the functional switch between the excitation energy transfer towards PSII RC, and the overexcitation quenching in the PSII core. Together with CP43, CP47 plays important roles in regulating excitation energy distribution in PSII core complexes.

Asymmetry in Charge Transfer Pathways Caused by Pigment–Protein Interactions in the Photosystem II Reaction Center Complex

This article discusses the photoinduced charge transfer (CT) kinetics within the reaction center complex of photosystem II (PSII RC). The PSII RC exhibits a structural symmetry in its arrangement of pigments forming two prominent branches, D1 and D2. Despite this symmetry, the CT has been observed to occur exclusively in the D1 branch. The mechanism to realize such functional asymmetry is yet to be understood. To approach this matter, we applied the theoretical tight-binding model of pigment excitations and simulated CT dynamics based upon the framework of an open quantum system. This simulation used a recently developed method of computation based on the quasi-adiabatic propagator path integral. A quantum CT state is found to be dynamically active when its site energy is resonant with the exciton energies of the PSII RC, regardless of the excitonic landscape we utilized. Through our investigation, it was found that the relative displacement between the local molecular energy levels of pigments can play a crucial role in realizing this resonance and therefore greatly affects the CT asymmetry in the PSII RC. Using this mechanism phenomenologically, we demonstrate that a near 100-to-1 ratio of reduction between the pheophytins in the D1 and D2 branches can be realized at both 77 K and 300 K. Our results indicate that the chlorophyll Chl D 1 is the most active precursor of the primary charge separation in the D1 branch and that the reduction of the pheophytins can occur within pico-seconds. Additionally, a broad resonance of the active CT state implies that a large static disorder observed in the CT state originates in the fluctuations of the relative displacements between the local molecular energy levels of the pigments in the PSII RC.

The primary donor of far-red photosystem II: ChlD1 or PD2?

Far-red light (FRL) Photosystem II (PSII) isolated from Chroococcidiopsis thermalis is studied using parallel analyses of low-temperature absorption, circular dichroism (CD) and magnetic circular dichroism (MCD) spectroscopies in conjunction with fluorescence measurements. This extends earlier studies (Nurnberg et al 2018 Science 360 (2018) 1210-1213). We confirm that the chlorophyll absorbing at 726 nm is the primary electron donor. At 1.8 K efficient photochemistry occurs when exciting at 726 nm and shorter wavelengths; but not at wavelengths longer than 726 nm. The 726 nm absorption peak exhibits a 21 ± 4 cm^-1 electrochromic shift due to formation of the semiquinone anion, Q_A^-. Modelling indicates that no other FRL pigment is located among the 6 central reaction center chlorins: P_D1, P_D2 Chl_D1, Chl_D2, Pheo_D1 and Pheo_D2. Two of these chlorins, Chl_D1 and P_D2, are located at a distance and orientation relative to Q_A^- so as to account for the observed electrochromic shift. Previously, Chl_D1 was taken as the most likely candidate for the primary donor based on spectroscopy, sequence analysis and mechanistic arguments. Here, a more detailed comparison of the spectroscopic data with exciton modelling of the electrochromic pattern indicates that P_D2 is at least as likely as Chl_D1 to be responsible for the 726 nm absorption. The correspondence in sign and magnitude of the CD observed at 726 nm with that predicted from modelling favors P_D2 as the primary donor. The pros and cons of P_D2 vs Chl_D1 as the location of the FRL-primary donor are discussed.

What can we still learn from the electrochromic band-shifts in Photosystem II?

Electrochromic band-shifts have been investigated in Photosystem II (PSII) from Thermosynechoccocus elongatus. Firstly, by using Mn-depleted PsbA1-PSII and PsbA3-PSII in which the Q_X absorption of Phe_D1 differs, a band-shift in the Q_X region of Phe_D2 centered at ~ 544 nm has been identified upon the oxidation, at pH 8.6, of Tyr_D. In contrast, a band-shift due to the formation of either Q_A^•- or Tyr_Z^• is observed in PsbA3-PSII at ~ 546 nm, as expected with E130 H-bonded to Phe_D1 and at ~ 544 nm as expected with Q130 H-bonded to Phe_D1. Secondly, electrochromic band-shifts in the Chla Soret region have been measured in O₂-evolving PSII in PsbA3-PSII, in the PsbA3/H198Q mutant in which the Soret band of P_D1 is blue shifted and in the PsbA3/T179H mutant. Upon Tyr_Z^•Q_A^•- formation the Soret band of P_D1 is red shifted and the Soret band of Chl_D1 is blue shifted. In contrast, only P_D1 undergoes a detectable S-state dependent electrochromism. Thirdly, the time resolved S-state dependent electrochromism attributed to P_D1 is biphasic for all the S-state transitions except for S₁ to S₂, and shows that: i) the proton release in S₀ to S₁ occurs after the electron transfer and ii) the proton release and the electron transfer kinetics in S₂ to S₃, in T. elongatus, are significantly faster than often considered. The nature of S₂Tyr_Z^• is discussed in view of the models in the literature involving intermediate states in the S₂ to S₃ transition.

Quantum biology revisited

Photosynthesis is a highly optimized process from which valuable lessons can be learned about the operating principles in nature. Its primary steps involve energy transport operating near theoretical quantum limits in efficiency. Recently, extensive research was motivated by the hypothesis that nature used quantum coherences to direct energy transfer. This body of work, a cornerstone for the field of quantum biology, rests on the interpretation of small-amplitude oscillations in two-dimensional electronic spectra of photosynthetic complexes. This Review discusses recent work reexamining these claims and demonstrates that interexciton coherences are too short lived to have any functional significance in photosynthetic energy transfer. Instead, the observed long-lived coherences originate from impulsively excited vibrations, generally observed in femtosecond spectroscopy. These efforts, collectively, lead to a more detailed understanding of the quantum aspects of dissipation. Nature, rather than trying to avoid dissipation, exploits it via engineering of exciton-bath interaction to create efficient energy flow.

Structural basis of light-harvesting in the photosystem II core complex

Photosystem II (PSII) is a membrane-spanning, multi-subunit pigment-protein complex responsible for the oxidation of water and the reduction of plastoquinone in oxygenic photosynthesis. In the present review, the recent explosive increase in available structural information about the PSII core complex based on X-ray crystallography and cryo-electron microscopy is described at a level of detail that is suitable for a future structure-based analysis of light-harvesting processes. This description includes a proposal for a consistent numbering scheme of protein-bound pigment cofactors across species. The structural survey is complemented by an overview of the state of affairs in structure-based modeling of excitation energy transfer in the PSII core complex with emphasis on electrostatic computations, optical properties of the reaction center, the assignment of long-wavelength chlorophylls, and energy trapping mechanisms.

New Perspectives on Photosystem II Reaction Centres

We apply the differential optical spectroscopy techniques of circular polarisation of luminescence (CPL) and magnetic CPL (MCPL) to the study of isolated reaction centres (RCs) of photosystem II (PS II). The data and subsequent analysis provide insights into aspects of the RC chromophore site energies, exciton couplings, and heterogeneities. CPL measurements are able to identify weak luminescence associated with the unbound chlorophyll-a (Chl-a) present in the sample. The overall sign and magnitude of the CPL observed relates well to the circular dichroism (CD) of the sample. Both CD and CPL are reasonably consistent with modelling of the RC exciton structure. The MCPL observed for the free Chl-a luminescence component in the RC samples is also easily understandable, but the MCPL seen near 680nm at 1.8K is anomalous, appearing to have a narrow, strongly negative component. A negative sign is inconsistent with MCPL of (exciton coupled) Qy states of either Chl-a or pheophytin-a (Pheo-a). We propose that this anomaly may arise as a result of the luminescence from a transient excited state species created following photo-induced charge separation within the RC. A comparison of CD spectra and modelling of RC preparations having a different number of pigments suggests that the non-conservative nature of the CD spectra observed is associated with the ‘special pair’ pigments PD1 and PD2.

Why are photosynthetic reaction centres dimeric?

All photosynthetic organisms convert solar energy into chemical energy through charge separation in dimeric reaction centres. It is unknown why early reaction centres dimerised and completely displaced their monomeric ancestors. Here, we discuss several proposed explanations for reaction-centre dimerism and conclude-with only weak assumptions about the primordial dimerisation event-that the most probable explanation for the dimerism is that it arose because it enhanced light-harvesting efficiency by deepening the excitonic trap, i.e., by enhancing the rate of exciton transfer from an antenna complex and decreasing the rate of back transfer. This effect would have outweighed the negative effect dimerisation would have had on charge transfer within the reaction centre. Our argument implies that dimerisation likely occurred after the evolution of the first antennas, and it explains why the lower-energy state of the special pair is bright.

Consequences of structural modifications in cytochrome b559 on the electron acceptor side of Photosystem II

Cytb₅₅₉ in Photosystem II is a heterodimeric b-type cytochrome. The subunits, PsbE and PsbF, consist each in a membrane α-helix. Mutants were previously designed and studied in Thermosynechococcus elongatus (Sugiura et al., Biochim Biophys Acta 1847:276-285, 2015) either in which an axial histidine ligand of the haem-iron was substituted for a methionine, the PsbE/H23M mutant in which the haem was lacking, or in which the haem environment was modified, the PsbE/Y19F and PsbE/T26P mutants. All these mutants remained active showing that the haem has no structural role provided that PsbE and PsbF subunits are present. Here, we have carried on the characterization of these mutants. The following results were obtained: (i) the Y19F mutation hardly affect the E_m of Cytb₅₅₉, whereas the T26P mutation converts the haem into a form with a E_m much below 0 mV (so low that it is likely not reducible by Q_B^-) even in an active enzyme; (ii) in the PsbE/H23M mutant, and to a less extent in PsbE/T26P mutant, the electron transfer efficiency from Q_A^- to Q_B is decreased; (iii) the lower E_m of the Q_A/Q_A^- couple in the PsbE/H23M mutant correlates with a higher production of singlet oxygen; (iv) the superoxide and/or hydroperoxide formation was not increased in the PsbE/H23M mutant lacking the haem, whereas it was significantly larger in the PsbE/T26P. These data are discussed in view of the literature to discriminate between structural and redox roles for the haem of Cytb₅₅₉ in the production of reactive oxygen species.

New insights on ChlD1 function in Photosystem II from site-directed mutants of D1/T179 in Thermosynechococcus elongatus

The monomeric chlorophyll, Chl_D1, which is located between the P_D1P_D2 chlorophyll pair and the pheophytin, Pheo_D1, is the longest wavelength chlorophyll in the heart of Photosystem II and is thought to be the primary electron donor. Its central Mg²⁺ is liganded to a water molecule that is H-bonded to D1/T179. Here, two site-directed mutants, D1/T179H and D1/T179V, were made in the thermophilic cyanobacterium, Thermosynechococcus elongatus, and characterized by a range of biophysical techniques. The Mn₄CaO₅ cluster in the water-splitting site is fully active in both mutants. Changes in thermoluminescence indicate that i) radiative recombination occurs via the repopulation of *Chl_D1 itself; ii) non-radiative charge recombination reactions appeared to be faster in the T179H-PSII; and iii) the properties of P_D1P_D2 were unaffected by this mutation, and consequently iv) the immediate precursor state of the radiative excited state is the Chl_D1⁺Pheo_D1^- radical pair. Chlorophyll bleaching due to high intensity illumination correlated with the amount of ¹O₂ generated. Comparison of the bleaching spectra with the electrochromic shifts attributed to Chl_D1 upon Q_A^- formation, indicates that in the T179H-PSII and in the WT*3-PSII, the Chl_D1 itself is the chlorophyll that is first damaged by ¹O₂, whereas in the T179V-PSII a more red chlorophyll is damaged, the identity of which is discussed. Thus, Chl_D1 appears to be one of the primary damage site in recombination-mediated photoinhibition. Finally, changes in the absorption of Chl_D1 very likely contribute to the well-known electrochromic shifts observed at ~430 nm during the S-state cycle.

A theoretical study on the dynamics of light harvesting in the dimeric photosystem II core complex: regulation and robustness of energy transfer pathways

Here we present our theoretical investigations into the light reaction in the dimeric photosystem II (PSII) core complex. An effective model for excitation energy transfer (EET) and primary charge separation (CS) in the PSII core complex was developed, with model parameters constructed based on molecular dynamics (MD) simulation data. Compared to experimental results, we demonstrated that this model faithfully reproduces the absorption spectra of the RC and core light-harvesting complexes (CP43 and CP47) as well as the full EET dynamics among the chromophores in the PSII core complex. We then applied master equation simulations and network analysis to investigate detailed EET plus CS dynamics in the system, allowing us to identify key EET pathways and produce a coarse-grained cluster model for the light reaction in the dimeric PSII core complex. We show that non-equilibrium energy transfer channels play important roles in the efficient light harvesting process and that multiple EET pathways exist between subunits of PSII to ensure the robustness of light harvesting in the system. Furthermore, we revealed that inter-monomer energy transfer dominated by the coupling between the two CLA625 molecules enables efficient energy exchange between two CP47s in the dimeric PSII core complex, which leads to significant energy pooling in the CP47 domain during the light reaction. Our study provides a blueprint for the design of light harvesting in the PSII core and show that a structure-based approach using molecular dynamics simulations and quantum chemistry calculations can be effectively utilized to elucidate the dynamics of light harvesting in complex photosynthetic systems.

Redox-state dependent blinking of single photosystem I trimers at around liquid-nitrogen temperature

Efficient light harvesting in a photosynthetic antenna system is disturbed by a ragged and fluctuating energy landscape of the antenna pigments in response to the conformation dynamics of the protein. This situation is especially pronounced in Photosystem I (PSI) containing red shifted chlorophylls (red Chls) with the excitation energy much lower than the primary donor. The present study was conducted to clarify light-harvesting dynamics of PSI isolated from Synechocystis sp. PCC6803 by using single-molecule spectroscopy at liquid‑nitrogen temperatures. Fluorescence emission at around 720 nm from the red Chls in single PSI trimers was monitored at 80-100 K. Intermittent variations in the emission intensities, so-called blinking, were frequently observed. Its time scale lay in several tens of seconds. The blinking amplitude depended on the redox state of the phylloquinone (A₁). Electrochromic shifts of Chls induced by the negative charge on A₁ were calculated based on the X-ray crystallographic structure. A Chl molecule, Chl-A839 (numbering according to PDB 5OY0), bound near A₁ was found to have a large electrochromic shift. This Chl has strong exciton coupling with neighboring Chl (A838) whose site energy was predicted to be determined by interaction with an arginine residue (ArgF84) [Adolphs et al., 2010]. A possible scenario of the blinking was proposed. Conformational fluctuations of ArgF84 seesaw the excitation-energy of Chl-A838, which perturbs the branching ratio of excitation-energy between the red Chl and the cationic form of P700 as a quencher. The electrochromic shift of Chl-A839 enhances the effect of the conformation dynamics of ArgF84.

Photochemistry beyond the red limit in chlorophyll f-containing photosystems

Photosystems I and II convert solar energy into the chemical energy that powers life. Chlorophyll a photochemistry, using red light (680 to 700 nm), is near universal and is considered to define the energy "red limit" of oxygenic photosynthesis. We present biophysical studies on the photosystems from a cyanobacterium grown in far-red light (750 nm). The few long-wavelength chlorophylls present are well resolved from each other and from the majority pigment, chlorophyll a. Charge separation in photosystem I and II uses chlorophyll f at 745 nm and chlorophyll f (or d) at 727 nm, respectively. Each photosystem has a few even longer-wavelength chlorophylls f that collect light and pass excitation energy uphill to the photochemically active pigments. These photosystems function beyond the red limit using far-red pigments in only a few key positions.

Wavelength-Dependent Exciton-Vibrational Coupling in the Water-Soluble Chlorophyll Binding Protein Revealed by Multilevel Theory of Difference Fluorescence Line-Narrowing

One of the most powerful line-narrowing techniques used to unravel the homogeneous lineshapes of inhomogeneously broadened systems is difference fluorescence line-narrowing spectroscopy. When this spectroscopy was applied to multichromophoric systems so far, the spectra were analyzed by an effective two-level system approach, composed of the electronic ground state and the lowest exciton state. An effective Huang-Rhys factor was assigned for the coupling of this state to the vibrations. Here, we extend this approach by including a multilevel line shape theory, which takes into account the excitonic coupling between pigments and thereby the effect of the delocalization of the excited states explicitly. In this way, it becomes possible to extract the spectral density of the local exciton-vibrational coupling. The theory is applied to the recombinant water-soluble chlorophyll binding protein reconstituted with chlorophyll a or b and reveals a significant decrease of the Huang-Rhys factor of the local exciton-vibrational coupling with decreasing transition energy of the chlorophylls. This decrease could be due to the increase in steric interactions reducing the flexibility of the environment and red-shifting the site energy of the pigments.

Strategy for Modeling the Electrostatic Responses of the Spectroscopic Properties of Proteins

For better understanding and more efficient use of the spectroscopic probes (vibrational and NMR) of the local electrostatic situations inside proteins, appropriate modeling of the properties of those probes is essential. The present study is devoted to examining the strategy for constructing such models. A more well-founded derivation than the ones in previous studies is given in constructing the models. Theoretical analyses are conducted on two representative example cases related to proteins, i.e., the peptide group of the main chains and the CO and NO ligands to the Fe²⁺ ion of heme, with careful treatment of the behavior of electrons in the electrostatic responses and with verification of consistency with observable quantities. It is shown that, for the stretching frequencies and NMR chemical shifts, it is possible to construct reasonable electrostatic interaction models that encompass the situations of hydration and uniform electric field environment and thus are applicable also to the cases of nonuniform electrostatic situations, which are highly expected for inside of proteins.

Exciton-Coupled Electron Transfer Process Controlled by Non-Markovian Environments

We theoretically investigate an exciton-coupled electron transfer (XCET) process that is conversion of an exciton into a charge transfer state. This conversion happens in an exciton transfer (XT) process, and the electron moves away in an electron transfer (ET) process in multiple environments (baths). This XCET process plays an essential role in the harvesting of solar energy in biological and photovoltaic materials. We develop a practical theoretical model to study the efficiency of the XCET process that occurs either in consecutive or concerted processes under the influence of non-Markovian baths. The role of quantum coherence in the XT-ET system and the baths is investigated using reduced hierarchal equations of motion (HEOM). This model includes independent baths for each XT and ET state, in addition to a XCET bath for the conversion process. We found that, while quantum system-bath coherence is important in the XT and ET processes, coherence between the XT and ET processes must be suppressed in order to realize that an efficient irreversible XCET process through the weak off-diagonal interaction between the XT and ET bridge sites arises from an XCET bath.