Structure of Dunaliella photosystem II reveals conformational flexibility of stacked and unstacked supercomplexes
Abstract
Photosystem II (PSII) generates an oxidant whose redox potential is high enough to enable water oxidation , a substrate so abundant that it assures a practically unlimited electron source for life on earth . Our knowledge on the mechanism of water photooxidation was greatly advanced by high-resolution structures of prokaryotic PSII . Here, we show high-resolution cryogenic electron microscopy (cryo-EM) structures of eukaryotic PSII from the green alga Dunaliella salina at two distinct conformations. The conformers are also present in stacked PSII, exhibiting flexibility that may be relevant to the grana formation in chloroplasts of the green lineage. CP29, one of PSII associated light-harvesting antennae, plays a major role in distinguishing the two conformations of the supercomplex. We also show that the stacked PSII dimer, a form suggested to support the organisation of thylakoid membranes , can appear in many different orientations providing a flexible stacking mechanism for the arrangement of grana stacks in thylakoids. Our findings provide a structural basis for the heterogenous nature of the eukaryotic PSII on multiple levels.
Editor's evaluation
The study provides exceptional new insights into the structural organization of the light energy capturing machinery of photosynthesis by resolving the structure of the water-splitting photosystem II and associated complexes isolated from the green alga Dunaliella using cryo-electron microscopy and modeling. The results indicate large-scale flexibility of photosystem II – light harvesting complex II supercomplex, which are likely to have functional consequences for the stacking of thylakoid membranes, interactions of complexes into mesoscale structures, which may control the funneling light energy into the photosystems or photoprotective mechanisms, controlling light-driven fluxes of electrons and protons.
https://doi.org/10.7554/eLife.81150.sa0Introduction
In eukaryotes, the light reaction of oxygenic photosynthesis occurs in chloroplasts. Four protein complexes essential for the light reactions reside in an elaborate membrane system of flattened sacs called thylakoids (Nelson and Ben-Shem, 2004). From these four complexes, the photosystem II (PSII) complex catalyses light-driven water oxidation and provides the electrons used for carbon fixation (Vinyard et al., 2013; Nelson and Junge, 2015).
Thylakoids form a physically continuous three-dimensional network, differentiated into two distinct physical domains: cylindrical stacked structures (called grana) and connecting single membrane regions (stroma lamellae). Photosystem I (PSI) is mainly located in the stroma lamellae, while PSII is found almost exclusively in the grana (Hankamer et al., 1997; Anderson, 2002; Kim et al., 2005). Grana stacking is a dynamic process dependent on the internal osmotic pressure, the luminal ion composition, and environmental cues and is thought to be supported by interactions among PSII complexes (Rubin et al., 1981; Barber et al., 1980; Kirchhoff, 2014; Liu et al., 2018; Dalal and Tripathy, 2018, Garab and Mustárdy, 2000; Kirchhoff et al., 2007).
PSII is a homodimer with a molecular mass of ~500 kDa, each monomer contains cofactors such as chlorophylls (Chls), quinones, carotenoids, and lipids which are coordinated by at least 20 protein subunits (Shen, 2015; Cox et al., 2020; Barber, 2004). In each PSII core, a cluster of four manganese (Mn) and one calcium (Ca) carries out H2O oxidation and O2 release (Nelson and Yocum, 2006; McEvoy and Brudvig, 2006). The eukaryotic reaction centre is a dimer surrounded by tightly bound monomeric light-harvesting complexes (LHCs) and trimeric LHCII complexes (Daskalakis, 2018; Mascoli et al., 2020; Barber et al., 2002). Two monomeric LHCs, CP26, and CP29 are located between LHCII trimers and PSII core subunits (Nield et al., 2000) additional LHCII trimers can bind PSII depending on light intensity and quality (Goldschmidt-Clermont and Bassi, 2015).
Although more than 3 billion years of evolution separate cyanobacteria, red algae, green algae and plants, and high-resolution PSII structures show that each PSII monomer along with its dimeric arrangement is highly conserved, especially in the membrane-bound regions of the PSII (Wei et al., 2016; Su et al., 2017; Pi et al., 2019; van Bezouwen et al., 2017; Sheng et al., 2019; Ago et al., 2016). Structural and spectroscopic investigations uncovered various aspects of PSII’s water splitting mechanism, but a complete model is still missing (Umena et al., 2011; Cox et al., 2014; Suga et al., 2019; Kupitz et al., 2014; Kern et al., 2018). Most of the mechanistic and structural studies of PSII were performed in thermophilic cyanobacteria, but structural studies of PSII from the eukaryotic lineage are lagging in terms of resolution and water molecules network (Ago et al., 2016; Sheng et al., 2019; Shen et al., 2019; Wei et al., 2016; Su et al., 2017; Pi et al., 2019; van Bezouwen et al., 2017; Nagao et al., 2019).
In this work, PSII was isolated from the halotolerant green alga Dunaliella salina. A high-resolution (2.43 Å) structure of PSII shows structural properties of the Dunaliella PSII supercomplex. The eukaryotic PSII appears to exist in two distinct core conformations that differ substantially in their inner dimer separation and the location of CP29, an important monomeric LHC. Structural analysis of stacked PSII dimers showed highly flexible interactions which can play a role in the dynamic organisation of chloroplast membranes. These findings introduce an additional, underlying, level of organisation, which can impact its excitation energy transfer (EET) properties and the overall organisation of the thylakoid membranes.
Results and discussion
Two distinct PSII conformations in green alga
Highly active PSII from D. salina cells was applied on glow-discharged holey carbon grids that were vitrified for cryo-EM structural determination (see Methods). Initial classification of the dataset showed that approximately 20% of the particle population were in a stacked PSII configuration, containing two PSII dimers facing each other on their stromal side (Figure 1—figure supplement 1). From the unstacked PSII dimers, approximately 20% were in the C2S configuration (two Cores, one Stable LHCII), which was previously identified by low-resolution cryo-EM (Drop et al., 2014; Figure 1—figure supplement 1c). The majority of the PSII particles contained two LHCII in the C2S2 configuration. The map of the C2S2 particles refined to a global resolution of 2.82 Å (Figure 1—figure supplement 2). Close examination of this map revealed that CP29, one of the monomeric LHC proteins, implicated as a junction for EET from LHCII trimers to PSII core, appeared to be in lower resolution than the rest of the supercomplex (Figure 1—figure supplement 2a). Indeed, when this particle set was further classified, two distinct PSII conformations of the PSII supercomplex became apparent. In these two conformations, the two PSII cores are shifted laterally with respect to each other (Figure 1). This lateral shift is accompanied by several other associated movements, most noticeably, a large movement of the CP29 subunit (Figure 1) in line with our initial observation. The two conformers were denoted compact and stretched PSII (C2S2COMP and C2S2STR, respectively), and the high-specific activity of 816 μmol O2/(mg Chl * hr) measured for the preparation prior to vitrification suggest that both are highly active. The final reconstruction of the compact orientation refined to an overall resolution of 2.43 Å, the highest of any eukaryotic PSII structures (PDB ID 7PI0; Figure 1a; Supplementary file 1 and Figure 1—figure supplements 1–3).
The C2S2COMP structure is similar to the previously determined C2S2 supercomplex from Chlamydomonas reinhardtii or higher plants (Shen et al., 2019; Sheng et al., 2019) and the cyanobacterial core structures (Umena et al., 2011; Kato et al., 2021). The second, stretched conformer was solved to 2.62 Å resolution and accounted for about 37% of the unstacked C2S2 PSII particles. Figure 1 and Videos 1–2 depict the superposition of the polypeptide chains of the two conformers, showing major differences in the location and orientation of PSII monomers. Superposition of the Dunaliella and Chlamydomonas C2S2 (PDB 6KAC) structures and maps suggests that the Chlamydomonas structure also contains these different conformers. This may explain the decreased local map resolution presented in the aforementioned subunits, compared to the rest of the cryo-EM map (Sheng et al., 2019).
Structures of D. salina unstacked PSII at high resolution
Thus far, available PSII structures suggested a single, highly conserved organisation of the two PSII cores (Su et al., 2017; Wei et al., 2016; Pi et al., 2019; Sheng et al., 2019; Umena et al., 2011; Ago et al., 2016). The high-resolution structures of Dunaliella C2S2COMP and C2S2STR provide a new perspective on the dynamic arrangement of eukaryotic PSII and the interaction of the core complex with its LHCs. To compare the C2S2COMP and C2S2STR, the core complexes were aligned (Figure 1 and Videos 1–2). Initial inspection showed that one of the major differences between the two conformations is the orientation of CP29 (Figure 1e). In C2S2STR CP29, helices A and C move towards the LHCII trimer of the opposite monomer and away from CP47, with helix B of CP29 serving as a rotation axis. Moreover, CP29COMP contained only 9 Chls, compared to 11 in CP29STR and 13 Chl in Chlamydomonas PSII CP29 (Sheng et al., 2019 and Supplementary file 2). Chls 605 and 616 were absent in both structures, and CP29COMP was also missing Chls 611 and 613. This might be attributed to the flexibility of CP29 C-terminus (which is proximal to 613), or some side chains and ligands rearrangement associated with the movement of CP29. The structure and b-factor of CP29COMP and CP29STR are similar (Figure 1—figure supplement 4); however, superposition of C2S2COMP, C2S2STR, and Chlamydomonas C2S2M2L2 (Sheng et al., 2019) suggests that CP29COMP may be an intermediate conformation between CP29STR, that is bound to the S-trimer through its C-terminus, and Chlamydomonas CP29, that binds the M- and L-trimers via the C-terminus (Figure 1—figure supplement 5).
In the C2S2STR conformation, the PSII monomers slid in the membrane plane along the central symmetry axis separating them (Figure 1b and Videos 1–2). The non-aligned core shows the extent of the shift in the core peptides together with the minor LHCs and LHCII trimer (Figure 1c and d). As a result, all the interactions at the core’s interface are modified, leading to local changes in chain orientations and the conformations of some loops. Core subunits at the centre of the monomer displayed a greater shift (D1, D2, CP47 CP43, and PsbO were displaced by 6–10 Å; Figure 1), and the peripheral subunits showed the largest shift and tilt compared to C2S2COMP (PsbE, PsbP, and CP26 moved by 13 Å, PsbZ showed the largest relocation of nearly 15 Å, and CP26 showed a maximal tilt of 16°; Figure 1). Multibody refinement (Nakane et al., 2018) of both C2S2COMP and C2S2STR demonstrated that the two PSII monomers in each conformation contain additional structural heterogeneity (Figure 1—figure supplements 6–7 and Videos 3–6).
Distinct CP29 conformations alter LHCII to PSII core connectivity
The observed conformational change of CP29 alters EET pathways from LHCII to the PSII core and may account for the differences between calculated and measured EET (Chmeliov et al., 2014; Chmeliov et al., 2016; Caffarri et al., 2011; Mascoli et al., 2020; Croce and van Amerongen, 2011; Croce and van Amerongen, 2020; van der Weij-de Wit et al., 2011). To assess changes in transfer rates between the stretched and compact orientations, we measured how the distances between the closest Chls of CP47 (PSII core), CP29, and LHCII change between the two PSII conformations (all reported Chl distances are measured from the central Mg atom). Overall, we find that CP29 and LHCII move closer to each other and away from CP47 in the stretched configuration. The average distances between CP29 Chls 603, 607, and 609 to the CP47 Chls 607 and 616, increased from 17 Å to 20 Å between the compact to stretched conformations, suggesting faster transfer rates from CP29 to CP47 in the compact conformation. In contrast to this, the average distances between CP29 Chls 604 and 612 to LHCII Chls 604 and 608 increased from 20 Å in the stretched conformation to 23 Å in the compact conformation, suggesting that transfer from LHCII to CP29 is slower in the compact orientation (Figure 1e). The missing CP29 Chls 611 and 613 form part of the interface to LHCII and are missing in the compact conformation, which should also contribute to slower transfer rates from LHCII to CP29 in the compact configuration (Figure 1—figure supplement 8). Altogether, transfer from LHCII to the PSII core should be considerably slower in the compact orientation from both distance and Chl occupancy considerations. Similar features of altered Chl conformations were identified in molecular dynamics (MDs) simulation of LHCII exploring its structural dynamics (Liguori et al., 2015) compared to its crystal structure (Liu et al., 2004). The analysis showed differences in the excitonic coupling of Chl clusters 606–607 and 611–612. MD suggested an increase in the interaction energies of 606–607 and a decrease in the interaction energies of similar proportion in 611–612 (Liguori et al., 2015). The 611–612 Chl pair was proposed as a light-harvesting regulator of EET from CP29 to CP47 and as a quenching site (Caffarri et al., 2011; Ruban et al., 2007; Novoderezhkin et al., 2005; Pascal et al., 2005), as its change in fluorescence yield was attributed to a protein conformational change that leads to a redistribution of the interpigment energetics (Valkunas et al., 2012).
The compact and stretched PSII conformations contain substantial levels of continuous structural heterogeneity
Using multibody refinement (Nakane et al., 2018), with each PSII monomer defined as a separate rigid body, significantly improved the resolution and map quality in both C2S2COMP and C2S2STR, showing that substantial structural heterogeneity exists in both datasets at the level of PSII monomers. Analysing the shape of the heterogeneity in C2S2COMP and C2S2STR, using principal component analysis (PCA) showed that the first six principal components (PCs) explain more than 85% of the variance in the data and consist of continuous heterogeneity (Figure 1—figure supplements 6–7; Videos 3–6). Substantial displacements of approximately 13 Å are observed between the two monomers in the compact conformation (Figure 1—figure supplement 6), and a larger range of displacements (up to 20 Å) exists in the stretched conformation (Figure 1—figure supplement 7). The direction of PCs describes translations perpendicular and parallel to the membrane plane. This suggests that both conformations are flexible and can respond to different membrane curvature (Videos 3–6). To examine the possible effects on energy transfer, we measured the change in intermonomer Chl distances across the different components. As expected, the PCs describing changes in the membrane plane markedly change some key distances between LHCII, CP29, and D1 across monomers (Figure 2). This means that within each PSII conformation, substantial levels of heterogeneity in transfer rates should be considered. Changes in Chl positions were observed in CP29 Chls linking CP29 to PSII core and those connecting CP29 with LHCII. These Chls moved by an average distance of more than 5 Å, in both conformations (Supplementary file 3). This implies that the association between PSII monomers and between PSII cores and LHCs contains a certain degree of freedom which can modulate EET; the entire assembly may be affected by changes in thylakoid membrane properties such as fluidity, composition, and curvature (Tardy and Havaux, 1997; Johnson and Wientjes, 2020).
Water channels and post-translational modifications in Dunaliella PSII
More than 1700 water molecules were detected in the C2S2COMP model (Figure 3a–b), the first detailed water molecules structure for a eukaryotic PSII. Overall, water molecules are clearly excluded from the membrane space in the PSII core, in contrast, the region occupied by LHC’s shows a relatively high number of water molecules in the membrane region. This stems from the presence of several conserved charged amino acids in these antennae and is probably important for the inclusion of such hydrophilic residues within the membrane. We used CAVER (Chovancova et al., 2012) to analyse the structure of internal cavities around the oxygen-evolving complex (OEC). As expected from the highly conserved environment around the OEC, the water channels identified previously in the high-resolution cyanobacterial core structure (Suga et al., 2019; Kaur et al., 2019) are clearly visible in the eukaryotic PSII, and overlap with the results of the internal cavity analysis, these are shown in Figure 3c and named ‘Large,’ ‘Narrow,’ and ‘Broad,’ following Kaur et al., 2019. When analysing the side chains lining the cavities around the OEC, a small hydrophobic patch, highly conserved in prokaryotes and eukaryotes (Figure 3—figure supplement 1), was identified at the beginning of the large channel (Figure 3d). This hydrophobic element may facilitate O2 release as part of the catalytic cycle (Figure 3d).
Several unique map densities were identified during model building, close to the OEC of both configurations a Na+ ion was modelled. This Na+ ion is coordinated by D1-His337, the backbone carbonyls of D1-Glu333, D1-Arg334, D2-Asn350, and a water molecule, in agreement with the recently identified (Wang et al., 2020) binding site (Figure 3—figure supplement 2a–b). This agrees with several studies showing that Na+ ions are required for optimal activity of PSII (Wang et al., 2020; Pogoryelov et al., 2003). Two additional densities, unique to C2S2COMP, were observed close to CP29-Ser84 and CP47-Cys218 in the stromal interface between CP29, CP47, and PsbH and within 10 Å of each other. These were modelled as post-translational modifications (PTMs) – Ser84 appears to be phosphorylated and Cys218 seems to be sulfinylated (Figure 3—figure supplement 2c–d). Thus far, PTMs were structurally seen in photosystems only as phosphorylated LHCII bound to PSI during state transition (Pan et al., 2018; Huang et al., 2021; Pan et al., 2021). Although they were not identified in-situ, several phosphorylation sites were shown to exist in CP29 large stromal loop (Chen et al., 2013; Liu et al., 2009; Poudyal et al., 2020; Hansson and Vener, 2003). CP29 phosphorylation was suggested to be linked with various stress responses, photosynthetic protein degradation, and state transition. Cysteine sulfinylation was shown to be linked to superoxide radical (O2.−) accumulation, which is subsequently converted by superoxide dismutase to hydrogen peroxide (H2O2) molecules (Sevilla et al., 2015; Rey et al., 2007; Matamoros and Becana, 2021). CP47-Cys218 is positioned on the outer edge of PSII, close to the stromal end of the thylakoid membrane, and thus is susceptible to oxidation by H2O2. The map density around Cys218 suggests two cysteine oxidation events which result in the formation of sulfinic acid (RS-O2H).
To summarize, the high-resolution structure of the eukaryotic PSII revealed two distinct states of the PSII complex, adding a new dimension to the known, large compositional heterogeneity of this important system (Croce and van Amerongen, 2011; Caffarri et al., 2009; Kouřil et al., 2020). The increased map resolution resulted in the identification of PTM’s, and several conserved hydrophobic residues near the OEC, which may serve as a pathway for the release of O2. In addition to the two distinct conformations, large levels of continuous structural heterogeneity were discovered within each individual state. Multibody analysis (Nakane et al., 2018) inherently treats the data as a collection of rigid bodies. This is a good approximation of the heterogeneity in photosynthetic systems but should be regarded as a conservative estimation to additional modes of heterogeneity which exist in this system within each body (Liguori et al., 2015).
The structure of D. salina stacked PSII at high resolution
The thylakoid membrane is made of two spatially distinct regions, stroma lamellae and grana stacks, each serving a different role in the photosynthetic process (Pribil et al., 2014; Koochak et al., 2019). Grana stacks size and numbers are affected by light intensity and ionic composition and can change rapidly (Wood et al., 2019). Membrane stacking depends on the presence of cations, mainly Mg2+, which is abundant in the thylakoid stroma (Ishijima et al., 2003), and between stacked PSII-LHCII (Wood et al., 2019). In vitro, suspending chloroplast membranes in low-salt medium cause grana unstacking, and addition of MgCl2 reverts the membranes back to their stacked organisation (Izawa and Good, 1966; Staehelin, 1976). Several low-resolution cryo-EM models of stacked PSII were obtained in recent years (Levitan et al., 2019; Grinzato et al., 2020; Albanese et al., 2017); these studies were also supplemented by mass spectrometry analysis detecting cross-linked regions across the stroma (Albanese et al., 2020), but a high-quality PSII structure that can shed light on the contribution of the supercomplex to thylakoid membrane stacking is missing.
The stacked PSII dataset refined to a 3.68 Å map after applying multibody refinement, with each dimer defined as a rigid body. Subsequently, the stacked particles were classified according to the higher quality PSII dimer, and two distinct populations of stacked PSII dimers were obtained, as observed for the unstacked PSII: one in the C2S2COMP conformation solved to 3.36 Å, and the other in the C2S2STR conformation solved to 3.84 Å (Figure 4; Figure 4—figure supplements 1 and 2). In both classes, the compact conformation exhibited the best fit for the second, lower resolution, PSII dimer.
Roughly 20 Å separate the two stacked dimers in both classes, as previously shown (Albanese et al., 2017). In several regions, this value decreases to approximately 10 Å (Figure 4d and l), owing to PSII stromal loops in core subunits and LHCs protruding into the space between the two dimers. Both PSII dimers are shifted by approximately 20 Å relative to each other rather than being perfectly aligned (Figure 4a; Figure 4—figure supplement 3a). Back projecting the stacked PSII onto an in-vivo observed stacked PSII shows that the dimensions of the purified stacked PSII closely match the intermembrane separation observed in vivo (Wietrzynski et al., 2020; Figure 4—figure supplement 4).
Altogether in the stacked PSII structure we do not observe any direct protein – protein interactions, this includes loops extending across the stromal gap, this contrasts with previous suggestions (Albanese et al., 2017) but can also stem from the absence of loosely associated PSII subunits (specifically PsbR). Below we discuss the extremely flexible nature of the stacked PSII dimer as revealed by multibody analysis. This is consistent with cross-linking results (Albanese et al., 2020) and strongly argues against direct protein – protein interactions across the stromal gap.
PCA showed extensive displacements and rotations across the population with stacked PSII dimers rotating relatively to their opposite dimer by as much as 30° and shifting by 80 Å in C2S2COMP, while in C2S2STR, the rotation is more restricted, showing a maximum of 19.8° and a shift of 57 Å (Figure 4 and Videos 7–10 ). The rotation axis of C2S2COMP appears to be broad region containing the N-termini stromal loops of D2 and CP29 on one dimer, and the stromal loop connecting D2 helices IV and V, CP43 N-terminus, and the C-termini of CP43, CP47, and PsbI on the opposite dimer (Figure 4—figure supplement 5). In the stacked C2S2STR, these stacking interactions also include a stromal loop from D1 which is pushed in the stromal gap by a change in the position of the PsbT C-terminus (green arrow in Figure 4k), this shift pushes this D1 loop (connecting helices IV and V) into the stromal space and closer to the adjacent dimer (Figure 4—figure supplement 3c). On the axis of rotation which consists of PSII core subunits, additional interactions between different LHCs seem to be essential to maintain stacking. All the rotation states include some degree of LHCs interaction across the stromal gap between opposite PSII dimers, and these seem to limit the extent of possible rotational states. In the stacked C2S2COMP particle set, the larger range of rotations means that at the extreme states CP26 and LHCII M2 are not involved in stacking interaction and can pair with additional complexes (Figure 4—figure supplement 5a), while in the stacked C2S2STR particle set, the smaller rotational range seem to be restricted by CP26 and LHCII M2 interactions (Figure 4—figure supplement 5b). These differences, when repeated over many stacked complexes (with additional LHCII complexes), can translate into substantial changes in thylakoid membrane stacking (Yakushevska et al., 2003).
The Mg to Mg distances between Chls in each of the stacked complexes (all above 50 Å) make EET between them inefficient. The closest protein contacts are found at the interface between core subunits from both PSII dimers and CP29, supported by peripheral interactions between LHCII trimer and CP26. Most of the PSII stromal surface is electronegative, and accordingly, most of the amino acids that seem to be involved in stacking interactions are either negatively charged or uncharged (Figure 4k). Interactions spanning 10 Å are probably insufficient to maintain PSII in its stacked arrangement; however, if mediated by a Mg2+ ion and two-to-four H2O molecules, stacking can be stabilized (Figure 4l). These interactions comply with the large degree of rotational freedom observed in the stacked dimers and with the strong dependance of stacked dimers on the presence of Mg2+ ions and may contribute to thylakoid membrane stacking (Staehelin, 1976). Indeed, a cation current counteracting the positive charges of the proton influx during light is known to occur in chloroplasts (Hind et al., 1974; Nami et al., 2021; Barber, 1980; Kaňa, 2016; Li et al., 2021; Kirchhoff et al., 2004; Puthiyaveetil et al., 2017). This has been suggested as the basis for some light-dependent alteration in the stromal spacing of thylakoid membranes (Puthiyaveetil et al., 2014; Kirchhoff et al., 2011). We suggest that the stacked PSII structure (which strongly depends on the presence of cations during purification) only relies on these weak interactions for its formation and is inherently extremely flexible in all dimensions but the PSII dimer separation distance. This flexibility may explain why stacked PSII structures are rarely detected in-vivo (Wietrzynski et al., 2020). However, when stacked PSII structures were detected using cryo-electron tomography, the identified configuration closely matched the stacked PSII dimer identified in this work (Wietrzynski et al., 2020); it is possible that in the native membrane state, range of motions in the stacked dimer or that the population of the extreme states increases, leading to the larger variability observed in vivo.
Summary
The structure of PSII from Dunaliella revealed an unexpected level of conformational flexibility in this highly conserved system. The two stable conformations appear to differ in their antennae connectivity and should be considered in PSII modelling attempts. Within each state, the large degree of structural heterogeneity also contributes to EET and may facilitate transitioning between the different states. In the stacked PSII dimer, we do not find any evidence for direct protein interaction connecting the two stacked systems, instead, long range electrostatic interaction between the core PSII subunits are flexible enough to allow for a wide range of motion, and their dependance on cation concentration provides a basis for light dependent regulation (Puthiyaveetil et al., 2017; Hind et al., 1974).
Methods
Dunaliella PSII sample preparation
D. salina (strain CONC-007) cells were cultured in a 10 l BG11 medium, supplemented with 1.5 M NaCl, 6 μg/ml ferric ammonium citrate, and 50 mM NaHCO3 at pH 8 (Caspy et al., 2020). The cells were grown with constant stirring and air bubbling under continuous white light (70 μE) at 25°C for 1 week. After reaching an OD730 of 0.4, the culture was harvested by centrifugation at 4000 g for 10 min and resuspended in a medium containing 50 mM HEPES pH 7.5, 300 mM sucrose, and 5 mM MgCl2. The cells were washed once in the same buffer and suspended in a buffer containing 25 mM MES, pH 6.5, 10 mM CaCl2, 10 mM MgCl2, 1 M betaine, 5 mM EDTA, and 12.5% glycerol. Protease-inhibitors cocktail was added to give final concentrations of 1 mM phenylmethylsulfonyl fluoride (PMSF), 1 µM pepstatin, 60 µM bestatin, and 1 mM benzamidine (Tokutsu et al., 2012). The cells were disrupted by an Avestin EmulsiFlex-C3 at 1500 psi. Unbroken cells and starch granules were removed by centrifugation at 5000 g for 5 min, and the membranes in the supernatant were precipitated by centrifugation in Ti70 rotor at 181,000 g for 1 hr. The pellet was suspended in a buffer containing 25 mM MES, pH 6.5, 10 mM CaCl2, 10 mM MgCl2, 1 M betaine, 5 mM EDTA, and 12.5% glycerol giving a Chl concentration of 0.4 mg/ml. n-Decyl-α-D-Maltopyranoside (α-DM) was added to a final concentration of 1%, and following stirring for 30 min at 4°C, the insoluble material was removed by centrifugation at 10,000 g for 5 min (Tokutsu et al., 2012). Supernatant was concentrated by centrifugation in TI-75 rotor at 377,000 g for 80 min. The pellet was suspended in the above buffer containing 0.3% α-DM at Chl concentration of about 1 mg/ml, loaded on sucrose gradients of 10 to 50% in SW-60 rotor and run at 336,000 g for 15 hr. Figure 1—figure supplement 9a shows the distribution of green bands in the tubes. The band containing PSII was concentrated by centrifugation at 550,000 g for 2 hr, and the pellet was suspended in a buffer containing 25 mM MES (pH 6.5), 1 mM CaCl2, 5 mM MgCl2, and 0.1% α-DM to give a Chl concentration of 2 mg Chl/ml. SDS-PAGE of the three bands is presented in Figure 1—figure supplement 9b. The final preparation exhibited oxygen evolution activity of 816 µmol O2/(mg Chl * hr) under 560 µmol photons * m−2 * s−1 illumination (Figure 1—figure supplement 9c).
Cryo-EM data collection and processing
Concentrated PSII solution (3 µl) was applied on glow-discharged holey carbon grids (Cu Quantifoil R1.2/1.3) that were vitrified for cryo-EM structural determination using a Vitrobot FEI (3 s blot at 4°C and 100% humidity). The images were collected using a 300 kV FEI Titan Krios electron microscope, with a slit width of 20 eV on a GIF-Quantum energy filter, at the EMBL cryo facility, Heidelberg, Germany. A Gatan Quantum K3-Summit detector was used in counting mode at a magnification of 130,000 (yielding a pixel size of 0.64 Å), with a total dose of 51.81 e Å−2. EPU was used to collect a total of 13,586 images, which were dose-fractionated into 40 movies frames, with defocus values of 0.8–1.9 μm at increments of 0.1 μm. The collected micrographs were motion-corrected and dose-weighted using MotionCor2 (Zheng et al., 2017). The contrast transfer function parameters were estimated using CtfFind v.4.1 (Rohou and Grigorieff, 2015). A total of 401,467 particles were picked using LoG reference-free picking in RELION3.1 (Zivanov et al., 2018). The picked particles were processed for reference-free two-dimensional (2D) averaging. After several rounds of 2D classification, which resulted in 253,804 particles, two initial models was generated using RELION3.1 (Zivanov et al., 2018), for the unstacked and stacked PSII.
3D classification of the unstacked PSII revealed two organisations of the LHCs surrounding the core complex – C2S and C2S2. C2S contained 21,066 particles were resampled at a pixel size of 0.896 Å, pooled together, and processed for 3D homogeneous refinement and multibody refinement (Nakane et al., 2018) using RELION3.1 (Zivanov et al., 2018), giving a final resolution of 3.61 Å. The C2S2 configuration was composed of 75,904 particles with a C2 symmetry, and these were resampled at a pixel size of 0.896 Å, pooled together, and processed for 3D homogeneous refinement and postprocessing using RELION3.1 (Zivanov et al., 2018), giving a final resolution of 2.82 Å. In an attempt to improve the map density of C2S2, mainly in the vicinity of CP29 and LHCII trimer, 3D classification without refinement was performed and revealed two distinct C2S2 conformations – compact (C2S2COMP) and stretched (C2S2STR). C2S2COMP was composed of 39,357 particles that undergone symmetry expansion, 3D homogeneous refinement, and multibody refinement (Nakane et al., 2018) in C1 symmetry to give a final resolution of 2.43 Å, and C2S2STR was composed of 23,014 particles that undergone symmetry expansion, 3D homogeneous refinement, and multibody refinement (Nakane et al., 2018) in C1 symmetry to give a final resolution of 2.62 Å.
23,874 particles that were assigned to the stacked PSII arrangement were resampled at a pixel size of 0.96 Å, pooled together, and processed for 3D homogeneous refinement and multibody refinement (Nakane et al., 2018) in C1 symmetry using RELION3.1 (Zivanov et al., 2018) and yielded a final resolution of 3.68 Å.
Focused refinement on each individual PSII complex yielded similar resolutions before multibody refinement (3.53 Å and 3.58 Å on each complex), showing both positions are occupied roughly by the same number of complexes. Focused classification was carried out on the upper dimer of the stacked PSII particles to determine if the compact and stretched conformations were also present in the stacked PSII arrangement. This analysis showed that the stacked PSII also contained a mixed population of the compact and stretched conformations. The compact set was composed of 9,567 particles, and these were pooled together and processed for 3D homogeneous refinement followed by multibody refinement (Nakane et al., 2018) to give a final resolution of 3.36 Å. The stretched set composed of 14,307 particles, these were pooled together and processed for 3D homogeneous refinement followed by multibody refinement (Nakane et al., 2018) to give a final resolution of 3.84 Å. Performing focused refinement on the lower PSII dimer of both conformations suggested a conformation mixture as well but was less conclusive due to the lower map quality of the lower PSII dimer, and both were fitted with the PSIICOMP model (using rigid body refinement) which gave the best overall fit to the map. All the reported resolutions were based on a gold-standard refinement, applying the 0.143 criterion on the Fourier shell correlation between the reconstructed half-maps. (Figure 1—figure supplement 2).
PCA was performed using relion_flex_analyse as detailed in Nakane et al., 2018. In short, the differences between each particle alignment parameters following the convergence of multibody refinement (at this point, each particle is aligned differently, optimally for each rigid body) are used to represent the heterogeneity in dataset. To generate state maps, the final, consensus, maps of each rigid body are translated along the PC axis to a position corresponding to the stated fraction of the particle population and then added together to generate the specific state map.
Maps
Focused maps obtained after multibody refinement were combined using the phenix_combined_focused_maps tool (Liebschner et al., 2019). Complete models were first refined into the consensus maps and then used to define the combined part of each map, per combined_focused_maps instructions.
Model building
To generate the C2S2 PSII, the cryo-EM structure of the C2S2 C. reinhardtii PSII model PDB 6KAC (Sheng et al., 2019) was selected. This model was fitted onto the cryo-EM density map using phenix.dock_in_map in the PHENIX suite (Liebschner et al., 2019) and manually rebuilt using Coot (Emsley et al., 2010). Stereochemical refinement was performed using phenix.real_space_refine in the PHENIX suite (Liebschner et al., 2019). The final model was validated using MolProbity (Chen et al., 2010). The refinement statistics are provided in Supplementary file 1. Local resolution was determined using ResMap (Kucukelbir et al., 2014), and the figures were generated using UCSF Chimera (Pettersen et al., 2004) and UCSF ChimeraX (Goddard et al., 2018). Representative cryo-EM densities are shown in Figure 1—figure supplement 3.
Data availability
The atomic coordinates have been deposited in the Protein Data Bank, with accession code 7PI0 (C2S2 COMP ), 7PI5 (C2S2 STR), 7PNK (C2S), 7PIN (stacked C2S2 COMP ) and 7PIW (stacked C2S2 STR ). The cryo-EM maps have been deposited in the Electron Microscopy Data Bank, with accession codes EMD-13429 (C2S2 COMP ), EMD-13430 (C2S2 STR ), EMD-13548 (C2S), EMD-13444 (stacked C2S2 COMP ) and EMD-13455 (stacked C2S2 STR ).
-
RCSB Protein Data BankID 7PI0. Unstacked compact Dunaliella PSII.
-
RCSB Protein Data BankID 7PI5. Unstacked stretched Dunaliella PSII.
-
RCSB Protein Data BankID 7PNK. Unstacked compact Dunaliella PSII.
-
RCSB Protein Data BankID 7PIN. Stacked compact Dunaliella PSII.
-
RCSB Protein Data BankID 7PIW. Stacked stretched Dunaliella PSII.
References
-
Novel features of eukaryotic photosystem II revealed by its crystal structure analysis from a red algaThe Journal of Biological Chemistry 291:5676–5687.https://doi.org/10.1074/jbc.M115.711689
-
The relationship between thylakoid stacking and salt induced chlorophyll fluorescence changesBiochimica et Biophysica Acta 591:92–103.https://doi.org/10.1016/0005-2728(80)90223-6
-
Organization of transmembrane helices in photosystem II: comparison of plants and cyanobacteriaPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 357:1329–1335.https://doi.org/10.1098/rstb.2002.1132
-
Water, water everywhere, and its remarkable chemistryBiochimica et Biophysica Acta 1655:123–132.https://doi.org/10.1016/j.bbabio.2003.10.011
-
Functional architecture of higher plant photosystem II supercomplexesThe EMBO Journal 28:3052–3063.https://doi.org/10.1038/emboj.2009.232
-
Structure and energy transfer pathways of the dunaliella salina photosystem I supercomplexBiochimica et Biophysica Acta. Bioenergetics 1861:148253.https://doi.org/10.1016/j.bbabio.2020.148253
-
MolProbity: all-atom structure validation for macromolecular crystallographyActa Crystallographica. Section D, Biological Crystallography 66:12–21.https://doi.org/10.1107/S0907444909042073
-
The significance of CP29 reversible phosphorylation in thylakoids of higher plants under environmental stressesJournal of Experimental Botany 64:1167–1178.https://doi.org/10.1093/jxb/ert002
-
Light harvesting in a fluctuating antennaJournal of the American Chemical Society 136:8963–8972.https://doi.org/10.1021/ja5027858
-
CAVER 3.0: a tool for the analysis of transport pathways in dynamic protein structuresPLOS Computational Biology 8:e1002708.https://doi.org/10.1371/journal.pcbi.1002708
-
Current understanding of the mechanism of water oxidation in photosystem II and its relation to XFEL dataAnnual Review of Biochemistry 89:795–820.https://doi.org/10.1146/annurev-biochem-011520-104801
-
Light-Harvesting and structural organization of photosystem II: from individual complexes to thylakoid membraneJournal of Photochemistry and Photobiology. B, Biology 104:142–153.https://doi.org/10.1016/j.jphotobiol.2011.02.015
-
Light-harvesting complex II (LHCII) and its supramolecular organization in Chlamydomonas reinhardtiiBiochimica et Biophysica Acta 1837:63–72.https://doi.org/10.1016/j.bbabio.2013.07.012
-
Features and development of cootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
Role of LHCII-containing macrodomains in the structure, function and dynamics of granaFunctional Plant Biology 27:723.https://doi.org/10.1071/PP99069_C1
-
Sharing light between two photosystems: mechanism of state transitionsCurrent Opinion in Plant Biology 25:71–78.https://doi.org/10.1016/j.pbi.2015.04.009
-
High-Light versus low-light: effects on paired photosystem II supercomplex structural rearrangement in pea plantsInternational Journal of Molecular Sciences 21:8643.https://doi.org/10.3390/ijms21228643
-
Structure and membrane organization of photosystem II in green plantsAnnual Review of Plant Physiology and Plant Molecular Biology 48:641–671.https://doi.org/10.1146/annurev.arplant.48.1.641
-
The relevance of dynamic thylakoid organisation to photosynthetic regulationBiochimica et Biophysica Acta. Bioenergetics 1861:148039.https://doi.org/10.1016/j.bbabio.2019.06.011
-
Role of ions in the regulation of light-harvestingFrontiers in Plant Science 7:1849.https://doi.org/10.3389/fpls.2016.01849
-
Entropy-assisted stacking of thylakoid membranesBiochimica et Biophysica Acta - Bioenergetics 1708:187–195.https://doi.org/10.1016/J.BBABIO.2005.03.011
-
Structural and functional self-organization of photosystem II in grana thylakoidsBiochimica et Biophysica Acta 1767:1180–1188.https://doi.org/10.1016/j.bbabio.2007.05.009
-
Structural changes of the thylakoid membrane network induced by high light stress in plant chloroplastsPhilosophical Transactions of the Royal Society of London. Series B, Biological Sciences 369:20130225.https://doi.org/10.1098/rstb.2013.0225
-
The structural and functional domains of plant thylakoid membranesThe Plant Journal 97:412–429.https://doi.org/10.1111/tpj.14127
-
Quantifying the local resolution of cryo-EM density mapsNature Methods 11:63–65.https://doi.org/10.1038/nmeth.2727
-
Macromolecular structure determination using X-rays, neutrons and electrons: recent developments in phenixActa Crystallographica. Section D, Structural Biology 75:861–877.https://doi.org/10.1107/S2059798319011471
-
Effects of chilling on the structure, function and development of chloroplastsFrontiers in Plant Science 9:1715.https://doi.org/10.3389/fpls.2018.01715
-
Design principles of solar light harvesting in plants: functional architecture of the monomeric antenna CP29Biochimica et Biophysica Acta. Bioenergetics 1861:148156.https://doi.org/10.1016/j.bbabio.2020.148156
-
Molecular responses of legumes to abiotic stress: post-translational modifications of proteins and redox signalingJournal of Experimental Botany 72:5876–5892.https://doi.org/10.1093/jxb/erab008
-
Water-Splitting chemistry of photosystem IIChemical Reviews 106:4455–4483.https://doi.org/10.1021/cr0204294
-
The complex architecture of oxygenic photosynthesisNature Reviews. Molecular Cell Biology 5:971–982.https://doi.org/10.1038/nrm1525
-
Structure and function of photosystems I and IIAnnual Review of Plant Biology 57:521–565.https://doi.org/10.1146/annurev.arplant.57.032905.105350
-
Structure and energy transfer in photosystems of oxygenic photosynthesisAnnual Review of Biochemistry 84:659–683.https://doi.org/10.1146/annurev-biochem-092914-041942
-
Excitation dynamics in the LHCII complex of higher plants: modeling based on the 2.72 angstrom crystal structureThe Journal of Physical Chemistry. B 109:10493–10504.https://doi.org/10.1021/jp044082f
-
UCSF chimera -- a visualization system for exploratory research and analysisJournal of Computational Chemistry 25:1605–1612.https://doi.org/10.1002/jcc.20084
-
Sodium dependency of the photosynthetic electron transport in the alkaliphilic cyanobacterium Arthrospira platensisJournal of Bioenergetics and Biomembranes 35:427–437.https://doi.org/10.1023/a:1027339814544
-
Structure and dynamics of thylakoids in land plantsJournal of Experimental Botany 65:1955–1972.https://doi.org/10.1093/jxb/eru090
-
CTFFIND4: fast and accurate defocus estimation from electron micrographsJournal of Structural Biology 192:216–221.https://doi.org/10.1016/j.jsb.2015.08.008
-
The structure of photosystem II and the mechanism of water oxidation in photosynthesisAnnual Review of Plant Biology 66:23–48.https://doi.org/10.1146/annurev-arplant-050312-120129
-
Reversible particle movements associated with unstacking and restacking of chloroplast membranes in vitroThe Journal of Cell Biology 71:136–158.https://doi.org/10.1083/jcb.71.1.136
-
Thylakoid membrane fluidity and thermostability during the operation of the xanthophyll cycle in higher-plant chloroplastsBiochimica et Biophysica Acta 1330:179–193.https://doi.org/10.1016/s0005-2736(97)00168-5
-
Revisiting the supramolecular organization of photosystem II in Chlamydomonas reinhardtiiThe Journal of Biological Chemistry 287:31574–31581.https://doi.org/10.1074/jbc.M111.331991
-
How photosynthetic proteins switchThe Journal of Physical Chemistry Letters 3:2779–2784.https://doi.org/10.1021/jz300983r
-
Charge separation is virtually irreversible in photosystem II core complexes with oxidized primary quinone acceptorThe Journal of Physical Chemistry. A 115:3947–3956.https://doi.org/10.1021/jp1083746
-
Photosystem II: the reaction center of oxygenic photosynthesisAnnual Review of Biochemistry 82:577–606.https://doi.org/10.1146/annurev-biochem-070511-100425
Article and author information
Author details
Funding
Israel Science Foundation (569/17)
- Nathan Nelson
Israel Science Foundation (199/21)
- Nathan Nelson
German-Israeli Foundation for Scientific Research and Development (G-1483-207/2018)
- Nathan Nelson
National Science Foundation (2034021)
- Yuval Mazor
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
Dr Yael Levi-Kalisman is gratefully acknowledged and thanked for vitrifying the samples. We also thank the Electron Microscopy Core Facility (EMCF) at the European Molecular Biology Laboratory (EMBL) for their support and Felix Weis for data collection and excellent technical support. Molecular graphics and analyses were performed with UCSF Chimera, developed by the Resource for Biocomputing, Visualisation, and Informatics at the University of California, San Francisco, with support from NIH P41-GM103311. Molecular graphics and analyses performed with UCSF ChimeraX, developed by the Resource for Biocomputing, Visualisation, and Informatics at the University of California, San Francisco, with support from National Institutes of Health R01-GM129325 and the Office of Cyber Infrastructure and Computational Biology, National Institute of Allergy and Infectious Diseases. This work was supported by The Israel Science Foundation (Grants No. 569/17 and 199/21), and by German-Israeli Foundation for Scientific Research and Development (GIF), Grant no. G-1483-207/2018. Y.M acknowledges the support by the National Science Foundation under Award No. 2034021.
Copyright
© 2023, Caspy et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,218
- views
-
- 179
- downloads
-
- 7
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Cell Biology
- Plant Biology
Plants distribute many nutrients to chloroplasts during leaf development and maturation. When leaves senesce or experience sugar starvation, the autophagy machinery degrades chloroplast proteins to facilitate efficient nutrient reuse. Here, we report on the intracellular dynamics of an autophagy pathway responsible for piecemeal degradation of chloroplast components. Through live-cell monitoring of chloroplast morphology, we observed the formation of chloroplast budding structures in sugar-starved leaves. These buds were then released and incorporated into the vacuolar lumen as an autophagic cargo termed a Rubisco-containing body. The budding structures did not accumulate in mutants of core autophagy machinery, suggesting that autophagosome creation is required for forming chloroplast buds. Simultaneous tracking of chloroplast morphology and autophagosome development revealed that the isolation membranes of autophagosomes interact closely with part of the chloroplast surface before forming chloroplast buds. Chloroplasts then protrude at the site associated with the isolation membranes, which divide synchronously with autophagosome maturation. This autophagy-related division does not require DYNAMIN-RELATED PROTEIN 5B, which constitutes the division ring for chloroplast proliferation in growing leaves. An unidentified division machinery may thus fragment chloroplasts for degradation in coordination with the development of the chloroplast-associated isolation membrane.
-
- Plant Biology
Photosynthetic organisms exhibit remarkable diversity in their light-harvesting complexes (LHCs). LHCs are associated with photosystem I (PSI), forming a PSI-LHCI supercomplex. The number of LHCI subunits, along with their protein sequences and pigment compositions, has been found to differ greatly among the PSI-LHCI structures. However, the mechanisms by which LHCIs recognize their specific binding sites within the PSI core remain unclear. In this study, we determined the cryo-electron microscopy structure of a PSI supercomplex incorporating fucoxanthin chlorophyll a/c-binding proteins (FCPs), designated as PSI-FCPI, isolated from the diatom Thalassiosira pseudonana CCMP1335. Structural analysis of PSI-FCPI revealed five FCPI subunits associated with a PSI monomer; these subunits were identified as RedCAP, Lhcr3, Lhcq10, Lhcf10, and Lhcq8. Through structural and sequence analyses, we identified specific protein–protein interactions at the interfaces between FCPI and PSI subunits, as well as among FCPI subunits themselves. Comparative structural analyses of PSI-FCPI supercomplexes, combined with phylogenetic analysis of FCPs from T. pseudonana and the diatom Chaetoceros gracilis, underscore the evolutionary conservation of protein motifs crucial for the selective binding of individual FCPI subunits. These findings provide significant insights into the molecular mechanisms underlying the assembly and selective binding of FCPIs in diatoms.