Low-Energy Chlorophyll States in the CP43 Antenna Protein Complex: Simulation of Various Optical Spectra. II

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

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

Energy dissipation efficiency in the CP43 assembly intermediate complex of photosystem II

Photosystem II in oxygenic organisms is a large membrane bound rapidly turning over pigment protein complex. During its biogenesis, multiple assembly intermediates are formed, including the CP43-preassembly complex (pCP43). To understand the energy transfer dynamics in pCP43, we first engineered a His-tagged version of the CP43 in a CP47-less strain of the cyanobacterium Synechocystis 6803. Isolated pCP43 from this engineered strain was subjected to advanced spectroscopic analysis to evaluate its excitation energy dissipation characteristics. These included measurements of steady-state absorption and fluorescence emission spectra for which correlation was tested with Stepanov relation. Comparison of fluorescence excitation and absorptance spectra determined that efficiency of energy transfer from β-carotene to chlorophyll a is 39 %. Time-resolved fluorescence images of pCP43-bound Chl a were recorded on streak camera, and fluorescence decay dynamics were evaluated with global fitting. These demonstrated that the decay kinetics strongly depends on temperature and buffer used to disperse the protein sample and fluorescence decay lifetime was estimated in 3.2-5.7 ns time range, depending on conditions. The pCP43 complex was also investigated with femtosecond and nanosecond time-resolved absorption spectroscopy upon excitation of Chl a and β-carotene to reveal pathways of singlet excitation relaxation/decay, Chl a triplet dynamics and Chl a → β-carotene triplet state sensitization process. The latter demonstrated that Chl a triplet in the pCP43 complex is not efficiently quenched by carotenoids. Finally, detailed kinetic analysis of the rise of the population of β-carotene triplets determined that the time constant of the carotenoid triplet sensitization is 40 ns.

Frequency-Domain Spectroscopic Study of the Photosystem I Supercomplexes, Isolated IsiA Monomers, and the Intact IsiA Ring

The PSI₃-IsiA₁₈ supercomplex is one of the largest and most complicated assemblies in photosynthesis. The IsiA ring, composed of 18 IsiA monomers (IsiA₁₈) surrounding the PSI trimer (PSI₃), forms under iron-deficient conditions in cyanobacteria and acts as a peripheral antenna. Based on the supercomplex structure recently determined via cryo-EM imaging, we model various optical spectra of the IsiA monomers and IsiA₁₈ ring. Comparison of the absorption and emission spectra of the isolated IsiA monomers and the full ring reveals that about 2.7 chlorophylls (Chls) are lost in the isolated IsiA monomers. The best fits for isolated monomers spectra are obtained assuming the absence of Chl 508 and Chl 517 and 70% loss of Chl 511. The best model describing all three hexamers and the entire ring suggests that the lowest energy pigments are Chls 511, 514, and 517. Based on the modeling results presented in this work, we conclude that there are most likely three entry points for EET from the IsiA₆ hexamer to the PSI core monomer, with two of these entry points likely being located next to each other (i.e., nine entry points from IsiA₁₈ to the PSI₃ trimer). Finally, we show that excitation energy transfer inside individual monomers is fast (<2 ps at T = 5 K) and at least 20 times faster than intermonomer energy transfer.

Response mechanisms to ocean warming exposure in Effrenium voratum (Symbiodiniaceae)

Ocean warming is an extreme environment event that has profound and lasting impacts on Symbiodiniaceae. However, their response mechanisms to elevated temperature exposure are poorly understood. In this study, the physiological and transcriptional responses of Effrenium voratum (Symbiodiniaceae) to ocean warming were examined. After exposure to 30 °C, no significant variations in growth, chlorophyll a, or photosynthetic and respiration rates were observed, while a higher temperature (34 °C) significantly reduced these physiological measurements. Meanwhile, lipid content and fatty acid composition were altered at high temperature (i.e., elevated degree of fatty acid saturation). Such biochemical constituents likely contributed to the mitigation of the negative effects of elevated temperatures. Furthermore, higher expression levels of genes related to the synthesis and elongation of fatty acids were detected at high temperature. The adjustment of lipids and fatty acid composition may be a potential mechanism by which E. voratum may survive under future global warming. ONE SENTENCE SUMMARY: The adjustment of lipids and fatty acid composition may be a potential mechanism by which E. voratum acclimate to future global warming.

