Introduction

Carotenoids are abundant natural pigments synthesized by higher plants, algae, and bacteria. These pigments are broadly utilized in photosynthesis in which they simultaneously play several roles. In various types of antenna complexes and photosystems, they supplement light harvesting of (bacterio)chlorophylls ((B)Chls) by enhancing absorption in the 400–500 nm range. Simultaneously in many complexes carotenoids also have photoprotective functions. Typically carotenoids quench potentially harmful triplet states of (B)Chls (Demmig-Adams and Adams 1993; Takaichi 1999) and in some cases (B)Chl excited singlet states (Dilbeck et al. 2016; Niedzwiedzki et al. 2016b; Staleva et al. 2015). In addition carotenoids are broadly used by various autotrophic and heterotrophic organisms for coloration (Britton et al. 1995). Most carotenoids are yellow or orange, which results from an electronic absorption in the spectral range between 400 and 550 nm but some rare examples of blue carotenoids with absorption expanded to 650 nm are also known (Polivka et al. 2010). The absorption band that defines carotenoids visible color is associated with the S0 → S2 electronic transition of the linear carbon–carbon double bond conjugation that is a characteristic structural motif of all carotenoids. This motif means that the spectroscopic properties of carotenoids are idealized by complementary linear polyenes. Polyenes are assigned to a C2h point group in which the symmetry of the electronic distribution will classify electronic states as rotationally symmetric (A) or asymmetric (B). Electronic transitions between those states are limited to those that change state inversion symmetry (g ↔ u) and Pariser alternancy symmetry (pseudoparity) (+ ↔ −). Consequently for the idealized polyenes (and ultimately for carotenoids) the first possible, S0 (11A g ) → S1 (21A g ), electronic transition would be strictly forbidden due to the same (Ag) symmetry and (−) pseudoparity of both electronic states (Christensen 1999; Hudson and Kohler 1974; Schulten and Karplus 1972). However, in a recent review on theoretical modeling of electronic excitations in polyenes, the authors emphasized that pseudoparity labels used to explain allowedness of electronic transitions for real polyenes and carotenoids should be treated as an exclusive feature of computational models and not be applied to real, typically asymmetric systems (Schmidt and Tavan 2012). That conclusion implies that for real polyenes and carotenoids transitions like S0 → S1 would be in principle allowed, although they have not yet been observed in absorption spectra of almost all of the hundreds of carotenoids and polyenes currently known (Britton et al. 2004). Very recently two carotenoids, synthetic deoxyperidinin (Greco et al. 2016) and geometric central-cis isomer of naturally occurring peridinin (Niedzwiedzki and Blankenship 2018), were proposed to reveal S0 → S1 in their absorption profiles. Presence of S0 → S1 absorption band was also proposed in cryogenic absorption spectra of short peridinin analogues (Niedzwiedzki et al. 2013). In order to explain complete absence of features associated with S0 → S1 transition in almost all other carotenoids, it was suggested that a displacement of both S0 and S1 electronic states plays a crucial role (Fiedor et al. 2016). It was assumed that both states are so significantly displaced that direct access from the ground state to S1 would require molecular deformation that simply cannot be achieved, and the transition is essentially silent due to lack of overlaps in Franck–Condon factors (Fiedor et al. 2016). These show that the understanding of excited state properties of carotenoids are far from complete (Hashimoto et al. 2018).

Numerous carotenoids are found in the light harvesting 2 (LH2) complexes produced by purple phototrophic bacteria. LH2 is a transmembrane pigment–protein complex playing the role of the peripheral light harvesting antenna expanding both absorption cross section and pigment density of the photosynthetic apparatus. It absorbs light and efficiently transfers excitation energy to the light harvesting complex 1 (LH1), which is bound to the reaction center where conversion of electronic excitation energy to chemical energy occurs (Cartron et al. 2014; Qian et al. 2013). The LH2 complex is a circular structure built from multiple heterodimers of α and β polypeptides. Each αβ heterodimer houses three BChls a and one carotenoid (which varies depending on bacterial species). Spectrally, LH2 complexes contain two distinct populations of BChl a designated B800 and B850, which is a consequence of their organization within the LH2 protein scaffold. The names of the BChl a populations are derived from the wavelength positions at which their Qy absorption bands appear. The B800 ring is a set of weakly coupled monomeric BChl a molecules with their macrocycles oriented parallel to the plane of photosynthetic membrane, located between adjacent pairs of αβ polypeptides. The B850 ring is an array of strongly coupled BChl a dimers located within each αβ subunit with macrocycles oriented perpendicular to the membrane plane (McDermott et al. 1995; Papiz et al. 2003; Prince et al. 1997). Two high-resolution X-ray structures of LH2 are available from two non-sulfur purple bacteria, Rhodopseudomonas acidophila strain 10050 and Phaeospirillum molischianum (Koepke et al. 1996; McDermott et al. 1995; Papiz et al. 2003; Prince et al. 1997) as well as low-resolution projection structures of LH2 from Rba. sphaeroides (Walz et al. 1998) and Allochromatium vinosum. These structures demonstrate that LH2 adopts circular form and has from 8 to 13 αβ subunits, depending on species (Kereiche et al. 2008; Koepke et al. 1996; McDermott et al. 1995; Papiz et al. 2003; Prince et al. 1997).

