Structural Factors That Alter the Redox Potential of Quinones in Cyanobacterial and Plant Photosystem I

Absence of electron-transfer-associated changes in the time-dependent X-ray free-electron laser structures of the photosynthetic reaction center

Using the X-ray free-electron laser (XFEL) structures of the photosynthetic reaction center from Blastochloris viridis that show light-induced time-dependent structural changes (Dods et al., (2021) Nature 589, 310-314), we investigated time-dependent changes in the energetics of the electron-transfer pathway, considering the entire protein environment of the protein structures and titrating the redox-active sites in the presence of all fully equilibrated titratable residues. In the dark and charge separation intermediate structures, the calculated redox potential (Em) values for the accessory bacteriochlorophyll and bacteriopheophytin in the electron-transfer-active branch (BL and HL) are higher than those in the electron-transfer-inactive branch (BM and HM). However, the stabilization of the charge-separated [PLPM]•+HL•- state owing to protein reorganization is not clearly observed in the Em(HL) values in the charge-separated 5 ps ([PLPM]•+HL•- state) structure. Furthermore, the expected chlorin ring deformation upon formation of HL•- (saddling mode) is absent in the HL geometry of the original 5 ps structure. These findings suggest that there is no clear link between the time-dependent structural changes and the electron-transfer events in the XFEL structures.

Redox Potentials of Iron-Sulfur Clusters in Type I Photosynthetic Reaction Centers

The electron transfer pathways in type I photosynthetic reaction centers, such as photosystem I (PSI) and reaction centers from green sulfur bacteria (GsbRC), are terminated by two Fe4S4 clusters, FA and FB. The protein structures are the basis of understanding how the protein electrostatic environment interacts with the Fe4S4 clusters and facilitates electron transfer. Using the protein structures, we calculated the redox potential (Em) values for FA and FB in PSI and GsbRC, solving the linear Poisson-Boltzmann equation. The FA-to-FB electron transfer is energetically downhill in the cyanobacterial PSI structure, while it is isoenergetic in the plant PSI structure. The discrepancy arises from differences in the electrostatic influences of conserved residues, including PsaC-Lys51 and PsaC-Arg52, located near FA. The FA-to-FB electron transfer is slightly downhill in the GsbRC structure. Em(FA) and Em(FB) exhibit similar levels upon isolation of the membrane-extrinsic PsaC and PscB subunits from the PSI and GsbRC reaction centers, respectively. The binding of the membrane-extrinsic subunit at the heterodimeric/homodimeric reaction center plays a key role in tuning Em(FA) and Em(FB).

Energetic Diversity in the Electron-Transfer Pathways of Type I Photosynthetic Reaction Centers

Photosynthetic reaction centers from heliobacteria (HbRC) and green sulfur bacteria (GsbRC) are homodimeric proteins and share a common ancestor with photosystem I (PSI), classified as type I reaction centers. Using the HbRC crystal structure, we calculated the redox potential (Em) values in the electron-transfer branches, solving the linear Poisson-Boltzmann equation and considering the protonation states of all titratable sites in the entire protein-pigment complex. Em(A-1) for bacteriochlorophyll g at the secondary site in HbRC (-1157 mV) is as low as Em(A-1) for chlorophyll a in PSI (-1173 mV). Em(A0/HbRC) is at the same level as Em(A0/GsbRC) and is 200 mV higher than Em(A0/PSI) due to the replacement of PsaA-Trp697/PsaB-Trp677 in PSI with PshA-Arg554 in HbRC. In contrast, Em(FX) for the Fe4S4 cluster in HbRC (-420 mV) is significantly higher than Em(FX) in GsbRC (-719 mV) and PSI (-705 mV) due to the absence of acidic residues that correspond to PscA-Asp634 in GsbRC and PsaB-Asp575 in PSI. It seems likely that type I reaction centers have evolved, adopting (bacterio)chlorophylls suitable for their light environments while maintaining electron-transfer cascades.

Energetics of the Electron Transfer Pathways in the Homodimeric Photosynthetic Reaction Center

