Multiple charge-separation pathways in photosystem II: modeling of transient absorption kinetics

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.

Perspective on Integrative Simulations of Bioenergetic Domains

Bioenergetic processes in cells, such as photosynthesis or respiration, integrate many time and length scales, which makes the simulation of energy conversion with a mere single level of theory impossible. Just like the myriad of experimental techniques required to examine each level of organization, an array of overlapping computational techniques is necessary to model energy conversion. Here, a perspective is presented on recent efforts for modeling bioenergetic phenomena with a focus on molecular dynamics simulations and its variants as a primary method. An overview of the various classical, quantum mechanical, enhanced sampling, coarse-grained, Brownian dynamics, and Monte Carlo methods is presented. Example applications discussed include multiscale simulations of membrane-wide electron transport, rate kinetics of ATP turnover from electrochemical gradients, and finally, integrative modeling of the chromatophore, a photosynthetic pseudo-organelle.

Unraveling quantum coherences mediating primary charge transfer processes in photosystem II reaction center

Photosystem II (PSII) reaction center (RC) is a unique complex that is capable of efficiently separating electronic charges across the membrane. The primary energy- and charge-transfer (CT) processes occur on comparable ultrafast timescales, which makes it extremely challenging to understand the fundamental mechanism responsible for the near-unity quantum efficiency of the transfer. Here, we elucidate the role of quantum coherences in the ultrafast energy and CT in the PSII RC by performing two-dimensional (2D) electronic spectroscopy at the cryogenic temperature of 20 kelvin, which captures the distinct underlying quantum coherences. Specifically, we uncover the electronic and vibrational coherences along with their lifetimes during the primary ultrafast processes of energy and CT. We construct an excitonic model that provides evidence for coherent energy and CT at low temperature in the 2D electronic spectra. The principles could provide valuable guidelines for creating artificial photosystems with exploitation of system-bath coupling and control of coherences to optimize the photon conversion efficiency to specific functions.

Far-red photosynthesis: Two charge separation pathways exist in plant Photosystem II reaction center

An alternative charge separation pathway in Photosystem II under the far-red light was proposed by us on the basis of electron transfer properties at 295 K and 5 K. Here we extend these studies to the temperature range of 77-295 K with help of electron paramagnetic resonance spectroscopy. Induction of the S2 state multiline signal, oxidation of Cytochrome b559 and ChlorophyllZ was studied in Photosystem II membrane preparations from spinach after application of a laser flashes in visible (532 nm) or far-red (730-750 nm) spectral regions. Temperature dependence of the S2 state signal induction after single flash at 730-750 nm (Tinhibition ~ 240 K) was found to be different than that at 532 nm (Tinhibition ~ 157 K). No contaminant oxidation of the secondary electron donors cytochrome b559 or chlorophyllZ was observed. Photoaccumulation experiments with extensive flashing at 77 K showed similar results, with no or very little induction of the secondary electron donors. Thus, the partition ratio defined as (yield of YZ/CaMn4O5-cluster oxidation):(yield of Cytb559/ChlZ/CarD2 oxidation) was found to be 0.4 at under visible light and 1.7 at under far-red light at 77 K. Our data indicate that different products of charge separation after far-red light exists in the wide temperature range which further support the model of the different primary photochemistry in Photosystem II with localization of hole on the ChlD1 molecule.

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

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.

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.

Photosynthetic Light-Harvesting (Antenna) Complexes-Structures and Functions

Chlorophylls and bacteriochlorophylls, together with carotenoids, serve, noncovalently bound to specific apoproteins, as principal light-harvesting and energy-transforming pigments in photosynthetic organisms. In recent years, enormous progress has been achieved in the elucidation of structures and functions of light-harvesting (antenna) complexes, photosynthetic reaction centers and even entire photosystems. It is becoming increasingly clear that light-harvesting complexes not only serve to enlarge the absorption cross sections of the respective reaction centers but are vitally important in short- and long-term adaptation of the photosynthetic apparatus and regulation of the energy-transforming processes in response to external and internal conditions. Thus, the wide variety of structural diversity in photosynthetic antenna “designs” becomes conceivable. It is, however, common for LHCs to form trimeric (or multiples thereof) structures. We propose a simple, tentative explanation of the trimer issue, based on the 2D world created by photosynthetic membrane systems.

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.

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.