Spheroidene (Sphe) is an open chain carotenoid produced as the major carotenoid in many species of purple phototrophic bacteria including the model organism Rba. sphaeroides (Takaichi 1999). Rba. sphaeroides LH2 complexes containing Sphe have been extensively studied with various static and time-resolved spectroscopic techniques (Angerhofer et al. 1995; Cong et al. 2008; Dilbeck et al. 2016; Fidler et al. 2013; Jimenez et al. 1996; Ostroumov et al. 2013; Polivka et al. 2002; Trinkunas et al. 2001). Sphe has nine conjugated carbon–carbon bonds in its backbone, a property that determines appearance of its absorption band at ~ 485 nm (0–0 vibronic peak) in low polarizable solvents like ethanol (Britton et al. 2004). Within the LH2 protein environment the 0–0 vibronic peak shifts to ~ 510 nm. Excited state properties of Sphe in various solvents have previously been studied at room temperature (RT) as well at 77 K. The lifetime of its S1 excited electronic singlet state is in the 7–9 ps range at RT and slows to ~ 12 ps at 77 K (Frank et al. 1997, 2000; Niedzwiedzki et al. 2007, 2009; Polivka et al. 1999, 2001, 2005; Rondonuwu et al. 2003; Zhang et al. 2000). For Sphe bound to the LH2 complex, the S1 excited state lifetime is substantially shorter (~ 1 ps) due to the fact that it efficiently transfers excitation energy to both B800 and B850 BChls (Cong et al. 2008; Polivka et al. 2007; Rondonuwu et al. 2004). As Sphe does not show a detectable absorption to the S1 electronic state, the exact value of the S1 state energy is still debatable. The S1 state energy level was previously investigated by spectroscopic methods including fluorescence emission (Fujii et al. 1998), resonance Raman scattering (Sashima et al. 1998) and S1 → S2 transient absorption (Niedzwiedzki et al. 2009; Polivka et al. 2001). The published S1 state energy levels were sensitive to the particular experimental techniques, yielding values spanning almost a 1000 cm−1 range (13,400–14,200 cm−1; 746–704 nm) in various solvents. Based on S1 → S2 TA spectra of LH2 containing Sphe, it was estimated that the S1 state energy of protein-bound Sphe is 13,400 cm−1 (Polivka et al. 2002) suggesting that the protein environment has little influence on the state energy. However, this value has not been confirmed with other techniques or by other research groups.

LH2 complexes from Rba. sphaeroides and Rhodoblastus acidophilus with alternative B800 pigments have previously been used to study energy transfer between the B800 and B850 (B)Chls and B850 LH2 complexes have previously been used to elucidate the role of B800 BChl a in carotenoid → B850 energy transfer (Bandilla et al. 1998; Herek et al. 2000; Linnanto and Korppi-Tommola 2002; Macpherson et al. 2001; Swainsbury et al. 2019). The primary focus of those studies was to elaborate excitation energy transfer between those (non-)native (B)Chls and B850 BChl a and to test the viability of in vivo pigment replacement to expand spectral coverage via protein engineering. However, a set of these LH2 complexes is also an excellent platform to study how energetic variation of the “B800” Qy band affects the ability of Sphe to transfer excitation from its S1 state to the (B)Chl within the B800 site, which is presented in this work. We demonstrate that the rate of energy transfer from the Sphe S1 state to “B800” is sensitive to the position of the (B)Chl Qy band, providing another indirect method to probe the energy level of S1 state of Sphe within LH2. The results support an S1 state energy of LH2-embeded Sphe in the vicinity of 13,400 cm−1 (746 nm).

Materials and methods

Bacterial growth and LH2 purification and reconstitution with non-native (B)Chls

Rba. sphaeroides cells lacking the CrtA gene [reported previously in (Chi et al. 2014)] were grown under 30 μmol photons s−1 m−2 illumination from OSRAM classic 116 W halogen bulbs in M22+ medium (Hunter and Turner 1988) supplemented with 0.1% w/v casamino acids. B800–B850 LH2 was purified by ion-exchange chromatography followed by size-exclusion as previously described (Swainsbury et al. 2019). To generate a set of LH2 complexes with alternative (B)Chls bound to the B800 site, purified LH2 was exchanged into pH 5.0 buffer containing lithium dodecyl sulfate (LDS). LDS selectively removes the B800 BChl a resulting in a loss of its characteristic 800 nm peak generating B850 LH2 named for its sole 850 nm band in the near-infra red region (Clayton and Clayton 1981; Fraser et al. 1999; Kramer et al. 1984; Robert and Frank 1988). Subsequently, B850 LH2 was incubated in excess of alternative (B)Chls (Chl a, Chl d, acChl a and BChl b) that readily bind to the vacant B800 site yielding complexes with new spectral features. For a detailed protocol see (Swainsbury et al. 2019).

Spectroscopic methods

Room temperature (RT) and 77 K steady-state absorption spectra were recorded using a UV-1800 spectrophotometer (Shimadzu, Japan). Fluorescence emission and excitation spectra of LH2 complexes were recorded at RT using an RF-6000 spectrofluorometer (Shimadzu, Japan). Fluorescence emission spectra were measured following excitation of the (0–0) vibronic peak of the carotenoid (Sphe) band at 510 nm with excitation and emission bandwidths of 5 nm, and with a 570-nm long pass glass filter at the detector entrance. Fluorescence excitation spectra were collected monitoring emission at 850 nm with excitation and emission bandwidths of 5 nm, with a 665-nm long pass glass filter at the detector entrance. Cryogenic measurements at 77 K were taken using a VNF-100 liquid nitrogen cryostat (Janis, USA). The samples were resuspended in 60/40 glycerol/buffer (v/v) mixture placed in 1-cm polymethyl methacrylate cuvettes and slowly frozen in nitrogen vapor.

