Chromophore−Chromophore and Chromophore−Protein Interactions in Monomeric Light-Harvesting Complex II of Green Plants Studied by Spectral Hole Burning and Fluorescence Line Narrowing

Mapping electronic decoherence pathways in molecules

Establishing the fundamental chemical principles that govern molecular electronic quantum decoherence has remained an outstanding challenge. Fundamental questions such as how solvent and intramolecular vibrations or chemical functionalization contribute to the decoherence remain unanswered and are beyond the reach of state-of-the-art theoretical and experimental approaches. Here we address this challenge by developing a strategy to isolate electronic decoherence pathways for molecular chromophores immersed in condensed phase environments that enables elucidating how electronic quantum coherence is lost. For this, we first identify resonance Raman spectroscopy as a general experimental method to reconstruct molecular spectral densities with full chemical complexity at room temperature, in solvent, and for fluorescent and non-fluorescent molecules. We then show how to quantitatively capture the decoherence dynamics from the spectral density and identify decoherence pathways by decomposing the overall coherence loss into contributions due to individual molecular vibrations and solvent modes. We illustrate the utility of the strategy by analyzing the electronic decoherence pathways of the DNA base thymine in water. Its electronic coherences decay in [Formula: see text]30 fs. The early-time decoherence is determined by intramolecular vibrations while the overall decay by solvent. Chemical substitution of thymine modulates the decoherence with hydrogen-bond interactions of the thymine ring with water leading to the fastest decoherence. Increasing temperature leads to faster decoherence as it enhances the importance of solvent contributions but leaves the early-time decoherence dynamics intact. The developed strategy opens key opportunities to establish the connection between molecular structure and quantum decoherence as needed to develop chemical strategies to rationally modulate it.

Recent Progress in Solution Structure Studies of Photosynthetic Proteins Using Small-Angle Scattering Methods

Utilized for gaining structural insights, small-angle neutron and X-ray scattering techniques (SANS and SAXS, respectively) enable an examination of biomolecules, including photosynthetic pigment-protein complexes, in solution at physiological temperatures. These methods can be seen as instrumental bridges between the high-resolution structural information achieved by crystallography or cryo-electron microscopy and functional explorations conducted in a solution state. The review starts with a comprehensive overview about the fundamental principles and applications of SANS and SAXS, with a particular focus on the recent advancements permitting to enhance the efficiency of these techniques in photosynthesis research. Among the recent developments discussed are: (i) the advent of novel modeling tools whereby a direct connection between SANS and SAXS data and high-resolution structures is created; (ii) the employment of selective deuteration, which is utilized to enhance spatial selectivity and contrast matching; (iii) the potential symbioses with molecular dynamics simulations; and (iv) the amalgamations with functional studies that are conducted to unearth structure-function relationships. Finally, reference is made to time-resolved SANS/SAXS experiments, which enable the monitoring of large-scale structural transformations of proteins in a real-time framework.

Excitation energy equilibration in a trimeric LHCII complex involves unusual pathways

We explore the energy equilibration within the LHCII trimer using various approaches, including the Redfield-Förster method (with different compartmentalization schemes) and the exact hierarchical equation of motion (HEOM). We demonstrate that the inter-monomeric migration in the trimeric LHCII complex is not limited to direct transfers between quasi-equilibrated chlorophylls (Chls) a, but also involves additional pathways with uphill transfers from Chls a to the stromal-side Chls b (connecting the Chls a clusters from different monomeric subunits). Although these uphill transfers are slow they still can increase the total rate of inter-monomeric transfers by a factor of 1.5. The same stromal-side Chls b also promote a depopulation of the Chl a604 long-lived state (blue-shifted and mixed with the lumenal-side Chls b). Due to the connection between the stromal- and lumenal-side Chls b clusters the intra- and inter-monomeric transfers from a604 to the main Chls a become faster by a factor of 1.6 and 1.75, respectively.

Force Decomposition Analysis: A Method to Decompose Intermolecular Forces into Physically Relevant Component Contributions