Two-dimensional electronic Stark spectroscopy of a photosynthetic dimer

Stark spectroscopy, which measures changes in the linear absorption of a sample in the presence of an external DC electric field, is a powerful experimental tool for probing the existence of charge-transfer (CT) states in photosynthetic systems. CT states often have small transition dipole moments, making them insensitive to other spectroscopic methods, but are particularly sensitive to Stark spectroscopy due to their large permanent dipole moment. In a previous study, we demonstrated a new experimental method, two-dimensional electronic Stark spectroscopy (2DESS), which combines two-dimensional electronic spectroscopy (2DES) and Stark spectroscopy. In order to understand how the presence of CT states manifest in 2DESS, here, we perform computational modeling and calculations of 2DESS as well as 2DES and Stark spectra, studying a photosynthetic dimer inspired by the photosystem II reaction center. We identify specific cases where qualitatively different sets of system parameters produce similar Stark and 2DES spectra but significantly different 2DESS spectra, showing the potential for 2DESS to aid in identifying CT states and their coupling to excitonic states.

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.

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.

Electron Transfer in Nitrogenase

Nitrogenase is the only enzyme capable of reducing N2 to NH3. This challenging reaction requires the coordinated transfer of multiple electrons from the reductase, Fe-protein, to the catalytic component, MoFe-protein, in an ATP-dependent fashion. In the last two decades, there have been significant advances in our understanding of how nitrogenase orchestrates electron transfer (ET) from the Fe-protein to the catalytic site of MoFe-protein and how energy from ATP hydrolysis transduces the ET processes. In this review, we summarize these advances, with focus on the structural and thermodynamic redox properties of nitrogenase component proteins and their complexes, as well as on new insights regarding the mechanism of ET reactions during catalysis and how they are coupled to ATP hydrolysis. We also discuss recently developed chemical, photochemical, and electrochemical methods for uncoupling substrate reduction from ATP hydrolysis, which may provide new avenues for studying the catalytic mechanism of nitrogenase.

Ultrafast Dynamics of Nonrigid Zinc-Porphyrin Arrays Mimicking the Photosynthetic "Special Pair"

Conjugated porphyrin arrays are heavily investigated as efficient molecular systems for photosynthesis and photocatalysis. Recently, a series of one-, two-, and six-zinc-porphyrin arrays, noncovalently linked through benzene-based hubs, have been synthesized with the aim of mimicking the structure and function of the bacteriochlorophyll "special pair" in photosynthetic reaction centers. The excitonically coupled porphyrin subunits are expected to activate additional excited state relaxation channels with respect to the monomer. Here, we unveil the appearance of such supramolecular electronic interactions using ultrafast transient absorption spectroscopy with sub-25 fs time resolution. Upon photoexcitation of the Soret band, we resolve energy trapping within ∼150 fs in a delocalized dark excitonic manifold. Moreover, excitonic interactions promote an additional fast internal conversion from the Q-band to the ground state with an efficiency of up to 60% in the hexamer. These relaxation pathways appear to be common loss channels that limit the lifetime of the exciton states in noncovalently bound molecular aggregates.

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.

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.

Quantum Chemical Modeling of the Photoinduced Activity of Multichromophoric Biosystems