Spectral densities and absorption spectra of the core antenna complex CP43 from photosystem II

Besides absorbing light, the core antenna complex CP43 of photosystem II is of great importance in transferring excitation energy from the antenna complexes to the reaction center. Excitation energies, spectral densities, and linear absorption spectra of the complex have been evaluated by a multiscale approach. In this scheme, quantum mechanics/molecular mechanics molecular dynamics simulations are performed employing the parameterized density functional tight binding (DFTB) while the time-dependent long-range-corrected DFTB scheme is applied for the excited state calculations. The obtained average spectral density of the CP43 complex shows a very good agreement with experimental results. Moreover, the excitonic Hamiltonian of the system along with the computed site-dependent spectral densities was used to determine the linear absorption. While a Redfield-like approximation has severe shortcomings in dealing with the CP43 complex due to quasi-degenerate states, the non-Markovian full second-order cumulant expansion formalism is able to overcome the drawbacks. Linear absorption spectra were obtained, which show a good agreement with the experimental counterparts at different temperatures. This study once more emphasizes that by combining diverse techniques from the areas of molecular dynamics simulations, quantum chemistry, and open quantum systems, it is possible to obtain first-principle results for photosynthetic complexes, which are in accord with experimental findings.

The structure of photosystem I from a high-light-tolerant cyanobacteria

Photosynthetic organisms have adapted to survive a myriad of extreme environments from the earth's deserts to its poles, yet the proteins that carry out the light reactions of photosynthesis are highly conserved from the cyanobacteria to modern day crops. To investigate adaptations of the photosynthetic machinery in cyanobacteria to excessive light stress, we isolated a new strain of cyanobacteria, Cyanobacterium aponinum 0216, from the extreme light environment of the Sonoran Desert. Here we report the biochemical characterization and the 2.7 Å resolution structure of trimeric photosystem I from this high-light-tolerant cyanobacterium. The structure shows a new conformation of the PsaL C-terminus that supports trimer formation of cyanobacterial photosystem I. The spectroscopic analysis of this photosystem I revealed a decrease in far-red absorption, which is attributed to a decrease in the number of long- wavelength chlorophylls. Using these findings, we constructed two chimeric PSIs in Synechocystis sp. PCC 6803 demonstrating how unique structural features in photosynthetic complexes can change spectroscopic properties, allowing organisms to thrive under different environmental stresses.

Evidence of Simultaneous Spectral Hole Burning Involving Two Tiers of the Protein Energy Landscape in Cytochrome b6f

Cytochrome b₆f, with one chlorophyll molecule per protein monomer, is a simple model system whose studies can help achieve a better understanding of nonphotochemical spectral hole burning (NPHB) and single-complex spectroscopy results obtained in more complicated photosynthetic chlorophyll-protein complexes. We are reporting new data and proposing an alternative explanation for spectral dynamics that was recently observed in cytochrome b₆f using NPHB. The relevant distribution of the tunneling parameter λ is a superposition of two components that are nearly degenerate in terms of the resultant NPHB yield and represent two tiers of the energy landscape responsible for NPHB. These two components likely burn competitively; we present the first demonstration of modeling a competitive NPHB process. Similar values of the NPHB yield result from distinctly different combinations of barrier heights, shifts along the generalized coordinate d, and/or masses of the entities involved in conformational changes m, with md² parameter different by a factor of 2.7. Consequently, in cytochrome b₆f, the first (at least) 10 h of fixed-temperature recovery preferentially probe different components of the barrier- and λ-distributions encoded into the spectral holes than thermocycling experiments. Both components most likely represent dynamics of the protein and not of the surrounding buffer/glycerol glass.