Transient absorption (TA) experiments were carried out using Helios, a femtosecond time-resolved pump-probe absorption spectrometer (UltrafastSystems LLC, USA) coupled to a Spectra-Physics femtosecond laser system described in detail previously (Niedzwiedzki et al. 2016b). All LH2 complexes were excited at 510 nm corresponding to (0–0) vibronic peak of Sphe absorption band. The energy of the excitation beam, focused on the sample in a spot approximately 1 mm diameter, was 200 nJ corresponding to 7 × 1013 photons cm−2 per pulse.

Data processing and global analysis

Prior to analysis, TA datasets were corrected for temporal dispersion using Surface Xplorer 4.0 (UltrafastSystems LLC, USA). The datasets were globally fitted with a kinetic model assuming a sequential population of excited states/species in cascade of nonreversible, decreasing rates (longer lifetimes). The fitting procedure of TA datasets gives EADS—evolution-associated difference spectra (van Stokkum et al. 2004). According to this model, the TA signal at any time delay and wavelength, ΔA(t,λ), can be reconstructed from superposition of nth Ci(t) and EADSi(λ) products.

$$\Delta A(t,\lambda ) = \sum\limits_{i = 1}^{n} {C_{i} (t){\text{EADS}}_{i} (\lambda )}$$
(1)

Ci(t) is the time-dependent concentration of ith EADS which is expressed as

$$\frac{{dC_{i} (t)}}{dt} = k_{i - 1} C_{i - 1} (t) - k_{i} C_{i} (t),i \ne 1,k_{i - 1} > k_{i}$$
(2)

and C1(t) is the populated by the excitation pulse that in spectrometer is represented as instrument response function, IRF:

$$\frac{{dC_{1} (t)}}{dt} = IRF(t) - k_{1} C_{1} (t)$$
(3)

The IRF was simulated by a Gaussian with a full-width at half-maximum (FWHM) of ~ 200 fs. Global analysis was performed using CarpetView 1.0 (Light Conversion Ltd., Lithuania). All plots were done in Origin 2019 (OriginLab Corp., USA).

Results

Steady-state absorption and fluorescence

Room temperature (RT) and 77 K absorption spectra (normalized at the Sphe (0–0) vibronic band at ~ 510 nm), and RT fluorescence emission spectra (normalized at their maxima) of the LH2 complexes are shown in Fig. 1. All LH2 complexes in this study contain the same Sphe carotenoid and native BChl a at the B850 site, the amplitudes of which are expected to remain consistent for all samples.

Fig. 1
figure 1

Absorption and fluorescence emission spectra of Rba. sphaeroides LH2 antenna complexes biding native and non-native (B)Chls in the B800 binding site. a RT and b 77 K absorption spectra of LH2 complexes: with a native BChl a in B800 binding site (B800–B850), B800-free (B850) and with the B800 site reconstituted with Chl a, acChl a, Chl d, and BChl b. c Normalized RT fluorescence emission spectra upon excitation at (0–0) vibronic band of Sphe (~ 510 nm)

The B800–B850 LH2 is used as the benchmark, naturally binding BChl a in the B800 site with an absorption maximum at 800 nm. In this sample, the B850 band is positioned at 849 nm and B800/B850 peak ratio is 0.72. Typically, for in-tact LH2 this ratio ranges between 0.7 and 0.8 at room temperature and variation depends on purification protocols and storage buffers (Cong et al. 2008; Dilbeck et al. 2016). The (0–0) vibronic band of Sphe absorption appears at 512 nm. Removal of BChl a from the B800 binding sites (B850) only marginally alters spectroscopic features associated with other pigments. The B850 band is slightly bathochromically shifted to 851 nm and Sphe absorption is hypsochromically shifted with the (0–0) vibronic band appearing at 507 nm, in agreement with a previous study on the same LH2 complex (Bandilla et al. 1998; Clayton and Clayton 1981; Polivka et al. 2007). These spectral shifts are the result of structural rearrangements and loss of interactions with the B800 BChl a, altering the microenvironments of the B850 BChl a and Sphe pigments (Kramer et al. 1984; Robert and Frank 1988; Saga et al. 2017). Incorporation of non-native (B)Chls (Chl a, Chl d, acChl a and BChl b) into B800 binding site has no meaningful effect on the positions of Sphe and B850 absorption bands which appear within a 2-nm range with respect to native B800-B850 LH2. The Qy band of the reconstituted (B)Chls appears at: 667 nm (Chl a), 691 nm (acChl a), 705 nm (Chl d), and ~ 830 nm (BChl b). Cryogenic temperatures (Fig. 1b) significantly improve the vibronic resolution of the absorption bands of Sphe and causes hypsochromic shifting of the Qy band of B800 (B)Chls to up to 5 nm. The Qy absorption band of BChl b is also clearly resolved at 77 K. Detailed information on the position of all spectral futures of these LH2 complexes at both RT and 77 K is listed in Table 1.