Multichromophoric biosystems represent a broad family with very diverse members, ranging from light-harvesting pigment-protein complexes to nucleic acids. The former are designed to capture, harvest, efficiently transport, and transform energy from sunlight for photosynthesis, while the latter should dissipate the absorbed radiation as quickly as possible to prevent photodamages and corruption of the carried genetic information. Because of the unique electronic and structural characteristics, the modeling of their photoinduced activity is a real challenge. Numerous approaches have been devised building on the theoretical development achieved for single chromophores and on model Hamiltonians that capture the essential features of the system. Still, a question remains: is a general strategy for the accurate modeling of multichromophoric systems possible? By using a quantum chemical point of view, here we review the advancements developed so far highlighting differences and similarities with the single chromophore treatment. Finally, we outline the important limitations and challenges that still need to be tackled to reach a complete and accurate picture of their photoinduced properties and dynamics.

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.

Intramolecular Vibrations Complement the Robustness of Primary Charge Separation in a Dimer Model of the Photosystem II Reaction Center

The energy conversion of oxygenic photosynthesis is triggered by primary charge separation in proteins at the photosystem II reaction center. Here, we investigate the impacts of the protein environment and intramolecular vibrations on primary charge separation at the photosystem II reaction center. This is accomplished by combining the quantum dynamic theories of condensed phase electron transfer with quantum chemical calculations to evaluate the vibrational Huang-Rhys factors of chlorophyll and pheophytin molecules. We report that individual vibrational modes play a minor role in promoting charge separation, contrary to the discussion in recent publications. Nevertheless, these small contributions accumulate to considerably influence the charge separation rate, resulting in subpicosecond charge separation almost independent of the driving force and temperature. We suggest that the intramolecular vibrations complement the robustness of the charge separation in the photosystem II reaction center against the inherently large static disorder of the involved electronic energies.

The wavelength of the incident light determines the primary charge separation pathway in Photosystem II

Charge separation is a key component of the reactions cascade of photosynthesis, by which solar energy is converted to chemical energy. From this photochemical reaction, two radicals of opposite charge are formed, a highly reducing anion and a highly oxidising cation. We have previously proposed that the cation after far-red light excitation is located on a component different from P_D1, which is the location of the primary electron hole after visible light excitation. Here, we attempt to provide further insight into the location of the primary charge separation upon far-red light excitation of PS II, using the EPR signal of the spin polarized ³P₆₈₀ as a probe. We demonstrate that, under far-red light illumination, the spin polarized ³P₆₈₀ is not formed, despite the primary charge separation still occurring at these conditions. We propose that this is because under far-red light excitation, the primary electron hole is localized on Chl_D1, rather than on P_D1. The fact that identical samples have demonstrated charge separation upon both far-red and visible light excitation supports our hypothesis that two pathways for primary charge separation exist in parallel in PS II reaction centres. These pathways are excited and activated dependent of the wavelength applied.

Exciton-vibrational resonance and dynamics of charge separation in the photosystem II reaction center

The dynamics of charge separation in the photosystem II reaction center (PSII-RC) in the presence of intramolecular vibrations with their frequency matching the energy gap between the exciton state acting as the primary electron donor and the first charge-transfer (CT) state are investigated. A reduced PSII-RC 4-state model explicitly including a CT state is analyzed within Redfield relaxation theory in the multidimensional exciton-vibrational (vibronic) basis. This model is used to study coherent energy/electron transfers and their spectral signatures obtained by two-dimensional electronic spectroscopy (2DES). Modeling of the time-resolved 2D frequency maps obtained by wavelet analysis reveals the origins of the coherences which produce the observed oscillating features in 2DES and allows comparing the lifetimes of the coherences. The results suggest faster excitonic decoherence as compared with longer-lived vibronic oscillations. The emerging picture of the dynamics unravels the role of resonant vibrations in sustaining the effective energy conversion in the PSII-RC. We demonstrate that the mixing of the exciton and CT states promoted by a resonant vibrational quantum allows faster penetration of excitation energy into the CT with subsequent dynamic localization at the bottom of the CT potential induced by the remaining non-resonant nuclear modes. The degree of vibration-assisted mixing and, correspondingly, the rate of primary charge separation, increases significantly in the case of electron-vibrational resonance. The observed features illustrate the principles of quantum design of the photosynthetic unit. These principles are connected with the phenomenon of coherent mixing within vibronic eigenstates, increasing the effectiveness of charge separation not only upon coherent and impulsive laser excitation utilized in the 2DES experiment, but also under natural conditions under non-coherent non-impulsive solar light illumination.

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.