Role of missing carotenoid in reducing the fluorescence of single monomeric photosystem II core complexes

The fluorescence of monomeric photosystem II core complexes (mPSIIcc) of the cyanobacterium Thermosynechococcus elongatus, originating from redissolved crystals, is investigated by using single-molecule spectroscopy (SMS) at 1.6 K. The emission spectra of individual mPSIIcc are dominated by sharp zero-phonon lines, showing the existence of different emitters compatible with the F685, F689, and F695 bands reported formerly. The intensity of F695 is reduced in single mPSIIcc as compared to single PSIIcc-dimers (dPSIIcc). Crystal structures show that one of the β-carotene (β-Car) cofactors located at the monomer-monomer interface in dPSIIcc is missing in mPSIIcc. This β-Car in dPSIIcc is in van der Waals distance to chlorophyll (Chl) 17 in the CP47 subunit. We suggest that this Chl contributes to the F695 emitter. A loss of β-Car cofactors in mPSIIcc preparations will lead to an increased lifetime of the triplet state of Chl 17, which can explain the reduced singlet emission of F695 as observed in SMS.

Ultrafast energy transfer within the photosystem II core complex

We report 2D electronic spectroscopy on the photosystem II core complex (PSII CC) at 77 K under different polarization conditions. A global analysis of the high time-resolution 2D data shows rapid, sub-100 fs energy transfer within the PSII CC. It also reveals the 2D spectral signatures of slower energy equilibration processes occurring on several to hundreds of picosecond time scales that are consistent with previous work. Using a recent structure-based model of the PSII CC [Y. Shibata, S. Nishi, K. Kawakami, J. R. Shen and T. Renger, J. Am. Chem. Soc., 2013, 135, 6903], we simulate the energy transfer in the PSII CC by calculating auxiliary time-resolved fluorescence spectra. We obtain the observed sub-100 fs evolution, even though the calculated electronic energy shows almost no dynamics at early times. On the other hand, the electronic-vibrational interaction energy increases considerably over the same time period. We conclude that interactions with vibrational degrees of freedom not only induce population transfer between the excitonic states in the PSII CC, but also reshape the energy landscape of the system. We suggest that the experimentally observed ultrafast energy transfer is a signature of excitonic-polaron formation.

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.

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.

Structure-Based Modeling of Fluorescence Kinetics of Photosystem II: Relation between Its Dimeric Form and Photoregulation

A photosystem II-enriched membrane (PSII-em) consists of the PSII core complex (PSII-cc) which is surrounded by peripheral antenna complexes. PSII-cc consists of two core antenna (CP43 and CP47) and the reaction center (RC) complex. Time-resolved fluorescence spectra of a PSII-em were measured at 77 K. The data were globally analyzed with a new compartment model, which has a minimum number of compartments and is consistent with the structure of PSII-cc. The reliability of the model was investigated by fitting the data of different experimental conditions. From the analysis, the energy-transfer time constants from the peripheral antenna to CP47 and CP43 were estimated to be 20 and 35 ps, respectively. With an exponential time constant of 320 ps, the excitation energy was estimated to accumulate in the reddest chlorophyll (Red Chl), giving a 692 nm fluorescence peak. The excited state on the Red Chl was confirmed to be quenched upon the addition of an oxidant, as reported previously. The calculations based on the Förster theory predicted that the excitation energy on Chl29 is quenched by ChlZD1(+), which is a redox active but not involved in the electron-transfer chain, located in the D1 subunit of RC, in the other monomer with an exponential time constant of 75 ps. This quenching pathway is consistent with our structure-based simulation of PSII-cc, which assigned Chl29 as the Red Chl. On the other hand, the alternative interpretation assigning Chl26 as the Red Chl was not excluded. The excited Chl26 was predicted to be quenched by another redox active ChlZD2(+) in the D2 subunit of RC in the same monomer unit with an exponential time constant of 88 ps.

