Friederike Allgöwer‡
a,
Abhishek Sirohiwal‡
a,
Ana P. Gamiz-Hernandez‡
a,
Maximilian C. Pöverlein
a,
Andrea Fantuzzi
b,
A. William Rutherford
b and
Ville R. I. Kaila
*a
aDepartment of Biochemistry and Biophysics, Stockholm University, 10691 Stockholm, Sweden. E-mail: ville.kaila@dbb.su.se
bDepartment of Life Sciences, Imperial College London, London SW7 2AZ, UK
First published on 12th May 2025
Photosystem II (PSII) is powered by the light-capturing properties of chlorophyll a pigments that define the spectral range of oxygenic photosynthesis. Some photosynthetic cyanobacteria can acclimate to growth in longer wavelength light by replacing five chlorophylls with long wavelength pigments in specific locations, including one in the reaction center (RC) (Science, 2018, 360, 1210–1213). However, the exact location and the nature of these long wavelength pigments still remain uncertain. Here we have addressed the color-tuning mechanism of the far-red light PSII (FRL-PSII) by excited state calculations at both the ab initio correlated (ADC2) and linear-response time-dependent density functional theory (LR-TDDFT) levels in combination with large-scale hybrid quantum/classical (QM/MM) simulations and atomistic molecular dynamics. We show that substitution of a single chlorophyll pigment (ChlD1) at the RC by chlorophyll d leads to a spectral shift beyond the far-red light limit, as a result of the protein electrostatic, polarization and electronic coupling effects that reproduce key structural and spectroscopic observations. Pigment substitution at the ChlD1 site further results in a low site energy within the RC that could function as a sink for the excitation energy and initiate the primary charge separation reaction, driving the water oxidation. Our findings provide a basis for understanding color-tuning mechanisms and bioenergetic principles of oxygenic photosynthesis at the far-red light limit.
Chlorophyll a (Chl a) pigments, abundantly present in the photosynthetic proteins, enable efficient light-harvesting of PSII with a maximum absorption around 670 nm (1.85 eV). It has been suggested that the specific Chl a involved in PSII photochemistry at 685 nm (1.81 eV) define the “red-limit” of oxygenic photosynthesis, where photons beyond 700 nm (energy <1.77 eV) do not efficiently trigger the charge separation responsible for water oxidation and release of O2. However, the cyanobacterium Chroococcidiopsis thermalis acclimatizes to far-red light, by replacing its conventional PSII (white-light PSII; WL-PSII) with a far-red isoform (FRL-PSII), which incorporates one Chl d and four Chl f pigments instead of Chl a into specific locations, while retaining all the other Chl a pigments in the antenna and reaction center (Fig. 1A).7,8 This site-specific ‘doping’ of the pigments in the FRL-PSII, results in a significant redshift of the absorption spectrum at room temperature with a peak absorption maximum around 715–725 nm (1.71–1.74 eV), partially resolved at cryo-temperatures.7
Locating the position and function of these far-red light pigments in PSII is crucial for understanding the mechanistic principles of far-red light driven oxygenic photosynthesis. Nürnberg et al.7 proposed that the far-red chlorophylls were not just acting as an extension to the antenna, but crucially, one of them was serving as the primary donor within the RC. This proposal7 was based on the efficient far-red light-driven charge separation in both room and low temperature, as well as hyper-luminescence compared to conventional cyanobacteria PSII. The finding that the Chl d and f pigments not only worked efficiently to absorb lower energy photons, but also acted as the primary electron donor in primary charge separation set both thermodynamic and kinetic boundaries on the chemically challenging photochemistry.7 Structural factors can be inferred by comparing the chlorophyll binding sites in WL-PSII and FRL-PSII, in particular by the identification of amino acid changes from WL to FRL that introduce new hydrogen-bonding partners for the modified formyl group, leading to a stabilization of the Chl d or f in the binding sites. Despite previous modelling efforts, spectroscopy, and structural studies, the exact location of all the far-red light absorbing pigments remain uncertain.7–13 For the primary donor, Nürnberg et al.7 proposed that the ChlD1 site was occupied by either a Chl f or the Chl d pigment, while the remaining four Chl f pigments were located in the CP43 and CP47 antennae (cf. ref. 7 for spectroscopic assignment). However, based on the assignment of low-temperature circular dichroism data Judd et al.9 discussed the possibility that the PD2 site could be occupied by Chl f, while accepting that the ChlD1 location was also a likely candidate. Ultrafast spectroscopic experiments10 were interpreted in a model consistent with a far-red pigment at ChlD1, the primary donor, but at that time PD2 was not considered as an option. Recent cryo-EM data (PDB ID: 7SA3)11 suggested a Chl d pigment at the ChlD1 site and four Chl f pigments in the integral CP43/CP47 antenna complexes, but these assignments were made at least in part based on the presence of hydrogen-bonding residues close to the putative formal group substitutions, rather than on direct detection of the formyl oxygens, which are structurally difficult to detect11,14 as the pigments differ only by the chemical substituents on C2 and C3 centers (Fig. 1C).
The changes in the pigments are accompanied by subtle alterations in the pigment binding pockets.7,8,12 In the FRL-PSII, the D1 core protein has D1-Tyr120 and D1-Thr155 near ChlD1, rather than D1-Phe119 and D1-Ala154 in WL-PSII (Fig. S1 and S2†).7 The positioning of D1-Tyr120 and D1-Thr155 are thought to facilitate hydrogen-bonding with the formyl group at C3 and C2 of Chl d and Chl f, respectively. In addition, an amino acid substitution was observed near PheoD1, which remains as a Pheo a in FRL-PSII, where the hydrogen-bonding D1-Gln130 in the WL-PSII is replaced by D1-Glu131 in the FRL-PSII. This substitution near PheoD1 is observed as a response to high-light conditions in Chl a-containing PSII (encoded by the psbA3 gene in cyanobacteria T. elongatus and in chloroplast PSII) tuning its reduction potential to higher values.15,16 However, the mechanisms by which these substitutions could facilitate pigment-binding and spectral tuning remain unclear.
To address the molecular basis of the spectral tuning mechanism to longer wavelengths and its bioenergetic consequences, we perform here large-scale multiscale quantum/classical (QM/MM) simulations and excited state calculations at the correlated ab initio quantum chemical level, using second-order algebraic diagrammatic construction (ADC(2)) theory,17 as well as linear-response time-dependent density functional theory (LR-TDDFT). To probe the functional dynamics of the substitutions, we also performed atomistic molecular dynamics simulation (MD) on models of both the FRL-PSII and WL-PSII (see Fig. S3 and Table S1 for computational details†) together with electrostatic calculations to probe redox tuning effects. Our atomistic models allow us to systematically probe the effect of Chl d and Chl f substitutions at distinct positions in the RC, and to derive spectral tuning principles. Our combined findings provide a basis for understanding how the RC traps excitation energy, and triggers charge-separation using the decreased energy available in far-red light.
![]() | ||
Fig. 2 Absorption spectrum of the FRL-PSII reaction center. (A) Computed absorption spectra based on the tetramer assembly computed at the RVS-LT-SOS-ADC(2)/MM level of theory and comparison with the experimental spectra of the FRL-PSII. The VEEs and the corresponding oscillator strengths were used to compute the spectrum, which was aligned with respect to the 674 nm peak to help the comparison (see ESI, Extended methods†). The FRL region (715–725 nm) is highlighted. Note that the experimental spectra contain the effect of substitutions also beyond the RC (see ref. 7). (B) Computed absorption spectra for substitutions at the ChlD1 site at the TDDFT/MM level (ωB97X-V/def2-TZVP, see Methods). The full spectra, including also the substitutions at the PD2 site, can be found in Fig. S4A.† (C) Effect of D1-Tyr120 and D1-Thr155 on the absorption spectrum computed at the TDDFT level, and comparison to experimental data. The absorption spectrum shows how selectively removing either D1-Tyr120 and D1-Thr155 blueshifts the far-red transition. (D) Structure of the QM region of the tetrameric pigment model, including the PD1, PD2, ChlD1, PheoD1, and surrounding amino acids and water molecules. Hydrogen bonds are shown in dashed orange lines, and hydrogen atoms are omitted for clarity. |
For the first-principles treatment of electron correlation effects, we applied the second-order algebraic diagrammatic construction theory, ADC(2), in combination with systematic virtual space reduction (RVS)18 and the Laplace transformed (LT) scaled opposite-spin (SOS) treatment.19–21 The LT-SOS-RVS approximations were carefully benchmarked, revealing a consistent, small shift (<0.1 eV; Table S13†) that was applied to all LT-SOS-RVS-ADC(2) values. All chlorophyll types showed low D1 diagnostics in the ADC(2) calculations, indicating a single-reference ground state that provides a reliable basis for the excited state calculations (Table S2†). The RVS-LT-SOS-ADC(2) approach allowed us to extend our model to account for the complete reaction center chlorophyll pigment core (6 pigments: 4 chlorophylls and 2 pheophytins) in their protein surroundings, providing a quantitative treatment of electronic excited states,22 an approach that has shown promising results for a range of photobiological systems.21,23–26 All QM/MM models were also treated at the LR-TDDFT level in combination with modern range-separated functionals (e.g. ωB97X-V;27 see Methods, and Table S3†), predicting results highly similar to those from our correlated QM/MM calculations. Our chlorophyll spectra, calculated as an ensemble average over QM/MM molecular dynamics (QM/MM-MD) simulations, reproduce the spectra computed using QM/MM optimized models (Fig. S4E†), thus indicating that our methodology is robust, despite the computationally highly challenging treatment of electronic excited states in extended multi-pigment models. However, as the absorption maxima of the FRL- and WL-PSII differ by only 0.1 eV and are thus within the error limit of the employed methodology, we rely here also on structural and energetic considerations by probing the influence of the environment on the Chl d and f substitutions with the RC.
At the QM/MM-ADC(2) level (Fig. 2A), when Chl d or Chl f are modeled at the ChlD1 site (with Chl a in the other RC pigment locations), we obtain long wavelength transitions at 1.769 eV/701 nm (Chl d @ ChlD1) and 1.749 eV/709 nm (Chl f @ ChlD1), respectively. Moreover, when Chl d is modeled at the PD2 site (Chl d @ PD2), a prominent absorption feature at 1.703 eV/728 nm is observed, while with Chl f at the PD2 site (Chl f @ PD2), an even longer wavelength absorption at 1.606 eV/772 nm arises from the local excitation at PD2 (Fig. S5†). Both models with PD2 substitutions are strongly redshifted, but the Chl d substitution is closer to the experimental absorption at 1.72 eV/720 nm.7
In all four FRL-PSII models as well as in the WL-PSII model, a shorter wavelength transition at 674 nm (1.840 eV) can be observed, corresponding to a local excitation at the PD1/PD2 pair (Fig. S5†). This spectral feature arises both at the ADC(2) and TDDFT level (Fig. 2), suggesting that overall our results are accurate, even when considering the intrinsic error sources of the methods.22,28
Conversely, the TDDFT-spectrum of the pigment clusters with Chl d at the ChlD1 site aligns well with the experimental spectrum for both the far-red band (709 nm/1.749 eV) and the PD1/PD2 band at 674 nm (1.840 eV), accurately reproducing of the relative intensities of the Q band (Fig. 2B, and Table S4†). Taken together, our findings suggest that Chl d at the ChlD1 site reproduces key optical features of the FRL-PSII, with further energetic and structural implications (see below).
Location | Chl d/f @ ChlD1 | Chl d/f @ PD2 | WL-PSII |
---|---|---|---|
PD1/PD2 | 1.936/1.943 | 1.867/1.786 | 1.946 |
ChlD1 | 1.904/1.897 | 2.007/2.006 | 1.993 |
The ChlD1 site shows the lowest Qy excitation energies when either Chl d (1.904 eV) or Chl f (1.897 eV) are modeled at the ChlD1 site (Table 1). We note that the excitation is localized on the ChlD1 (Fig. S6†), supporting that the long wavelength pigment could act as a trap for the excitation energy. The differences between the Qy energies of ChlD1 and PD1/PD2 are rather small (0.032 eV for Chl d @ ChlD1, 0.046 eV for Chl f @ChlD1), suggesting that the exciton could be shared amongst multiple sites (see Table S5†).
For the WL-PSII model, the lowest Qy excitation comprises the PD1/PD2 pair together with the ChlD1 site. In comparison, the site energy of the Chl a pigment at ChlD1 in the WL-PSII model (1.993 eV) is at shorter wavelengths, as expected, relative to Chl d (1.904 eV) and Chl f (1.897 eV). The excitonic couplings between PD1 and PD2 are significantly larger when Chl d (97 cm−1) or Chl f (105 cm−1) occupies the ChlD1 site, whereas for Chl a in the WL-PSII model we obtain an excitonic coupling of 77 cm−1 that is slightly lower than the experimental estimates of 80–150 cm−1 (see Extended methods, and Table S12†).29–31
When Chl d or Chl f occupy the PD2 site, the lowest Qy excitation energies are predicted at the PD1/PD2 pair (1.867 eV for Chl d, 1.786 eV for Chl f; Table 1). Relative to the WL-PSII model (1.946 eV) these Qy energies are significantly red shifted. The first excited state (S1) of the PD1/PD2 pair comprises a local excitation on PD2 with a small contribution from PD1 (Fig. S5†) that can be considered a “delocalized” state, whilst we obtain stronger excitonic coupling of 180–192 cm−1 amongst the Chl a/d and a/f pairs at the PD1/PD2 site.
Overall, the site energies support that the ChlD1 substitutions lead to trapping of the excitation energy at the ChlD1 site, similar to what is known about the charge separation process in the WL-PSII, whilst the PD2 substitutions could result in trapping of the excitation energy at the PD1/PD2 site and a non-dominant charge transfer pathway.32,33 Taken together, the predicted site energies thus support that the ChlD1 substitutions lead to a productive charge separation pathway in the FRL-PSII.
Interestingly, the protein electrostatics, including hydrogen-bonding effects, also have a redshifting effect for Chl d (−0.003 eV), but blueshifts Chl f (+0.056 eV) for the ChlD1 site, further suggesting that protein electrostatics favors Chl d binding. The differences in electrostatic tuning effects can be attributed to the distinct charge density difference upon the photoexcitation (Fig. S11†) that is favored by the proximity of the C3 formyl group on Chl d to the hydrogen-bonding D1-Tyr120/D1-Thr155 (Fig. S7†), resulting in an overall redshift for Chl d, but no hydrogen-bond and an overall blueshift for Chl f (see also Fig. 2C and S8†). When D1-Tyr120 and D1-Thr155 are modeled into the WL-PSII (replacing D1-Phe119 and D1-Ala154), a blueshift occurs of only 5 nm (0.015 eV) in the Qy band of the Chl a at the ChlD1 site (Fig. 3 and Table S7†).
Excitonic coupling with surrounding pigments and quantum mechanical polarization effects further redshift the excitation energy (0.032 eV) relative to the monomeric pigment for Chl d (Fig. 3). Overall, these findings suggest that the intrinsic light-capturing properties of the Chl d pigment are enhanced relative to Chl f at the ChlD1 site by a combination of strain, electrostatics, and coupling/polarization effects, which result in tuning its color and redox properties as the primary electron donor (see below).
![]() | ||
Fig. 4 Analysis of MD simulations of the FRL-PSII and WL-PSII. (A) ChlD1 and its potential hydrogen-bonding partners near the formyl group, namely D1-Tyr120 and D1-Thr155 (D1-Phe119 and D1-Ala154 in the WL-PSII). (B and C) Hydrogen-bonding distance between vinyl (Chl a)/formyl (Chl d/f) group of ChlD1 and the (B) D1-Tyr120/D1-Phe119 or (C) D1-Thr155/D1-Ala154 (FRL-PSII/WL-PSII) from the MD simulations. The dotted black line shows the average distance from the WL-PSII MD simulation as a reference. MD distances are shown as a box plot, while QM/MM optimized distances are indicated as colored dashed lines. Boxes indicate the second and third quartiles of the distribution. (D) Water molecules around ChlD1. The water molecule W2 is bound between W1 and the C133 carbonyl group of ChlD1. The cryo-EM density of the FRL-PSII (PDB ID: 7SA3, EMD-24943)12 is consistent with the water positions obtained from MD simulations (see also main text). (E) Occupancy of the W2 position in each of the MD models. In the Chl d/f @ PD2 simulations W2 escapes the binding pocket. |
The PD2 pigment is stabilized through a hydrogen-bonding interaction to D2-Tyr191 via the C132 group in the FRL-PSII model (D2-Trp191 in WL-PSII), as well as stacking interactions with PD1. To this end, no hydrogen-bonding interactions are seen between the protein and the formyl groups Chl d or f when they occupy the PD2 site (Fig. S10A†). The formyl group of Chl d, when in the PD2 location, is highly dynamic in the MD simulation, sampling multiple rotameric states. In contrast, the formyl group of Chl f remains in a single rotameric conformation during the MD simulations (Fig. S10B†).
In addition to these local interactions, we observe differences in the dynamics of the hydrogen-bonding network around the ChlD1, especially for a specific water molecule, W2. In the MD simulations W2 forms a hydrogen-bonding link (67% occupancy) between W1, the water that is the axial ligand to Mg2+ of ChlD1, and its carbonyl (C132–COOCH3) substituent (Fig. 4D and E). W2 is structurally conserved at the ChlD1 site, and its density is well-resolved in both the WL-PSII and FRL-PSII.5,12,34 Modifications of the coordinating water molecule and its hydrogen-bonding network can affect both optical and redox properties of the chlorophyll (Table S3 and Fig. S10†).35–37
In contrast, for our MD simulations of Chl d and Chl f at the PD2 site, W2 near ChlD1 escapes from the pigment pocket, possibly due to steric effects, repositioning D2-Leu205 and the propionate group of ChlD1 that could affect the local hydration (Fig. 4). These changes are further supported by the smaller sidechain of D1-Leu173 in the FRL-PSII relative to D1-Met172 in the WL-PSII (Fig. S10D†).
Our MD simulations suggest that both Chl d and Chl f at ChlD1 lead to subtle structural changes within the RC. We observe a small increase in the PD1/PD2 inter-pigment distance (3.65 Å average, edge-to-edge) with Chl d at the ChlD1 site relative to the WL-PSII (3.53 Å) (Fig. S10A†) that could influence the electronic coupling (Fig. S13†), charge separation, and electron–hole localization following the charge-separation process. For Chl f at the ChlD1 site, the PD1/PD2 inter-pigment distance increases to 3.73 Å, most likely due to steric effects.
In contrast, our MD simulations with Chl d and Chl f at PD2, revealed only minor perturbation in the inter-pigment spacing of PD1/PD2 (3.55 Å for both Chl a/f and Chl a/d composition). Nevertheless, our calculations suggest that the presence of the long wavelength chlorophyll would affect both the excitonic coupling (Fig. S13†) and redox properties of the pigments (see section below; Table S9†). Although coarse-grained MD simulations have been used to model PSII systems on longer timescales,38 the atomic detail provided by atomistic MD simulations are central for describing the detailed pigment surroundings.
To probe how the Chl d substitution at the ChlD1 site affects the charge separation in the FRL-PSII, we estimated redox potentials based on Poisson–Boltzmann electrostatic calculations coupled with Monte Carlo sampling (PBE/MC), and probed the charge separation kinetics using a kinetic master equation model based on electron transfer theory (see Extended methods†).
Our PBE/MC calculations suggested that the oxidation potential of the ChlD1 site increases by 145 mV in the FRL-PSII model with Chl d and by 198 mV with Chl f relative to the WL-PSII model (Fig. S16†). Following the excitation forming the state at 1.72 eV, the charge separation in the FRL-PSII is thus expected to result in a ChlD1+/PheoD1− state at 1.57 eV with Chl d at the ChlD1 site, computed based on the shifts in redox potentials relative to the WL-PSII. To this end, the PBE/MC calculations suggest that the Chl d modifications at ChlD1, stabilize the charge-separated states due to the hydrogen-bonding effects to surrounding residues by ca. 0.15 eV relative to the WL-PSII system (Fig. 5A), and could thus compensate for the 0.10 eV lower photon energy in the FRL-PSII, while maintaining a significant driving force needed for the water oxidation and quinone reduction.
![]() | ||
Fig. 5 Probing the charge separation mechanism of FRL-PSII. (A) Effect of the chlorophyll substitution at the ChlD1 site on the charge separation energetics based on redox calculations and experimental data. (B) CT kinetics and state populations along the ChlD1+PheoD1−-mediated pathway for the FRL-PSII (Chl d @ ChlD1; bright colors) and WL-PSII (transparent colors). (C) Kinetics of charge transfer: light-driven CT at ChlD1 leads to a ChlD1+/PheoD1− state, whilst the subsequent electron transfer leads to reduction of QB and oxidation of the Mn4O5Ca cluster. The arrows indicate the order of events predicted by our kinetic model. The charge transfer rates reported correspond to the FRL system (with Chl d at the ChlD1 site), whilst rates that are unique for the WL-PSII model, are shown in grey. See Fig. S15, S16 and Table S10 and ESI, Extended methods for detailed models.† |
Interestingly, our PBE/MC models suggest that with Chl f at ChlD1, the ChlD1+Pheo− charge separated-state becomes more stabilized possibly due to its intrinsically higher oxidation potential relative to Chl d and a (Fig. 5A and S16†). The shifted redox properties result in less driving force to stabilize the subsequent electron transfer steps and could favor re-combination and harmful singlet oxygen production, as it is energetically closer to the Mnox/QB− state (1.25 eV).
Our kinetic models based on the determined energetics (Fig. 5 and S15†), predict that the ChlD1+Pheo−-mediated charge separation pathway is fast, with a halftime reaching QB reduction/Mn oxidation within ca. 300 μs (cf. ref. 16 and 39), which is rate-limited by the slow QA → QB reaction.40 The shifted redox properties at ChlD1 thus support a charge separation kinetics comparable to that in the WL-PSII, despite the lower photon energy and shifted redox properties.
When Chl f occupies the ChlD1 site, our predictions also showed redshifted features that could be consistent with the experimental spectra. However, Chl f in the ChlD1 site is not capable of binding to the far-red specific, conserved hydrogen-bonding residues, and the presence of which resulted in a small blueshift of the Chl f Qy absorption and a lowered binding constant.
When Chl d and Chl f were modeled at the PD2 site, it led to a large redshift of the longest wavelength absorption significantly beyond the experimentally observed range (Fig. 2A and S4A†). Moreover, these models did not show similar structural stabilization as for the chemical substitution at the ChlD1 site, as there were no obvious conserved changes in the FRL-PSII that could accommodate, either sterically or by the provision of hydrogen-bonds, the additional formyl groups present on Chl d or f. In addition, the presence of Chl f at the PD2 site perturbs the structure of the pigment cluster resulting in the loss of the W2 water that plays a role in the hydrogen-bond network associated with the W1 water that coordinates the ChlD1. Changes in this hydrogen-bonding network can result in spectroscopic and functional changes making the PSII susceptible to photodamage.36,37 Importantly, the charge transfer kinetics in PSII have been extensivity studied (cf. ref. 6, 16 and 39 and ref. therein). Our kinetic models are consistent with the observed kinetics for the WT-PSII, and reproduce also key observations of the far-red light system (see ref. 13), with redox tuning effects favoring the initial charge separation.
Taken together, our combined findings thus indicate that the Chroococcidiopsis thermalis FRL-PSII has a long wavelength chlorophyll at the ChlD1 site, with Chl d reproducing key electronic and structural features, and support initial arguments for the ChlD1 location for the long wavelength pigments (see Introduction).7,9,12 The employed multiscale treatment provides an accurate description of the electronic structure of the pigments, but at a high computational cost. In this regard, recently developed machine learning-based methods 41,42 could provide cost-effective alternatives for the prediction of optical spectra and site energies (cf. also ref. 43).
Our excited-state calculations and analysis of transition orbitals (Fig. S6 and S7†) showed that a long wavelength chlorophyll at the ChlD1 site gave rise to a low-energy Qy excitation that could act as the primary donor and result in a ChlD1+/PheoD1− state (Fig. 5A). This primary radical pair state transfers the cation radical (the electron hole) to PD1/PD2 pair and transfers the electron to QA, creating a transient charge-separated (PD1PD2)+/QA− state, with a t1/2 of 0.2 ns (Fig. 5B), as suggested by our kinetic model, and with subsequent electron transfer steps forming and eventually MnoxQB− as the final state (Fig. 5C). Interestingly, the chemical substitution of ChlD1 compensates for the lower photon energy by the tuned redox properties that could further favor the charge transfer towards the QB sites (Fig. 5A). In this regard, recent 2D-electronic spectroscopy (2D-ES) studies of the FRL-PSII44 found experimental support for a similar (PD1PD2)+/ChlD1− mediated pathway, although this charge separation was assumed to be initiated from a state in which PheoD1 is pre-reduced. Our kinetic models indicate that Chl a at the ChlD1 site (and Chl d/f in PD2) lacks the driving force necessary for initial charge separation in the FRL-PSII, while for Chl f at this site, the ChlD1+/PheoD1− state is energetically further stabilized, and could potentially promote electron back transfer and charge recombination, which in turn leads to the formation of triplet chlorophyll and the production of harmful singlet oxygen. Chl d is unique among the pigments as it supports reasonable charge separation energetics at the ChlD1 site, and facilitates fast charge transfer kinetics, which are comparable to the WL-PSII according to our kinetic model. In this regard, our electrostatic calculations, suggest that the charge-separated ChlD1+/PheoD1− pair is stabilized by ca. 0.15 eV to compensate for the ca. 0.1 eV lower photon energy. Although our current models reveal distinct differences between the FRL- and WL-PSII, elucidating the exact mechanistic details of the charge transfer kinetics will require further detailed structure-based electron transfer calculations.
In contrast, when either Chl d or f occupies the PD2 site, our calculations suggest that the excitation energy is likely to be trapped within the RC (Fig. S15, Table S5 and S6†). This retention implies two possible scenarios for functional PSII: either an energetically uphill and kinetically slow energy transfer to ChlD1 – the putative site of primary charge separation in the WL-PSII (see above; and ref. 31, 32 and 45–47) – or a direct charge separation event involving PD1/PD2. While modification of the pigment at the PD2 site would indeed expand the spectral range of PSII (Fig. 2 and S4†), we speculate that the putative (PD1PD2)+/ChlD1− state with Chl f at PD2 would result in a charge recombination.
Our calculations suggest that the specific FRL-PSII surroundings around ChlD1 support the redshifting effect of Chl d that could enhance its intrinsic light-capturing properties at the far-red-light limit, while Chl f at the same position led to blueshifting effects on the intrinsic low-energy Qy state. We could link these spectral shifts to protein-induced electrostatic and hydrogen-bonding effects, leading to an energetic preference of the Chl d at the ChlD1 site. This effect arises specifically via the FRL-PSII specific D1-Tyr120 and D1-Thr155 residues that form hydrogen-bonding interactions with Chl d. In summary, considering both the pigment–protein interactions at the ChlD1 site, spectral tuning effects, as well as possible bioenergetic advantages in the charge separation process (Fig. 5), our calculations strongly support that Chl d at the ChlD1 site enables the FRL-PSII to capture low energy photons beyond the far-red light limit.
Footnotes |
† Electronic supplementary information (ESI) available. See DOI: https://doi.org/10.1039/d4sc07473d |
‡ These authors contributed equally. |
This journal is © The Royal Society of Chemistry 2025 |