Table 1 Positions of spectral features in RT and 77 K absorption spectra of LH2 complexes consisting various (B)Chls at the B800 site

Fluorescence emission spectra (normalized at their maxima) of all LH2 complexes recorded at RT are provided in Fig. 1c. Fluorescence emission peaks were at 857 nm for all except B850-LH2, which was at 860 nm. The weak emission band visible at 805 nm in the spectrum of B800–B850 LH2 (805/857 ratio of 0.03) is associated with direct fluorescence emission from B800 BChl a. Similarly, the expansion of the emission band between 810 and 840 nm in the spectrum of BChl b-B850 LH2 is associated with direct fluorescence emission from BChl b. Analogous, weak emission bands were recorded also for Chl-containing LH2s but are outside of the shown scale and were omitted for clarity.

Fluorescence excitation spectra (monitoring emission from B850 BChl a at 860 nm) overlaid with absorptance (1-transmittance (1-T)) spectra are shown in Fig. 2. Energy transfer between non-native (B)Chls and B850 BChl a was extensively elaborated in past study (Swainsbury et al. 2019) and will not be discussed here. The main focus of this investigation is to determine how the overall Sphe-to-B850 BChl a energy transfer is affected by incorporating alternative (B)Chls into the B800 site. These pigments serve the role of intermediary in energy transfer between Sphe and the B850 BChl a, as has been well established for the native B800 BChl a (Macpherson et al. 2001; Polivka et al. 2007). All spectra were normalized at maximum of Qx band of BChl a at 590 nm, assuming 100% energy transfer efficiency (EET) to B850 Qy providing an internal standard to which the Sphe signals can be compared. In the benchmark B800-B850 LH2 complex, Sphe-to-B850 BChl a (Sphe → B850) EET was 85%, comparable to results from previous studies (Cong et al. 2008; Dilbeck et al. 2016; Rondonuwu et al. 2004). Removal of BChl a from B800 site (B850) causes a substantial reduction of the Sphe → B850 EET to 67%, clearly demonstrating that a substantial portion of excitation energy from the carotenoid is intermediately routed via B800 in agreement with previous studies (Macpherson et al. 2001; Polivka et al. 2007). Incorporation of BChl b into the B800 site restores Sphe → B850 EET to range comparable with the benchmark sample. Recovery of Sphe → B850 EET is smaller if Chls are reconstituted into B800 site, for Chl d, Sphe → B850 EET increases only marginally to 70% and to 74% for (ac)Chl a. Due to a substantial energy gap between the hypothetical fluorescence emission from the S2 electronic state of Sphe (hypothetically a mirror image of the absorption around the (0–0) vibronic band) and absorption spectrum of Qy band of reconstituted (B)Chls, Sphe S2 → (B)Chls Soret → B850 excitation energy transfer is negligible. It is more complicated for BChl a and b because those pigments have strong absorption associated with the Qx band at ~ 600 nm that substantially overlaps the theoretical S2 emission spectrum. This spectral overlap means that the Sphe S2 → BChls Qx → B850 channel is also viable. However, for Chls with negligible strength of the Qx transition, the Sphe S1 → Chls Qy → B850 route will provide most of the extra contribution to the overall Sphe → B850 EET.

Fig. 2
figure 2

Fluorescence excitation spectra (Exc) and absorptance (1-transmittance (1-T)) spectra of Rba. sphaeroides LH2 complexes biding native and non-native (B)Chls in the B800 binding site recorded at RT. Fluorescence excitation profiles were recorded monitoring emission at 860 nm from B850 band. The spectral profiles are normalized to the same value at the Qx band of BChl a (590 nm) assuming 100% EET to B850 Qy. Sphe-to-B850 EET was calculated at the indicated wavelength (vertical dashed line)

As mentioned in the introduction, the energy of the first excited singlet state (S1) of Sphe has previously been measured via various spectroscopic methods. Fluorescence emission spectroscopy revealed an extremely small quantum yield of fluorescence from the S1 state of 10−7 and a state energy of 14,200 cm−1 (704 nm) (Fujii et al. 1998), which was later supported by resonance Raman scattering (Sashima et al. 1998). On the other hand, transient absorption studies probing the S1 → S2 transition of Sphe in solvents placed the S1 state energy at 13,400 cm−1 (746 nm) (Niedzwiedzki et al. 2009; Polivka et al. 2001) and similar S1 state energy was proposed for Sphe embedded in LH2 based on results from TA (Polivka et al. 2002). Having in hand a set of LH2 complexes substituted with various (B)Chls in the B800 site allowed us to probe Sphe S1 state energy with TA using a very different approach. Steady-state absorption spectra of the LH2 complexes demonstrate that the Qy band of reconstituted (B)Chls spans the 666–830 nm range and partially overlaps with the proposed 704–746 nm range of the Sphe S1 state. Therefore, one might expect that at some point both energies may coincide. When the Sphe S1 → (B)Chls Qy excitation transfer route becomes permissible, it will manifest as a decrease in lifetime of the Sphe S1 excited state, which was examined by TA upon excitation of the Sphe (0–0) vibronic band.

Transient absorption