On the performance of a photosystem II reaction centre-based photocell

The photosystem II reaction centre is the photosynthetic complex responsible for oxygen production on Earth. Its water splitting function is particularly favoured by the formation of a stable charge separated state via a pathway that starts at an accessory chlorophyll. Here we envision a photovoltaic device that places one of these complexes between electrodes and investigate how the mean current and its fluctuations depend on the microscopic interactions underlying charge separation in the pathway considered. Our results indicate that coupling to well resolved vibrational modes does not necessarily offer an advantage in terms of power output but can lead to photo-currents with suppressed noise levels characterizing a multi-step ordered transport process. Besides giving insight into the suitability of these complexes for molecular-scale photovoltaics, our work suggests a new possible biological function for the vibrational environment of photosynthetic reaction centres, namely, to reduce the intrinsic current noise for regulatory processes.

Primary Charge Separation in the Photosystem II Reaction Center Revealed by a Global Analysis of the Two-dimensional Electronic Spectra

The transfer of electronic charge in the reaction center of Photosystem II is one of the key building blocks of the conversion of sunlight energy into chemical energy within the cascade of the photosynthetic reactions. Since the charge transfer dynamics is mixed with the energy transfer dynamics, an effective tool for the direct resolution of charge separation in the reaction center is still missing. Here, we use experimental two-dimensional optical photon echo spectroscopy in combination with the theoretical calculation to resolve its signature. A global fitting analysis allows us to clearly and directly identify a decay pathway associated to the primary charge separation. In particular, it can be distinguished from regular energy transfer and occurs on a time scale of 1.5 ps under ambient conditions. This technique provides a general tool to identify charge separation signatures from the energy transport in two-dimensional optical spectroscopy.

Quantum design of photosynthesis for bio-inspired solar-energy conversion

Photosynthesis is the natural process that converts solar photons into energy-rich products that are needed to drive the biochemistry of life. Two ultrafast processes form the basis of photosynthesis: excitation energy transfer and charge separation. Under optimal conditions, every photon that is absorbed is used by the photosynthetic organism. Fundamental quantum mechanics phenomena, including delocalization, underlie the speed, efficiency and directionality of the charge-separation process. At least four design principles are active in natural photosynthesis, and these can be applied practically to stimulate the development of bio-inspired, human-made energy conversion systems.

Quantum - coherent dynamics in photosynthetic charge separation revealed by wavelet analysis

Experimental/theoretical evidence for sustained vibration-assisted electronic (vibronic) coherence in the Photosystem II Reaction Center (PSII RC) indicates that photosynthetic solar-energy conversion might be optimized through the interplay of electronic and vibrational quantum dynamics. This evidence has been obtained by investigating the primary charge separation process in the PSII RC by two-dimensional electronic spectroscopy (2DES) and Redfield modeling of the experimental data. However, while conventional Fourier transform analysis of the 2DES data allows oscillatory signatures of vibronic coherence to be identified in the frequency domain in the form of static 2D frequency maps, the real-time evolution of the coherences is lost. Here we apply for the first time wavelet analysis to the PSII RC 2DES data to obtain time-resolved 2D frequency maps. These maps allow us to demonstrate that (i) coherence between the excitons initiating the two different charge separation pathways is active for more than 500 fs, and (ii) coherence between exciton and charge-transfer states, the reactant and product of the charge separation reaction, respectively; is active for at least 1 ps. These findings imply that the PSII RC employs coherence (i) to sample competing electron transfer pathways, and ii) to perform directed, ultrafast and efficient electron transfer.