Circularly polarized luminescence spectroscopy reveals low-energy excited states and dynamic localization of vibronic transitions in CP43

Circularly polarized luminescence (CPL) spectroscopy is an established but relatively little-used technique that monitors the chirality of an emission. When applied to photosynthetic pigment assemblies, we find that CPL provides sensitive and detailed information on low-energy exciton states, reflecting the interactions, site energies and geometries of interacting pigments. CPL is the emission analog of circular dichroism (CD) and thus spectra explore the optical activity only of fluorescent states of the pigment-protein complex and consequently the nature of the lowest-energy excited states (trap states), whose study is a critical area of photosynthesis research. In this work, we develop the new approach of temperature-dependent CPL spectroscopy, over the 2-120 K temperature range, and apply it to the CP43 proximal antenna protein of photosystem II. Our results confirm strong excitonic interactions for at least one of the two well-established emitting states of CP43 named "A" and "B". Previous structure-based models of CP43 spectra are evaluated in the light of the new CPL data. Our analysis supports the assignments of Shibata et al. [Shibata et al. J. Am. Chem. Soc. 135 (2013) 6903-6914], particularly for the highly-delocalized B-state. This state dominates CPL spectra and is attributed predominantly to chlorophyll a's labeled Chl 634 and Chl 636 (alternatively labeled Chl 43 and 45 by Shibata et al.). The absence of any CPL intensity in intramolecular vibrational sidebands associated with the delocalized "B" excited state is attributed to the dynamic localization of intramolecular vibronic transitions.

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.

Variation of exciton-vibrational coupling in photosystem II core complexes from Thermosynechococcus elongatus as revealed by single-molecule spectroscopy

The spectral properties and dynamics of the fluorescence emission of photosystem II core complexes are investigated by single-molecule spectroscopy at 1.6 K. The emission spectra are dominated by sharp zero-phonon lines (ZPLs). The sharp ZPLs are the result of weak to intermediate exciton-vibrational coupling and slow spectral diffusion. For several data sets, it is possible to surpass the effect of spectral diffusion by applying a shifting algorithm. The increased signal-to-noise ratio enables us to determine the exciton-vibrational coupling strength (Huang-Rhys factor) with high precision. The Huang-Rhys factors vary between 0.03 and 0.8. The values of the Huang-Rhys factors show no obvious correlation between coupling strength and wavelength position. From this result, we conclude that electrostatic rather than exchange or dispersive interactions are the main contributors to the exciton-vibrational coupling in this system.

Carotenoid triplet states in photosystem II: coupling with low-energy states of the core complex

The photo-excited triplet states of carotenoids, sensitised by triplet-triplet energy transfer from the chlorophyll triplet states, have been investigated in the isolated Photosystem II (PSII) core complex and PSII-LHCII (Light Harvesting Complex II) supercomplex by Optically Detected Magnetic Resonance techniques, using both fluorescence (FDMR) and absorption (ADMR) detection. The absence of Photosystem I allows us to reach the full assignment of the carotenoid triplet states populated in PSII under steady state illumination at low temperature. Five carotenoid triplet ((3)Car) populations were identified in PSII-LHCII, and four in the PSII core complex. Thus, four (3)Car populations are attributed to β-carotene molecules bound to the core complex. All of them show associated fluorescence emission maxima which are relatively red-shifted with respect to the bulk emission of both the PSII-LHCII and the isolated core complexes. In particular the two populations characterised by Zero Field Splitting parameters |D|=0.0370-0.0373 cm(-1)/|E|=0.00373-0.00375 cm(-1) and |D|=0.0381-0.0385 cm(-1)/|E|=0.00393-0.00389 cm(-1), are coupled by singlet energy transfer with chlorophylls which have a red-shifted emission peaking at 705 nm. This observation supports previous suggestions that pointed towards the presence of long-wavelength chlorophyll spectral forms in the PSII core complex. The fifth (3)Car component is observed only in the PSII-LHCII supercomplex and is then assigned to the peripheral light harvesting system.

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.

