Primary Photophysical Processes in Photosystem II: Bridging the Gap between Crystal Structure and Optical Spectra

Excitation landscape of the CP43 photosynthetic antenna complex from multiscale simulations

Photosystem II (PSII), the principal enzyme of oxygenic photosynthesis, contains two integral light harvesting proteins (CP43 and CP47) that bind chlorophylls and carotenoids. The two intrinsic antennae play crucial roles in excitation energy transfer and photoprotection. CP43 interacts most closely with the reaction center of PSII, specifically with the branch of the reaction center (D1) that is responsible for primary charge separation and electron transfer. Deciphering the function of CP43 requires detailed atomic-level insights into the properties of the embedded pigments. To advance this goal, we employ a range of multiscale computational approaches to determine the site energies and excitonic profile of CP43 chlorophylls, using large all-atom models of a membrane-bound PSII monomer. In addition to time-dependent density functional theory (TD-DFT) used in the context of a quantum-mechanics/molecular-mechanics setup (QM/MM), we present a thorough analysis using the perturbed matrix method (PMM), which enables us to utilize information from long-timescale molecular dynamics simulations of native PSII-complexed CP43. The excited state energetics and excitonic couplings have both similarities and differences compared with previous experimental fits and theoretical calculations. Both static TD-DFT and dynamic PMM results indicate a layered distribution of site energies and reveal specific groups of chlorophylls that have shared contributions to low-energy excitations. Importantly, the contribution to the lowest energy exciton does not arise from the same chlorophylls at each system configuration, but rather changes as a function of conformational dynamics. An unexpected finding is the identification of a low-energy charge-transfer excited state within CP43 that involves a lumenal (C2) and the central (C10) chlorophyll of the complex. The results provide a refined basis for structure-based interpretation of spectroscopic observations and for further deciphering excitation energy transfer in oxygenic photosynthesis.

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.

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.

Quasiparticle Self-Consistent GW-Bethe-Salpeter Equation Calculations for Large Chromophoric Systems

The GW-Bethe–Salpeter equation (BSE) method is promising for calculating the low-lying excitonic states of molecular systems. However, so far it has only been applied to rather small molecules and in the commonly implemented diagonal approximations to the electronic self-energy, it depends on a mean-field starting point. We describe here an implementation of the self-consistent and starting-point-independent quasiparticle self-consistent (qsGW)-BSE approach, which is suitable for calculations on large molecules. We herein show that eigenvalue-only self-consistency can lead to an unfaithful description of some excitonic states for chlorophyll dimers while the qsGW-BSE vertical excitation energies (VEEs) are in excellent agreement with spectroscopic experiments for chlorophyll monomers and dimers measured in the gas phase. Furthermore, VEEs from time-dependent density functional theory calculations tend to disagree with experimental values and using different range-separated hybrid (RSH) kernels does change the VEEs by up to 0.5 eV. We use the new qsGW-BSE implementation to calculate the lowest excitation energies of the six chromophores of the photosystem II (PSII) reaction center (RC) with nearly 2000 correlated electrons. Using more than 11,000 (6000) basis functions, the calculation could be completed in less than 5 (2) days on a single modern compute node. In agreement with previous TD-DFT calculations using RSH kernels on models that also do not include environmental effects, our qsGW-BSE calculations only yield states with local characters in the low-energy spectrum of the hexameric complex. Earlier works with RSH kernels have demonstrated that the protein environment facilitates the experimentally observed interchromophoric charge transfer. Therefore, future research will need to combine correlation effects beyond TD-DFT with an explicit treatment of environmental electrostatics.

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.

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.

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.

In Vitro Enzymatic Activity Assays Implicate the Existence of the Chlorophyll Cycle in Chlorophyll b-Containing Cyanobacteria

In plants, chlorophyll (Chl) a and b are interconvertible by the action of three enzymes-chlorophyllide a oxygenase, Chl b reductase (CBR) and 7-hydroxymethyl chlorophyll a reductase (HCAR). These reactions are collectively referred to as the Chl cycle. In plants, this cyclic pathway ubiquitously exists and plays essential roles in acclimation to different light conditions at various developmental stages. By contrast, only a limited number of cyanobacteria species produce Chl b, and these include Prochlorococcus, Prochloron, Prochlorothrix and Acaryochloris. In this study, we investigated a possible existence of the Chl cycle in Chl b synthesizing cyanobacteria by testing in vitro enzymatic activities of CBR and HCAR homologs from Prochlorothrix hollandica and Acaryochloris RCC1774. All of these proteins show respective CBR and HCAR activity in vitro, indicating that both cyanobacteria possess the potential to complete the Chl cycle. It is also found that CBR and HCAR orthologs are distributed only in the Chl b-containing cyanobacteria that habitat shallow seas or freshwater, where light conditions change dynamically, whereas they are not found in Prochlorococcus species that usually habitat environments with fixed lighting. Taken together, our results implicate a possibility that the Chl cycle functions for light acclimation in Chl b-containing cyanobacteria.

On the interface of light-harvesting antenna complexes and reaction centers in oxygenic photosynthesis

