Two-dimensional electronic spectra of the photosynthetic apparatus of green sulfur bacteria

TL;DR

This paper presents a theoretical model for the two-dimensional electronic spectra of green sulfur bacteria's photosynthetic apparatus, elucidating energy transfer mechanisms and environmental interactions.

Contribution

It introduces a novel model calculation that accurately reproduces experimental spectra and transfer times, revealing environment-assisted transport as key to energy transfer.

Findings

Model matches experimental spectral features and transfer times

Identifies environment-assisted transport as crucial for energy transfer

Infers coupling strengths between chlorosome, FMO complex, and environment

Abstract

Advances in time resolved spectroscopy have provided new insight into the energy transmission in natural photosynthetic complexes. Novel theoretical tools and models are being developed in order to explain the experimental results. We provide a model calculation for the two-dimensional electronic spectra of Cholorobaculum tepidum which correctly describes the main features and transfer time scales found in recent experiments. From our calculation one can infer the coupling of the antenna chlorosome with the environment and the coupling between the chlorosome and the Fenna-Matthews-Olson complex. We show that environment assisted transport between the subunits is the required mechanism to reproduce the experimental two-dimensional electronic spectra.

Figures

Tables

Table 1. Table 1: Parameter of the coupled antenna chlorosome / FMO complex. S A 0 superscript subscript 𝑆 𝐴 0 S_{A}^{0} , λ A 0 superscript subscript 𝜆 𝐴 0 \lambda_{A}^{0} and μ A 0 superscript subscript 𝜇 𝐴 0 \mu_{A}^{0} have been used for absorption spectra calculations in Fig.  1 . The dipole orientations are given with respect to the protein data bank 3ENI FMO structure.

reorganization energy ${\lambda }_A$ (${\lambda }_A^0$)	$100\left(300\right)$ cm-1
coupling to BChl 1 and 6    $J_A$	$3$ cm-1
peak absorption $S_A$ ($S_A^0$)	$12940\left(13300\right)$ cm-1
dipole strength ${\mu }_A$(${\mu }_A^0$)	$2.1\left(8.4\right)$ $\times {\mu }_{\text{FMO}}$
antenna orientation (parallel) ${\hat{\mu }}_{A,1}$, ${\hat{\mu }}_{A,2}$	$(0.208,-0.788,0.580)$, $(-0.790,0.215,0.575)$
perpendicular ${\hat{\mu }}_{A,⟂}$	$(1,1,1)/\sqrt{3}$

Full text

Introduction

Recent experiments, for the first time, track the energy flow through the whole photosynthetic apparatus of Cholorobaculum tepidum (C. tepidum) [1]. C. tepidum are green sulfur bacteria that perform anoxygenic photosynthesis under the lowest light intensities known [2]. Their photosynthetic apparatus consists of specialized antenna structures known as chlorosome, which are linked via the Fenna-Matthews-Olson (FMO) complexes to the reaction centers (RC) where oxidation takes place. The chlorosome is a large ellipsoidal structure formed by self-assembled aggregates of $\sim 10^{5}$ bacteriochlorophyll (BChl) c pigments, organized into large oligomers with relatively little protein involvement [3]. The chlorosome rests upon a paracrystalline aggregate of BChl a pigments and CsmA proteins, the baseplate. The homotrimeric BChl a pigment-protein FMO complex sits between the chlorosome baseplate and the RC [4] (see Fig. 1). The FMO complex consists of three monomeric units, each composed of seven BChl a pigments [5] and an additional satellite BChl a pigment later discovered [6] surrounded by proteins. On a larger scale, recent electron microscopy studies propose the existence of several supercomplexes formed by two RC and four FMO trimers attached to each chlorosome antenna [7].

Optical measurements and atomistic-based calculations together with mass spectroscopy data have determined an energy funnel structure resulting in directed energy flow in C. tepidum from the antenna to the RC [4]. Light absorbed at high energy in the antenna system moves within the photosynthetic unit as electronic excitation of the individual pigments, delocalised into exciton states. The coupling of the exciton to the molecular motion and vibration leads to thermal dissipation, which eventually directs the exciton to the RC located spatially next to the energetically lowest exciton state [2]. Already within a single pigment, reorganization processes lower the energy of the excited state due to the shift of the nuclear coordinates to the equilibrium positions [8]. The transfer time from the antenna through the FMO complex to the reaction center depends critically on the strength of the inter-pigment dipole-dipole couplings, in addition to the vibrational dissipation. For achieving efficient transport, the transfer time must be short in comparison with the loss rate due to exciton decay [9, 10, 11].