Computational quantum chemistry can be more than just numerical experiments when methods are specifically adapted to investigate chemical concepts. One important example is the development of energy decomposition analysis (EDA) to reveal the physical driving forces behind intermolecular interactions. In EDA, typically the interaction energy from a good-quality density functional theory (DFT) calculation is decomposed into multiple additive components that unveil permanent and induced electrostatics, Pauli repulsion, dispersion, and charge-transfer contributions to noncovalent interactions. Herein, we formulate, implement, and investigate decomposing the forces associated with intermolecular interactions into the same components. The resulting force decomposition analysis (FDA) is potentially useful as a complement to the EDA to understand chemistry, while also providing far more information than an EDA for data analysis purposes such as training physics-based force fields. We apply the FDA based on absolutely localized molecular orbitals (ALMOs) to analyze interactions of water with sodium and chloride ions as well as in the water dimer. We also analyze the forces responsible for geometric changes in carbon dioxide upon adsorption onto (and activation by) gold and silver anions. We also investigate how the force components of an EDA-based force field for water clusters, namely MB-UCB, compare to those from force decomposition analysis.

Light-Harvesting Complex II Adopts Different Quaternary Structures in Solution as Observed Using Small-Angle Scattering

The high-resolution crystal structure of the trimeric major light-harvesting complex of photosystem II (LHCII) is often perceived as the basis for understanding its light-harvesting and photoprotective functions. However, the LHCII solution structure and its oligomerization or aggregation state may generally differ from the crystal structure and, moreover, also depend on its functional state. In this regard, small-angle scattering experiments provide the missing link by offering structural information in aqueous solution at physiological temperatures. Herein, we use small-angle scattering to investigate the solution structures of two different preparations of solubilized LHCII employing the nonionic detergents n-octyl-β-d-glucoside (OG) and n-dodecyl-β-D-maltoside (β-DM). The data reveal that the LHCII-OG complex is equivalent to the trimeric crystal structure. Remarkably, however, we observe─for the first time─a stable oligomer composed of three LHCII trimers in the case of the LHCII-β-DM preparation, implying additional pigment-pigment interactions. The latter complex is assumed to mimic trimer-trimer interactions which play an important role in the context of photoprotective nonphotochemical quenching.

Nature of low-energy exciton levels in light-harvesting complex II of green plants as revealed by satellite hole structure

Persistent non-photochemical hole burning at 4.2 K is an efficient experimental tool to unravel position and nature of low-energy excitonic states in pigment-protein complexes. This is demonstrated here for the case of the trimeric chlorophyll (Chl) a/b light-harvesting complexes of Photosystem II (LHC II) of green plants, where previous work (Pieper et al. J Phys Chem B 103:2412, 1999a) reported a highly localized lowest energy state at 680 nm. At that time, this finding appeared to be consistent with the contemporary knowledge about the LHC II structure, which mainly suggested the presence of weakly coupled Chl heterodimers. Currently, however, it is widely accepted that the lowest state is associated with an excitonically coupled trimer of Chl molecules at physiological temperatures. This raises the question, why an excitonically coupled state has not been identified by spectral hole burning. A re-inspection of the hole burning data reveals a remarkable dependence of satellite hole structure on burn fluence, which is indicative of the excitonic coupling of the low-energy states of trimeric LHC II. At low fluence, the satellite hole structure of the lowest/fluorescing  ~ 680 nm state is weak with only one shallow satellite hole at 649 nm in the Chl b spectral range. These findings suggest that the lowest energy state at  ~ 680 nm is essentially localized on a Chl a molecule, which may belong to a Chl a/b heterodimer. At high fluence, however, the lowest energy hole shifts blue to  ~ 677 nm and is accompanied by two satellite holes at  ~ 673 and 663 nm, respectively, indicating that this state is excitonically coupled to other Chl a molecules. In conclusion, LHC II seems to possess two different, but very closely spaced lowest energy states at cryogenic temperatures of 4.2 K.

Exciton properties and optical spectra of light harvesting complex II from a fully atomistic description

We present a fully atomistic simulation of linear optical spectra (absorption, fluorescence and circular dichroism) of the Light Harvesting Complex II (LHCII) trimer using a hybrid approach, which couples a quantum chemical description of the chlorophylls with a classical model for the protein and the external environment (membrane and water). The classical model uses a polarizable Molecular Mechanics force field, thus allowing mutual polarization effects in the calculations of the excitonic properties. The investigation is performed both on the crystal structure and on structures generated by a μs long classical molecular dynamics simulation of the complex within a solvated membrane. The results show that this integrated approach not only provides a good description of the excitonic properties and optical spectra without the need for additional refinements of the excitonic parameters, but it also allows an atomistic investigation of the relative importance of electronic, structural and environment effects in determining the optical spectra.