Transient absorption measurements were taken at RT and 77 K. Figure 3 shows TA spectra of Sphe recorded at time delays corresponding to maximum signal of S1 → Sn excited state absorption (ESA) (0.5–0.9 ps) (Fig. 3a) and at delay times after complete decay of S1 state (12–17 ps) (Fig. 3b). For better comparability, all TA spectra were normalized at the (0-1) vibronic band of Sphe (~ 477 nm). To aid visualization, inverted steady-state absorption spectra of the B800-B850 and B850 LH2 are also provided (dashed lines) allowing better comparison of bleaching bands in the TA spectra with the shapes and positions of the Sphe vibronic bands. The TA spectra comprise several features, some of which are not associated with the carotenoid and correspond to bleaching of Qx band of BChl a and b, and bleaching of Qy bands of the reconstituted Chls. The bleaching of Qx for B850 BChl a (595 nm—RT, 597 nm—77 K) and Qy of Chl a (670 nm—RT, 665 nm—77 K), acChl a (696 nm—RT, 693 nm—77 K) and Chl d (709 nm—RT, 705 nm—77 K) were visible in the data. Bleaching of the Qx band of BChl b is not clearly noticeable. It should be noted that bleaching of the Chl Qy bands may not correspond well with their position in the steady-state absorption spectrum due to the overlapping contribution of a typically red-shifted stimulated emission.

Fig. 3
figure 3

Transient absorption spectra and dynamics of Sphe S1 state in LH2 complexes with various (B)Chls in the B800 site. a, b RT and c, d 77 K TA spectra of Sphe recorded at time delays corresponding to maximum of S1 → Sn ESA (0.5–0.9 ps, a, c) and at delay times immediately after complete decay of the S1 state (12–17 ps, b, d). For better comparability, all TA spectra were normalized at Sphe (0–1) vibronic band. Inverted absorption spectra of B800–B850 and B850 LH2 are also shown. e, f Sphe S1 state dynamics probed at maximum of S1 → Sn ESA. The traces were recorded at RT and 77 K and are normalized at their maxima

The main transient feature associated with the carotenoid is the S1 → Sn ESA band that peaks at ~ 556 nm. Note that after normalization of the spectra to the ground state bleaching at 477 nm, the amplitude of the S1 → Sn ESA of Sphe differs for LH2s containing different (B)Chls in the B800 site. The bleaching of the Sphe ground state absorption should be the same in all cases because it corresponds to the initial pool of excited carotenoids. However, the amplitude of S1 → Sn ESA may differ in each LH2 as it will depend on what fraction of excited molecules (S2 state) will internally convert to the S1 state. There are various channels competing with S2 → S1 internal conversion (IC) including possible energy transfer to the Qy bands of the (B)Chls directly or via Qx, or S2 mediated singlet-to-triplet fission. Given that the (B)Chls in the B800 site differently participate in these processes, there is no guarantee that in all LH2 complexes the same portion of excited carotenoids will internally convert to the S1 state. Data show that the decay of the spectral feature at ~ 530 nm takes substantially longer than that of S1 → Sn ESA, and that its true spectral shape became apparent at delay times after S1 → Sn ESA had completely decayed. Those spectra are shown in Fig. 3b and d. This spectral feature was previously observed and abbreviated as the S* transient spectrum, where S* corresponds to some vague excited singlet state (Cong et al. 2008; Papagiannakis et al. 2002, 2006). However, our recent TA investigations on various carotenoids in LH2 and LH1 demonstrated that S* may have different origins depending on the complex. In native B800-B850 LH2 S* could be explained as a combination of ground state bleaching of carotenoid within a Car·+-B800 BChl a pair, which forms with a reasonable yield upon direct excitation of carotenoid (Cong et al. 2008; Niedzwiedzki et al. 2016a), and a temporary electrochromic response of neutral carotenoid in proximity to a Car·+-B800 BChl a pair. In LH1 the S* feature is associated with quickly decaying carotenoid triplets populated via S2 singlet fission as it is reminiscent of T1Tn ESA (Kakitani et al. 2007; Niedzwiedzki et al. 2017). Therefore, the broad and flat positive band in the spectrum of the B800–B850 LH2 suggests that its TA spectrum is dominated by temporal electroabsorption of neutral Sphe, combined with bleaching of ground state absorption of Sphe involved in Sphe·+-B800 BChl a pair as was previously demonstrated for LH2 containing the carotenoid neurosporene (Niedzwiedzki et al. 2016a). In the B850 and Chl a-B850 LH2, the sharp positive band at 527 nm (RT) and 529 nm (77 K), which dominated TA spectra after decay of the S1 state of Sphe, suggests that the TA signals are mostly contributed by T1Tn ESA for Sphe bound to LH2 (Kakitani et al. 2007). Figure 3e and f shows Sphe S1 state dynamics (rise and decay) at RT and at 77 K in all studied LH2 systems. The dynamics were recorded at maximum of S1 → Sn ESA band and for better comparability are normalized at the maximum. Differences in Sphe S1 excited state lifetime upon incorporation of various (B)Chls into the B800 site were apparent. In order to get more precise information on dynamics of Sphe S1 state, kinetic traces from Fig. 3e and f were subjected to fitting with sum of exponentially decaying components convoluted by the IRF function of the spectrometer. Additionally, the TA datasets were also independently subjected to global fitting.

Global fitting