The composition of the light harvesting complexes is determined using X-ray, mass spectroscopy and electron microscopy experiments. The structure of the FMO protein [12, 13], the chlorosome [3] and more recently the RC [14] and the baseplate [15] have been characterised using these methods. Fitting of optical spectra and direct structure-based quantum-chemistry calculations have been used to obtain the coupling and relaxation rates within the FMO [16] and antenna [17, 18] subsystems. Recent optical experiments analyse the time-dependent tranfer between the chlorosome and the baseplate [19, 20] and from the FMO to the RC [21]. Due to the high computational complexity and missing parameters, theoretical attempts to describe the transfer within the whole photosynthetic system are limited [22]. Using scattering resistant two-dimensional electronic spectroscopy (2DES), Dostal et al show that across the entire apparatus of C. tepidum the energy flows at a timescale of tens of picoseconds [1]. This relatively slow transfer is augmented by a faster transfer within the subunits on a timescale of sub picoseconds to a few picoseconds. 2DES is a very powerful experimental technique with high temporal and spectral resolution that has been successfully applied to study the energy-transfer pathways in photosynthetic complexes, in particular for the FMO system [23, 24]. In contrast to absorption or fluorescence measurements, 2DES provides the full correlation map between the excitations and the probing wavelengths as a function of time after initial light absorption [25, 26, 27]. This allows one to infer the couplings between the electronic transitions and additionally to probe the finite time scale of the transfer and reorganization processes (Fig. 1).

In this paper we present a theoretical modelling of the joint chlorosome and FMO system, that allows us to reproduce the experimental 2DES [1]. We use the exact hierarchical equations of motion (HEOM) formalism [28, 29, 11] that considers the exciton and the vibrational environment on the same footing. Moreover, HEOM is capable of treating the different reorganisation energies in the different subunits found in the photosynthetic apparatus. In particular the inclusion of the tightly packed antenna chlorosome leads to a larger reorganisation energy [30] ($\sim 300$ cm*-1*), which is almost tenfold increased compared to the FMO subunit [16]. By varying the model parameters we derive the effective coupling and connectivity between the antenna chlorosome and the FMO protein which matches the experimental data. In addition we extract the relative orientation of the antenna chlorosome with respect to the FMO trimer.

Results

The energy funnel structure of the C. tepidum is revealed in the calculated absorption spectra shown in Fig. 1. It presents a strong absorbing chlorosome peak at high energy followed by a sequence of lower peaks located in the FMO region. Both, the baseplate and the RC, are low absorbing subunits and are not distinguishable in the experimental absorption spectra and 2DES signal of the whole photosynthetic unit [1]. We have thus excluded them from our model calculations. From the line shapes and positions of the FMO complex and chlorosome peaks in the linear absorption we infer that the antenna chlorosome can be modeled as an effective exciton with strong environmental coupling (parameters in Table 1).

The linear absorption signal yields the dipole strengths and approximate site energies, but does not provide insight into the exciton dynamics and reorganization shifts. These processes are unfolded with the 2DES (Fig. 2). In 2DES a sequence of three laser pulses creates coherences between the ground state and between the exciton states, which are read out by a fourth pulse. The first and last time intervals between the pulses are converted by a Fourier transform to the frequency domain, with $\omega_{1}$ and $\omega_{3}$ denoting the excitation and emission frequencies, respectively. The remaining interval between second and third pulse sets the delay time $\tau$. In 2DES the energy transfer is studied by comparing the peaked signals in the frequency domain at different delay times. To avoid the overexposure of the FMO complex by the more absorptive chlorosome subunit, the experimental 2DES are measured with a different laser spectrum compared to the linear absorption spectra [1]. We do account for the modified laser spectrum by shifting the position of the chlorosome from $S_{A}^{0}$ to the peak position $S_{A}$ resulting from the convolution with the 2DES laser spectrum. The weaker laser intensity in the chlorosome region is taken into account by an effectively reduced dipole strength, Tab. 1.