Photosynthetic pigment-protein complexes (PPCs) accomplish light-energy capture and photochemistry in natural photosynthesis. In this review, we examine three pigment protein complexes in oxygenic photosynthesis: light-harvesting antenna complexes and two reaction centers: Photosystem II (PSII), and Photosystem I (PSI). Recent technological developments promise unprecedented insights into how these multi-component protein complexes are assembled into higher order structures and thereby execute their function. Furthermore, the interfacial domain between light-harvesting antenna complexes and PSII, especially the potential roles of the structural loops from CP29 and the PB-loop of ApcE in higher plant and cyanobacteria, respectively, are discussed. It is emphasized that the structural nuances are required for the structural dynamics and consequently for functional regulation in response to an ever-changing and challenging environment.

Large Synthetic Molecule that either Folds or Aggregates through Weak Supramolecular Interactions Determined by Solvent

Weak noncovalent interactions between large disclike molecules in poorly solvating media generally lead to the formation of fibers where the molecules stack atop one another. Here, we show that a particular chiral spacing group between large aromatic moieties, which usually lead to columnar stacks, in this case gives rise to an intramolecularly folded structure in relatively polar solvents, but in very apolar solvents forms finite aggregates. The molecule that displays this behavior has a C 3 symmetric benzene-1,3,5-tris(3,3'-diamido-2,2'-bipyridine) (BTAB) core with three metalloporphyrin units appended to it through short chiral spacers. Quite well-defined chromophore arrangements are evident by circular dichroism (CD) spectroscopy of this compound in solution, where clear exciton coupled bands of porphyrins are observed. In more polar solvents where the molecules are dispersed, a relatively weak CD signal is observed as a result of intramolecular folding, a feature confirmed by molecular modeling. The intramolecular folding was confirmed by measuring the CD of a C 2 symmetric analogue. The C 3 symmetric BTAB cores that would normally be expected to stack in a chiral arrangement in apolar solvents show no indication of CD, suggesting that there is no transfer of chirality through it (although the expected planar conformation of the 2,2'-bipyridine unit is confirmed by NMR spectroscopy). The incorporation of the porphyrins on the 3,3'-diamino-2,2'-bipyridine moiety spaced by a chiral unit leaves the latter incapable of assembling through supramolecular π-π stacking. Rather, modeling indicates that the three metalloporphyrin units interact, thanks to van der Waals interactions, favoring their close interactions over that of the BTAB units. Atomic force microscopy shows that, in contrast to other examples of molecules with the same core, disclike aggregates (rather than fibrillar one dimensional aggregates) are favored by the C 3 symmetric molecule. The closed structures are formed through nondirectional interlocking of porphyrin rings. The chiral spacer between the rigid core and the porphyrin moieties is undoubtedly important in determining the outcome in polar or less polar solvents, as modeling shows that this joint in the molecule has two favored conformations that render the molecule relatively flat or convex.

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.

Intramolecular strong electronic coupling in a discretely H-aggregated phthalocyanine dimer connected with a rigid linker

This is the first example of a discrete dimer of phthalocyanines in an H-aggregate form. Intramolecular strong electronic coupling was observed between the phthalocyanines electrochemically.

Two-Dimensional Spectroscopy of Chlorophyll a Excited-State Equilibration in Light-Harvesting Complex II

Excited-state relaxation dynamics and energy-transfer processes in the chlorophyll a (Chl a) manifold of the light-harvesting complex II (LHCII) were examined at physiological temperature using femtosecond two-dimensional electronic spectroscopy (2DES). The experiments were done under conditions free from singlet-singlet annihilation and anisotropic decay. Energy transfer between the different domains of the Chl a manifold was found to proceed on time scales from hundreds of femtoseconds to five picoseconds, before reaching equilibration. No component slower than 10 ps was observed in the spectral equilibration dynamics. We clearly observe the bidirectional (uphill and downhill) energy transfer of the equilibration process between excited states. This bidirectional energy flow, although implicit in the modeling and simulation of the EET processes, has not been observed in any prior transient absorption studies. Furthermore, we identified the spectral forms associated with the different energy transfer lifetimes in the equilibration process.

Thylakoid membrane model of the Chl a fluorescence transient and P700 induction kinetics in plant leaves

A new Thylakoid model is presented, which describes in detail the electron/proton transfer reactions between membrane protein complexes including photosystems II and I (PSII, PSI), cytochrome (Cyt) b 6 f, mobile plastoquinone PQ pool in the thylakoid membrane, plastocyanin in lumen and ferredoxin in stroma, reduction of NADP via FNR and cyclic electron transfer. The Thylakoid model parameters were fitted both to Chl fluorescence induction data (FI) and oxido-reductions of P700 (ΔA 810) measured from 20 μs up to 20 s in pea leaves. The two-wave kinetics of FI and ΔA 810 (O(JI)PSM and OABCDE) were described quantitatively, provided that the values of membrane electrochemical potential components ΔΨ(t), pHL(t)/pHS(t) are in physiologically relevant ranges. The time courses on the time scale from nanoseconds to tens of seconds of oxido-reduction changes of ET components as well as concentrations of proton/ions (K+, Cl-) were calculated. We assume a low constant FNR activity over this period. Charge movements across the thylakoid membrane by passive leakage and active ATPase transport and proton buffer reactions are simulated. The dynamics of charge fluxes during photosynthetic induction under low light (PFD 200 μmol photons m-2 s-1) were analyzed. The initial wave of P700 oxidation within 20 ms during independent operation of PSI and PSII was followed after 50 ms by PSI donor-side reduction from reduced PQ pool via Cyt b 6 f site. The Cyt b 6 f reactions contribute to the stabilization of fluxes in the time range 1 s < t < 10 s. The detailed analysis of Chl a fluorescence at the PSM stage (t > 10 s) would need the investigation of FNR activation effect in order to explain the transitions between cyclic and linear electron transport.

Multichromophoric π-Conjugation: Modular Design for Gated and Cascade Energy Transfer

Multichromophore arrays allow for cascade energy transfer. As an isoelectronic analogue of indacenyl, bis(triazolo)benzene features a fused tricyclic skeleton that rigidly places two π-extended triazoles in close proximity. Such triazole-based fluorophores behave as electronically independent modules in the ground states, but become tightly coupled upon photoexcitation for highly efficient excitation energy transfer (EET) that can be gated by external stimuli. Taking this donor-acceptor fluorophore system a step further, we have designed and implemented a cascade EET. Here, the initial excitation takes part in a circular relay to arrive at the longest-wavelength emitting site as the final destination. Modularly constructed triazoloarenes should serve as versatile platforms for chemically controlled optical signaling.

Quantum Chemical Studies of Light Harvesting

The design of optimal light-harvesting (supra)molecular systems and materials is one of the most challenging frontiers of science. Theoretical methods and computational models play a fundamental role in this difficult task, as they allow the establishment of structural blueprints inspired by natural photosynthetic organisms that can be applied to the design of novel artificial light-harvesting devices. Among theoretical strategies, the application of quantum chemical tools represents an important reality that has already reached an evident degree of maturity, although it still has to show its real potentials. This Review presents an overview of the state of the art of this strategy, showing the actual fields of applicability but also indicating its current limitations, which need to be solved in future developments.

Multiscale model of light harvesting by photosystem II in plants

The first step of photosynthesis in plants is the absorption of sunlight by pigments in the antenna complexes of photosystem II (PSII), followed by transfer of the nascent excitation energy to the reaction centers, where long-term storage as chemical energy is initiated. Quantum mechanical mechanisms must be invoked to explain the transport of excitation within individual antenna. However, it is unclear how these mechanisms influence transfer across assemblies of antenna and thus the photochemical yield at reaction centers in the functional thylakoid membrane. Here, we model light harvesting at the several-hundred-nanometer scale of the PSII membrane, while preserving the dominant quantum effects previously observed in individual complexes. We show that excitation moves diffusively through the antenna with a diffusion length of 50 nm until it reaches a reaction center, where charge separation serves as an energetic trap. The diffusion length is a single parameter that incorporates the enhancing effect of excited state delocalization on individual rates of energy transfer as well as the complex kinetics that arise due to energy transfer and loss by decay to the ground state. The diffusion length determines PSII's high quantum efficiency in ideal conditions, as well as how it is altered by the membrane morphology and the closure of reaction centers. We anticipate that the model will be useful in resolving the nonphotochemical quenching mechanisms that PSII employs in conditions of high light stress.

Bottom-Up Hierarchical Self-Assembly of Chiral Porphyrins through Coordination and Hydrogen Bonds

A series of chiral synthetic compounds is reported that shows intricate but specific hierarchical assembly because of varying positions of coordination and hydrogen bonds. The evolution of the aggregates (followed by absorption spectroscopy and temperature-dependent circular dichroism studies in solution) reveal the influence of the proportion of stereogenic centers in the side groups connected to the chromophore ring in their optical activity and the important role of pyridyl groups in the self-assembly of these chiral macrocycles. The optical activity spans 2 orders of magnitude depending on composition and constitution. Two of the aggregates show very high optical activity even though the isolated chromophores barely give a circular dichroism signal. Molecular modeling of the aggregates, starting from the pyridine-zinc(II) porphyrin interaction and working up, and calculation of the circular dichroism signal confirm the origin of this optical activity as the chiral supramolecular organization of the molecules. The aggregates show a broad absorption range, between approximately 390 and 475 nm for the transitions associated with the Soret region alone, that spans wavelengths far more than the isolated chromophore. The supramolecular assemblies of the metalloporphyrins in solution were deposited onto highly oriented pyrolitic graphite in order to study their hierarchy in assembly by atomic force microscopy. Zero and one-dimensional aggregates were observed, and a clear dependence on deposition temperature was shown, indicating that the hierarchical assembly took place largely in solution. Moreover, scanning electron microscopy images of porphyrins and metalloporphyrins precipitated under out-of-equilibrium conditions showed the dependence of the number and position of chiral amide groups in the formation of a fibrillar nanomaterial. The combination of coordination and hydrogen bonding in the complicated assembly of these molecules-where there is a clear hierarchy for zinc(II)-pyridyl interaction followed by hydrogen-bonding between amide groups, and then van der Waals interactions-paves the way for the preparation of molecular materials with multiple chromophore environments.

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.

Critical assessment of the emission spectra of various photosystem II core complexes

We evaluate low-temperature (low-T) emission spectra of photosystem II core complexes (PSII-cc) previously reported in the literature, which are compared with emission spectra of PSII-cc obtained in this work from spinach and for dissolved PSII crystals from Thermosynechococcus (T.) elongatus. This new spectral dataset is used to interpret data published on membrane PSII (PSII-m) fragments from spinach and Chlamydomonas reinhardtii, as well as PSII-cc from T. vulcanus and intentionally damaged PSII-cc from spinach. This study offers new insight into the assignment of emission spectra reported on PSII-cc from different organisms. Previously reported spectra are also compared with data obtained at different saturation levels of the lowest energy state(s) of spinach and T. elongatus PSII-cc via hole burning in order to provide more insight into emission from bleached and/or photodamaged complexes. We show that typical low-T emission spectra of PSII-cc (with closed RCs), in addition to the 695 nm fluorescence band assigned to the intact CP47 complex (Reppert et al. J Phys Chem B 114:11884-11898, 2010), can be contributed to by several emission bands, depending on sample quality. Possible contributions include (i) a band near 690-691 nm that is largely reversible upon temperature annealing, proving that the band originates from CP47 with a bleached low-energy state near 693 nm (Neupane et al. J Am Chem Soc 132:4214-4229, 2010; Reppert et al. J Phys Chem B 114:11884-11898, 2010); (ii) CP43 emission at 683.3 nm (not at 685 nm, i.e., the F685 band, as reported in the literature) (Dang et al. J Phys Chem B 112:9921-9933, 2008; Reppert et al. J Phys Chem B 112:9934-9947, 2008); (iii) trap emission from destabilized CP47 complexes near 691 nm (FT1) and 685 nm (FT2) (Neupane et al. J Am Chem Soc 132:4214-4229, 2010); and (iv) emission from the RC pigments near 686-687 nm. We suggest that recently reported emission of single PSII-cc complexes from T. elongatus may not represent intact complexes, while those obtained for T. elongatus presented in this work most likely represent intact PSII-cc, since they are nearly indistinguishable from emission spectra obtained for various PSII-m fragments.

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

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

Primary radical ion pairs in photosystem II core complexes

Ultrafast absorption spectroscopy with 20-fs resolution was applied to study primary charge separation in spinach photosystem II (PSII) reaction center (RC) and PSII core complex (RC complex with integral antenna) upon excitation at maximum wavelength 700-710 nm at 278 K. It was found that the initial charge separation between P680* and ChlD1 (Chl-670) takes place with a time constant of ~1 ps with the formation of the primary charge-separated state P680* with an admixture of: P680*((1-δ)) (P680(δ+)ChlD1(δ-)), where δ ~ 0.5. The subsequent electron transfer from P680(δ+)ChlD1(δ-) to pheophytin (Pheo) occurs within 13 ps and is accompanied by a relaxation of the absorption band at 670 nm (ChlD1(δ-)) and bleaching of the PheoD1 bands at 420, 545, and 680 nm with development of the Pheo(-) band at 460 nm. Further electron transfer to QA occurs within 250 ps in accordance with earlier data. The spectra of P680(+) and Pheo(-) formation include a bleaching band at 670 nm; this indicates that Chl-670 is an intermediate between P680 and Pheo. Stimulated emission kinetics at 685 nm demonstrate the existence of two decaying components with time constants of ~1 and ~13 ps due to the formation of P680(δ+)ChlD1(δ-) and P680(+)PheoD1(-), respectively.

Modeling of fluorescence line-narrowed spectra in weakly coupled dimers in the presence of excitation energy transfer

This work describes simple analytical formulas to describe the fluorescence line-narrowed (FLN) spectra of weakly coupled chromophores in the presence of excitation energy transfer (EET). Modeling studies for dimer systems (assuming low fluence and weak coupling) show that the FLN spectra (including absorption and emission spectra) calculated for various dimers using our model are in good agreement with spectra calculated by: (i) the simple convolution method and (ii) the more rigorous treatment using the Redfield approach [T. Renger and R. A. Marcus, J. Chem. Phys. 116, 9997 (2002)]. The calculated FLN spectra in the presence of EET of all three approaches are very similar. We argue that our approach provides a simplified and computationally more efficient description of FLN spectra in the presence of EET. This method also has been applied to FLN spectra obtained for the CP47 antenna complex of Photosystem II reported by Neupane et al. [J. Am. Chem. Soc. 132, 4214 (2010)], which indicated the presence of uncorrelated EET between pigments contributing to the two lowest energy (overlapping) exciton states, each mostly localized on a single chromophore. Calculated and experimental FLN spectra for CP47 complex show very good qualitative agreement.

Ultrafast infrared spectroscopy in photosynthesis

In recent years visible pump/mid-infrared (IR) probe spectroscopy has established itself as a key technology to unravel structure-function relationships underlying the photo-dynamics of complex molecular systems. In this contribution we review the most important applications of mid-infrared absorption difference spectroscopy with sub-picosecond time-resolution to photosynthetic complexes. Considering several examples, such as energy transfer in photosynthetic antennas and electron transfer in reaction centers and even more intact structures, we show that the acquisition of ultrafast time resolved mid-IR spectra has led to new insights into the photo-dynamics of the considered systems and allows establishing a direct link between dynamics and structure, further strengthened by the possibility of investigating the protein response signal to the energy or electron transfer processes. This article is part of a Special Issue entitled: Vibrational spectroscopies and bioenergetic systems.

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

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

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

Formation and decay of P680 (P(D1)-P(D2))⁺PheoD1⁻ radical ion pair in photosystem II core complexes

Under physiological conditions (278 K) femtosecond pump-probe laser spectroscopy with 20-fs time resolution was applied to study primary charge separation in spinach photosystem II (PSII) core complexes excited at 710 nm. It was shown that initial formation of anion radical band of pheophytin molecule (Pheo⁻) at 460 nm is observed with rise time of ~11ps. The kinetics of the observed rise was ascribed to charge separation between Chl (chlorophyll a) dimer, primary electron donor in PSII (P680*) and Pheo located in D1 protein subunit (PheoD1) absorbing at 420 nm, 545 nm and 680 nm with formation of the ion-radical pair P680⁺PheoDI⁻. The subsequent electron transfer from Pheo(D1)⁻ to primary plastoquinone electron acceptor (Q(A)) was accompanied by relaxation of the 460-nm band and occurred within ~250 ps in good agreement with previous measurements in Photosystem II-enriched particles and bacterial reaction centers. The subtraction of the P680⁺ spectrum measured at 455 ps delay from the spectra at 23 ps or 44 ps delay reveals the spectrum of Pheo(DI)⁻, which is very similar to that measured earlier by accumulation method. The spectrum of Pheo(DI)⁻ formation includes a bleaching (or red shift) of the 670 nm band indicating that Chl-670 is close to Pheo(D1). According to previous measurements in the femtosecond-picosecond time range this Chl-670 was ascribed to Chl(D1) [Shelaev, Gostev, Vishnev, Shkuropatov, Ptushenko, Mamedov, Sarkisov, Nadtochenko, Semenov and Shuvalov, J. Photochemistry and Photobiology, B: Biology 104 (2011) 45-50]. Stimulated emission at 685 nm was found to have two decaying components with time constants of ~1ps and ~14ps. These components appear to reflect formation of P680⁺Chl(D1)⁻ and P680⁺Pheo(D1)⁻, respectively, as found earlier. This article is part of a special issue entitled: photosynthesis research for sustainability: keys to produce clean energy.

Light harvesting in photosystem II

Water oxidation in photosynthesis takes place in photosystem II (PSII). This photosystem is built around a reaction center (RC) where sunlight-induced charge separation occurs. This RC consists of various polypeptides that bind only a few chromophores or pigments, next to several other cofactors. It can handle far more photons than the ones absorbed by its own pigments and therefore, additional excitations are provided by the surrounding light-harvesting complexes or antennae. The RC is located in the PSII core that also contains the inner light-harvesting complexes CP43 and CP47, harboring 13 and 16 chlorophyll pigments, respectively. The core is surrounded by outer light-harvesting complexes (Lhcs), together forming the so-called supercomplexes, at least in plants. These PSII supercomplexes are complemented by some "extra" Lhcs, but their exact location in the thylakoid membrane is unknown. The whole system consists of many subunits and appears to be modular, i.e., both its composition and organization depend on environmental conditions, especially on the quality and intensity of the light. In this review, we will provide a short overview of the relation between the structure and organization of pigment-protein complexes in PSII, ranging from individual complexes to entire membranes and experimental and theoretical results on excitation energy transfer and charge separation. It will become clear that time-resolved fluorescence data can provide invaluable information about the organization and functioning of thylakoid membranes. At the end, an overview will be given of unanswered questions that should be addressed in the near future.

Selective and differential optical spectroscopies in photosynthesis

Photosynthetic pigments are inherently intense optical absorbers and have strong polarisation characteristics. They can also luminesce strongly. These properties have led optical spectroscopies to be, quite naturally, key techniques in photosynthesis. However, there are typically many pigments in a photosynthetic assembly, which when combined with the very significant inhomogeneous and homogeneous linewidths characteristic of optical transitions, leads to spectral congestion. This in turn has made it difficult to provide a definitive and detailed electronic structure for many photosynthetic assemblies. An electronic structure is, however, necessary to provide a foundation for any complete description of fundamental processes in photosynthesis, particularly those in reaction centres. A wide range of selective and differential spectral techniques have been developed to help overcome the problems of spectral complexity and congestion. The techniques can serve to either reduce spectral linewidths and/or extract chromophore specific information from unresolved spectral features. Complementary spectral datasets, generated by a number of techniques, may then be combined in a 'multi-dimensional' theoretical analysis so as to constrain and define effective models of photosynthetic assemblies and their fundamental processes. A key example is the work of Renger and his group (Raszewski, Biophys J 88(2):986-998, 2005) on PS II reaction centre assemblies. This article looks to provide an overview of some of these techniques and indicate where their strengths and weaknesses may lie. It highlights some of our own contributions and indicates areas where progress may be possible.

Structure-based modeling of energy transfer in photosynthesis

We provide a minimal model for a structure-based simulation of excitation energy transfer in pigment-protein complexes (PPCs). In our treatment, the PPC is assembled from its building blocks. The latter are defined such that electron exchange occurs only within, but not between these units. The variational principle is applied to investigate how the Coulomb interaction between building blocks changes the character of the electronic states of the PPC. In this way, the standard exciton Hamiltonian is obtained from first principles and a hierarchy of calculation schemes for the parameters of this Hamiltonian arises. Possible extensions of this approach are discussed concerning (i) the inclusion of dispersive site energy shifts and (ii) the inclusion of electron exchange between pigments. First results on electron exchange within the special pair of photosystem II of cyanobacteria and higher plants are presented and compared with earlier results on purple bacteria. In the last part of this mini-review, the coupling of electronic and nuclear degrees of freedom is considered. First, the standard exciton-vibrational Hamiltonian is parameterized with the help of a normal mode analysis of the PPC. Second, dynamical theories are discussed that exploit this Hamiltonian in the study of dissipative exciton motion.

Functional compartmental modeling of the photosystems in the thylakoid membrane at 77 K

Time-resolved fluorescence spectroscopy measurements at 77 K on thylakoid membrane preparations and isolated photosynthetic complexes thereof were investigated using target analysis with the aim of building functional compartmental models for the photosystems in the thylakoid membrane. Combining kinetic schemes with different spectral constraints enabled us to resolve the energy transfer pathways and decay characteristics of the different emissive species. We determined the spectral and energetic properties of the red Chl pools in both photosystems and quantified the formation of LHCII-LHCI-PSI supercomplexes in the transition from native to unstacked thylakoid membranes.

A structure-based model of energy transfer reveals the principles of light harvesting in photosystem II supercomplexes

Photosystem II (PSII) initiates photosynthesis in plants through the absorption of light and subsequent conversion of excitation energy to chemical energy via charge separation. The pigment binding proteins associated with PSII assemble in the grana membrane into PSII supercomplexes and surrounding light harvesting complex II trimers. To understand the high efficiency of light harvesting in PSII requires quantitative insight into energy transfer and charge separation in PSII supercomplexes. We have constructed the first structure-based model of energy transfer in PSII supercomplexes. This model shows that the kinetics of light harvesting cannot be simplified to a single rate limiting step. Instead, substantial contributions arise from both excitation diffusion through the antenna pigments and transfer from the antenna to the reaction center (RC), where charge separation occurs. Because of the lack of a rate-limiting step, fitting kinetic models to fluorescence lifetime data cannot be used to derive mechanistic insight on light harvesting in PSII. This model will clarify the interpretation of chlorophyll fluorescence data from PSII supercomplexes, grana membranes, and leaves.

Evolutionary changes in chlorophyllide a oxygenase (CAO) structure contribute to the acquisition of a new light-harvesting complex in micromonas

Chlorophyll b is found in photosynthetic prokaryotes and primary and secondary endosymbionts, although their light-harvesting systems are quite different. Chlorophyll b is synthesized from chlorophyll a by chlorophyllide a oxygenase (CAO), which is a Rieske-mononuclear iron oxygenase. Comparison of the amino acid sequences of CAO among photosynthetic organisms elucidated changes in the domain structures of CAO during evolution. However, the evolutionary relationship between the light-harvesting system and the domain structure of CAO remains unclear. To elucidate this relationship, we investigated the CAO structure and the pigment composition of chlorophyll-protein complexes in the prasinophyte Micromonas. The Micromonas CAO is composed of two genes, MpCAO1 and MpCAO2, that possess Rieske and mononuclear iron-binding motifs, respectively. Only when both genes were introduced into the chlorophyll b-less Arabidopsis mutant (ch1-1) was chlorophyll b accumulated, indicating that cooperation between the two subunits is required to synthesize chlorophyll b. Although Micromonas has a characteristic light-harvesting system in which chlorophyll b is incorporated into the core antennas of reaction centers, chlorophyll b was also incorporated into the core antennas of reaction centers of the Arabidopsis transformants that contained the two Micromonas CAO proteins. Based on these results, we discuss the evolutionary relationship between the structures of CAO and light-harvesting systems.

Photosystem II does not possess a simple excitation energy funnel: time-resolved fluorescence spectroscopy meets theory

The experimentally obtained time-resolved fluorescence spectra of photosystem II (PS II) core complexes, purified from a thermophilic cyanobacterium Thermosynechococcus vulcanus, at 5–180 K are compared with simulations. Dynamic localization effects of excitons are treated implicitly by introducing exciton domains of strongly coupled pigments. Exciton relaxations within a domain and exciton transfers between domains are treated on the basis of Redfield theory and generalized Förster theory, respectively. The excitonic couplings between the pigments are calculated by a quantum chemical/electrostatic method (Poisson-TrEsp). Starting with previously published values, a refined set of site energies of the pigments is obtained through optimization cycles of the fits of stationary optical spectra of PS II. Satisfactorily agreement between the experimental and simulated spectra is obtained for the absorption spectrum including its temperature dependence and the linear dichroism spectrum of PS II core complexes (PS II-CC). Furthermore, the refined site energies well reproduce the temperature dependence of the time-resolved fluorescence spectrum of PS II-CC, which is characterized by the emergence of a 695 nm fluorescence peak upon cooling down to 77 K and the decrease of its relative intensity upon further cooling below 77 K. The blue shift of the fluorescence band upon cooling below 77 K is explained by the existence of two red-shifted chlorophyll pools emitting at around 685 and 695 nm. The former pool is assigned to Chl45 or Chl43 in CP43 (Chl numbering according to the nomenclature of Loll et al. Nature2005, 438, 1040) while the latter is assigned to Chl29 in CP47. The 695 nm emitting chlorophyll is suggested to attract excitations from the peripheral light-harvesting complexes and might also be involved in photoprotection.

Photosynthetic reaction center as a quantum heat engine

Two seemingly unrelated effects attributed to quantum coherence have been reported recently in natural and artificial light-harvesting systems. First, an enhanced solar cell efficiency was predicted and second, population oscillations were measured in photosynthetic antennae excited by sequences of coherent ultrashort laser pulses. Because both systems operate as quantum heat engines (QHEs) that convert the solar photon energy to useful work (electric currents or chemical energy, respectively), the question arises whether coherence could also enhance the photosynthetic yield. Here, we show that both effects arise from the same population-coherence coupling term which is induced by noise, does not require coherent light, and will therefore work for incoherent excitation under natural conditions of solar excitation. Charge separation in light-harvesting complexes occurs in a pair of tightly coupled chlorophylls (the special pair) at the heart of photosynthetic reaction centers of both plants and bacteria. We show the analogy between the energy level schemes of the special pair and of the laser/photocell QHEs, and that both population oscillations and enhanced yield have a common origin and are expected to coexist for typical parameters. We predict an enhanced yield of 27% in a QHE motivated by the reaction center. This suggests nature-mimicking architectures for artificial solar energy devices.

Understanding photosynthetic light-harvesting: a bottom up theoretical approach

We discuss a bottom up approach for modeling photosynthetic light-harvesting. Methods are reviewed for a full structure-based parameterization of the Hamiltonian of pigment-protein complexes (PPCs). These parameters comprise (i) the local transition energies of the pigments in their binding sites in the protein, the site energies; (ii) the couplings between optical transitions of the pigments, the excitonic couplings; and (iii) the spectral density characterizing the dynamic modulation of pigment transition energies and excitonic couplings by protein vibrations. Starting with quantum mechanics perturbation theory, we provide a microscopic foundation for the standard PPC Hamiltonian and relate the expressions obtained for its matrix elements to quantities that can be calculated with classical molecular mechanics/electrostatics approaches including the whole PPC in atomic detail and using charge and transition densities obtained with quantum chemical calculations on the isolated building blocks of the PPC. In the second part of this perspective, the Hamiltonian is utilized to describe the quantum dynamics of excitons. Situations are discussed that differ in the relative strength of excitonic and exciton-vibrational coupling. The predictive power of the approaches is demonstrated in application to different PPCs, and challenges for future work are outlined.

Simulations of the two-dimensional electronic spectroscopy of the photosystem II reaction center

We report simulations of the two-dimensional electronic spectroscopy of the Q(y) band of the D1-D2-Cyt b559 photosystem II reaction center at 77 K. We base the simulations on an existing Hamiltonian that was derived by simultaneous fitting to a wide range of linear spectroscopic measurements and described within modified Redfield theory. The model obtains reasonable agreement with most aspects of the two-dimensional spectra, including the overall peak shapes and excited state absorption features. It does not reproduce the rapid equilibration from high energy to low energy excitonic states evident by a strong cross-peak below the diagonal. We explore modifications to the model to incorporate new structural data and improve agreement with the two-dimensional spectra. We find that strengthening the system-bath coupling and lowering the degree of disorder significantly improves agreement with the cross-peak feature, while lessening agreement with the relative diagonal/antidiagonal width of the 2D spectra. We conclude that two-dimensional electronic spectroscopy provides a sensitive test of excitonic models of the photosystem II reaction center and discuss avenues for further refinement of such models.

Chlorophyll a fluorescence induction: a personal perspective of the thermal phase, the J-I-P rise

The fast (up to 1 s) chlorophyll (Chl) a fluorescence induction (FI) curve, measured under saturating continuous light, has a photochemical phase, the O-J rise, related mainly to the reduction of Q(A), the primary electron acceptor plastoquinone of Photosystem II (PSII); here, the fluorescence rise depends strongly on the number of photons absorbed. This is followed by a thermal phase, the J-I-P rise, which disappears at subfreezing temperatures. According to the mainstream interpretation of the fast FI, the variable fluorescence originates from PSII antenna, and the oxidized Q(A) is the most important quencher influencing the O-J-I-P curve. As the reaction centers of PSII are gradually closed by the photochemical reduction of Q(A), Chl fluorescence, F, rises from the O level (the minimal level) to the P level (the peak); yet, the relationship between F and [Q(A) (-)] is not linear, due to the presence of other quenchers and modifiers. Several alternative theories have been proposed, which give different interpretations of the O-J-I-P transient. The main idea in these alternative theories is that in saturating light, Q(A) is almost completely reduced already at the end of the photochemical phase O-J, but the fluorescence yield is lower than its maximum value due to the presence of either a second quencher besides Q(A), or there is an another process quenching the fluorescence; in the second quencher hypothesis, this quencher is consumed (or the process of quenching the fluorescence is reversed) during the thermal phase J-I-P. In this review, we discuss these theories. Based on our critical examination, that includes pros and cons of each theory, as well mathematical modeling, we conclude that the mainstream interpretation of the O-J-I-P transient is the most credible one, as none of the alternative ideas provide adequate explanation or experimental proof for the almost complete reduction of Q(A) at the end of the O-J phase, and for the origin of the fluorescence rise during the thermal phase. However, we suggest that some of the factors influencing the fluorescence yield that have been proposed in these newer theories, as e.g., the membrane potential ΔΨ, as suggested by Vredenberg and his associates, can potentially contribute to modulate the O-J-I-P transient in parallel with the reduction of Q(A), through changes at the PSII antenna and/or at the reaction center, or, possibly, through the control of the oxidation-reduction of the PQ-pool, including proton transfer into the lumen, as suggested by Rubin and his associates. We present in this review our personal perspective mainly on our understanding of the thermal phase, the J-I-P rise during Chl a FI in plants and algae.

Probing Electron-Transfer Times in Photosynthetic Reaction Centers by Hole-Burning Spectroscopy

A brief discussion is presented of transient hole-burned (HB) spectra (and the information that they provide) obtained for isolated reaction centers (RCs) from wild-type (WT) Rhodobacter sphaeroides, RCs containing zinc-bacteriochlorophylls (Zn-BChls), and RCs of Photosystem II (PSII) from spinach and Chlamydomonas reinhardtii . The shape of the spectral density and the strength of electron-phonon coupling in bacterial RCs are discussed. We focus, however, on heterogeneity of isolated PS II RCs from spinach and, in particular, Chlamydomonas reinhardtii , site energies of active (electron acceptor) and inactive pheophytins, the nature of the primary electron donor(s), and the possibility of multiple charge-separation (CS) pathways in the isolated PSII RC. We conclude with comments on current efforts in HB spectroscopy in the area of photosynthesis and future directions in HB spectroscopy.

Primary electron donor(s) in isolated reaction center of photosystem II from Chlamydomonas reinhardtii

Isolated reaction centers (RCs) from wild-type Chlamydomonas (C.) reinhardtii of Photosystem II (PSII), at different levels of intactness, were studied to provide more insight into the nature of the charge-separation (CS) pathway(s). We argue that previously studied D1/D2/Cytb559 complexes (referred to as RC680), with ChlD1 serving as the primary electron donor, contain destabilized D1 and D2 polypeptides and, as a result, do not provide a representative model system for the intact RC within the PSII core. The shapes of nonresonant transient hole-burned (HB) spectra obtained for more intact RCs (referred to as RC684) are very similar to P(+)QA(-) - PQA absorbance difference and triplet minus singlet spectra measured in PSII core complexes from Synechocystis PCC 6803 [Schlodder et al. Philos. Trans. R. Soc. London, Ser. B2008, 363, 1197]. We show that in the RC684 complexes, both PD1 and ChlD1 may serve as primary electron donors, leading to two different charge separation pathways. Resonant HB spectra cannot distinguish the CS times corresponding to different paths, but it is likely that the zero-phonon holes (ZPHs) observed in the 680-685 nm region (corresponding to CS times of ∼1.4-4.4 ps) reveal the ChlD1 pathway; conversely, the observation of charge-transfer (CT) state(s) in RC684 (in the 686-695 nm range) and the absence of ZPHs at λB > 685 nm likely stem from the PD1 pathway, for which CS could be faster than 1 ps. This is consistent with the finding of Krausz et al. [Photochem. Photobiol. Sci.2005, 4, 744] that CS in intact PSII core complexes can be initiated at low temperatures with fairly long-wavelength excitation. The lack of a clear shift of HB spectra as a function of excitation wavelength within the red-tail of the absorption (i.e., 686-695 nm) and the absence of ZPHs suggest that the lowest-energy CT state is largely homogeneously broadened. On the other hand, in usually studied destabilized RCs, that is, RC680, for which CT states have never been experimentally observed, ChlD1 is the most likely electron donor.

Site energies of active and inactive pheophytins in the reaction center of Photosystem II from Chlamydomonas reinhardtii

It is widely accepted that the primary electron acceptor in various Photosystem II (PSII) reaction center (RC) preparations is pheophytin a (Pheo a) within the D1 protein (Pheo(D1)), while Pheo(D2) (within the D2 protein) is photochemically inactive. The Pheo site energies, however, have remained elusive, due to inherent spectral congestion. While most researchers over the past two decades placed the Q(y)-states of Pheo(D1) and Pheo(D2) bands near 678-684 and 668-672 nm, respectively, recent modeling [Raszewski et al. Biophys. J. 2005, 88, 986 - 998; Cox et al. J. Phys. Chem. B 2009, 113, 12364 - 12374] of the electronic structure of the PSII RC reversed the assignment of the active and inactive Pheos, suggesting that the mean site energy of Pheo(D1) is near 672 nm, whereas Pheo(D2) (~677.5 nm) and Chl(D1) (~680 nm) have the lowest energies (i.e., the Pheo(D2)-dominated exciton is the lowest excited state). In contrast, chemical pigment exchange experiments on isolated RCs suggested that both pheophytins have their Q(y) absorption maxima at 676-680 nm [Germano et al. Biochemistry 2001, 40, 11472 - 11482; Germano et al. Biophys. J. 2004, 86, 1664 - 1672]. To provide more insight into the site energies of both Pheo(D1) and Pheo(D2) (including the corresponding Q(x) transitions, which are often claimed to be degenerate at 543 nm) and to attest that the above two assignments are most likely incorrect, we studied a large number of isolated RC preparations from spinach and wild-type Chlamydomonas reinhardtii (at different levels of intactness) as well as the Chlamydomonas reinhardtii mutant (D2-L209H), in which the active branch Pheo(D1) is genetically replaced with chlorophyll a (Chl a). We show that the Q(x)-/Q(y)-region site energies of Pheo(D1) and Pheo(D2) are ~545/680 nm and ~541.5/670 nm, respectively, in good agreement with our previous assignment [Jankowiak et al. J. Phys. Chem. B 2002, 106, 8803 - 8814]. The latter values should be used to model excitonic structure and excitation energy transfer dynamics of the PSII RCs.