The results of global fitting are shown in Fig. 4. The data were fitted according to unidirectional sequential model assuming irreversible and slower rates in each subsequent decay step. This kind of fitting of TA data gives EADS—evolution-associated difference spectra (van Stokkum et al. 2004). These are given for each LH2 complex at both RT and 77 K. The EADS were normalized to the amplitudes of their time-dependent concentrations C(t) to represent their contributions to the TA spectra. For reference, each panel also includes an inverted absorption spectrum of the LH2 that was adjusted to closely match bleaching of the Sphe absorption band. Each EDAS is color coded and their lifetimes are provided in the panel legend. In addition, legends contain short descriptions for each EADS. If the EADS lifetime was shorter than FWHM of the IRF (~200 fs) of the TA spectrometer, the lifetime was denoted as < IRF. It should be noted that EADS do not necessary represent the true transient spectra of a particular molecular species (e.g., specific excited state of Sphe or (B)Chl). However, if EADS are spectrally and temporally distinct their interpretation and assignment is rather straightforward. The first EADS (black solid line, present in all LH2s) is associated with decay of Sphe S2 excited state, known to have lifetime in the range of 150–250 fs depending on polarizability of the surrounding medium (Ricci et al. 1996). This EADS also has a clear S2 manifestation—bleaching of the state absorption and adjacent stimulated emission visible as negative bands at the long wavelength side of the absorption bleaching. The EADS with lifetime between 0.4 and 0.6 ps (dash-dot-dash, red), that is present in all fits except for 77 K B800–B850 LH2 is associated with vibrational relaxation of Sphe S1 (hotS1) excited electronic state. The EADS with characteristic sharp transient absorption band (S1 → Sn ESA band) at ~ 555 nm (solid red) is associated with decay of vibrationally equilibrated S1 state. Note that the S1 lifetime varies with (B)Chl incorporated to the B800 site from 1.6 ps (1.22 ps at 77 K) in the B800–B850 LH2 to 3.2 ps (3.93 ps at 77 K) in the B850 LH2 with various intermediate values for other (B)Chls. The EADS shown in solid blue with lifetime of ~ 12 ps (~ 6 ps at 77 K) for the B800–B850 was originally described as the “S*” feature associated with transient band of vague excited singlet state (Cong et al. 2008; Papagiannakis et al. 2002, 2006). As mentioned earlier, the spectral signature of this feature can be described as a combination electrochromic response of non-excited neutral carotenoid placed next to Sphe·+-B800 BChl a pair, always populated with significant yield in native LH2 complexes as previously demonstrated experimentally (Cong et al. 2008; Niedzwiedzki et al. 2016a; Polivka et al. 2004) and theoretically (Wormit and Dreuw 2006), and ground state bleaching of Sphe·+. This “S*” EADS is present only in global fitting results for the native B800–B850 LH2 as the presence of B800 BChl a is necessary for formation of the cation–anion pair, which appears to be unable to form with alternative B800 pigments. The blue color-coded EADS in the fitting results of the remaining LH2 complexes differs from the EADS present in fittings of the B800–B850 LH2. EADS associated with the electrochromic response have a much broader positive band while for all other “blue” EADS the bands are narrower, and these characteristics are enhanced at 77 K. Another EADS with lifetime of 77 ps has the same shape as EADS with an infinite lifetime; therefore, both are associated with transient band of the same excited state. Previously, these spectro-kinetic components were observed in neurosporene-containing LH2 with B800 BChl a removed and in LH1 complexes (Gradinaru et al. 2001; Niedzwiedzki et al. 2016a). These EADS are associated with fast and slow decay of the carotenoid triplet state, T1. It was demonstrated that in the absence of B800 BChl a in the antenna complex (B850 LH2 or any LH1) upon carotenoid direct excitation, T1 state very efficiently formed (Niedzwiedzki et al. 2016a, 2017) most likely via singlet fission (Tavan and Schulten 1987), presumably while carotenoid molecule is in the excited S2 state. This process apparently occurs for Sphe in all LH2 complexes reconstituted with non-native pigments in the B800 site—the S − T (singlet-minus-triplet) band of Sphe (ground state bleaching with adjunct T1Tn ESA band) is very evident. It is worth noticing that lifetime of the Sphe triplet formed via singlet fission was not uniform in all cases, which would suggest it is associated with triplet–triplet annihilation within the carotenoid and is not affected by (B)Chl reconstitution. In addition one would expect instantaneous triplet–triplet annihilation, which is linkly to occur within the temporal resolution of the spectrometer. Therefore, it is difficult to explain why the lifetime of the Sphe triplet pool varies substantially from complex to complex. Note that long-lived (essentially infinite in this time scale) Sphe triplets in the B800–B850 LH2 are formed via sensitization of the BChl a triplet and this process is slow and may not be as efficient as singlet fission. The cyan EADS with a lifetime of 1.2 ns (1.4 ns at 77 K) is associated with the recovery of excited B850 BChl a (B850*) and the major signature associated with this pigment is bleaching of the Qx band at 600 nm. Note that in this case the EADS provides valuable information for the effective lifetimes of the molecular species that are formed during excitation decay; however, due to the complexity of the real pathway of excitation decay, their spectral shapes are only approximate and in some cases may be far away from real ones. This is due to mixing of simultaneously appearing/decaying signals. In order to visualize it, a more realistic, branching model that breaks the effective lifetimes into microscopic rates was also performed for the RT TA data of B800–B850 LH2 and is shown in Figure S1 for comparative purposes. A summary of the lifetimes for the excited Sphe S1 state are gathered in Table 2 along with results from single wavelength fits.

