|Home | About | Journals | Submit | Contact Us | Français|
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.
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 known2. 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 ~105 bacteriochlorophyll (BChl) c pigments, organized into large oligomers with relatively little protein involvement3. 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 RC4 (see Fig. 1). The FMO complex consists of three monomeric units, each composed of seven BChl a pigments5 and an additional satellite BChl a pigment later discovered6 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 antenna7.
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 RC4. 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 state2. 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 positions8. 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 decay9,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 protein12,13, the chlorosome3 and more recently the RC14 and the baseplate15 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 FMO16 and antenna17,18 subsystems. Recent optical experiments analyse the time-dependent tranfer between the chlorosome and the baseplate19,20 and from the FMO to the RC21. Due to the high computational complexity and missing parameters, theoretical attempts to describe the transfer within the whole photosynthetic system are limited22. 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 picoseconds1. 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 system23,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 absorption25,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 2DES1. We use the exact hierarchical equations of motion (HEOM) formalism11,28,29 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 energy30 (~300cm−1), which is almost tenfold increased compared to the FMO subunit16. 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.
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 unit1. 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 ω1 and ω3 denoting the excitation and emission frequencies, respectively. The remaining interval between second and third pulse sets the delay time τ. 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 spectra1. We do account for the modified laser spectrum by shifting the position of the chlorosome from to the peak position SA 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, Table 1.
The panels in Fig. 2 show the results of our 2DES calculation at T=150K for delay times ranging from 30 fs to 150ps. The signals are peaked around (ω1, ω3)=(i, j), where denote the energies of the FMO exciton and antenna chlorosome states. For the shortest delay time τ=30fs the intensity is distributed along the diagonal line ω1=ω3. The diagonal peaks indicate the energies where the light was absorbed. With increasing delay time 1–150ps 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 ps1. At 1ps, 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 ω3 frequency due to the chlorosome reorganization processes. At 5ps 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 50ps, while cross-peaks intensities (A, 4) and (A, 5) remain almost constant. Starting from 100ps, when the delay time τ 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 temperature31, we have computed the 2DES at temperature 150K, while the experimental 2DES are recorded at 77K. 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 JA=3cm−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.
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 ω1 and emission frequency ω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 experiments1. Note that inclusion of the RC1 into our model would also mainly affect the relative weights of the cross-peaks in the SE spectra, which are determined by the thermal population weights. Additionally we observe how the 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 peak1, we estimate the reorganizational shift of the chlorosome antenna to be λA~300–400cm−1, which is in line with theoretical predictions30 and inline with our absorption spectra calculation in Fig. 1.
Next we establish the effective coupling JA between the FMO unit and the chlorosome antenna with its larger reorganization energy λA. The chosen coupling strength JA is determined from a comparison with the experimentally recorded time scales of the exciton relaxation, namely ~70ps scale for the antenna decay into the reaction center1. 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 60ps for two different laser spectra, the redshifted one used in the 2DES calculations and the unshifted one, in Fig. 4. Increasing the coupling JA leads to monotonous decrease of the chlorosome antenna population, while changing the reorganization energy for fixed coupling gives the optimal value of λA for an efficient depletion11. For the broad band laser covering all chlorosome energies, an optimal value 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 cm−1 in Fig. 4a. In both cases, the coupling JA=3cm−1 matches the experimentally observed depletion of the chlorosome antenna population to 0.3 after 60ps. Here, we have chosen JA=3cm−1 and λA=100cm−1 for the 2DES simulation. The decay time of the chlorosome antenna population in Fig. 4 is only weakly affected by temperature (ps, 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 NB, ND nitrogen atoms of the BChl PDB structure33,34. Chemical labelling and mass spectrometry data have established the orientation of the FMO protein on the baseplate membrane4 and the chlorosome3 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 (Table 1).
The energy transport in the photosynthetic system is modeled using the Frenkel exciton description35, including the external electromagnetic field from the laser pulses in the impulsive limit25. 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 matrix26,36. The excitonic site energies and couplings are taken from Adolphs and Renger16, (Table 4, col. 2 and Table 1 col. 4 respectively). The antenna is added as an additional pigment with energy in the chlorosome region SA=12940cm−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 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 chlorosome4. 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 FMO32,37. Both, the baseplate and the eight pigment have not been observed in 2DES1. Therefore we do not include them in the model, but directly couple the chlorosome to the FMO BChl 1 and BChl6. 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 mth pigment is given by Drude-Lorentz shape with a correlation time γ−1=50fs. We set for the seven FMO BChls cm−1.
2DES signals are obtained from a sequence of three ultra-short laser pulses and from recording the resulting signal at a fixed time delay τ between the second and third pulse. The resulting signal is proportional to the third order optical response function , where the excitation frequency is denoted by ω1 and the emission frequency by ω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 units11,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 150K and truncating 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 complex36 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 polarization36. In addition we average over two orthogonal orientations of the chlorosome antenna in the baseplate plane.
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-scales39,40. We show excellent agreement with the main features of the experimental 2DES1 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 chlorosomes2 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.
How to cite this article: Kramer, T. and Rodŕıguez, M. Two-dimensional electronic spectra of the photosynthetic apparatus of green sulfur bacteria. Sci. Rep. 7, 45245; doi: 10.1038/srep45245 (2017).
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
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.
The authors declare no competing financial interests.
Author Contributions T.K. and M.R. contributed equally to the theoretical analysis and writing of the manuscript.