Evidence for coherent mixing of excited and charge-transfer states in the major plant light-harvesting antenna, LHCII

LHCII, the major light harvesting antenna from plants, plays a dual role in photosynthesis. In low light it is a light-harvester, while in high light it is a quencher that protects the organism from photodamage. The switching mechanism between these two orthogonal conditions is mediated by protein dynamic disorder and photoprotective energy dissipation. The latter in particular is thought to occur in part via spectroscopically 'dark' states. We searched for such states in LHCII trimers from spinach, at both room temperature and at 77 K. Using 2D electronic spectroscopy, we explored coherent interactions between chlorophylls absorbing on the low-energy side of LHCII, which is the region that is responsible for both light-harvesting and photoprotection. 2D beating frequency maps allow us to identify four frequencies with strong excitonic character. In particular, our results show the presence of a low-lying state that is coupled to a low-energy excitonic state. We assign this to a mixed excitonic-charge transfer state involving the state with charge separation within the Chl a603-b609 heterodimer, borrowing some dipole strength from the Chl a602-a603 excited states. Such a state may play a role in photoprotection, in conjunction with specific and environmentally controlled realizations of protein dynamic disorder. Our identification and assignment of the coherences observed in the 2D frequency maps suggests that the structure of exciton states as well as a mixing of the excited and charge-transfer states is affected by coupling of these states to resonant vibrations in LHCII.

Excitation energy transfer in phycobiliproteins of the cyanobacterium Acaryochloris marina investigated by spectral hole burning

The cyanobacterium Acaryochloris marina developed two types of antenna complexes, which contain chlorophyll-d (Chl d) and phycocyanobilin (PCB) as light-harvesting pigment molecules, respectively. The latter membrane-extrinsic complexes are denoted as phycobiliproteins (PBPs). Spectral hole burning was employed to study excitation energy transfer and electron-phonon coupling in PBPs. The data reveal a rich spectral substructure with a total of four low-energy electronic states whose absorption bands peak at 633, 644, 654, and at about 673 nm. The electronic states at ~633 and 644 nm can be tentatively attributed to phycocyanin (PC) and allophycocyanin (APC), respectively. The remaining low-energy electronic states including the terminal emitter at 673 nm may be associated with different isoforms of PC, APC, or the linker protein. Furthermore, the hole burning data reveal a large number of excited state vibrational frequencies, which are characteristic for the chromophore PCB. In summary, the results are in good agreement with the low-energy level structure of PBPs and electron-phonon coupling parameters reported by Gryliuk et al. (BBA 1837:1490-1499, 2014) based on difference fluorescence line-narrowing experiments.

Observation of Electronic Excitation Transfer Through Light Harvesting Complex II Using Two-Dimensional Electronic-Vibrational Spectroscopy

Light-harvesting complex II (LHCII) serves a central role in light harvesting for oxygenic photosynthesis and is arguably the most important photosynthetic antenna complex. In this work, we present two-dimensional electronic-vibrational (2DEV) spectra of LHCII isolated from spinach, demonstrating the possibility of using this technique to track the transfer of electronic excitation energy between specific pigments within the complex. We assign the spectral bands via comparison with the 2DEV spectra of the isolated chromophores, chlorophyll a and b, and present evidence that excitation energy between the pigments of the complex are observed in these spectra. Finally, we analyze the essential components of the 2DEV spectra using singular value decomposition, which makes it possible to reveal the relaxation pathways within this complex.

Subtle spectral effects accompanying the assembly of bacteriochlorophylls into cyclic light harvesting complexes revealed by high-resolution fluorescence spectroscopy