Fig. 4
figure 4

Global fitting results (EADS) of TA datasets obtained after excitation at (0–0) vibronic band of Sphe. The EADS were multiplied by maximal value of their time-dependent concentration (Ci) and they represent relative contributions in the TA datasets. After the lifetime value of each EADS, a short description of the process associated with it is provided. For comparative purposes an inverted absorption spectrum of each LH2 complex is also shown. EADS—evolution-associated difference spectra, < IRF—lifetime shorter than FWHM of IRF function

Table 2 Sphe excited S1 state dynamics obtained from global and single wavelength fittings

Discussion

Estimate of S1 state energy for LH2-bound spheroidene

Global fitting of energy transfer from Sphe to substituted (B)Chls in the LH2 B800 site, particularly the trend in the lifetime of the Sphe S1 state, can be used to indirectly estimate the S1 state energy level of carotenoid bound within the LH2 complex. Currently for Sphe bound to LH2, only one experimental estimate of 13,400 cm−1 (746 nm) is available (Polivka et al. 2002). More classical spectroscopic measurements via recording fluorescence emission positioned the Sphe S1 state energy at 14,200 cm−1 (704 nm) in solvent (Fujii et al. 1998). Currently it is uncertain which value is more accurate as it is assumed that the solvent environment will not alter S1 state energy, meaning it may span the 13,400–14,200 cm−1 range. Our data show that the Sphe S1 → (B)Chls Qy energy transfer rate depends upon the (B)Chl reconstructed in the B800 site; therefore, the energy of Qy band, and by consequence its overlap with the hypothetical Sphe S1 state, influences the energy transfer process. This is an indication that energy transfer between Sphe and (B)Chls occurs via resonance transfer, characteristic for weakly coupled molecules (Scholes et al. 1997), which in principle requires good spectral overlap between emission of a donor (carotenoid) and absorption of an acceptor ((B)Chl Qy band):

$$k_{\text{ET}} = \frac{1}{c\hbar }\left| T \right|^{2} J$$
(4)
$$J = \int\limits_{{ - \infty }}^{\infty } {F_{D} \left( \nu \right)\varepsilon _{A} \left( \nu \right)d\nu }$$
(5)

where kET is the rate constant for Sphe S1 → (B)Chl Qy energy transfer, T is the electronic coupling term, J is the spectral overlap integral, FD is a donor emission (S1, Sphe), and \(\varepsilon_{\text{A}}\) is an acceptor absorption (Qy, (B)Chl).

Figure 5a shows spectra for the hypothetical emission from the S1 state of LH2-bound Sphe and Qy absorption band of (B)Chls reconstructed to LH2, illustrating the overlap between them. Because we cannot acquire fluorescence emission spectra from the Sphe S1 state, it was mimicked by taking RT S1 emission of 3,4,5,6-tetrahydrospheroidene, a spheroidene analogue with shorter conjugation (N = 8), recorded in petroleum ether and shifting it to the expected S1 state energy of Sphe. The (0–0) vibronic band of the fluorescence spectrum was positioned at the hypothetical energies of the S1 state (13,400 cm−1—solid red, 14,200 cm−1—dash-dot red). The electronic coupling term (T) between Sphe and the B800-bound (B)Chl was assumed to be a constant as most of physical properties of Sphe and (B)Chl (distance, orientation, etc.) should remain similar. What clearly changes is the spectral overlap integral (J). The dependence of J on the energy/position of (B)Chl Qy band, and its overlap with the Sphe S1 emission spectrum, is expected to relate directly to experimentally obtained Sphe S1 → (B)Chl QykET rates. The kET rates were obtained as follows. The effective/observed Sphe S1 state decay rate kobs (kobs = 1/τS1) obtained from fitting of data from native and reconstituted LH2s consists of kobs = kint + kET850 + kET800, where kint is intrinsic decay of the S1 state back to the ground state, and kET850, kET800 rates are associated with excitation energy transfer to both the B800 and B850 sites of LH2. For the LH2 without B800 (B850 LH2, rate denoted kobs): kobs = kint + kET850. Clearly, it follows that the kET800 rate can be obtained as: kET800 = kobs − kobs. The Sphe S1 → (B)Chl QykET rates obtained from fitting results of TA datasets are plotted in Fig. 5b in black symbols (open circles, RT, solid circles, 77 K). On the same graph, the spectral overlaps for both Sphe S1 energies are also plotted in blue symbols. As we do not have 77 K fluorescence emission of 3,4,5,6-tetrahydrospheroidene, spectral overlaps for cryogenic temperature were impossible to calculate. The overlaps have unitless values and the data were scaled to get the best coincidence with corresponding set of experimentally calculated kET800 rates for comparative purposes. It is noticeable that kET800 for Sphe S1 → BChl b Qy is an apparent outlier from the overall trend. Our past studies demonstrated that this pigment correctly binds within the B800 site (Swainsbury et al. 2019) so this result is rather surprising. However, in one characteristic BChl b is unique with respect to all other (B)Chls. Once bound, its Qy absorption band substantially overlaps (even at cryogenic temperature) with Qy absorption band of the B850 BChl a. Consequently, B850 BChl a → B800 BChl b excitation transfer is possible as it will be thermodynamically allowed. In this scenario, thermal equilibration of excited B850 BChl a and B800 BChl b will occur and the system will collectively decay with a common rate. This scheme has consequences for observed lifetime (kob)−1 of the Sphe S1 state in BChl b LH2; statistically, upon excitation of the Sphe S2 state a substantial pool of excitations will be transferred to B850 BChl a by the S2 → Qx route. After instantaneous IC, B850 BChl a → B800 BChl b feedback could populate BChl b Qy before it can accept excitation from the Sphe S1 state. Additionally, for B800 BChl b, Sphe S2 → BChl b Qx excitation energy transfer is possible, which quickly populates BChl b Qy via IC. With a few possible excitation transfer routes competing for the same population of molecules, the Sphe S1 → BChl b Qy route is probably the most disadvantageous as it must first go through vibrational equilibration (as demonstrated in the fitting results) before transferring excitation energy. As the Sphe S1 → BChl b Qy process may not be fast enough to compete with the other two processes, it is likely to be less efficient than one might expect; thus, BChl b was excluded from drawing conclusions on the overall trend. For the remaining (B)Chls, a remarkable correlation occurs between the kET800 and spectral overlaps for Sphe S1 state energy at 13,400 cm−1. On the other hand, the trend in spectral overlaps for Sphe S1 state energy of 14,200 cm−1 is flatter and does not correlate well with the experimental kET800. This is an indication that the previously suggested energy of 13,400 cm−1 is a good estimate of the actual value. This has additional confirmation in kET800 plot at 77 K, which is still valuable despite our inability to compare it with spectral overlaps. It is expected that at 77 K vibronic bands of the fluorescence spectrum of Sphe will substantially narrow and spectral overlap will be more sensitive to position of the (B)Chl Qy band. If the Sphe S1 state energy is indeed 14,200 cm−1 (704 nm) there should be a substantial difference in kET rates between Chl a, acChl a and Chl d because it is likely that only the Qy band of Chl d (705 nm) will overlap the narrow emission from the Sphe (0–0) vibronic band. However, in the plot of experimental kET rates all three Chls are essentially the same, meaning the spectral overlap is similar for all of them. This suggests that the hypothetical Sphe fluorescence spectrum is shifted to lower energies (longer wavelengths) than the proposed 14,200 cm−1 level, again suggesting it resides closer to the 13,400 cm−1 level.