Photosynthetic reaction centers from a green sulfur bacterium (GsbRC), the PscA/PscA proteins, and photosystem I (PSI), PsaA/PsaB proteins, share structural similarities. Here, we report the redox potential (Em) values of GsbRC by solving the linear Poisson-Boltzmann equation and considering the protonation states of all titratable sites in the entire GsbRC protein and identify the factors that shift the Em values with respect to PSI. The Em values for one-electron reduction of the accessory (A-1) and adjacent (A0) chlorophylls in GsbRC are 100-250 mV higher than those in PSI, whereas the Em values for the Fe4S4 cluster (FX) are at the same level. The PsaA-Trp697/PsaB-Trp677 pair in PSI, which forms the A1-quinone binding site, is replaced with PscA-Arg638 in GsbRC. PsaB-Asp575 in PSI, which is responsible for the Em difference between A1A and A1B quinones in PSI, is absent in GsbRC. These discrepancies also contribute to the upshift in Em(A-1) and Em(A0) in GsbRC with respect to PSI. It seems likely that the upshifted Em for chlorophylls in GsbRC ultimately originates from the characteristics of the electrostatic environment that corresponds to the A1 site of PSI.

Conformational Changes and H-Bond Rearrangements during Quinone Release in Photosystem II

In photosystem II (PSII) and photosynthetic reaction centers from purple bacteria (PbRC), the electron released from the electronically excited chlorophyll is transferred to the terminal electron acceptor quinone, Q_B. Q_B accepts two electrons and two protons before leaving the protein. We investigated the molecular mechanism of quinone exchange in PSII, conducting molecular dynamics (MD) simulations and quantum mechanical/molecular mechanical (QM/MM) calculations. MD simulations suggest that the release of Q_B leads to the transformation of the short helix (D1-Phe260 to D1-Ser264), which is adjacent to the stromal helix de (D1-Asn247 to D1-Ile259), into a loop and to the formation of a water-intake channel. Water molecules enter the Q_B binding pocket via the channel and form an H-bond network. QM/MM calculations indicate that the H-bond network serves as a proton-transfer pathway for the reprotonation of D1-His215, the proton donor during Q_BH^-/Q_BH₂ conversion. Together with the absence of the corresponding short helix but the presence of Glu-L212 in PbRC, it seems likely that the two type-II reaction centers undergo quinone exchange via different mechanisms.

Mechanism of Mixed-Valence Fe2.5+···Fe2.5+ Formation in Fe4S4 Clusters in the Ferredoxin Binding Motif

Most low-potential Fe₄S₄ clusters exist in the conserved binding sequence CxxCxxC (CⁿCⁿ⁺³Cⁿ⁺⁶). Fe(II) and Fe(III) at the first (Cⁿ) and third (Cⁿ⁺⁶) cysteine ligand sites form a mixed-valence Fe^2.5+···Fe^2.5+ pair in the reduced Fe(II)₃Fe(III) cluster. Here, we investigate the mechanism of how the conserved protein environment induces mixed-valence pair formation in the Fe₄S₄ clusters, F_X, F_A, and F_B in photosystem I, using a quantum mechanical/molecular mechanical approach. Exchange coupling between Fe sites is predominantly determined by the shape of the Fe₄S₄ cluster, which is stabilized by the preorganized protein electrostatic environment. The backbone NH and CO groups in the conserved CxxCxxC and adjacent helix regions orient along the Fe^Cn···Fe^C(n+6) axis, generating an electric field and stabilizing the Fe^Cn(II)Fe^C(n+6)(III) state in F_A and F_B. The overlap of the d orbitals via -S- (superexchange) is observed for the single Fe^Cn(II)···Fe^C(n+6)(III) pair, leading to the formation of the mixed-valence Fe^2.5+···Fe^2.5+ pair. In contrast, several superexchange Fe(II)···Fe(III) pairs are observed in F_X due to the highly symmetric pair of the CDGPGRGGTC sequences. This is likely the origin of F_X serving as an electron acceptor in the two electron transfer branches.

Room temperature XFEL crystallography reveals asymmetry in the vicinity of the two phylloquinones in photosystem I

Photosystem I (PS I) has a symmetric structure with two highly similar branches of pigments at the center that are involved in electron transfer, but shows very different efficiency along the two branches. We have determined the structure of cyanobacterial PS I at room temperature (RT) using femtosecond X-ray pulses from an X-ray free electron laser (XFEL) that shows a clear expansion of the entire protein complex in the direction of the membrane plane, when compared to previous cryogenic structures. This trend was observed by complementary datasets taken at multiple XFEL beamlines. In the RT structure of PS I, we also observe conformational differences between the two branches in the reaction center around the secondary electron acceptors A_1A and A_1B. The π-stacked Phe residues are rotated with a more parallel orientation in the A-branch and an almost perpendicular confirmation in the B-branch, and the symmetry breaking PsaB-Trp673 is tilted and further away from A_1A. These changes increase the asymmetry between the branches and may provide insights into the preferential directionality of electron transfer.

Electron Acceptor-Donor Iron Sites in the Iron-Sulfur Cluster of Photosynthetic Electron-Transfer Pathways