The panels in Fig. 2 show the results of our 2DES calculation at $T=150$ K for delay times ranging from 30 fs to 150 ps. The signals are peaked around $(\omega_{1},\omega_{3})=(i,j)$ , where $i=1,\ldots,7,A$ denote the energies of the FMO exciton and antenna chlorosome states. For the shortest delay time $\tau=30$ fs the intensity is distributed along the diagonal line $\omega_{1}=\omega_{3}$. The diagonal peaks indicate the energies where the light was absorbed. With increasing delay time $1-150$ ps the chlorosome peak $(A,A)$ weakens in amplitude and off-diagonal cross-peaks $(i,j)$ appear below it. The redistribution of the diagonal peaks into cross peaks is the indication of the excitation flow from the antenna unit to the FMO complex. This process continues until the chlorosome peak fades around $\tau\sim 200$ ps [1]. At $1$ ps, the diagonal peaks in the FMO transitions $(i,i)$, $i=2$-$5$, decay into off-diagonal peaks, indicating energy transfer to the energetically lower FMO states. The antenna peak becomes broader and shifts to a lower $\omega_{3}$ frequency due to the chlorosome reorganization processes. At $5$ ps the intra-FMO complex flows have reached a steady state, while the chlorosome peak $(A,A)$ still decays and transfers energy into the FMO. Cross-peaks $(A,1)$ and $(A,2)$ start to appear and intensify at 50 ps, while cross-peaks intensities $(A,4)$ and $(A,5)$ remain almost constant. Starting from $100$ ps, when the delay time $\tau$ exceeds the chlorosome thermalization time, the cross peak at the intersection of the antenna and the energetically lowest FMO state $1$ (situated close to the RC on FMO BChl 3 and 4) dominates the $(A,i)$ signals. In this regime, we only observe variations in the 2DES due to the slow decay into the thermal equilibrium state.

To ensure numerical convergence of the HEOM method, which increases with increasing temperature [31], we have computed the 2DES at temperature $150$ K, while the experimental 2DES are recorded at $77$ K. Compared to the experimental values, we obtain a lower $(A,1)$ peak intensity due to the different thermal populations of the lowest eigenstate of the FMO complex ($0.6$ vs. $0.9$). This is in agreement with the temperature dependent experimental results in Fig. S1 of Ref. 32. The time scales inferred from the 2DES are only minimally affected by the change in temperature, while the reorganization energy has a more dramatic effect (see below). In our calculations we have used a small dipole-dipole coupling $J_{A}=3$ cm*-1* between the antenna and the BChl $1$ and $6$. This coupling rate, together with the reorganisation energy and the onsite energy landscape are responsible for the transfer time of the energy flow between the chlorosome antenna and the FMO complex subunits.

Model parameters

The 2DES excitation spectrum superimposes three different excitation mechanisms, namely stimulated emission (SE), ground state bleaching (GSB) and excited state absorption (ESA), see Fig. 3. The interpretation of time-scales seen in 2DES is complicated by the different peak positions and dynamics underlying the different contributions. Only the SE part tracks directly the excitonic energy transfer towards the equilibrium state, while GSB contributes with a static background signal resembling the linear absorption spectra. The ESA process is also unique to 2DES, where the pulse sequence leads to a further excitation of a second exciton. This results in a negative signal, compared to the GSB and SE contributions. Grey areas at each panel show the integrated 2D spectra over the excitation $\omega_{1}$ and emission frequency $\omega_{3}$. At large delay times, the top line of the combined signal (Fig. 3) represents the absorption spectra and the lateral dashed line represents the transient absorption spectra. We observe in Fig. 3 that the SE contribution, which represents the dynamical energy transfer towards the thermal equilibrium state, shows a similar contribution of $(A,2)$ with respect to $(A,1)$. This is not expected at the lower cryogenic temperatures in the experiments [1]. Note that inclusion of the RC [1] into our model would also mainly affect the relative weights of the cross-peaks $(A,i=1,2)$ in the SE spectra, which are determined by the thermal population weights. Additionally we observe how the $(A,i=4,5)$ cross-peaks in the combined signal (Fig. 3) are absent in the individual SE contribution. They are a feature of the GSB and ESA contributions in Fig. 3 and remain constant for long delay times, see Fig. 2. The most prominent feature in systems with larger reorganization energies, such as the chlorosome antenna, is the shift of the SE $(A,A)$ peak from the diagonal of the 2DES to the bare exciton energy (Fig. 3). With the addition of the GB and ESA parts, which are not undergoing this shift, the combined 2DES (Fig. 3) features a shift of about half the reorganization energy. From the difference of the excitation and emission frequencies in the 2DES observed at the chlorosome peak [1], we estimate the reorganizational shift of the chlorosome antenna to be $\lambda_{A}\sim 300-400$ cm*-1*, which is in line with theoretical predictions [30] and inline with our absorption spectra calculation in Fig. 1.