Spectral shift mechanisms of chlorophylls in liquids and proteins

Origins of non-excitonic spectral shifts of chlorophylls that can reach -1,000 cm(-1) in pigment-protein complexes are actively debated in literature. We investigate possible shift mechanisms, basing on absorption and fluorescence measurements in large number of liquids. Transition wavelength in solvent-free state was estimated (±2 nm) for chlorophyll a (Chl a, 647 nm), Chl b (624 nm), bacteriochlorophyll a (BChl a, 752 nm), and pheophytines. The dispersive-repulsive shift is a predominating mechanism. It depends on polarizability difference between the ground and the excited state Δα and the Lorenz-Lorentz function of refractive index of solvent (n). The approximate (± 2Å(3)) increase of polarizability Δα is close to 15Å(3) for S(1) bands of Chl a, BChl a, and BPheo a, slightly larger for Chl b (18Å(3)), and less for Pheo a (11Å(3)). The effect of solvent polarity, expressed in terms of static dielectric permittivity (ε) is relatively minor, but characteristic for different pigments and transitions. Remarkably, maximum influence of ε on S(1) band of BChl a is less (-20 ± 10 cm(-1)) than that for Chl a (-50 ± 10 cm(-1)), and not correlated with dipole moment changes on excitation Δμ (∼2D and 0.1 ± 0.1D, respectively). Hydrogen bonding in protic solvents produces red shifts in Chl a (-60 cm(-1)) and BChl a (-100 cm(-1)), but not in Chl b. Second axial ligand of BChl a has no influence on the S(1) band, whereas the S(2) transition suffers a -400 to -600 cm(-1) down shift. Aromatic character of solvent is responsible for a ∼-100 cm(-1) red shift of both Q transitions in BChl a. The S(1) bands in chlorophylls are relatively insensitive with respect to dielectric properties and specific solvation. Therefore, nontrivial mechanisms, yielding large site-energy shifts are expected in photosynthetic chlorophyll-proteins.

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.

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.

Structure-based simulation of linear optical spectra of the CP43 core antenna of photosystem II