Electrostatic Asymmetry in the Reaction Center of Photosystem II

The exciton Hamiltonian of the chlorophyll (Chl) and pheophytin (Pheo) pigments in the reaction center (RC) of photosystem II is computed based on recent crystal structures by using the Poisson-Boltzmann/quantum-chemical method. Computed site energies largely confirm a previous model inferred from fits of optical spectra, in which Chl_D1 has the lowest site energy, while that of Pheo_D1 is higher than that of Pheo_D2. The latter assignment has been challenged recently under reference to mutagenesis experiments. We argue that these data are not in contradiction to our results. We conclude that Chl_D1 is the primary electron donor in both isolated RCs and intact core complexes at least at cryogenic temperatures. The main source of asymmetry in site energies is the charge distribution in the protein. Because many small contributions from various structural elements have to be taken into account, it can be assumed that this asymmetry was established in evolution by global optimization of the RC protein.

Long-range energy transport in photosystem II

We simulate the long-range inter-complex electronic energy transfer in photosystem II-from the antenna complex, via a core complex, to the reaction center-using a non-Markovian (ZOFE) quantum master equation description that allows the electronic coherence involved in the energy transfer to be explicitly included at all length scales. This allows us to identify all locations where coherence is manifested and to further identify the pathways of the energy transfer in the full network of coupled chromophores using a description based on excitation probability currents. We investigate how the energy transfer depends on the initial excitation-localized, coherent initial excitation versus delocalized, incoherent initial excitation-and find that the overall energy transfer is remarkably robust with respect to such strong variations of the initial condition. To explore the importance of vibrationally enhanced transfer and to address the question of optimization in the system parameters, we systematically vary the strength of the coupling between the electronic and the vibrational degrees of freedom. We find that the natural parameters lie in a (broad) region that enables optimal transfer efficiency and that the overall long-range energy transfer on a ns time scale appears to be very robust with respect to variations in the vibronic coupling of up to an order of magnitude. Nevertheless, vibrationally enhanced transfer appears to be crucial to obtain a high transfer efficiency, with the latter falling sharply for couplings outside the optimal range. Comparison of our full quantum simulations to results obtained with a "classical" rate equation based on a modified-Redfield/generalized-Förster description previously used to simulate energy transfer dynamics in the entire photosystem II complex shows good agreement for the overall time scales of excitation energy transport.

How exciton-vibrational coherences control charge separation in the photosystem II reaction center

In photosynthesis absorbed sun light produces collective excitations (excitons) that form a coherent superposition of electronic and vibrational states of the individual pigments. Two-dimensional (2D) electronic spectroscopy allows a visualization of how these coherences are involved in the primary processes of energy and charge transfer. Based on quantitative modeling we identify the exciton-vibrational coherences observed in 2D photon echo of the photosystem II reaction center (PSII-RC). We find that the vibrations resonant with the exciton splittings can modify the delocalization of the exciton states and produce additional states, thus promoting directed energy transfer and allowing a switch between the two charge separation pathways. We conclude that the coincidence of the frequencies of the most intense vibrations with the splittings within the manifold of exciton and charge-transfer states in the PSII-RC is not occurring by chance, but reflects a fundamental principle of how energy conversion in photosynthesis was optimized.

Light Absorption and Energy Transfer in the Antenna Complexes of Photosynthetic Organisms

