Modeling study of non-line-narrowed hole-burned spectra in weakly coupled dimers and multi-chromophoric molecular assemblies

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.

How Well Does the Hole-Burning Action Spectrum Represent the Site-Distribution Function of the Lowest-Energy State in Photosynthetic Pigment-Protein Complexes?

For the first time, we combined Monte Carlo and nonphotochemical hole burning (NPHB) master equation approaches to allow for ultrahigh-resolution (<0.005 cm^-1, smaller than the typical homogeneous line widths at 5 K) simulations of the NPHB spectra of dimers and trimers of interacting pigments. These simulations reveal significant differences between the zero-phonon hole (ZPH) action spectrum and the site-distribution function (SDF) of the lowest-energy state. The NPHB of the lowest-energy pigment, following the excitation energy transfer (EET) from the higher-energy pigments which are excited directly, results in the shifts of all excited states. These shifts affect the ZPH action spectra and EET times derived from the widths of the spectral holes burned in the donor-dominated regions. The effect is present for a broad variety of realistic antihole functions, and it is maximal at relatively low values of interpigment coupling (V ≤ 5 cm^-1) where the use of the Förster approximation is justified. These findings need to be considered in interpreting various optical spectra of photosynthetic pigment-protein complexes for which SDFs (describing the inhomogeneous broadening) are often obtained directly from the ZPH action spectra. Water-soluble chlorophyll-binding protein (WSCP) was considered as an example.

The lowest-energy chlorophyll of photosystem II is adjacent to the peripheral antenna: Emitting states of CP47 assigned via circularly polarized luminescence

The identification of low-energy chlorophyll pigments in photosystem II (PSII) is critical to our understanding of the kinetics and mechanism of this important enzyme. We report parallel circular dichroism (CD) and circularly polarized luminescence (CPL) measurements at liquid helium temperatures of the proximal antenna protein CP47. This assembly hosts the lowest-energy chlorophylls in PSII, responsible for the well-known "F695" fluorescence band of thylakoids and PSII core complexes. Our new spectra enable a clear identification of the lowest-energy exciton state of CP47. This state exhibits a small but measurable excitonic delocalization, as predicated by its CD and CPL. Using structure-based simulations incorporating the new spectra, we propose a revised set of site energies for the 16 chlorophylls of CP47. The significant difference from previous analyses is that the lowest-energy pigment is assigned as Chl 612 (alternately numbered Chl 11). The new assignment is readily reconciled with the large number of experimental observations in the literature, while the most common previous assignment for the lowest energy pigment, Chl 627(29), is shown to be inconsistent with CD and CPL results. Chl 612(11) is near the peripheral light-harvesting system in higher plants, in a lumen-exposed region of the thylakoid membrane. The low-energy pigment is also near a recently proposed binding site of the PsbS protein. This result consequently has significant implications for our understanding of the kinetics and regulation of energy transfer in PSII.

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

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

Selective and differential optical spectroscopies in photosynthesis

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