The linear optical spectra (absorbance, linear dichroism, circular dichroism, fluorescence) of the CP43 (PsbC) antenna of the photosystem II core complex (PSIIcc) pertaining to the S(0) → S(1) (Q(Y)) transitions of the chlorophyll (Chl) a pigments are simulated by applying a combined quantum chemical/electrostatic method to obtain excitonic couplings and local transition energies (site energies) on the basis of the 2.9 Å resolution crystal structure (Guskov et al., Nat Struct Mol Biol 16:334-342, 2009). The electrostatic calculations identify three Chls with low site energies (Chls 35, 37, and 45 in the nomenclature of Loll et al. (Nature 438:1040-1044, 2005). A refined simulation of experimental spectra of isolated CP43 suggests a modified set of site energies within 143 cm(-1) of the directly calculated values (root mean square deviation: 80 cm(-1)). In the refined set, energy sinks are at Chls 37, 43, and 45 in agreement with earlier fitting results (Raszewski and Renger, J Am Chem Soc 130:4431-4446, 2008). The present structure-based simulations reveal that a large part of the redshift of Chl 37 is due to a digalactosyldiacylglycerol lipid. This finding suggests a new role for lipids in PSIIcc, namely the tuning of optical spectra and the creation of an excitation energy funnel towards the reaction center. The analysis of electrostatic pigment-protein interactions is used to identify amino acid residues that are of potential interest for an experimental approach to an assignment of site energies and energy sinks by site-directed mutagenesis.

Calculation of chromophore excited state energy shifts in response to molecular dynamics of pigment-protein complexes

The absorption and energy transfer properties of photosynthetic pigments are strongly influenced by their local environment or "site." Local electrostatic fields vary in time with protein and chromophore molecular movement and thus transiently influence the excited state transition properties of individual chromophores. Site-specific information is experimentally inaccessible in many light-harvesting pigment-proteins due to multiple chromophores with overlapping spectra. Full quantum mechanical calculations of each chromophores excited state properties are too computationally demanding to efficiently calculate the changing excitation energies along a molecular dynamics trajectory in a pigment-protein complex. A simplified calculation of electrostatic interactions with each chromophores ground to excited state transition, the so-called charge density coupling (CDC) for site energy, CDC, has previously been developed to address this problem. We compared CDC to more rigorous quantum chemical calculations to determine its accuracy in computing excited state energy shifts and their fluctuations within a molecular dynamics simulation of the bacteriochlorophyll containing light-harvesting Fenna-Mathews-Olson (FMO) protein. In most cases CDC calculations differed from quantum mechanical (QM) calculations in predicting both excited state energy and its fluctuations. The discrepancies arose from the inability of CDC to account for the differing effects of charge on ground and excited state electron orbitals. Results of our study show that QM calculations are indispensible for site energy computations and the quantification of contributions from different parts of the system to the overall site energy shift. We suggest an extension of QM/MM methodology of site energy shift calculations capable of accounting for long-range electrostatic potential contributions from the whole system, including solvent and ions.

Optical properties, excitation energy and primary charge transfer in photosystem II: theory meets experiment

In this review we discuss structure-function relationships of the core complex of photosystem II, as uncovered from analysis of optical spectra of the complex and its subunits. Based on descriptions of optical difference spectra including site directed mutagenesis we propose a revision of the multimer model of the symmetrically arranged reaction center pigments, described by an asymmetric exciton Hamiltonian. Evidence is provided for the location of the triplet state, the identity of the primary electron donor, the localization of the cation and the secondary electron transfer pathway in the reaction center. We also discuss the stationary and time-dependent optical properties of the CP43 and CP47 subunits and the excitation energy transfer and trapping-by-charge-transfer kinetics in the core complex.

Parameters of the protein energy landscapes of several light-harvesting complexes probed via spectral hole growth kinetics measurements

The parameters of barrier distributions on the protein energy landscape in the excited electronic state of the pigment/protein system have been determined by means of spectral hole burning for the lowest-energy pigments of CP43 core antenna complex and CP29 minor antenna complex of spinach Photosystem II (PS II) as well as of trimeric and monomeric LHCII complexes transiently associated with the pea Photosystem I (PS I) pool. All of these complexes exhibit sixty to several hundred times lower spectral hole burning yields as compared with molecular glassy solids previously probed by means of the hole growth kinetics measurements. Therefore, the entities (groups of atoms), which participate in conformational changes in protein, appear to be significantly larger and heavier than those in molecular glasses. No evidence of a small (∼1 cm(-1)) spectral shift tier of the spectral diffusion dynamics has been observed. Therefore, our data most likely reflect the true barrier distributions of the intact protein and not those related to the interface or surrounding host. Possible applications of the barrier distributions as well as the assignments of low-energy states of CP29 and LHCII are discussed in light of the above results.

Investigating the early stages of photosystem II assembly in Synechocystis sp. PCC 6803: isolation of CP47 and CP43 complexes

Biochemical characterization of intermediates involved in the assembly of the oxygen-evolving Photosystem II (PSII) complex is hampered by their low abundance in the membrane. Using the cyanobacterium Synechocystis sp. PCC 6803, we describe here the isolation of the CP47 and CP43 subunits, which, during biogenesis, attach to a reaction center assembly complex containing D1, D2, and cytochrome b(559), with CP47 binding first. Our experimental approach involved a combination of His tagging, the use of a D1 deletion mutant that blocks PSII assembly at an early stage, and, in the case of CP47, the additional inactivation of the FtsH2 protease involved in degrading unassembled PSII proteins. Absorption spectroscopy and pigment analyses revealed that both CP47-His and CP43-His bind chlorophyll a and β-carotene. A comparison of the low temperature absorption and fluorescence spectra in the Q(Y) region for CP47-His and CP43-His with those for CP47 and CP43 isolated by fragmentation of spinach PSII core complexes confirmed that the spectroscopic properties are similar but not identical. The measured fluorescence quantum yield was generally lower for the proteins isolated from Synechocystis sp. PCC 6803, and a 1-3-nm blue shift and a 2-nm red shift of the 77 K emission maximum could be observed for CP47-His and CP43-His, respectively. Immunoblotting and mass spectrometry revealed the co-purification of PsbH, PsbL, and PsbT with CP47-His and of PsbK and Psb30/Ycf12 with CP43-His. Overall, our data support the view that CP47 and CP43 form preassembled pigment-protein complexes in vivo before their incorporation into the PSII complex.

Lowest electronic states of the CP47 antenna protein complex of photosystem II: simulation of optical spectra and revised structural assignments

In this work, we present simulated steady-state absorption, emission, and nonresonant hole burning (HB) spectra for the CP47 antenna complex of photosystem II (PS II) based on fits to recently refined experimental data (Neupane et al. J. Am. Chem. Soc. 2010, 132, 4214). Excitonic simulations are based on the 2.9 Å resolution structure of the PS II core from cyanobacteria (Guskov et al. Nat. Struct. Mol. Biol. 2009, 16, 334), and allow for preliminary assignment of the chlorophylls (Chls) contributing to the lowest excitonic states. The search for realistic site energies was guided by experimental constraints and aided by simple fitting algorithms. The following experimental constraints were used: (i) the oscillator strength of the lowest-energy state should be approximately ≤0.5 Chl equivalents; (ii) the excitonic structure must explain the experimentally observed red-shifted (∼695 nm) emission maximum; and (iii) the excitonic interactions of all states must properly describe the broad (non-line-narrowed, NLN) HB spectrum (including its antihole) whose shape is extremely sensitive to the excitonic structure of the complex, especially the lowest excitonic states. Importantly, our assignments differ significantly from those previously reported by Raszewski and Renger (J. Am. Chem. Soc. 2008, 130, 4431), due primarily to differences in the experimental data simulated. In particular, we find that the lowest state localized on Chl 526 possesses too high of an oscillator strength to fit low-temperature experimental data. Instead, we suggest that Chl 523 most strongly contributes to the lowest excitonic state, with Chl 526 contributing to the second excitonic state. Since the fits of nonresonant holes are more restrictive (in terms of possible site energies) than those of absorption and emission spectra, we suggest that fits of linear optical spectra along with HB spectra provide more realistic site energies.

Insight into the electronic structure of the CP47 antenna protein complex of photosystem II: hole burning and fluorescence study

We report low temperature (T) optical spectra of the isolated CP47 antenna complex from Photosystem II (PSII) with a low-T fluorescence emission maximum near 695 nm and not, as previously reported, at 690-693 nm. The latter emission is suggested to result from three distinct bands: a lowest-state emission band near 695 nm (labeled F1) originating from the lowest-energy excitonic state A1 of intact complexes (located near 693 nm and characterized by very weak oscillator strength) as well as emission peaks near 691 nm (FT1) and 685 nm (FT2) originating from subpopulations of partly destabilized complexes. The observation of the F1 emission is in excellent agreement with the 695 nm emission observed in intact PSII cores and thylakoid membranes. We argue that the band near 684 nm previously observed in singlet-minus-triplet spectra originates from a subpopulation of partially destabilized complexes with lowest-energy traps located near 684 nm in absorption (referred to as AT2) giving rise to FT2 emission. It is demonstrated that varying contributions from the F1, FT1, and FT2 emission bands led to different maxima of fluorescence spectra reported in the literature. The fluorescence spectra are consistent with the zero-phonon hole action spectra obtained in absorption mode, the profiles of the nonresonantly burned holes as a function of fluence, as well as the fluorescence line-narrowed spectra obtained for the Q(y) band. The lowest Q(y) state in absorption band (A1) is characterized by an electron-phonon coupling with the Huang-Rhys factor S of approximately 1 and an inhomogeneous width of approximately 180 cm(-1). The mean phonon frequency of the A1 band is 20 cm(-1). In contrast to previous observations, intact isolated CP47 reveals negligible contribution from the triplet-bottleneck hole, i.e., the AT2 trap. It has been shown that Chls in intact CP47 are connected via efficient excitation energy transfer to the A1 trap near 693 nm and that the position of the fluorescence maximum depends on the burn fluence. That is, the 695 nm fluorescence maximum shifts blue with increasing fluence, in agreement with nonresonant hole burned spectra. The above findings provide important constraints and parameters for future excitonic calculations, which in turn should offer new insight into the excitonic structure and composition of low-energy absorption traps.

Low-temperature protein dynamics of the B800 molecules in the LH2 light-harvesting complex: spectral hole burning study and comparison with single photosynthetic complex spectroscopy

Previously published and new spectral hole burning (SHB) data on the B800 band of LH2 light-harvesting antenna complex of Rps. acidophila are analyzed in light of recent single photosynthetic complex spectroscopy (SPCS) results (for a review, see Berlin et al. Phys. Life Rev. 2007, 4, 64.). It is demonstrated that, in general, SHB-related phenomena observed for the B800 band are in qualitative agreement with the SPCS data and the protein models involving multiwell multitier protein energy landscapes. Regarding the quantitative agreement, we argue that the single-molecule behavior associated with the fastest spectral diffusion (smallest barrier) tier of the protein energy landscape is inconsistent with the SHB data. The latter discrepancy can be attributed to SPCS probing not only the dynamics of of the protein complex per se, but also that of the surrounding amorphous host and/or of the host-protein interface. It is argued that SHB (once improved models are developed) should also be able to provide the average magnitudes and probability distributions of light-induced spectral shifts and could be used to determine whether SPCS probes a set of protein complexes that are both intact and statistically relevant. SHB results are consistent with the B800 --> B850 energy-transfer models including consideration of the whole B850 density of states.

The CP43 proximal antenna complex of higher plant photosystem II revisited: modeling and hole burning study. I

The CP43 core antenna complex of photosystem II is known to possess two quasi-degenerate "red"-trap states (Jankowiak, R. et al. J. Phys. Chem. B 2000, 104, 11805). It has been suggested recently (Zazubovich, V.; Jankowiak, R. J. Lumin. 2007, 127, 245) that the site distribution functions of the red states (A and B) are uncorrelated and that narrow holes are burned in the subpopulations of chlorophylls (Chls) from states A and B that are the lowest-energy Chl in their complex and previously thought not to transfer energy. This model of uncorrelated excitation energy transfer (EET) between the quasidegenerate bands is expanded by taking into account both electron-phonon and vibrational coupling. The model is applied to fit simultaneously absorption, emission, zero-phonon action, and transient hole burned (HB) spectra obtained for the CP43 complex with minimized contribution from aggregation. It is demonstrated that the above listed spectra can be well-fitted using the uncorrelated EET model, providing strong evidence for the existence of efficient energy transfer between the two lowest energy states, A and B (either from A to B or from B to A), in CP43. Possible candidate Chls for the low-energy A and B states are discussed, providing a link between CP43 structure and spectroscopy. Finally, we propose that persistent holes originate from regular NPHB accompanied by the redistribution of oscillator strength due to excitonic interactions, rather than photoconversion involving Chl-protein hydrogen bonding, as suggested before ( Hughes J. L. et al. Biochemistry 2006, 45, 12345 ). In the accompanying paper (Reppert, M.; Zazubovich, V.; Dang, N. C.; Seibert, M.; Jankowiak, R. J. Phys. Chem. B 2008, 9934), it is demonstrated that the model discussed in this manuscript is consistent with excitonic calculations, which also provide very good fits to both transient and persistent HB spectra obtained under non-line-narrowing conditions.