Exact simulation of pigment-protein complexes unveils vibronic renormalization of electronic parameters in ultrafast spectroscopy.

The primary steps of photosynthesis rely on the generation, transport, and trapping of excitons in pigment-protein complexes (PPCs). Generically, PPCs possess highly structured vibrational spectra, combining many discrete intra-pigment modes and a quasi-continuous of protein modes, with vibrational and electronic couplings of comparable strength. The intricacy of the resulting vibronic dynamics poses significant challenges in establishing a quantitative connection between spectroscopic data and underlying microscopic models. Here we show how to address this challenge using numerically exact simulation methods by considering two model systems, namely the water-soluble chlorophyll-binding protein of cauliflower and the special pair of bacterial reaction centers. We demonstrate that the inclusion of the full multi-mode vibronic dynamics in numerical calculations of linear spectra leads to systematic and quantitatively significant corrections to electronic parameter estimation. These multi-mode vibronic effects are shown to be relevant in the longstanding discussion regarding the origin of long-lived oscillations in multidimensional nonlinear spectra.

How the molecular structure determines the flow of excitation energy in plant light-harvesting complex II.

Excitation energy transfer in the light-harvesting complex II of higher plants is modeled using excitonic couplings and local transition energies determined from structure-based calculations recently (Müh et al., 2010). A theory is introduced that implicitly takes into account protein induced dynamic localization effects of the exciton wavefunction between weakly coupled optical and vibronic transitions of different pigments. Linear and non-linear optical spectra are calculated and compared with experimental data reaching qualitative agreement. High-frequency intramolecular vibrational degrees of freedom are found important for ultrafast subpicosecond excitation energy transfer between chlorophyll (Chl) b and Chla, since they allow for fast dissipation of the excess energy. The slower ps component of this transfer is due to the monomeric excited state of Chlb 605. The majority of exciton relaxation in the Chla spectral region is characterized by slow ps exciton equilibration between the Chla domains within one layer and between the lumenal and stromal layers in the 10-20ps time range. Subpicosecond exciton relaxation in the Chla region is only found within the terminal emitter domain (Chls a 610/611/612) and within the Chla 613/614 dimer. Deviations between measured and calculated exciton state life times are obtained for the intermediate spectral region between the main absorbance bands of Chla and Chlb that indicate that besides Chlb 608 another pigment should absorb there. Possible candidates, so far not identified by structure-based calculations, but by fitting of optical spectra and mutagenesis studies, are discussed. Additional mutagenesis studies are suggested to resolve this issue.

Water soluble chlorophyll binding protein of higher plants: a most suitable model system for basic analyses of pigment-pigment and pigment-protein interactions in chlorophyll protein complexes.

This short review paper describes spectroscopic studies on pigment-pigment and pigment-protein interactions of chlorophyll (Chl) a and b bound to the recombinant protein of class IIa water soluble chlorophyll protein (WSCP) from cauliflower. Two Chls form a strongly excitonically coupled open sandwich dimer within the tetrameric protein matrix. In marked contrast to the mode of excitonic coupling of Chl and bacterio-Chl molecules in light harvesting complexes and reaction centers of all photosynthetic organisms, the unique structural pigment array in the Chl dimer of WSCP gives rise to an upper excitonic state with a large oscillator strength. This property opens the way for thorough investigations on exciton relaxation processes in Chl-protein complexes. Lifetime measurements of excited singlet states show that the unusual stability towards photodamage of Chls bound to WSCP, which lack any protective carotenoid molecule, originates from a high diffusion barrier to interaction of molecular dioxygen with Chl triplets. Site selective spectroscopic methods provide a wealth of information on the interactions of the Chls with the protein matrix and on the vibronic structure of the pigments. The presented data and discussions illustrate the great potential of WSCP as a model system for systematic experimental and theoretical studies on the functionalizing of Chls by the protein matrix. It opens the way for further detailed analyses and a deeper understanding of the properties of pigment protein complexes.

Thermally activated superradiance and intersystem crossing in the water-soluble chlorophyll binding protein.

The crystal structure of the class IIb water-soluble chlorophyll binding protein (WSCP) from Lepidium virginicum is used to model linear absorption and circular dichroism spectra as well as excited state decay times of class IIa WSCP from cauliflower reconstituted with chlorophyll (Chl) a and Chl b. The close agreement between theory and experiment suggests that both types of WSCP share a common Chl binding motif, where the opening angle between pigment planes in class IIa WSCP should not differ by more than 10 degrees from that in class IIb. The experimentally observed (Schmitt et al. J. Phys. Chem. B 2008, 112, 13951) decrease in excited state lifetime of Chl a homodimers with increasing temperature is fully explained by thermally activated superradiance via the upper exciton state of the dimer. Whereas a temperature-independent intersystem crossing (ISC) rate is inferred for WSCP containing Chl a homodimers, that of WSCP with Chl b homodimers is found to increase above 100 K. Our quantum chemical/electrostatic calculations suggest that a thermally activated ISC via an excited triplet state T4 is responsible for the latter temperature dependence.