In photosystem I, two electron-transfer pathways via quinones (A_1A and A_1B) are merged at the iron-sulfur Fe₄S₄ cluster F_X into a single pathway toward the other two Fe₄S₄ clusters F_A and F_B. Using a quantum mechanical/molecular mechanical approach, we identify the redox-active Fe sites in the clusters. In F_A and F_B, the Fe site, which does not belong to the CxxCxxCxxxCP motif, serves as an electron acceptor/donor. F_X has two independent electron acceptor Fe sites for A- and B-branch electron transfers, depending on the Asp-B575 protonation state, which causes the A_1A-to-F_X electron transfer to be uphill and the A_1B-to-F_X electron transfer to be downhill. The two asymmetric electron-transfer pathways from A₁ to F_X and the separation of the electron acceptor and donor Fe sites are likely associated with the specific role of F_X in merging the two electron transfer pathways into the single pathway.

Phylloquinone is the principal Mehler reaction site within photosystem I in high light

Photosynthesis is a vital process, responsible for fixing carbon dioxide, and producing most of the organic matter on the planet. However, photosynthesis has some inherent limitations in utilizing solar energy, and a part of the energy absorbed is lost in the reduction of O2 to produce the superoxide radical (O2•-) via the Mehler reaction, which occurs principally within photosystem I (PSI). For decades, O2 reduction within PSI was assumed to take place solely in the distal iron-sulfur clusters rather than within the two asymmetrical cofactor branches. Here, we demonstrate that under high irradiance, O2 photoreduction by PSI primarily takes place at the phylloquinone of one of the branches (the A-branch). This conclusion derives from the light dependency of the O2 photoreduction rate constant in fully mature wild-type PSI from Chlamydomonas reinhardtii, complexes lacking iron-sulfur clusters, and a mutant PSI, in which phyllosemiquinone at the A-branch has a significantly longer lifetime. We suggest that the Mehler reaction at the phylloquinone site serves as a release valve under conditions where both the iron-sulfur clusters of PSI and the mobile ferredoxin pool are highly reduced.

Conserved residue PsaB-Trp673 is essential for high-efficiency electron transfer between the phylloquinones and the iron-sulfur clusters in Photosystem I

Despite the high level of symmetry between the PsaA and PsaB polypeptides in Photosystem I, some amino acids pairs are strikingly different, such as PsaA-Gly693 and PsaB-Trp673, which are located near a cluster of 11 water molecules between the A_1A and A_1B quinones and the F_X iron-sulfur cluster. In this work, we changed PsaB-Trp673 to PsaB-Phe673 in Synechocystis sp. PCC 6803. The variant contains ~ 85% of wild-type (WT) levels of Photosystem I but is unable to grow photoautotrophically. Both time-resolved and steady-state optical measurements show that in the PsaB-W673F variant less than 50% of the electrons reach the terminal iron-sulfur clusters F_A and F_B; the majority of the electrons recombine from A_1A^- and A_1B^-. However, in those reaction centers which pass electrons forward the transfer is heterogeneous: a minor population shows electron transfer rates from A_1A^- and A_1B^- to F_X slightly slower than that of the WT, whereas a major population shows forward electron transfer rates to F_X slowed to the ~ 10 µs time range. Competition between relatively similar forward and backward rates of electron transfer from the quinones to the F_X cluster account for the relatively low yield of long-lived charge separation in the PsaB-W673F variant. A higher water content and its increased mobility observed in MD simulations in the interquinone cavity of the PsaB-W673F variant shifts the pK of PsaB-Asp575 and allows its deprotonation in situ. The heterogeneity found may be rooted in protonation state of PsaB-Asp575, which controls whether electron transfer can proceed beyond the phylloquinone cofactors.

Assessment of the orientation and conformation of pigments in protein binding sites from infrared difference spectra

Time resolved FTIR difference spectroscopy (DS) has been used to study photosystem I (PSI) with the disubstituted 1,4-naphthoquinones acequinocyl (AcQ) and lapachol (Lpc) incorporated into the A₁ binding site. AcQ is a 2-acetoxy-3-dodecyl-1,4-naphthoquinone, Lpc is a 2-hydroxy-3-(3-methyl-2-butenyl)-1,4-naphthoquinone. To assess whether the experimental spectra are specific to different orientations of the quinone and their substitutions ONIOM-type QM/MM vibrational frequency calculations were undertaken for various orientations of the pigments and side-chain conformations in the A₁ binding site. Comparison of calculated and experimental spectra for the reduced species (semiquinone anion) suggests that the orientation for the naphthoquinone ring in the binding site and specific side-chain conformations can be identified based on the spectra. In native PSI phylloquinone (PhQ) in the A₁ binding site binds with its phytyl chain ortho to the hydrogen bonded carbonyl group. This is not found to be the case for the hydrocarbon tail of AcQ, which is meta to the H-bonded carbonyl group. In contrast, Lpc in PSI binds with its hydrocarbon tail also ortho to the H-bonded carbonyl group. Furthermore, comparison of calculated and experimental spectra indicates which conformations the acetoxy group of AcQ and the hydroxy group of Lpc adopt in the A₁ binding site.

