New insights Into the evolution of the electron transfer from cytochrome f to photosystem I in the green and red branches of photosynthetic eukaryotes

New Insights into the Evolution of the Electron Transfer from Cytochrome fto Photosystem I in the Green and Red Branches of Photosynthetic Eukaryotes Carmen Castell 1 , Luis A. Rodríguez-Lumbreras 2 , Manuel Hervás 1 , Juan Fernández-Recio 2 and José A. Navarro 1, * 1 Instituto de Bioquímica Vegetal y Fotosíntesis, CSIC and Universidad de Sevilla, cicCartuja, Sevilla, Spain 2 Instituto de Ciencias de la Vid y del Vino (ICVV), CSIC—Universidad de La Rioja—Gobierno de La Rioja, Logro~ no, Spain *Corresponding author: Email, [email protected]; Fax, þ34 954 46 01 65. (Received 10 February 2021; Accepted 15 March 2021) In cyanobacteria and most green algae of the eukaryotic green lineage, the copper-protein plastocyanin (Pc) alternatively replaces the heme-protein cytochrome c 6 (Cc 6 ) as the soluble electron carrier from cytochrome f(Cf) to photosystem I (PSI). The functional and structural equivalence of ‘green’ Pc and Cc 6 has been well established, representing an example of convergent evolution of two unrelated proteins. However, plants only produce Pc, despite having evolved from green algae. On the other hand, Cc 6 is the only soluble donor available in most species of the red lineage of photosynthetic organisms, which includes, among others, red algae and diatoms. Interestingly, Pc genes have been identified in oceanic diatoms, probably acquired by horizontal gene transfer from green algae. However, the mechanisms that regulate the expression of a functional Pc in diatoms are still unclear. In the green eukaryotic lineage, the transfer of electrons from Cf to PSI has been characterized in depth. The conclusion is that in the green lineage, this process involves strong electrostatic interactions between partners, which ensure a high affinity and an efficient electron transfer (ET) at the cost of limiting the turnover of the process. In the red lineage, recent kinetic and structural modeling data suggest a different strategy, based on weaker electrostatic interactions between partners, with lower affinity and less efficient ET, but favoring instead the protein exchange and the turnover of the process. Finally, in diatoms the interaction of the acquired green-type Pc with both Cf and PSI may not yet be optimized. Keywords: Cytochrome c 6 •Cytochrome f•Electron transfer •Photosynthetic green and red lineages •Photosystem I • Plastocyanin. Introduction Oxygenic photosynthesis is a crucial process in the acquisition and storage of solar energy in the biosphere. Globally, it is estimated that photosynthesis captures 310 21 J year 1 and is responsible for the annual assimilation of 100 Gt of carbon (Hall 1976,Field et al. 1998). Photosynthetic productivity is shared roughly equally between terrestrial systems and the oceans, making photosynthesis a global process spread across the entire Earth's surface. At the molecular level, the productivity of photosynthetic organisms is related to the efficiency of their photosynthetic electron transfer (ET) chain, in which photo-induced electron transport is coupled to the generation of both energy-rich compounds and a proton motive force. The main features of the sophisticated machinery that performs oxygenic photosynthesis are preserved from prokaryotic cyanobacteria to plants. Thus, in the photosynthetic chain, three large membrane complexes coexist: photosystem II (PSII), cytochrome b 6 f(Cb 6 f) and photosystem I (PSI), which are connected by mobile electron transporters (Hervás et al. 2003,Blankenship 2014)(Fig. 1). In plants, ET between the Cb 6 f complex and PSI is carried out in the thylakoid lumen of the chloroplast solely by the copperprotein plastocyanin (Pc), whereas the output of electrons from PSI is directed to the iron-sulfur protein ferredoxin (Fd) (Fig. 1A). However, this situation must be taken as an end point in the evolutionary path of the green lineage that, starting from cyanobacteria (considered the ancestors of the eukaryotic chloroplast) and green algae, leads to plants (Fig. 2A) (Falkowski et al. 2004,De la Rosa et al. 2006,Sétif 2006). Conversely, in aquatic systems, photosynthetic microorganisms have maintained alternative redox protein pairs, which can be understood as an adaptation to fluctuations and limitations in the bioavailability of metals in these ecosystems (specifically iron (Fe) and copper (Cu)) (De la Rosa et al. 2006,Moore and Braucher 2008,Nouet et al. 2011,Marchetti et al. 2012). In particular, in most cyanobacteria and unicellular green algae, Pc acts by alternatively replacing the heme-protein cytochrome c 6 (Cc 6 ) (considered the initial ancestral carrier) when copper is available. Furthermore, the flavoprotein flavodoxin (Fld) replaces Fd under iron-limiting conditions (Hervás et al. 2003, Sancho 2006,Sétif 2006)(Fig. 1A). Thus, plants have lost Cc 6 and Fld as alternative proteins along their evolution within the green lineage, probably due to the particularities of the process Plant Cell Physiol. – ©The Author(s) 2021. Published by Oxford University Press on behalf of Japanese Society of Plant Physiologists. This is an Open Access article distributed under the terms of the Creative Commons Attribution Non-Commercial License (https:// creativecommons. org/ licenses/ bync/ 4. 0/ ), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the original work is properly cited. For commercial reuse, please contact [email protected] Special Issue - Review 62(7): 1082 1093 doi:https://doi.org/10.1093/pcp/pcab044, Advance Access publication on 27 March 2021, available online at https://academic.oup.com/pcp Downloaded from https://academic.oup.com/pcp/article/62/7/1082/6194660 by Universidad de Sevilla user on 10 December 2021 of land colonization with respect to the bioavailability of iron and copper (De la Rosa et al. 2006,Pierella Karlusich et al. 2014). The functional and structural equivalences of Pc and Cc 6 , and of Fd and Fld, have been well established (Hervás et al. 1995, Hervás et al. 2003,Hippler and Drepper 2006,Sétif 2006, Pierella Karlusich et al. 2014,Bendall and Howe 2016), and therefore, both couples represent an example of evolutionary functional convergence of unrelated proteins (De la Rosa et al. 2006). Interestingly, the existence of alternative non-iron proteins is only partially maintained in the red lineage of photosynthetic organisms that diverged from the green lineage along evolution. The red lineage includes, among others, red algae and diatoms (Falkowski et al. 2004)(Fig. 2A). Thus, Fld is present as an alternative to Fd in most algal taxa of the red lineage, including red algae, cryptophytes, stramenopiles (that include diatoms) and haptophytes (Pierella Karlusich et al. 2014). However, Cc 6 is the only soluble donor available in red algae and most algae of the red-plastid lineage, as evidenced by direct protein characterization, genome sequences and transcriptomic analysis of different red algae and stramenopiles, including diatoms (Akazaki et al. 2009,Bowler et al. 2010,Groussman et al. 2015). These works have confirmed a widespread predominance of the petJ gene encoding Cc 6 and have validated this protein as the electron carrier from Cb 6 f to PSI in the majority of organisms within the red-plastid lineage (Blaby-Haas and Merchant 2012,Groussman et al. 2015). However, although all the stramenopiles whose genomes have been sequenced have genes for Cc 6 , genome sequencing and metatranscriptomics analysis have also identified petE genes, encoding Pc, in several species of oceanic diatoms. This has been explained as gene acquisition from green alga by horizontal gene transfer, probably as an adaptation to the limitation of iron in their natural habitats (Peers and Price 2006,Marchetti et al. 2012, Groussman et al. 2015,Hippmann et al. 2017). In fact, a particular strain of the open ocean diatom Thalassiosira oceanica appears to constitutively replace Cc 6 by Pc, as the two genes that encode Cc 6 in its genome are weakly expressed, presumably reducing its overall need for iron (Peers and Price 2006, Lommer et al. 2012,Kong and Price 2020). Finally, it is interesting to note that although the Fd/Fld pair has well-conserved common functional characteristics in the vast majority of photosynthetic organisms (Pierella Karlusich et al. 2014), the Cc 6 /Pc pair presents significant variations in different types of organisms, both in the mechanism of interaction with its partners and in their interaction areas (DíazQuintana et al. 2008). Although the functional and structural properties of the ‘green’Cc 6 /Pc pair—and the comparison with equivalent cyanobacterial systems—have been characterized in depth (Hippler and Drepper 2006,Díaz-Quintana et al. 2008,Bendall and Howe 2016), recent works have provided new data on the ET mechanism and protein features in the red-plastid lineage (Bernal-Bayard et al. 2013,Bernal-Bayard et al. 2015,Roncel et al. 2016,Antoshvili et al. 2019,Castell et al. 2021). These data suggest that different strategies have beendevelopedintheevolutionoftheETfromCb 6 ftoPSIin the green and red branches of eukaryotic photosynthetic organisms. Interaction of the Soluble Carriers with PSI The interaction of Pc and Cc 6 with PSI has been analyzed in vitro and in vivo in several green algae and plants (reviewed in De la Rosa et al. 2006,Hippler and Drepper 2006), and more Fig. 1 (A) Photosynthetic chain and electron transfer pathways (shown by arrows) based on the existence of alternative soluble electron carrier proteins. PQ, plastoquinone; PsaF, PSI subunit involved—in green algae and plants—in the electrostatic interaction with Pc and Cc 6 , the two alternative soluble electron carriers between Cb 6 f and PSI. Backbone representation of photosynthetic proteins of the green alga Chlamydomonas reinhardtii:(B)luminalexposedpartofCf intheCb 6 fcomplex(ProteinDataBank, PDBcode, 1cfm);(C)Pc(PDBcode, 2plt);(D) Cc 6 (PDBcode,1cyj); and (E) luminal exposed part of PsaF in PSI (PDB code, 6ijo). The copper (in Pc) and heme (in Cf and Cc 6 ) cofactors are highlighted. 1083 Plant Cell Physiol. –62(7): 1082 1093 (2021) doi:https://doi.org/10.1093/pcp/pcab044 Downloaded from https://academic.oup.com/pcp/article/62/7/1082/6194660 by Universidad de Sevilla user on 10 December 2021 recently in diatoms (Bernal-Bayard et al. 2013,Bernal-Bayard et al. 2015). In addition to an extensive kinetic characterization of the intermolecular ET processes, structural determinations or models are already available for PSI, Cc 6 and Pc of plants or green algae, as well as for Cc 6 of red algae and diatoms. Structural features Pc is a small (ca. 100 amino acids) protein with a b-barrel tertiary structure that binds a type-1 blue copper cofactor (Fig. 1C). On the other hand, Cc 6 (ca. 90 amino acids) is a typical class I c-type cytochrome with a-helix-based tertiary structure and a heme group (reviewed in De la Rosa et al. 2006)(Fig. 1D). Because anoxygenic photosynthesis uses cytochromes as electron carriers rather than Pc-like proteins, it is considered that Cc 6 evolved as the ancestral carrier in primitive oxygenic organisms. Both Pc and Cc 6 show a well-defined ET area—determined by the solvent accessibility of the heme group in Cc 6 or the copper ligands in Pc—which is surrounded by a hydrophobic surface patch (Guss and Freeman 1983,Fraz~ ao et al. 1995,Inoue et al. 1999,Crowley et al. 2002,Díaz-Moreno et al. 2005; and see Figs. 1C, D,2B). In addition, Cc 6 and Pc show an electrostatically charged surface area at an equivalent location relative to the ET sites and the hydrophobic patches (Fig. 2B)(De la Rosa et al. 2006). The electrostatic nature of this patch defines the type of interaction of Pc and Cc 6 with their partners. Thus, in cyanobacteria, the electrostatic patch can be relatively neutral or even positive, but in Pc and Cc 6 of the eukaryotic green lineage, this electrostatic patch always has a strong negative charge (De la Rosa et al. 2006,Hippler and Drepper 2006) (Fig. 2B). One intriguing exception to this general description is fern Pc, in which the acidic region relocates and surrounds the hydrophobic patch, resulting in distinctive electrostatic properties (Kohzuma et al. 1999). The negative patch of green Cc 6 and Pc has a counterpart in PSI. Biochemical and structural analyses conclude that, in addition to a hydrophobic binding pocket in the PSI core (see below), in green algae and higher plants, there is a strongly positively charged N-terminal region of the PsaF subunit of PSI exposed to the lumen. This N-terminal extension (ca. 70 amino acids; Figs. 1E,2B) is responsible for the bulk electrostatic interaction between PSI and the negatively charged Pc and Cc 6 donors (Farah et al. 1995,Hippler et al. 1996,Hippler et al. 1998,Sommer et al. 2006). In the red lineage, the solved crystal structures of Cc 6 from red algae and diatoms show that this soluble carrier maintains the global fold of the ‘green’Cc 6 (Yamada et al. 2000,Akazaki Fig. 2 (A) Schematic representation of the evolution of the main lineages of photosynthetic eukaryotes and the presence of the electron carriers Pc andCc 6 .In thegreenlineage, plants(thatonlyproducePc) haveevolvedfrom greenalgae,whichcanproducePcandCc 6 .In theredlineage,redalgae only produce Cc 6 and are the origin, through secondary endosymbiosis events, of the red-type chloroplast of other secondary algae groups. (B) Electrostatic potential of the charged surface area—involved in the electrostatic interaction with its partners—of Cc 6 , Pc and the PsaF subunit from the green alga C. reinhardtii, and of the Cc 6 (PDB code, 3dmi) and modeled PsaF (Bernal-Bayard et al. 2015) from the diatom P. tricornutum. The modeled acquired ‘green-type’Pc of the diatom T.oceanica is also shown (in brackets). The proteins are oriented as shown in Fig. 1, with Pc and Cc 6 displayingtheirelectrostaticareas infront andwith thehydrophobicpatchesatthe topof bothproteins. Electrostaticpotentialvaluesareshown on a scale from red to blue, corresponding to 10.0 and þ10.0 kcal mol 1 , respectively, at 298 K. K A and K D , association and dissociation constants, respectively; k ET , electron transfer rate constant. See the text for more details. 1084 C. Castell et al. jElectron transfer from cytochrome fto photosystem I in the green and red branches Downloaded from https://academic.oup.com/pcp/article/62/7/1082/6194660 by Universidad de Sevilla user on 10 December 2021 et al. 2009). However, with respect to protein surface and functional areas, the Cc 6 of red algae and diatoms still preserves the hydrophobic ET area, but the negative character of the electrostatic area is sensibly reduced (Fig. 2B)(Bernal-Bayard et al. 2013). Interestingly, this decrease in the electrostatic properties of Cc 6 corresponds to equivalent changes in the PSI partner as compared with ‘green’systems (Fig. 2B). The recently determined structure of PSI of the red alga Cyanidioschyzon merolae (at 4 Å resolution) (Antoshvili et al. 2019) shows both the conservation of the PSI hydrophobic interaction site and the Nterminal region in PsaF. However, the luminal PsaF extension of this red alga displays a drastic decrease in the number of positive groups, which is also found in other algae of the red lineage, such as diatoms (Fig. 2B). This evidences the existence of significant electrostatic differences in this specific area of PSI between the organisms of the green and red lineages. In conclusion, although the evolution of ET to PSI in red-type organisms has also led to complementary electrostatic interactions between acidic and basic patches in Cc 6 and PsaF, respectively, the electrostatic character of both partners is similarly reduced (Bernal-Bayard et al. 2013,Bernal-Bayard et al. 2015)(Fig. 2B). The red Pc A particular case is the existence of petE genes, encoding Pc, that have been annotated in the genomes of several strains of oceanic diatom species, as T. oceanica and Fragilariopsis cylindrus, as well as in the haptophyte Emiliania huxleyi and the dinoflagellate Karenia brevis (Nosenko et al. 2006,Peers and Price 2006,Blaby-Haas and Merchant 2012). Furthermore, evidence for a functional Pc has been reported in the case of T. oceanica, where the holoprotein was detected and partially purified and sequenced (Peers and Price 2006). However, the intracellular levels of Pc reported in T. oceanica (3 mM Pc, Peers and Price 2006,Kong and Price 2020) are remarkably low compared to those reported for Pc in cyanobacteria, green algae and plants (Navarro et al. 2011), or to the Cc 6 content described in other diatoms under iron sufficient conditions (Cc 6 ca. 28 mM in the coastal diatom Phaeodactylum tricornutum)(Bernal-Bayard et al. 2013,Roncel et al. 2016). On the other hand, the Pc level in T. oceanica is comparable to the Cc 6 levels described in P. tricornutum under iron-limiting conditions (4–7mM) (Roncel et al. 2016,Castell et al. 2021). In any case, the Pctranslated region in diatoms is closely related to the Pc from green algae. Consequently, the ‘red’acquired Pc actually shows the typical strong negative electrostatics of a ‘green-type’Pc (Fig. 2B; and see below) (Castell et al. 2021). The expression of a functional Pc in diatoms raises, however, two questions that remain to be clarified. First, it requires a specific copper-transporting system to the thylakoid that has not yet been found. High-affinity copper uptake systems have been described in T. oceanica (Kong and Price 2019), but a specific thylakoid transport system is also expected to include at least a P-type ATPase transporter located in the thylakoid membrane, which delivers Cu to the lumen to be incorporated into Pc (Guo et al. 2010,Guo et al. 2015,Nouet et al. 2011,BlabyHaas and Merchant 2012,Kong and Price 2019). This is an important point since, in the absence of specific mechanisms of copper import into the thylakoid lumen, the copper available to be incorporated into Pc would probably be limited (Ho et al. 2003,Levy et al. 2008,Twining and Baines 2013,Castell et al. 2021). The second question lies in the regulatory mechanism by which a putative Pc could act as an alternative to Cc 6 in diatoms. In cyanobacteria and green algae, copper levels radically determine the expression of one or another protein through precisely regulated mechanisms, including transcription regulators, metal sensors and proteases that have not been found in diatoms (Merchant et al. 2020,García-Ca~ nas et al. 2021). In addition, it has been reported in diatoms that not only low levels of copper but also sufficient levels of iron induce a drastic reduction in Pc transcripts, which, on the contrary, increase under iron-limiting conditions (Lommer et al. 2012,Marchetti et al. 2012,Hippmann et al. 2017,Rizkallah et al. 2020). Recently, a moderate 3-fold decrease in the Pc content (from 3 to 1 mM) has been reported in T. oceanica under copper deficiency (Kong and Price 2020). On the other hand, Cc 6 levels in diatoms seem to depend mainly on the availability of iron (Roncel et al. 2016, Castell et al. 2021). Kinetic analysis In addition to the knowledge of the electrostatic properties of green-type Pc, Cc 6 and psaF, a large set of fast time-resolved kinetic data on PSI reduction is available. All these data indicate that the long-range electrostatic attraction between strongly charged complementary patches—acidic in Pc and Cc 6 and basic in PsaF—drives the formation of an initial transient complex in green algae and plants. Further short-range hydrophobic interactions would also be involved in the reorganization of the complex and its fine tuning, to reach an effective configuration leading to an efficient ET (Hippler et al. 1996,Hippler et al. 1998, Sommer et al. 2006). In fact, the structural properties of PSI and the electron donors in the green lineage ensure both efficient partner association (K A up to 10 5 M 1 ) and ET rates (k ET up to 10 5 s 1 )(Bottin and Mathis 1985,Haehnel et al. 1994, Hervás et al. 1995,Drepper et al. 1996,Sigfridsson et al. 1996, Hippler et al. 1998,Molina-Heredia et al. 2003,Klughammer and Schreiber 2016). Again, the exception to this general description is fern Pc: its distinctive electrostatic properties still allow this protein to interact with PSI but with a lower affinity (K A 10 4 M 1 ) and ET rate (k ET 10 3 s 1 )(Navarro et al. 2004). The high efficiency of the green lineage in the formation of the [donor:PSI] complex and in the ET step has, however, a derived disadvantage. Obviously, PSI reduction must be followed by the dissociation of the reactants, but the too strong electrostatic interaction in the green complexes implies a limitation to this dissociation and, therefore, for a new subsequent redox cycle (Finazzi et al. 2005,Kuhlgert et al. 2012). In fact, the release of oxidized Pc (and probably Cc 6 ) from PSI represents the kinetic limiting step in the ET from Cb 6 f to PSI in eukaryotic ‘green’organisms (Drepper et al. 1996,Finazzi et al. 2005). In the red lineage, the kinetic analysis of PSI reduction by Cc 6 in the diatom P. tricornutum indicated the occurrence of a mechanism similar to that previously described in the green lineage systems. Thus, ET takes place according to 1085 Plant Cell Physiol. –62(7): 1082 1093 (2021) doi:https://doi.org/10.1093/pcp/pcab044 Downloaded from https://academic.oup.com/pcp/article/62/7/1082/6194660 by Universidad de Sevilla user on 10 December 2021 an initial protein–protein encounter complex, followed by a reorganization process leading to an optimized final configuration (Bernal-Bayard et al. 2013). However, the kinetic data make evident the weaker electrostatic nature of the interaction, in agreement with the also weaker surface charges of both Cc 6 and PsaF. Consequently, and in comparison with the green-type systems, the red-type system of P. tricornutum showed lower values both for the binding affinity of Cc 6 to PSI (K A 710 3 M 1 ) and ET (k ET 210 4 s 1 ), in addition to a lower efficiency in the formation of the properly arranged [Cc 6 :PSI] complex (Bernal-Bayard et al. 2013), all of which intriguingly reminds the case of the fern [Pc:PSI] system. In this respect, ET kinetics determined by Dual-KLAS/NIR in the cyanobacterium Synechocystis sp. PCC 6803 (with a Pc also having decreased electrostatic properties) showed a much tighter coupling of the redox state of Pc to PSI that can be explained as a higher complex turnover (Theune et al. 2021). Interestingly, an in vitro kinetic analysis showed that the interaction of P. tricornutum PSI with green algae Cc 6 is more efficient kinetically than the native diatom system, with three times and two times higher K A and k ET values, respectively (Bernal-Bayard et al. 2013). This suggests that the native diatom Cc 6 /PSI couple (and maybe the fern complex) represents a compromise between a lower efficiency in ET and a facilitated turnover, due to a lower binding affinity. Finally, there are no functional data available for the interaction of red PSI with the native ‘green-type’ acquired Pc (as in T. oceanica). However, kinetic data indicated that P. tricornutum PSI is able to react with the acidic Pc of green algae and plants, although it does much less efficiently than with P. tricornutum Cc 6 (k ET 310 2 s 1 vs. 2.2 10 4 s 1 , respectively), demonstrating the formation of less-optimal functional complexes (Bernal-Bayard et al. 2013,Bernal-Bayard et al. 2015). However, this efficiency was even doubled in Pc mutants of Chlamydomonas reinhardtii (green alga) in which negative charges were replaced for positive ones, trying to mimic the electrostatics of diatom Cc 6 (Bernal-Bayard et al. 2015). Therefore, while an alternative Pc in diatoms could represent a favorable adaptation to iron limitation, its ‘green-type’ character may limit its ET efficiency (Bernal-Bayard et al. 2015, Antoshvili et al. 2019). Even so, a recent work (Castell et al. 2021) has shown that the heterologous expression in P. tricornutum of the C. reinhardtii Pc mutant more effective in reducing diatom PSI (Pc E85K; Bernal-Bayard et al. 2015) seems to promote an increased growth under iron-limiting conditions, and so to act as a functional alternative to native Cc 6 (Castell et al. 2021). 3D structural data on the complexes Until very recently there were not solved 3D structures at atomic resolution for intermolecular complexes of PSI with the electron donors Pc or Cc 6 . However, from kinetic data, mutagenesis studies, mass spectrometry analysis and computational modeling, it has been possible to propose structural models for the binding of the two soluble electron donors to PSI in green organisms. These models show that this interaction comprises both a global electrostatic attraction and short-range hydrophobic contacts (Hippler et al. 1996,Sommer et al. 2004,Sommer et al. 2006). The latter involve the hydrophobic patches of Pc/Cc 6 and the luminal i/j loops of the PsaA/B subunits in PSI, including the PsaA W651/PsaB W627 residues near P700 (photosystem I primary donor) as part of the ET pathway (Sommer et al. 2004,Busch and Hippler 2011). Electrostatic interactions mainly involve complementary charges in the Pc/ Cc 6 negative areas and common positive groups in the PsaF subunit but also include negative groups in PsaB that are required to favor the unbinding of the oxidized donor from PSI (Sommer et al. 2006,Busch and Hippler 2011). The main features of the interaction between green PSI and Pc have been recently corroborated in the high-resolution structure of a triple complex of plant PSI with Fd and Pc, solved by cryo-electron microscopy (Caspy et al. 2020). The analysis of the PSI–Pc binding confirms that the association of Pc depends principally on hydrophobic interactions, mainly—but not only—around the two tryptophan residues of PSI near P700, with the two copper-ligand His residues in Pc located adjacent to these tryptophan groups. In addition, stabilization of Pc binding occurs mostly through interactions between the positive patch in PsaF and negatively charged Pc residues (D42, D44, E43 and E45), although there are also other additional—but weaker—interactions of PsaF with other negative groups in Pc (Caspy et al. 2020). In the red lineage, the structure of the [Cc 6 :PSI] transient ET complex of the diatom P. tricornutum has been analyzed by computational docking and compared to the equivalent model complex in the green lineage (Fig. 3)(Bernal-Bayard et al. 2015). The P. tricornutum model nicely explains why the Cc 6 /PSI redox couple shows a lower efficiency in diatoms compared to green systems, regarding both the formation of the ET complex and the ET reaction. Thus, in the most energetically favorable docking orientations of the diatom [Cc 6 :PSI] complex, the two tryptophan residues of PSI (PsaA W652/PsaB W624), equivalent to those forming the hydrophobic binding site and the ET path to P700 in green organisms, were located at a relatively long Fig. 3 Best-energy docking models for efficient ET (PsaA W652/PsaB W624and Cc 6 heme groups at <3.0 Ådistance) between P. tricornutum PSI and P. tricornutum or M. braunii Cc 6 . Membrane lipids and PsaA/B subunits of PSI (in light grey), the PsaF subunit of PSI (in dark grey), P700 (in green) and Cc 6 of P. tricornutum (in blue) or M. braunii (in light brown)aredepicted.PSIW652/W624groups,locatedintheETpathway between P700 and the cytochromes, are also shown. PSI:Cc 6 models correspond to the coordinates described in Bernal-Bayard et al. (2015). 1086 C. Castell et al. jElectron transfer from cytochrome fto photosystem I in the green and red branches Downloaded from https://academic.oup.com/pcp/article/62/7/1082/6194660 by Universidad de Sevilla user on 10 December 2021 distance (10–20 Å) from the Cc 6 heme. On the contrary, the [Cc 6 :PSI] docking models that seemed to have the most efficient orientations for ET (PsaA W652/PsaB W624 within 3.0 Å distance from the Cc 6 heme) did not show electrostatic interactions involving the PsaF positive patch or any other relevant favorable interactions (Fig. 3, and see below) (Bernal-Bayard et al. 2015). Altogether, this suggests that the most efficient ET orientations in the P. tricornutum [Cc 6 :PSI] complex are not energetically optimized (Bernal-Bayard et al. 2015). Interestingly, it has been recently reported in the red algae C. merolae that the reduction in isolated PSI by a cyanobacterial Cc 6 was not affected by the removal of the PsaF subunit (Antoshvili et al. 2019). Therefore, it is possible that in the red lineage the PsaF subunit may have a minor role in the interaction with the Cc 6 donor compared to the green lineage. It should be noted that the low-energy docking orientations between the PSI of P. tricornutum and the Cc 6 of the green alga Monoraphidium braunii, placed the two tryptophans of PSI and the Cc 6 heme at shorter distances, in accordance with the more efficient kinetic data of their cross-reaction (Bernal-Bayard et al. 2013,Bernal-Bayard et al. 2015). The docking model of the P. tricornutum [Cc 6 :PSI] complex (Fig. 3) shows a different orientation with respect to the complex previously described in the C. reinhardtii green system (Hippler et al. 1996,Sommer et al. 2004,Sommer et al. 2006, Busch and Hippler 2011). In the latter, there are specific interactions involving the negative patch of Cc 6 and positive residues of PsaF, as well as a strong salt bridge between the adjacent R66 and D65 groups in Cc 6 and the R623/D624 pair of PsaB in PSI, which results in a distance between redox cofactors of 14 Å (Sommer et al. 2006,Busch and Hippler 2011). However, the D65 group is not conserved in P. tricornutum, and thus, it cannot stabilize this orientation, resulting in a less favorable overall binding energy. Interestingly, the docking model of the complex between M. braunii Cc 6 and P. tricornutum PSI shows strong interactions between the negative patch of Cc 6 and positive residues of PsaF, as described in the equivalent native complex of C. reinhardtii (Fig. 3)(Sommer et al. 2006). Finally, in diatoms a 3D structure of the interaction of the native ‘green-type’acquired Pc and PSI is not available. However, docking simulations performed for the interaction between P. tricornutum PSI and C. reinhardtii Pc (BernalBayard et al. 2015) showed that the best-energy model for an efficient docking displays a strong electrostatic interaction between the PsaF positive patch and Pc negative residues (D42, E43 and D44), which are equivalent to those described in the Pc/PSI complex of plants (Bernal-Bayard et al. 2015,Caspy et al. 2020). In agreement with this model, it was recently reported that the reduction reaction of PSI of the red alga C. merolae by a plant Pc was strongly inhibited when depleting the PsaF subunit (Antoshvili et al. 2019). Interaction of Cb 6 f with the Soluble Carriers Structural and kinetic analyses After their oxidation by PSI, Pc and Cc 6 have to be reduced again by cytochrome f(Cf), a subunit of the Cb 6 f complex and the exit point of electrons to the luminal soluble electron acceptors (Cramer 2019). Cf (a c-type cytochrome with ca. 290 amino acids) have been structurally characterized in cyanobacteria and organisms of the green lineage, both as the entire protein within the Cb 6 f complex and as a soluble truncated form (Martinez et al. 1994,Chi et al. 2000,Kurisu et al. 2003, Stroebel et al. 2003). Unique features of Cf include the N-terminal tyrosine Y1 as an unusual heme-binding residue (Martinez et al. 1994,Chi et al. 2000). Moreover, the overall protein fold is also unusual for a c-type cytochrome, as it displays a large and elongated hydrophilic domain (ca. 250 amino acids) with a b-barrel-based structure, exposed to the thylakoid lumen (Kurisu et al. 2003,Stroebel et al. 2003). This hydrophilic domain is anchored to the membrane by a single transmembrane a-helix (Fig. 1B). In addition, the lumen-exposed region exhibits a large and a small, flat domain, connected by a bend, with the larger domain containing the heme group and the Y1 ligand, which are solvent exposed (Martinez et al. 1994) (Figs. 1B,4). The green-type Cf shows a positively charged surface patch at the interface of the large and small domains, complementary to the negative region of Pc or Cc 6 (Martinez et al. 1994,Soriano et al. 1997)(Fig. 4). As will be discussed later, in organisms from the green lineage both the large and the small domain are involved in the docking of Pc or Cc 6 . Regarding the red lineage, the surface electrostatic potential computed on structural models of P. tricornutum and T. oceanica Cf (Fig. 4) shows the conservation of the ET exposed area around the Y1 and the heme, as well as a positive area placed in an equivalent position relative to the green-type Cf (Fig. 4). However, due to the substitution of positive groups in the interface of the large and small domains, in the red lineage the positive electrostatic potential of this area in the Cf surface is sensibly reduced (Fig. 4). Therefore, similar and parallel changes occur in Cf and psaF, also in accordance with the decreased electrostatic character of Cc 6 . However, in the case of T. oceanica Cf, it is interesting to note the presence of positive groups outside the usual ET region around the heme-exposed area (Fig. 4). Kinetic data for native proteins are only available for the Pc/ Cf interaction (reviewed in Bendall and Howe 2016), since the overlap between the absorbance spectra of Cc 6 and Cf prevents to accurately follow the kinetics of the ET reaction between both native cytochromes (Grove and Kostic 2003). In eukaryotes, the Pc/Cf interaction has been studied in green algae and plants (Qin and Kostic 1993,Meyer et al. 1993,Kannt et al. 1996, Soriano et al. 1997,Gong et al. 2000,Schlarb-Ridley et al. 2003). The final conclusion, obtained not only from kinetic data but also from computational and structural NMR studies (see below), is that ET from Cf to Pc occurs following a similar mechanism (binding, reorganization and ET) to that described for PSI reduction (Cruz-Gallardo et al. 2012,Schilder and Ubbink 2013). Again, this mechanism involves long-range electrostatic attraction for initial binding and preorientation of the two partners, as well as short-range hydrophobic contacts to achieve a productive configuration for ET. Given that previous structural and functional data have established the general functional equivalence of Pc and Cc 6 (De la Rosa et al. 2006, 1087 Plant Cell Physiol. –62(7): 1082 1093 (2021) doi:https://doi.org/10.1093/pcp/pcab044 Downloaded from https://academic.oup.com/pcp/article/62/7/1082/6194660 by Universidad de Sevilla user on 10 December 2021 Hippler and Drepper 2006), the hypothesis that a similar mechanism controls the ET reaction of Cf with Cc 6 in green-type organisms seems reasonable. Structure of the complexes in the green lineage Several NMR-based structures of plant [Cf:Pc] complexes have been described (Ubbink et al. 1998,Lange et al. 2005)(Fig. 5A). In addition, Brownian Dynamics (BD) simulations have studied the Cf interaction with the soluble electron acceptors in plants and in the green alga C. reinhardtii (Gross and Pearson 2003, Musiani et al. 2005,Haddadian and Gross 2006,Fedorov et al. 2019). However, these types of studies have been absent until now in organisms of the red lineage. The NMR-based structure of the plant [Cf:Pc] complex (Fig. 5A) shows a relatively large complex interface, involving both the large and small domains of Cf, forming the so-called ‘side-on’relative orientation (in contraposition to the ‘head-on’ orientation that appears in some complexes of cyanobacteria) (Ubbink et al. 1998,Cruz-Gallardo et al. 2012). As expected, the complex shows electrostatic interactions between positively charged residues in Cf and negatively charged groups in Pc, as well as interactions between both hydrophobic patches, thus providing an efficient ET pathway (Fe–Cu distance of 10.9 Å) (Ubbink et al. 1998)(Fig. 5A). Although a detailed analysis on NMR-based docking models of the [Cf:Cc 6 ] complex has been carried out in the cyanobacterium Nostoc (Díaz-Moreno et al. 2014), in eukaryotic organisms, the available data are limited to BD simulations applied to the green [Cf:Cc 6 ] (and [Cf:Pc]) complex of C. reinhardtii (Gross and Pearson 2003,Haddadian and Gross 2005,Haddadian and Gross 2006; and see below). These studies demonstrated the functional equivalence of the hydrophobic and electrostatic areas in both green Pc and Cc 6 . As described previously for the Cf:Pc interaction, hydrophobic residues around the heme of Cf (including Y1) interact with hydrophobic residues surrounding the Cc 6 heme. On the other hand, electrostatic attractive interactions involve positive groups in Cf—located in both the large and small domains—and negative groups in Cc 6 . This results in a ‘side-on’orientation, which also shows a relatively large interface (Haddadian and Gross 2005; and see Fig. 5B). In conclusion, the comparison of the green [Cf:Pc] and [Cf:Cc 6 ] complexes confirms the key role of the electrostatic and hydrophobic areas in protein binding and, therefore, in achieving an efficient ET, but also the role of the small Cf domain in stabilizing the proper docking orientations of Pc and Cc 6 . Structural modeling of the complexes in the red lineage As mentioned above, there are no structural data available about the [Cf:Cc 6 ] (and [Cf:Pc]) complexes in organisms of the red lineage. However, ET proteins, and in particular, ctype cytochromes, have been shown to be structurally quite rigid, and thus, BD and docking approaches treating Cf and Cc 6 as rigid bodies have been widely used (Haddadian and Gross 2005,Haddadian and Gross 2006,Cruz-Gallardo et al. 2012, Bernal-Bayard et al. 2015). Here, we have carried out docking simulations to study the [Cf:Cc 6 ] complex in the red lineage by using pyDock, a protein–protein rigid-body docking protocol that uses electrostatics and desolvation energy to score the generated docking poses (Jiménez-García et al. 2013; and see Supplementary Appendix S1). To validate our analysis, the same approach has also been used to explore the [Cf:Cc 6 ] complex of the green alga C. reinhardtii, previously studied by BD simulations (Gross and Pearson 2003,Haddadian and Gross 2005,Haddadian and Gross 2006). Our results indicate that Fig. 4 (Left) Backbone representation of the exposed luminal domain of Cf of the green alga C. reinhardtii. The view displays in front both the heme group and the bound Tyr1 (in green). (Right) Surface electrostatic potential distributions of C. reinhardtii Cf and the structural models of P. tricornutum and T. oceanica Cf. The proteins are oriented as shown in the backbone representation on the left, and electrostatic potential values are shown on the same scale as in Fig. 2. See the text and Supplementary Appendix S1 for more details. 1088 C. Castell et al. jElectron transfer from cytochrome fto photosystem I in the green and red branches Downloaded from https://academic.oup.com/pcp/article/62/7/1082/6194660 by Universidad de Sevilla user on 10 December 2021 the best-energy docking models for the C. reinhardtii complex converge toward an orientation with the shortest cofactors distance (Fig. 5B). This orientation is basically identical to the models previously described (Haddadian and Gross 2005, Haddadian and Gross 2006). In the case of P. tricornutum, the energy-distance landscape of the docking models of the native [Cf:Cc 6 ] complex also shows that the lowest-energy models converge toward orientations with a short distance between the heme iron in Cf and the heme group in Cc 6 (Supplementary Fig. S1). As it was the case in the [Cc 6 :PSI] complex (Bernal-Bayard et al. 2015; and see above), the P. tricornutum [Cf:Cc 6 ] docking complex has a different orientation compared with the previously described model of the equivalent green complex in C. reinhardtii (Haddadian and Gross 2005)(Fig. 5B, C). Thus, the docking models of the P. tricornutum [Cf:Cc 6 ] complex with the lowest—most favorable—energies and optimal distances between the Fe in Cf and the heme cofactor in Cc 6 (energies from 25 to 22 a.u.; distances between 8 and 10 Å) have a ‘head-on’orientation (Fig. 5C), more similar to that described previously in cyanobacterial systems. These orientations display smaller interfaces and weaker electrostatic interactions, with hydrophobic interactions being more relevant (Crowley et al. 2001, Cruz-Gallardo et al. 2012)(Fig. 5C). In the ‘side-on’configuration observed in the complex of C. reinhardtii, negatively charged residues of Cc 6 (D2, E54, K57, D65, R66, E70 and E71) have significant contacts with positive residues in Cf, in particular with the K65, K188 and K189 groups. However, in P. tricornutum, the E54 and D65 groups in Cc 6 and K189 in Cf are not conserved and cannot stabilize this ‘side-on’orientation. On the contrary, the orientation in P. tricornutum involves relevant electrostatic interactions between E28 in Cc 6 and residues K65 and K187 in Cf, as well as between residues R66 and K57 in Cc 6 and residues E108 and E164 in Cf. As a consequence, the small domain of P. tricornutum Cf appears to have only a minor role in the interaction with Cc 6 (Fig. 5C). In fact, the simulated deletion of the Cf small domain did not significantly affect the docking to Cc 6 (Fig. 5D). This latter result is in contrast with the docking simulation in C. reinhardtii, in which the deletion of the Cf small domain affected the binding positions on Cf of both Pc and Cc 6 (Haddadian and Gross 2006). Interestingly, ‘side-on’ docking orientations can be found in P. tricornutum that are relatively similar to the best-energy models obtained in C. reinhardtii. However, although these orientations have short distances between the Fe in Cf and the heme cofactor in Cc 6 (between 9.0 and 10.5 Å), they have unfavorable higher energies (from 15 to 9 a.u.) (Supplementary Fig. S1). Finally, it is interesting to analyze the modeled [Cf:Pc] complex of T. oceanica in comparison with the above-described [Cf: Cc 6 ] complex of P. tricornutum, as well as with the equivalent [Cf:Pc] complex of the green lineage. Remarkably, the docking between T. oceanica Cf and Pc produced a much larger population of low-energy docking orientations than the P. tricornutum [Cf:Cc 6 ] complex. This indicates higher binding affinities, possibly as a result of the stronger electrostatics of the acquired ‘green-type’Pc (Fig. 2B,Supplementary Fig. S1). However, in T. oceanica, the energy-distance landscape does not converge toward a single prevalent low-energy orientation with a short distance between Cu (in Pc) and Fe (in the Cf heme). On the contrary, it results in a series of different orientations in a similar range of energies and distances (Fig. 6A–C,Supplementary Fig. S1). This suggests that the Cf:Pc interaction in T. oceanica is not optimized in a single group of solutions with a favorable combination of interactions and distances for ET. Therefore, different protein–protein orientations might be functionally possible and could coexist. The best-energy/distance balanced dockings of the T. oceanica Cf:Pc modeled complex include: (i) a ‘headon’configuration, similar to some cyanobacterial complexes, involving the hydrophobic patches of both Cf and Pc with no interactions between the electrostatic areas or the Cf small domain (energy of 24.2 a.u. and the shortest Fe–Cu distance Fig. 5 (A) Representative structure for the plant Cf:Pc complex (turnip CfandspinachPc;obtainedfromPDBcode2pcf)(Ubbinketal.1998).Pc is colored in blue and the copper-bound His87 is shown. Best-energy docking models for efficient ET between Cf and Cc 6 (in red) of (B) the green alga C. reinhardtii (rank 1, docking energy 32.0 a.u., distance between Fe in Cf to heme in Cc 6 of 8.2 Å) and (C) the diatom P. tricornutum (rank 6, docking energy 24.7 a.u., distance between Fe in Cf to heme in Cc 6 of 8.0 Å, the shortest distance model). (D) Superimposed docking models of the P. tricornutum [Cf:Cc 6 ] complex (in C) and the model (in blue) corresponding to a truncated Cf without the small domain (rank 2, docking energy 31.0 a.u., distance between Fe in Cf to heme in Cc 6 of 8.4 Å). See Supplementary Appendix S1 for more details. 1089 Plant Cell Physiol. –62(7): 1082 1093 (2021) doi:https://doi.org/10.1093/pcp/pcab044 Downloaded from https://academic.oup.com/pcp/article/62/7/1082/6194660 by Universidad de Sevilla user on 10 December 2021 of 11.1 Å) (Fig. 6A); (ii) a ‘side-on’configuration, more similar to the green complexes (Figs. 5A,6B), involving the electrostatic and hydrophobic patches of both proteins and the small domain of Cf (energy of 31.1 a.u.; Fe–Cu distance of 12.7 Å) (Fig. 6B); and (iii) an intermediate configuration, which includes the hydrophobic patches and some residues of the electrostatic patches (energy of 30.1 a.u.; Fe–Cu distance of 12.2 Å) (Fig. 6C). The fact that Pc is not present in many diatom species suggests that its heterologous acquisition may be relatively recent, and therefore, the [Cf:Pc] complex in T. oceanica would not have yet achieved a prevalent configuration optimized for ET. Surprisingly, the clearly more favorable energies (energy of 38 to 36 a.u.) correspond to a cluster of very similar orientations, in which the electrostatic patches of both proteins establish strong interactions, also involving the additional positive groups outside the usual region of ET in Cf. However, in these orientations, the Cu and the Fe are located at a long distance from each other (Fe–Cu distance of 20 Å) (Fig. 6D, Supplementary Fig. S1). Obviously, these orientations also show a long distance from the Fe in Cf to the groups that form the hydrophobic patch of Pc, and in particular to the Cu-binding H87 residue (H88 in the secuence of T. oceanica Pc) (Fig. 6D). The imidazole of H87 defines the shortest path from the copper atom to the protein surface (6 Å) and forms the ET pathway in Pc from Cf and to PSI in cyanobacteria, green algae and plants (Fig. 5A), as clearly demonstrated by kinetic and structural data (Ubbink et al. 1998,Hippler and Drepper 2006,Caspy et al. 2020). However, in the apparently unproductive orientations of lowest energy in T. oceanica, the highly conserved Y84 residue in Pc (typically referred to as Y83 in cyanobacterial and eukaryotic Pc) points directly toward the heme-binding Y1 of Cf (Y1–Y84 distance of 5.1 Å; Fe–Y84 distance of 9.9 Å) (Fig. 6D). Based on the reactivity of Pc with inorganic complexes, as well as on theoretical studies, it was previously proposed that Y83 could act as another potential ET pathway in Pc (reviewed in Redinbo et al. 1994). However, this hypothesis was later discarded in the interaction of Pc with Cf and PSI (Ubbink et al. 1998,Hippler and Drepper 2006,Caspy et al. 2020). Therefore, it could be a matter of discussion whether this alternative ET pathway involving Y83—not active in green Pcs—is active or not in the heterologously acquired diatom Pc. Final Remarks One of the most fascinating facts in the evolution of the ET from Cb 6 f to PSI is that while the overall backbone structures of the involved proteins are well conserved from prokaryotic to eukaryotic organisms, there are however significant differences, mainly affecting the electrostatic surface potential, which imply the appearance of functional areas with their own distinctive characteristics in the different types of photosynthetic organisms. The green and red lineages comprise the two main lines in the evolution of eukaryotic photosynthetic organisms from a common ancestor. This ancestral organism resulted from an endosymbiotic event, in which a cyanobacterium ancestor gave origin to the eukaryotic chloroplast. Thus, in both the greenand red-type chloroplasts, the exit point is cyanobacterial photosynthesis, in which a shorter PsaF protein does not playacrucialroleinthebindingofthetwosolublealternative electron donors, Pc and Cc 6 . Moreover, in cyanobacteria, the two soluble carriers do not have electrostatic characteristics widely shared among these organisms, predominating proteins of either almost neutral or strongly positive character, that have a complementary counterpart on the surface of the Cf partners. In eukaryotic organisms, the evolution of the ET from Cb 6 fto PSI has developed conserved negative electrostatic patches in the soluble carriers, parallel to positive complementary areas of interaction in the two membrane complexes (in Cf and PsaFPSI). However, within this general framework, the green line has developed areas of interaction with a stronger electrostatic charge, with a higher affinity between partners and a more Fig.6 (A–C)Representativebest-energydocking modelsforefficientET between the modeled Cf and Pc (in blue) of the diatom T. oceanica (selected by the shorter distances between Fein Cf and Cu in Pc). The Cf heme, the iron-bound Tyr1 and the copper-bound His88 in Pc are shown. (A) Rank 62, docking energy 24.2 a.u., distance between Fe in Cf to Cu in Pc of 11.1Å (the shortest distance model). (B) Rank 8, docking energy 31.1 a.u., distance between Fe in Cf to Cu in Pc of 12.7 Å.(C) Rank 16,docking energy30.1a.u.,distancebetweenFeinCf to Cu in Pc of 12.2 Å. (D) Best-energy docking model showing the Tyr84 groupinPc(rank1,docking energy38.4 a.u.,distance betweenFe inCf to Cu in Pc of 20.5 Å). See the text andSupplementary Appendix S1 for more details. 1090 C. Castell et al. jElectron transfer from cytochrome fto photosystem I in the green and red branches Downloaded from https://academic.oup.com/pcp/article/62/7/1082/6194660 by Universidad de Sevilla user on 10 December 2021