Excited state dynamics in recombinant water-soluble chlorophyll proteins (WSCP) from cauliflower investigated by transient fluorescence spectroscopy.

The present study describes the fluorescence emission properties of recombinant water-soluble chlorophyll (Chl) protein (WSCP) complexes reconstituted with either Chl a or Chl b alone (Chl a only or Chl b only WSCP, respectively) or mixtures of both pigments at different stoichiometrical ratios. Detailed investigations were performed with time and space correlated ps fluorescence spectroscopy within the temperature range from 10 to 295 K. The following points were found: (a) The emission spectra at room temperature (295 K) are well characterized by bands with a dominating Lorentzian profile broadened due to phonon scattering and peak positions located at 677, 684 and 693 nm in the case of Chl a only WSCP and at 665, 675 and 689 nm for Chl b only WSCP. In addition, all spectra contain minor bands in the longer wavelength region. (b) The emission spectra at 10 K of samples suspended in buffer containing 50% glycerol are dominated by bands peaking at 668 nm for Chl b only WSCP and at 685 nm for Chl a only WSCP and samples reconstituted with mixtures of Chl a and Chl b. (c) At 10 K and in buffer with 50% glycerol the decay kinetics of WSCP samples with Chl a only are dominated by a component with a time constant of 6.2 (+/-0.2) ns at 685 nm while those of WSCP containing mixtures of Chl a and Chl b are characterized by a slightly shorter value of 6.0 (+/-0.2) ns. WSCP containing Chl b only exhibits a distinctly longer value of 7.0 (+/-0.3) ns at an emission wavelength of 668 nm. (d) The decay associated emission spectra at 10 K of all samples exhibit at least 3 decay components with time constants of 80-120 ps, 2-4 ns and 6-7 ns in 50% glycerol. These results are consistently described within the framework of our previously presented model (J. Phys. Chem. B 2007, 111, No. 46, 13325; J. Phys. Chem. B 2007, 111, No. 35, 10487) , for the structural motifs of chlorophyll binding to the tetrameric protein matrix of WSCP. It is shown that formation of strongly coupled open sandwich dimers does not lead to quenching of 1Chl a* or 1Chl b*.

The primary electron donor of photosystem II of the cyanobacterium Acaryochloris marina is a chlorophyll d and the water oxidation is driven by a chlorophyll a/chlorophyll d heterodimer.

We present a theoretical analysis of the flash-induced absorbance difference spectrum assigned to the formation of the secondary radical pair P(+)QA(-) in photosystem II of the chlorophyll d-containing cyanobacterium Acaryochloris marina. An exciton Hamiltonian determined previously for chlorophyll a-containing photosystem II complexes is modified to take into account the occupancy of certain binding sites by chlorophyll d instead of chlorophyll a. Different assignments of the reaction center pigments to chlorophyll a or d from the literature are investigated in the calculation of the absorbance difference spectrum. A quantitative explanation of the experimental data requires one chlorophyll a molecule per reaction center, located at the site of P(D1). The remaining sites are occupied by chlorophyll d and pheophytin a. By far, the lowest site energy is found for the accessory chlorophyll of the D1 branch, Chl(D1), which represents the sink of excitation energy and therefore the primary electron donor. The cationic state is stabilized at the chlorophyll a, which drives the oxidation of water.

Pigment-pigment and pigment-protein interactions in recombinant water-soluble chlorophyll proteins (WSCP) from cauliflower.

Plants contain water-soluble chlorophyll-binding proteins (WSCPs) that function neither as antennas nor as components of light-induced electron transfer of photosynthesis but are likely constituents of regulatory protective pathways in particular under stress conditions. This study presents results on the spectroscopic properties of recombinant WSCP from cauliflower reconstituted with chlorophyll b (Chl b) alone or with mixtures of Chl a and Chl b. Two types of experiments were performed: (a) measurements of stationary absorption spectra at 77 and 298 K and CD spectra at 298 K and (b) monitoring of laser flash-induced transient absorption changes with a resolution of 200 fs in the time domain of up to 100 ps. On the basis of a theoretical analysis outlined by Renger et al. (J. Phys. Chem. B 2007, 111, 10487) the data obtained in part (a) are interpreted within a model where tetrameric WSCP binds predominantly two Chl molecules in the form of an excitonically coupled "open sandwich" dimer with a tilt angle of about 30 degrees between the chlorin planes. The time-resolved measurements on Chl a/Chl b heterodimers are described by two exponential kinetics with time constants of 400 fs and 7 ps. These kinetics are assumed to reflect a heterogeneous population of WSCPs with Chl dimers either in excitonic coupled "open sandwich" or weakly coupled geometric arrays. The 400 fs component is assigned to excited-state relaxations from the upper to the lower excitonic level of the strongly coupled "open sandwich" dimer, while the 7-8 ps component probably indicates excitation energy transfer from 1Chl b* to Chl a in a dimer array with weak coupling due to significantly longer mutual distances between the chlorin rings.