Proton transfer pathway from the oxygen-evolving complex in photosystem II substantiated by extensive mutagenesis

We report a structure-based biological approach to identify the proton-transfer pathway in photosystem II. First, molecular dynamics (MD) simulations were conducted to analyze the H-bond network that may serve as a Grotthuss-like proton conduit. MD simulations show that D1-Asp61, the H-bond acceptor of H₂O at the Mn₄CaO₅ cluster (W1), forms an H-bond via one water molecule with D1-Glu65 but not with D2-Glu312. Then, D1-Asp61, D1-Glu65, D2-Glu312, and the adjacent residues, D1-Arg334, D2-Glu302, and D2-Glu323, were thoroughly mutated to the other 19 residues, i.e., 114 Chlamydomonas chloroplast mutant cells were generated. Mutation of D1-Asp61 was most crucial. Only the D61E and D61C cells grew photoautotrophically and exhibit O₂-evolving activity. Mutations of D2-Glu312 were less crucial to photosynthetic growth than mutations of D1-Glu65. Quantum mechanical/molecular mechanical calculations indicated that in the PSII crystal structure, the proton is predominantly localized at D1-Glu65 along the H-bond with D2-Glu312, i.e., pK_a(D1-Glu65) > pK_a(D2-Glu312). The potential-energy profile shows that the release of the proton from D1-Glu65 leads to the formation of the two short H-bonds between D1-Asp61 and D1-Glu65, which facilitates downhill proton transfer along the Grotthuss-like proton conduit in the S₂ to S₃ transition. It seems possible that D1-Glu65 is involved in the dominant pathway that proceeds from W1 via D1-Asp61 toward the thylakoid lumen, whereas D2-Glu312 and D1-Arg334 may be involved in alternative pathways in some mutants.

Reversible inhibition and reactivation of electron transfer in photosystem I

In photosystem I (PSI) complexes at room temperature electron transfer from A₁^- to F_X is an order of magnitude faster on the B-branch compared to the A-branch. One factor that might contribute to this branch asymmetry in time constants is TrpB673 (Thermosynechococcus elongatus numbering), which is located between A_1B and F_X. The corresponding residue on the A-branch, between A_1A and F_X, is GlyA693. Here, microsecond time-resolved step-scan FTIR difference spectroscopy at 77 K has been used to study isolated PSI complexes from wild type and TrpB673Phe mutant (WB673F mutant) cells from Synechocystis sp. PCC 6803. WB673F mutant cells require glucose for growth and are light sensitive. Photoaccumulated FTIR difference spectra indicate changes in amide I and II protein vibrations upon mutation of TrpB673 to Phe, indicating the protein environment near F_X is altered upon mutation. In the WB673F mutant PSI samples, but not in WT PSI samples, the phylloquinone molecule that occupies the A₁ binding site is likely doubly protonated following long periods of repetitive flash illumination at room temperature. PSI with (doubly) protonated quinone in the A₁ binding site are not functional in electron transfer. However, electron transfer functionality can be restored by incubating the light-treated mutant PSI samples in the presence of added phylloquinone.

Kinetics and Energetics of Phylloquinone Reduction in Photosystem I: Insight From Modeling of the Site Directed Mutants

Two phylloquinone molecules (A ₁), one being predominantly coordinated by PsaA subunit residues (A _1A) the other by those of PsaB (A _1B), act as intermediates in the two parallel electron transfer chains of Photosystem I. The oxidation kinetics of the two phyllosemiquinones by the iron-sulfur cluster F_X differ by approximately one order of magnitude, with A 1 A - being oxidized in about 200 ns and A 1 B - in about 20 ns. These differences are generally explained in terms of asymmetries in the driving force for F_X reduction on the two electron transfer chains. Site directed mutations of conserved amino acids composing the A ₁ binding site have been engineered on both reaction center subunits, and proved to affect selectively the oxidation lifetime of either A 1 A - , for PsaA mutants, or A 1 B - , for PsaB mutants. The mutation effects are here critically reviewed, also by novel modeling simulations employing the tunneling formalism to estimate the electron transfer rates. Three main classes of mutation effects are in particular addressed: (i) those leading to an acceleration, (ii) those leading to a moderated slowing (~5-folds), and (iii) those leading to a severe slowing (>20-folds) of the kinetics. The effect of specific amino acid perturbations contributing to the poising of the phylloquinones redox potential and, in turn, to PSI functionality, is discussed.