Next we establish the effective coupling $J_{A}$ between the FMO unit and the chlorosome antenna with its larger reorganization energy $\lambda_{A}$ . The chosen coupling strength $J_{A}$ is determined from a comparison with the experimentally recorded time scales of the exciton relaxation, namely $\sim 70$ ps scale for the antenna decay into the reaction center [1]. The transfer time, seen as population decrease of the antenna, becomes shortest for specific conjunctions of reorganization energies and coupling strengths. We show the population transfer after $60$ ps for two different laser spectra, the redshifted one used in the 2DES calculations and the unshifted one, in Fig. 4. Increasing the coupling $J_{A}$ leads to monotonous decrease of the chlorosome antenna population, while changing the reorganization energy for fixed coupling gives the optimal value of $\lambda_{A}$ for an efficient depletion [11]. For the broad band laser covering all chlorosome energies, an optimal value $\lambda^{0}_{A}\sim 300-400$ cm*-1* leads to the fastest inter-unit transfer in Fig. 4b. For the red-shifted laser and antenna peak, a similar transfer efficiency is achieved for $\lambda_{A}\sim 100-200$ cm*-1* in Fig. 4 a. In both cases, the coupling $J_{A}=3$ cm*-1* matches the experimentally observed depletion of the chlorosome antenna population to $0.3$ after $60$ ps. Here, we have chosen $J_{A}=3$ cm*-1* and $\lambda_{A}=100$ cm*-1* for the 2DES simulation. The decay time of the chlorosome antenna population $p(A)={\rm e}^{-t/t_{\rm decay}}$ in Fig. 4 is only weakly affected by temperature ($t_{\rm decay,300K}=52.5$ ps, $t_{\rm decay,100K}=57$ ps).

In addition to the position of the spectral peaks and the excitonic couplings, 2DES is sensitive to the dipole orientation. The FMO complex dipoles are oriented along the $N_{B}$, $N_{D}$ nitrogen atoms of the BChl PDB structure [33, 34]. Chemical labelling and mass spectrometry data have established the orientation of the FMO protein on the baseplate membrane [4] and the chlorosome [3] as shown in Fig. 1. For the tubular antenna complex we orient the effective dipole parallel to the baseplate. If instead a perpendicular antenna-baseplate orientation is chosen, the 2DES calculations differ and the cross peak between the chlorosome antenna and lowest exciton state $1$ is diminished, Fig. 5. The lowest exciton state $1$, located mainly on BChl 3, is in particular diminished by the perpendicular antenna orientation with respect to the baseplate due to the BChl 3 dipole orientation parallel to the baseplate plane (see Fig. 1). To account for different rotations of the antenna within the baseplate plane, we average over two orthogonal configurations (Tab. 1).

Methods

Exciton parameters