Refinement of a structural model of a pigment-protein complex by accurate optical line shape theory and experiments.

Time-local and time-nonlocal theories are used in combination with optical spectroscopy to characterize the water-soluble chlorophyll binding protein complex (WSCP) from cauliflower. The recombinant cauliflower WSCP complexes reconstituted with either chlorophyll b (Chl b) or Chl a/Chl b mixtures are characterized by absorption spectroscopy at 77 and 298 K and circular dichroism at 298 K. On the basis of the analysis of these spectra and spectra reported for recombinant WSCP reconstituted with Chl a only (Hughes, J. L.; Razeghifard, R.; Logue, M.; Oakley, A.; Wydrzynski, T.; Krausz, E. J. Am. Chem. Soc. U.S.A. 2006, 128, 3649), the "open-sandwich" model proposed for the structure of the pigment dimer is refined. Our calculations show that, for a reasonable description of the data, a reduction of the angle between pigment planes from 60 degrees of the original model to about 30 degrees is required when exciton relaxation-induced lifetime broadening is included in the analysis of optical spectra. The temperature dependence of the absorption spectrum is found to provide a unique test for the two non-Markovian theories of optical spectra. Based on our data and the 1.7 K spectra of Hughes et al. (2006), the time-local partial ordering prescription theory is shown to describe the experimental results over the whole temperature range between 1.7 K and room temperature, whereas the alternative time-nonlocal chronological ordering prescription theory fails at high temperatures. Modified-Redfield theory predicts sub-100 fs exciton relaxation times for the homodimers and a 450 fs time constant in the heterodimers. Whereas the simpler Redfield theory gives a similar time constant for the homodimers, the one for the heterodimers deviates strongly in the two theories. The difference is explained by multivibrational quanta transitions in the protein which are neglected in Redfield theory.

Intermolecular coulomb couplings from ab initio electrostatic potentials: application to optical transitions of strongly coupled pigments in photosynthetic antennae and reaction centers.

An accurate and numerically efficient method for the calculation of intermolecular Coulomb couplings between charge densities of electronic states and between transition densities of electronic excitations is presented. The coupling of transition densities yields the Förster type excitation energy transfer coupling, and from the charge density coupling, a shift in molecular excitation energies results. Starting from an ab initio calculation of the charge and transition densities, atomic partial charges are determined such as to fit the resulting electrostatic potentials of the different states and the transition. The different intermolecular couplings are then obtained from the Coulomb couplings between the respective atomic partial charges. The excitation energy transfer couplings obtained in the present TrEsp (transition charge from electrostatic potential) method are compared with couplings obtained from the simple point-dipole and extended dipole approximations and with those from the ab initio transition density cube method of Krüger, Scholes, and Fleming. The present method is of the same accuracy as the latter but computationally more efficient. The method is applied to study strongly coupled pigments in the light-harvesting complexes of green sulfur bacteria (FMO), purple bacteria (LH2), and higher plants (LHC-II) and the "special pairs" of bacterial reaction centers and reaction centers of photosystems I and II. For the pigment dimers in the antennae, it is found that the mutual orientation of the pigments is optimized for maximum excitonic coupling. A driving force for this orientation is the Coulomb coupling between ground-state charge densities. In the case of excitonic couplings in the "special pairs", a breakdown of the point-dipole approximation is found for all three reaction centers, but the extended dipole approximation works surprisingly well, if the extent of the transition dipole is chosen larger than assumed previously. For the "special pairs", a large shift in local transition energies is found due to charge density coupling.

Simulations of frequency-domain spectra: structure-function relationships in photosynthetic pigment-protein complexes.

The theory of dissipative exciton motion in chromophore complexes is applied to develop an approximate scheme for the simulation of frequency-domain linear absorption and circular dichroism. Besides lifetime broadening of the exciton lines and the inclusion of vibrational satellites in the spectra, the computations also account for static disorder. In applying the theory to a pigment protein complex of the photosynthetic light harvesting complex LHC-II of green plants the temperature dependence of linear absorption can be well reproduced.