Modelling electron transfer in photosystem I: limits and perspectives

Uncovering the parameters underlying the electron transfer (ET) in photosynthetic reaction centres is of importance for understanding the molecular mechanisms underpinning their functionality. The reductive nature of most cofactors involved in photosynthetic ET makes the direct estimation of their properties difficult. Photosystem I (PSI) operates in a highly reducing regime, making the assessment of cofactor properties even more difficult. Kinetic modelling coupled to a non-adiabatic description of ET is a useful approach in overcoming this hindrance. Here we review the theory and modelling approaches that have been used in assessing parameters associated with ET reactions in PSI, with particular attention to ET reactions involving the phylloquinones and the iron-sulphur clusters. In most modelling studies, the goal is to estimate the driving force of ET, which is usually associated with the cofactor midpoint potentials. The driving force is sensitive to many factors, which define the ET rate, i.e. the reorganisation energy, the coupling with nuclear modes and the electronic matrix elements, which are explored and discussed here. The importance of an inclusive modelling of both forward and reverse ET processes is discussed and highlighted. It is shown that although estimates are indeed sensitive to the exact parameter sets employed in the modelling, a general consensus is still attained, pointing to a scenario where Δ G A 1 A → F X 0 / Δ G A 1 B → F X 0 is weakly endergonic/exergonic, respectively. It is emphasised that to further refine those estimates, it will require a joint effort between computational modelling and more wide-ranging experimental studies.

Energetic insights into two electron transfer pathways in light-driven energy-converting enzymes

We report E_m values of (bacterio-)chlorophylls for one-electron reduction in both electron-transfer branches of PbRC, PSI, and PSII.

We report redox potentials (E_m) for one-electron reduction for all chlorophylls in the two electron-transfer branches of water-oxidizing enzyme photosystem II (PSII), photosystem I (PSI), and purple bacterial photosynthetic reaction centers (PbRC). In PSI, E_m values for the accessory chlorophylls were similar in both electron-transfer branches. In PbRC, the corresponding E_m value was 170 mV less negative in the active L-branch (B_L) than in the inactive M-branch (B_M), favoring B_L˙^– formation. This contrasted with the corresponding chlorophylls, Chl_D1 and Chl_D2, in PSII, where E_m(Chl_D1) was 120 mV more negative than E_m(Chl_D2), implying that to rationalize electron transfer in the D1-branch, Chl_D1 would need to serve as the primary electron donor. Residues that contributed to E_m(Chl_D1) < E_m(Chl_D2) simultaneously played a key role in (i) releasing protons from the substrate water molecules and (ii) contributing to the larger cationic population on the chlorophyll closest to the Mn₄CaO₅ cluster (P_D1), favoring electron transfer from water molecules. These features seem to be the nature of PSII, which needs to possess the proton-exit pathway to use a protonated electron source—water molecules.

Electron Transfer through the Acceptor Side of Photosystem I: Interaction with Exogenous Acceptors and Molecular Oxygen

This review considers the state-of-the-art on mechanisms and alternative pathways of electron transfer in photosynthetic electron transport chains of chloroplasts and cyanobacteria. The mechanisms of electron transport control between photosystems (PS) I and II and the Calvin-Benson cycle are considered. The redistribution of electron fluxes between the noncyclic, cyclic, and pseudocyclic pathways plays an important role in the regulation of photosynthesis. Mathematical modeling of light-induced electron transport processes is considered. Particular attention is given to the electron transfer reactions on the acceptor side of PS I and to interactions of PS I with exogenous acceptors, including molecular oxygen. A kinetic model of PS I and its interaction with exogenous electron acceptors has been developed. This model is based on experimental kinetics of charge recombination in isolated PS I. Kinetic and thermodynamic parameters of the electron transfer reactions in PS I are scrutinized. The free energies of electron transfer between quinone acceptors A_1A/A_1B in the symmetric redox cofactor branches of PS I and iron-sulfur clusters F_X, F_A, and F_B have been estimated. The second-order rate constants of electron transfer from PS I to external acceptors have been determined. The data suggest that byproduct formation of superoxide radical in PS I due to the reduction of molecular oxygen in the A₁ site (Mehler reaction) can exceed 0.3% of the total electron flux in PS I.