The energy transport in the photosynthetic system is modeled using the Frenkel exciton description [35], including the external electromagnetic field from the laser pulses $H_{\rm field}(t)$ in the impulsive limit [25]. Besides the seven single exciton states of the FMO system and the antenna chlorosome, additionally 28 two exciton states are explicitly included in the Hamiltonian and in the dipole matrix [26, 36]. The excitonic site energies and couplings are taken from Adolphs and Renger [16], (Tab. 4, col. 2 and Tab. 1 col. 4 respectively). The antenna is added as an additional pigment with energy in the chlorosome region $S_{A}=12940$ cm*-1* (see Fig. 2 in Ref. 1). This artificially red-shifted energy increases the FMO signal compared to the otherwise dominant absorbance of the chlorosome. Only for the linear absorption calculation (measured with a broader frequency laser covering the complete chlorosome spectrum, see Fig. 1 in Ref. 1) we take $S_{A}^{0}=13300$ cm*-1*. The chlorosome is linked to the FMO via the baseplate located on top of BChl 1, BChl 6 and BChl 8 (Fig. 1a). The baseplate is set parallel to the plane of the Mg center atoms of the BChl 8 pigments of the trimeric FMO complex. BChl 6 and BChl 1 have been considered before as entry channel for energy flow from the chlorosome [4]. The more recent discovery of an eight chlorophyll in proximity of the baseplate and BChl 1 also suggest that this pigment links the chlorosome to the FMO [37, 32]. Both, the baseplate and the eight pigment have not been observed in 2DES [1]. Therefore we do not include them in the model, but directly couple the chlorosome to the FMO BChl 1 and BChl 6. For the vibrational couplings of the BChls and the chlorosome, we consider eight independent sets of harmonic oscillators coupled linearly to each pigment. The spectral density of the vibrations of the $m$th pigment is given by Drude-Lorentz shape $J_{m}(\omega)=2\lambda_{m}\frac{\omega\gamma}{\omega^{2}+\gamma^{2}}$ with a correlation time $\gamma^{-1}=50$ fs. We set for the seven FMO BChls $\lambda_{\rm FMO}=35$ cm*-1*.

2DES calculation

2DES signals are obtained from a sequence of three ultra-short laser pulses and from recording the resulting signal at a fixed time delay $\tau$ between the second and third pulse. The resulting signal is proportional to the third order optical response function $S^{(3)}(\omega_{3},\tau,\omega_{1})$, where the excitation frequency is denoted by $\omega_{1}$ and the emission frequency by $\omega_{3}$. The frequency resolved 2DES are computed from a time propagation of the density matrix using the numerically exact hierarchical equation of motion (HEOM) method on massively parallel graphic processing units [11, 36, 38]. In contrast to perturbative treatments using Redfield or Förster theory, HEOM accounts for the different reorganization process of the FMO complex compared to the chlorosome antenna. With increasing reorganization energy, the computational and memory requirements for HEOM increase, in particular for computing the two exciton states. By using the high temperature approximation at $150$ K and truncating $\lambda_{A}=100$ cm*-1* we ensure converged results of all simulations and capture the larger reorganizational shift of the antenna. The effect of disorder is not directly considered here, but has been discussed before for various disorder parameters of the FMO complex [36] and does not affect the presented analysis. 2DES are often measured for two different laser setups, stemming from the rephasing and non-rephasing contributions. All spectra are rotationally averaged over 20 different orientations of the joint complex with respect to the laser polarization [36]. In addition we average over two orthogonal orientations of the chlorosome antenna in the baseplate plane.

Conclusion

By modeling 2DES of the photosynthetic unit of C. tepidum, we have studied energy transfer processes between the chlorosome and the FMO complex. We show how the energy funnels from higher to lower energies not only within the FMO complex but also within subunits of the photosynthetic apparatus. The structural variety of the subunits in the photosynthetic apparatus of C. tepidum leads to different reorganisational strengths. To account for these differences, we use the HEOM formalism which is required to properly model the varying environmental couplings and time-scales [39, 40]. We show excellent agreement with the main features of the experimental 2DES [1] and demonstrate that energy transfer processes occur on a slow scale of tens of picoseconds between the subsystems. The transfer time is still fast compared to the nanoseconds lifetime of excitons in isolated chlorosomes [2] and allows the photosynthetic apparatus to function with high efficiency.

From comparison with the experimental data, we estimate the effective reorganization energy of the chlorosome antenna, the inter unit coupling, and the relative orientation between subsystems. We find that the energy is transferred from the antenna chlorosome to the FMO complex via the strong environmental coupling, but with only weak dipole-dipole interactions between the subunits. The strong reorganisational energies in the chlorosome result in ranges of couplings which optimise the transfer time between the subunits.

Acknowledgements

The work was supported by the North-German Supercomputing Alliance (HLRN) and by the German Research Foundation (DFG) grant KR 2889. M.R. has received funding from the European Union’s Horizon 2020 research and innovation programme under the Marie Sklodowska-Curie grant agreement No 707636.

Author contributions statement

T.K. and M.R. contributed equally to the theoretical analysis and writing of the manuscript.

Additional information

The authors declare no competing financial interests.