The process of photosynthesis is initiated by the capture of sunlight by a network of light-absorbing molecules (chromophores), which are also responsible for the subsequent funneling of the excitation energy to the reaction centers. Through evolution, genetic drift, and speciation, photosynthetic organisms have discovered many solutions for light harvesting. In this review, we describe the underlying photophysical principles by which this energy is absorbed, as well as the mechanisms of electronic excitation energy transfer (EET). First, optical properties of the individual pigment chromophores present in light-harvesting antenna complexes are introduced, and then we examine the collective behavior of pigment-pigment and pigment-protein interactions. The description of energy transfer, in particular multichromophoric antenna structures, is shown to vary depending on the spatial and energetic landscape, which dictates the relative coupling strength between constituent pigment molecules. In the latter half of the article, we focus on the light-harvesting complexes of purple bacteria as a model to illustrate the present understanding of the synergetic effects leading to EET optimization of light-harvesting antenna systems while exploring the structure and function of the integral chromophores. We end this review with a brief overview of the energy-transfer dynamics and pathways in the light-harvesting antennas of various photosynthetic organisms.

Efficiency of energy funneling in the photosystem II supercomplex of higher plants

The investigation of energy transfer properties in photosynthetic multi-protein networks gives insight into their underlying design principles. Here, we discuss the excitonic energy transfer mechanisms of the photosystem II (PS-II) C2S2M2 supercomplex, which is the largest isolated functional unit of the photosynthetic apparatus of higher plants. Despite the lack of a definite energy gradient in C2S2M2, we show that the energy transfer is directed by relaxation to low energy states. C2S2M2 is not organized to form pathways with strict energetically downhill transfer, which has direct consequences for the transfer efficiency, transfer pathways and transfer limiting steps. The exciton dynamics is sensitive to small changes in the energetic layout which, for instance, are induced by the reorganization of vibrational coordinates. In order to incorporate the reorganization process in our numerical simulations, we go beyond rate equations and use the hierarchically coupled equation of motion approach (HEOM). While transfer from the peripheral antenna to the proteins in proximity to the reaction center occurs on a faster time scale, the final step of the energy transfer to the RC core is rather slow, and thus the limiting step in the transfer chain. Our findings suggest that the structure of the PS-II supercomplex guarantees photoprotection rather than optimized efficiency.

Challenges facing an understanding of the nature of low-energy excited states in photosynthesis

While the majority of the photochemical states and pathways related to the biological capture of solar energy are now well understood and provide paradigms for artificial device design, additional low-energy states have been discovered in many systems with obscure origins and significance. However, as low-energy states are naively expected to be critical to function, these observations pose important challenges. A review of known properties of low energy states covering eight photochemical systems, and options for their interpretation, are presented. A concerted experimental and theoretical research strategy is suggested and outlined, this being aimed at providing a fully comprehensive understanding.

A multi-pathway model for photosynthetic reaction center

Charge separation occurs in a pair of tightly coupled chlorophylls at the heart of photosynthetic reaction centers of both plants and bacteria. Recently it has been shown that quantum coherence can, in principle, enhance the efficiency of a solar cell, working like a quantum heat engine. Here, we propose a biological quantum heat engine (BQHE) motivated by Photosystem II reaction center (PSII RC) to describe the charge separation. Our model mainly considers two charge-separation pathways which is more than that typically considered in the published literature. We explore how these cross-couplings increase the current and power of the charge separation and discuss the effects of multiple pathways in terms of current and power. The robustness of the BQHE against the charge recombination in natural PSII RC and dephasing induced by environments is also explored, and extension from two pathways to multiple pathways is made. These results suggest that noise-induced quantum coherence helps to suppress the influence of acceptor-to-donor charge recombination, and besides, nature-mimicking architectures with engineered multiple pathways for charge separations might be better for artificial solar energy devices considering the influence of environments.

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.

Role of coherent vibrations in energy transfer and conversion in photosynthetic pigment-protein complexes

Oscillatory features of two-dimensional spectra of photosynthetic pigment-protein complexes during few picoseconds after electronic excitations of chlorophylls in various pigment-proteins were recently related to the coherent nuclear vibrations. It has been also speculated that the vibrations may assist the excitonic energy transfer and charge separation, hence contributing to energy transport and energy conversion efficiency. Here, we consider three theoretical approaches usually used for characterization of the excitation dynamics and charge separation, namely Redfield, Förster, and Marcus model descriptions, regarding this question. We show that two out of the three mechanisms require explicit resonances of excitonic splittings and the nuclear vibration frequencies. However, the third one related to the electron transfer is in principle off resonant.

