Fluorescence emission spectra
Depending on the light intensity, the total NPQ implementation time in cyanobacterial cells is typically a few minutes (Rakhimberdieva et al. 2004; Wilson et al. 2006). PAM fluorometry and/or steady-state fluorescence spectra measured at the end of the quenching process do not reveal subtle acquisition spectral changes during NPQ evolution. Time-resolved spectrofluorimetry requires prior mathematical modelling and therefore has some shortcomings. We have taken advantage of tracking changes in the emitted fluorescence spectra of Nostoc over tens of seconds.
The initial steady-state fluorescence emission spectrum of a Nostoc sample pre-aged in the dark for 5 min contains the well-known bands of phycobiliproteins and chlorophyll (Fig. 1a). The peak at 574 nm belongs to PE and the band at 660 nm is attributed to bulk APC, while the separate band of PC with the known position at 640 nm is not visible. The separate band at 682 nm of PBLcm is not resolved as it is superimposed by the most intense band at 685 nm of PS II (Brecht et al. 2014) due to their close superposition. The bright shoulder at about 730 nm is attributed to the sum of long wavelength satellites of the enumerated bands plus some Photosystem I emission.
Under our conditions, the maximum NPQ level was reached within 180 s (3 min) against the background of a slow NPQ process and an even slower reverse process controlled by fluorescence recovery protein (FRP) (Kirilovsky and Kerfeld 2016). This time has been divided into shorter segments to show the successive decrease in intensity of the fluorescence bands. NPQ does not affect PE, as the well-resolved 574 nm peak does not change in height and shape during the whole 180 s and could serve as an additional convenient benchmark to compare the intensity decrease of other bands (Fig. 1А). This result corresponds to the fact that the excitation energy transfer from the PE to the PBS core remains at its maximum level and is therefore not affected by the NPQ, and is consistent with the data that the OCPR does not come into direct anchoring contact with the lateral cylinders of the PBS (Scott et al. 2006; Stadnichuk et al. 2009; Wilson et al. 2008).
Changes in the amplitude of other fluorescence spectral bands have been observed by tracking the NPQ formation over the first 30 s. The highest degree of quenching was registered for the main fluorescence emission peak around 685 nm. However, as in the original spectrum of the dark-adapted sample (0 s), it was not possible to separate the ApcE band from the fluorescence superposition with the 685 nm band of PS II in the spectra measured after 30, 60, 120 or 180 s of OCP-dependent phototransformation (Fig. 1a). The resolution of these closely overlapping spectral bands depends on their relative heights during the quenching process. To reveal the problem in more detail, and to be able to see the NPQ dynamics on the scale of tens of seconds we resort to difference NPQ spectra (Rakhimberdieva et al. 2010) normalised to their peak positions. It was shown that after the first 30 s of high light irradiation, the main peak in the NPQ spectrum shifted from 685 to 682.5 nm, indicating the presence of the ApcE band (Fig. 1b). This is the first direct in vivo evidence that quenching is associated with a change in ApcE fluorescence emission. The much smaller degree of NPQ intensity change in the 660 nm region (Fig. 1b) compared to the 682.5 and 685 nm bands means that the bulk APC in the PBS core is very unlikely to serve as a site for OCPR docking in vivo. Again, no bands shorter than 682 nm were observed in these repeated spectral experiments, indicating that ApcD (679 nm) is unlikely to be involved in contact with OCPR. This result provides additional evidence for the anchoring role of the PBLcm domain of ApcE to previously estimated data (Elanskaya et al. 2018; Stadnichuk et al. 2012; Zhang et al. 2014).
Rate Of Excitation Energy Transfer From Pbs To Ps Ii
It is well established that PBSs are located on the surface of the PS II dimers with their two core basal cylinders in the plane of the thylakoid membrane (Adir et al. 2009; Arteni et al. 2009; Bald et al. 1996; Bryant and Canniffe 2018). Several details of the possible relative positions of the PBS and PS II dimer with their two-fold symmetry have been revealed by cryo-EM and cross-linking/mass spectrometry techniques (Arteni et al. 2009; Chang et al. 2015). To date, a reliable way to study the formation of water-soluble PBS and membranous PS II supercomplexes remains the use of spatial computer modelling based on the known crystallographic and EM data (Bryant and Canniffe 2018; Chang et al. 2015; Zlenko et al. 2017).
A distance from the nearest phycobilin chromophores within the two basal PBS core cylinders to the antennal chlorophyll molecule(s) of the PS II pigment-protein dimer is estimated to be in the range of 40–48 Å (Chang et al. 2015; Krasilnikov et al. 2020). PBLcm/ApcF-containing APC disk of each basal cylinder is responsible for the protrusion of this disk at the bottom of the PBS core. This protrusion fits well with the hole on the cytoplasmic side of PSII and forms a tight interaction between PBS and PSII (Chang et al. 2015; Krasilnikov et al. 2020). Considering the thickness of the apoprotein layer over the near-surface chlorophylls of PS II and the gap created by this protrusion and the amorphous PBLcm loop exposed from the PBS core towards the thylakoid membrane, the distance of 42 Å most likely offered by the model (Krasilnikov et al. 2020) was used here to provide an opportunity for energy transfer from PBS to PS II. The criterion for the sufficiency of the energy transfer only from PBLcm was the agreement of the transfer time determined according to the calculations and established experimentally.
It is implicitly assumed that PBS transfers the absorbed energy to PS II via the FRET mechanism, which is characterized by the notion of the Förster radius (R0) ([Förster 1948; Lacowicz 2006):
\(Ro=\sqrt[6]{A{{\Phi }}_{f}{\chi }^{2}J\left(\lambda \right)/{n}^{4}}\) (Eq. 1),
where A = 9000ln10/128(π)5NA ≈ 8.8 10− 25 mol is a constant. Φf is the fluorescence quantum yield of phycocyanobilin chromophores in APC trimers equal to 0.6 (Matamala et al. 2007), n = 1.33 is a refractive index of the surrounding medium (Grabowski and Gantt 1978); the overlap integral J(λ) is a function of the normalised PBLcm fluorescence emission spectrum, peaked at 682 nm, taken as the donor and the chlorophyll extinction coefficient in PS II as the acceptor. According to the calculations in (Krasilnikov et al. 2020), it is equal to 5.74 10− 13 cm6 Mol− 1 for PBLcm/PS II spectral overlap. The transition dipole moments of the PBLcm phycocyanobilin chromophore and the nearest antennal chlorophyll of CP43 appear to be parallel in favorable relative orientation for FRET (Krasilnikov et al. 2020), and thus the orientation factor χ2 reaches its maximum equal to 4. Based on these values, R0 was found to be 86 Å. The characteristic time τ of the FRET rate was then calculated as follows:
\(\tau ={\tau }_{d}{\left(\frac{R}{{R}_{0}}\right)}^{6}\) (Eq. 2),
where τd = 1.5–1.6 ns (Matamala et al. 2007) is the excited state lifetime of the donor phycocyanobilin chromophore, and R equals 42 Å above. The characteristic time obtained was found to be around 20 ps, which agrees very well with the value determined experimentally for a final step of the PBS → PS II plausible energy transfer pathway (Acuña et al. 2018). According to our present considerations, other phycobilins including the ApcD and ApcF chromophores form an extensive network for energy transfer within the PBS, which routes (Zhao et al. 1992) end in the PBLcm, minimising direct energy transfer from other phycobilins of the PBS core to PS II. The principles of the PBS core architecture are the same in both the hemidiscoidal and bulkier block-type and hemiellipsoidal PBSs of red algae. In particular, the distance between the PBLcm and ApcD chromophores in the PBS cores of the red algae Griffithsia pacifica (Zhang et al. 2017), Porphyridium purpureum (Ma et al. 2020), the cyanobacteria Anabaena 7120, Synechococcus 7002 (Zheng et al. 2021) and Thermosynechococcus vulcanus (Kawakami et al. 2022) is about 30 Å in all cases determined by cryo-EM microscopy. This implies that the ApcD → PBLcm energy transfer reveals photophysical heterogeneity of PBS transfer pathways (Squires et al. 2019) and must be more efficient than the independent transfer from ApcD to PS II chlorophyll (Krasilnikov et al. 2020).
Excitonic Coupling Of Ocp And Apce
The interchromophore distance between the hECN of OCPR and potential neighbouring phycobilin chromophores, due to equally curved lateral surfaces of different APC trimer disks in the PBS core, has to be about 25 Å (Dominguez-Martin et al. 2021; Krasilnikov et al. 2020; Stadnichuk et al. 2015). This conclusion is based on the docking model of APC and OCP crystal structures and the calculated temporal thermodynamic stability of this assembled protein supercomplex (Stadnichuk et al. 2015; Zlenko et al. 2016), as well as the in vitro obtaining of one of the possible complexes of OCPR and bulk APC as part of the PBS core (Zhang et al. 2014). A priori, the relatively close distance of 25 Å gives reason to consider FRET or excitonic energy transfer (Dominguez-Martin et al. 2021; Krasilnikov et al. 2020; Stadnichuk et al. 2015) from PBS to OCP with subsequent heat dissipation from the ketocarotenoid excited state as the basis for the NPQ mechanism.
The spectral overlap of the S0 → S2 absorption spectrum of OCPR with the fluorescence emission spectrum of PBLcm was found to be minimal (Fig. 2a), and therefore the time of the estimated FRET calculations in this proposed case should be equal to ~ 1450 ps (Stadnichuk et al. 2015). A low–lying S1 state of carotenoids is dipole forbidden for absorption via one-photon transition and is therefore absent from linear absorption spectra measurements. The S1 state energy of hECN with the 0–0 band lying for OCPR at 14,000 ± 200 cm− 1 (713 nm) has been determined by time-resolved spectroscopy methods (Polivka et al. 2013). Due to a very small value of the S0 → S1 transition density of carotenoids, the extinction coefficient of hECN used for FRET calculations was assumed to be ≤ 1000 M− 1cm− 1 (Stadnichuk et al. 2015); as a result of the strong 1/R6 time dependence of the transfer efficiency and the small extinction coefficient, the estimated time of energy migration in this also proposed state should finally be 1650–1700 ps (Stadnichuk et al. 2015). The results of both types of calculations are far from the experimentally determined 20 ps when compared with the step of energy migration from PBLcm to PS II (Acuña et al. 2018), and therefore exclude the S1 and S2 states of hECN from the OCPR-dependent energy quenching in the framework of FRET theory.
A positive theory for the NPQ implementation of excitonic energy transfer from PBLcm to OCPR was established in 2015 (Stadnichuk et al. 2015). An absolute prerequisite for excitonic coupling is the complete coincidence of the donor fluorescence and acceptor absorption energy states (Krueger et al. 1998). The much broader S0-S1 absorption spectrum of OCPR, constructed by red batochromic shifting of the known S0 → S2 transition spectrum to 0–0 band energy value of S1 compared to the fluorescence spectrum of PBLcm (peak position at 682 nm corresponds to 14,665 cm− 1), here shows such a complete spectral overlap (Fig. 2b). The original 3D model of OCP docking to APC in the PBS core (Stadnichuk et al. 2015) was improved (Krasilnikov et al. 2020), to take into account the separation of OCP N- and C-polypeptide domains upon photoactivation of OCPO (Bondanza et al. 2020; Kerfeld et al. 2017) and the quenching activity of only the N-terminal domain of the formed OCPR, which is geometrically equivalent to RCP (Lewerenz et al. 2014). The theory of the exciton interaction between the OCP ketocarotenoid and the sterically neighbouring phycocyanobilin chromophore of APC was developed and described in detail in (Stadnichuk et al. 2015). The characteristic time τq of the excitonic energy transfer from phycocyanobilin of PBLcm to hECN of OCPR was written as follows (Krasilnikov et al. 2020; Stadnichuk et al. 2015):
\({\tau }_{q}=\frac{\pi \hslash {R}^{3}{n}^{2}}{\left|\chi \right|{\mu }_{a}{\mu }_{d}}\) (Eq. 3),
where, according to the geometric model (Krasilnikov et al. 2020), R = 24.5 Å is a distance between the centres of the dipoles of phycobilin and hECN obtained by the docking procedure; n = 1.33 is the refractive index (unitless) as in Eq. 1; µd = 12 D is the known transition dipole moment of the phycocyanobilin chromophore in Debeyes (Matamala et al. 2007; Ren et al. 2013) and µa = 1.8 D was previously obtained from quantum chemical calculation of the virtual S0 → S1 transition dipole moment of hECN in OCPR (Stadnichuk et al. 2015).
In the combined tripartite RCP-PBLcm-PS II model created (Krasilnikov et al. 2020), the main condition was considered above the plausible energy transfer from PBS to PS II. The transition dipole moment of phycocyanobilin of PBLcm in the model is optimally directed for photosynthetic energy transfer (|χ| = 2) to the nearest chlorophyll of PS II (#47 of CP43) (Loll et al. 2005). Relatedly, the direction of the hECN transient dipole moment vector is secondary. The combined geometric positions of the PBLcm-PSII and PBLcm-RCP pairs then only allow the placement of the RCP such that the carotenoid hECN molecule would be nearly perpendicular to the axis of the PBS core cylinders. Any other orientation of the RCP with the S0 → S1 transition dipole moment of hECN directed along the conjugated chain of the hECN carotenoid molecule (Hashimoto et al. 2018) is limited by the strong steric clashes. As a result, the orientation factor |χ| between the transition dipole moments of hECN and the chromophore of ApcE appeared to be small, only 0.24–0.42. The duration of energy transfer time and quenching by the excitonic mechanism is inversely proportional to the values of the orientation factor and the transition dipole moments of the interacting chromophores (Eq. 3). Despite the unfavourably low value of the factor |χ|, the probable characteristic time obtained here was in the range of 1.5–2.6 ps, which cannot yet be more precisely determined. The larger time previously calculated in (Stadnichuk et al. 2015) did not exceed 10 ps. Nevertheless, the obtained values of the characteristic time were two orders of magnitude smaller than purely speculative for the OCPR calculations within FRET theory and 2–5 times smaller than the final step time of energy migration from ApcE to PS II. (Acuña et al. 2018).