Fig. 5
figure 5

Spectral overlaps between hypothetical Sphe S1 fluorescence emission and (B)Chl Qy band absorption and Sphe S1 → B800 (B)Chl Qy rates (kET800) as function of (B)Chl Qy energy. (a) Hypothetical fluorescence emission spectrum of Sphe from the S1 state for two different energy levels (13,400 cm−1 and 14,200 cm−1) along with Qy band absorption spectra of all (B)Chls present in studied LH2 complexes. These spectra were used to calculate spectral overlaps (J). (b) Experimental Sphe S1 → B800 (B)Chl QykET rates and calculated spectral overlaps (unitless, multiplied by scaling factor for the best coincidence) as a function of energy of Qy band of (B)Chl present in the B800 site

In LH2 with non-native (B)Chls, the overall Sphe → B850 EET not clearly linked with the efficiency of the Sphe S1 → (B)Chl Qy transfer route

The static fluorescence excitation experiments shown in Fig. 2 demonstrate that the trend in the Sphe → B850 EET does not follow with Sphe S1 → B800 (B)Chl QykET800. This is due to several factors that may vary from complex to complex. One is the yield of singlet fission for Sphe in LH2s reconstructed with non-native (B)Chls. As this process practically eliminates substantial pool of excited carotenoids from participation in the energy transfer process (decreasing the overall Sphe → B850 EET), its presence is highly undesirable. It seems that the process is limited in the native B800-B850 LH2 but is always present at higher but different degree in LH2s with altered B800 sites. The second factor is less than 100% EET between reconstituted B800 (B)Chls and B850, as can be seen for Chl d in Fig. 2. Therefore, excitation energy transfer from the carotenoid to the intermediate B800 (B)Chl a may experience additional losses before it is accepted by B850. The last factor that should be considered is the Sphe S2 → B800 (B)Chl Qx excitation energy transfer route. It seems that for BChls a and b in the B800 site this process may be considerably more efficient than the three Chls. This is because both BChls have a distinct Qx band ~ 600 nm that allows energy transfer via S2 state while the Chls lack significant absorption in this spectral region.

Summary

In this work we performed spectroscopic analysis of six variants of the LH2 antenna complex from Rba. sphaeroides, including native B800–B850, B800-free LH2 (B850) and four LH2 s with non-native (B)Chls reconstituted into the B800 site. We particularly focused how these modifications affect excitation energy transfer between naturally bound carotenoid Sphe and the B800 (B)Chls. Interestingly there is no apparent trend in the overall energy transfer from Sphe to B850 BChl a. This could be explained due to contributions of few, sometimes competing processes, that become active (or inactive) in LH2 s with altered B800 sites. On the other hand, the trend in the Sphe S1 → B800 (B)Chl Qy energy transfer is noticeable and experimental outcomes could be used to check the validity of previously proposed energy values of the S1 state of Sphe embedded in LH2. The results support the LH2-bound Sphe S1 state energy in the vicinity of 13,400 cm−1 (746 nm).