We have observed that an assembly of the bacteriochloropyll a molecules into B850 and B875 groups of cyclic bacterial light-harvesting complexes LH2 and LH1, respectively, results an almost total loss of the intra-molecular vibronic structure in the fluorescence spectrum, and simultaneously, an essential enhancement of its phonon sideband due to electron-phonon coupling. While the suppression of the vibronic coupling in delocalized (excitonic) molecular systems is predictable, as also confirmed by our model calculations, a boost of the electron-phonon coupling is rather unexpected. The latter phenomenon is explained by exciton self-trapping, promoted by mixing the molecular exciton states with charge transfer states between the adjacent chromophores in the tightly packed B850 and B875 arrangements. Similar, although less dramatic trends were noted for the light-harvesting complexes containing chlorophyll pigments.

Two-Dimensional Electronic Spectroscopy of Light-Harvesting Complex II at Ambient Temperature: A Joint Experimental and Theoretical Study

We have performed broad-band two-dimensional (2D) electronic spectroscopy of light-harvesting complex II (LHCII) at ambient temperature. We found that electronic dephasing occurs within ∼60 fs and inhomogeneous broadening is approximately 120 cm(-1). A three-dimensional global fit analysis allows us to identify several time scales in the dynamics of the 2D spectra ranging from 100 fs to ∼10 ps and to uncover the energy-transfer pathways in LHCII. In particular, the energy transfer between the chlorophyll b and chlorophyll a pools occurs within ∼1.1 ps. Retrieved 2D decay-associated spectra also uncover the spectral positions of corresponding diagonal peaks in the 2D spectra. Residuals in the decay traces exhibit periodic modulations with different oscillation periods. However, only one of them can be associated with the excitonic cross-peaks in the 2D spectrum, while the remaining ones are presumably of vibrational origin. For the interpretation of the spectroscopic data, we have applied a refined exciton model for LHCII. It reproduces the linear absorption, circular dichroism, and 2D spectra at different waiting times. Several components of the energy transport are revealed from theoretical simulations that agree well with the experimental observations.

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.

Protein dynamics tunes excited state positions in light-harvesting complex II

Light harvesting and excitation energy transfer in photosynthesis are relatively well understood at cryogenic temperatures up to ∼100 K, where crystal structures of several photosynthetic complexes including the major antenna complex of green plants (LHC II) are available at nearly atomic resolution. The situation is much more complex at higher or even physiological temperatures, because the spectroscopic properties of antenna complexes typically undergo drastic changes above ∼100 K. We have addressed this problem using a combination of quasielastic neutron scattering (QENS) and optical spectroscopy on native LHC II and mutant samples lacking the Chl 2/Chl a 612 pigment molecule. Absorption difference spectra of the Chl 2/Chl a 612 mutant of LHC II reveal pronounced changes of spectral position and their widths above temperatures as low as ∼80 K. The complementary QENS data indicate an onset of conformational protein motions at about the same temperature. This finding suggests that excited state positions in LHC II are affected by protein dynamics on the picosecond time scale. In more detail, this means that at cryogenic temperatures the antenna complex is trapped in certain protein conformations. At higher temperature, however, a variety of conformational substates with different spectral position may be thermally accessible. At the same time, an analysis of the widths of the absorption difference spectra of Chl 2/Chl a 612 reveals three different reorganization energies or Huang-Rhys factors in different temperature ranges, respectively. These findings imply that (dynamic) pigment-protein interactions fine-tune electronic energy levels and electron-phonon coupling of LHC II for efficient excitation energy transfer at physiological temperatures.

Excitation energy transfer and electron-vibrational coupling in phycobiliproteins of the cyanobacterium Acaryochloris marina investigated by site-selective spectroscopy

In adaption to its specific environmental conditions, the cyanobacterium Acaryochloris marina developed two different types of light-harvesting complexes: chlorophyll-d-containing membrane-intrinsic complexes and phycocyanobilin (PCB) - containing phycobiliprotein (PBP) complexes. The latter complexes are believed to form a rod-shaped structure comprising three homo-hexamers of phycocyanin (PC), one hetero-hexamer of phycocyanin and allophycocyanin (APC) and probably a linker protein connecting the PBPs to the reaction centre. Excitation energy transfer and electron-vibrational coupling in PBPs have been investigated by selectively excited fluorescence spectra. The data reveal a rich spectral substructure with a total of five low-energy electronic states with fluorescence bands at 635nm, 645nm, 654nm, 659nm and a terminal emitter at about 673 nm. The electronic states at ~635 and 645 nm are tentatively attributed to PC and APC, respectively, while an apparent heterogeneity among PC subunits may also play a role. The other fluorescence bands may be associated with three different isoforms of the linker protein. Furthermore, a large number of vibrational features can be identified for each electronic state with intense phonon sidebands peaking at about 31 to 37cm⁻¹, which are among the highest phonon frequencies observed for photosynthetic antenna complexes. The corresponding Huang-Rhys factors S fall in the range between 0.98 (terminal emitter), 1.15 (APC), and 1.42 (PC). Two characteristic vibronic lines at about 1580 and 1634cm⁻¹ appear to reflect CNH⁺ and CC stretching modes of the PCB chromophore, respectively. The exact phonon and vibrational frequencies vary with electronic state implying that the respective PCB chromophores are bound to different protein environments. This article is part of a special issue entitled: photosynthesis research for sustainability: keys to produce clean energy.