Mixing of exciton and charge-transfer states in light-harvesting complex Lhca4

Structure-based modeling of spectra of the wild-type Lhca4 and NH mutant enables us to build the exciton model of the complex that includes a charge-transfer state mixed with the excited-state manifold.

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.

Analysis of electron donors in photosystems in oxygenic photosynthesis by photo-CIDNP MAS NMR

Both photosystem I and photosystem II are considerably similar in molecular architecture but they operate at very different electrochemical potentials. The origin of the different redox properties of these RCs is not yet clear. In recent years, insight was gained into the electronic structure of photosynthetic cofactors through the application of photochemically induced dynamic nuclear polarization (photo-CIDNP) with magic-angle spinning NMR (MAS NMR). Non-Boltzmann populated nuclear spin states of the radical pair lead to strongly enhanced signal intensities that allow one to observe the solid-state photo-CIDNP effect from both photosystem I and II from isolated reaction center of spinach (Spinacia oleracea) and duckweed (Spirodela oligorrhiza) and from the intact cells of the cyanobacterium Synechocystis by (13)C and (15)N MAS NMR. This review provides an overview on the photo-CIDNP MAS NMR studies performed on PSI and PSII that provide important ingredients toward reconstruction of the electronic structures of the donors in PSI and PSII.

Primary electron transfer processes in photosynthetic reaction centers from oxygenic organisms

This minireview is written in honor of Vladimir A. Shuvalov, a pioneer in the area of primary photochemistry of both oxygenic and anoxygenic photosyntheses (See a News Report: Allakhverdiev et al. 2014). In the present paper, we describe the current state of the formation of the primary and secondary ion-radical pairs within photosystems (PS) II and I in oxygenic organisms. Spectral-kinetic studies of primary events in PS II and PS I, upon excitation by ~20 fs laser pulses, are now available and reviewed here; for PS II, excitation was centered at 710 nm, and for PS I, it was at 720 nm. In PS I, conditions were chosen to maximally increase the relative contribution of the direct excitation of the reaction center (RC) in order to separate the kinetics of the primary steps of charge separation in the RC from that of the excitation energy transfer in the antenna. Our results suggest that the sequence of the primary electron transfer reactions is P680 → ChlD1 → PheD1 → QA (PS II) and P700 → A 0A/A 0B → A 1A/A 1B (PS I). However, alternate routes of charge separation in PS II, under different excitation conditions, are not ruled out.

Spectral exhibition of electron-vibrational relaxation in P* state of Rhodobacter sphaeroides reaction centers

Electron-vibrational relaxation in the excited state of the primary electron donor, bacteriochlorophyll dimer P, in the reaction centers (RCs) of purple photosynthetic bacteria Rhodobacter sphaeroides is modeled. A multimode model of three states (i.e., the ground state Pg, initially excited P1*, and relaxed excited P2*) is used to calculate the incoherent dynamics of the difference (ΔA) spectra on a femtosecond timescale for the YM210 W mutant RCs. The relaxation processes are described by the step-ladder model. The model shows that the electron-vibrational relaxation in the excited state of P is visualized by the transient red shift of the stimulated emission from P*. The dynamics of this shift is observed as a change in the ΔA spectrum shape in its red-most part, within a few hundreds of femtoseconds after excitation. As a result, an initial rise in the red-side ΔA kinetics is delayed with respect to the blue-side kinetics. The time constant of the P1* → P2* electronic relaxation (54 fs) and the Pg, P1*, and P2* vibrational relaxations (120 fs), used in the model, provided the best fit of the experimental time-resolved ΔA spectra and kinetics at 90 and 293 K. The possible nature of the P1* → P2* electronic relaxation is discussed.