Refined structure-based simulation of plant light-harvesting complex II: linear optical spectra of trimers and aggregates

Linear optical spectra of solubilized trimers and small lamellar aggregates of the major light-harvesting complex II (LHCII) of higher plants are simulated employing excitonic couplings and site energies of chlorophylls (Chls) computed on the basis of the two crystal structures by a combined quantum chemical/electrostatic approach. A good agreement between simulation and experiment is achieved (except for the circular dichroism in the Chl b region), if vibronic transitions of Chls are taken into account. Site energies are further optimized by refinement fits of optical spectra. The differences between refined and directly calculated values are not significant enough to decide, whether the crystal structures are closer to trimers or aggregates. Changes in the linear dichroism spectrum upon aggregation are related to site energy shifts of Chls b601, b607, a603, a610, and a613, and are interpreted in terms of conformational changes of violaxanthin and the two luteins involving their ionone rings. Chl a610 is the energy sink at 77K in both conformations. An analysis of absorption spectra of trimers perpendicular and parallel to the C(3)-axis (van Amerongen et al. Biophys. J. 67 (1994) 837-847) shows that only Chl a604 close to neoxanthin is significantly reoriented in trimers compared to the crystal structures. Whether this pigment is orientated in aggregates as in the crystal structures, can presently not be determined faithfully. To finally decide about pigment reorientations that could be of relevance for non-photochemical quenching, further polarized absorption and fluorescence measurements of aggregates or detergent-depleted LHCII would be helpful. This article is part of a Special Issue entitled: Photosynthesis Research for Sustainability: from Natural to Artificial.

Quantum chemical description of absorption properties and excited-state processes in photosynthetic systems

The theoretical description of the initial steps in photosynthesis has gained increasing importance over the past few years. This is caused by more and more structural data becoming available for light-harvesting complexes and reaction centers which form the basis for atomistic calculations and by the progress made in the development of first-principles methods for excited electronic states of large molecules. In this Review, we discuss the advantages and pitfalls of theoretical methods applicable to photosynthetic pigments. Besides methodological aspects of excited-state electronic-structure methods, studies on chlorophyll-type and carotenoid-like molecules are discussed. We also address the concepts of exciton coupling and excitation-energy transfer (EET) and compare the different theoretical methods for the calculation of EET coupling constants. Applications to photosynthetic light-harvesting complexes and reaction centers based on such models are also analyzed.

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.

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.

Pathways of energy flow in LHCII from two-dimensional electronic spectroscopy

Photosynthetic light-harvesting complexes absorb energy and guide photoexcitations to reaction centers with speed and efficacy that produce near-perfect efficiency. Light harvesting complex II (LHCII) is the most abundant light-harvesting complex and is responsible for absorbing the majority of light energy in plants. We apply two-dimensional electronic spectroscopy to examine energy flow in LHCII. This technique allows for direct mapping of excitation energy pathways as a function of absorption and emission wavelength. The experimental and theoretical results reveal that excitation energy transfers through the complex on three time scales: previously unobserved sub-100 fs relaxation through spatially overlapping states, several hundred femtosecond transfer between nearby chlorophylls, and picosecond energy transfer steps between layers of pigments. All energy is observed to collect into the energetically lowest and most delocalized states, which serve as exit sites. We examine the angular distribution of optimal energy transfer produced by this delocalized electronic structure and discuss how it facilitates the exit step in which the energy moves from LHCII to other complexes toward the reaction center.