# Calredoxin represents a novel type of calcium-dependent sensor-responder connected to redox regulation in the chloroplast

Ana Karina Hochmal, Karen Zinzius, Ratana Charoenwattanasatien, Philipp Gäbelein, Risa Mutoh, Hideaki Tanaka, Stefan Schulze, Gai Liu, Martin Scholz, André Nordhues, et al.

<figure><img src="image_1.png" /></figure>

## ► To cite this version:

Ana Karina Hochmal, Karen Zinzius, Ratana Charoenwattanasatien, Philipp Gäbelein, Risa Mutoh, et al.. Calredoxin represents a novel type of calcium-dependent sensor-responder connected to redox regulation in the chloroplast. Nature Communications, 2016, 7, pp.11847. 10.1038/ncomms11847. hal-01412644

HAL Id: hal-01412644

https://hal.science/hal-01412644

Submitted on 27 May 2020

HAL is a multi-disciplinary open access archive for the deposit and dissemination of scientific research documents, whether they are published or not. The documents may come from teaching and research institutions in France or abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est destinée au dépôt et à la diffusion de documents scientifiques de niveau recherche, publiés ou non, émanant des établissements d'enseignement et de recherche français ou étrangers, des laboratoires publics ou privés.

<figure><img src="image_2.png" /></figure>
<figure><img src="image_3.png" /></figure>

# Calredoxin represents a novel type of calcium-dependent sensor-responder connected to redox regulation in the chloroplast

Ana Karina Hochmal<fnref n="1" />,*, Karen Zinzius<fnref n="1" />,*, Ratana Charoenwattanasatien²,*, Philipp Gäbelein<fnref n="1" />, Risa Mutoh²,³, Hideaki Tanaka²,³, Stefan Schulze<fnref n="1" />, Gai Liu<fnref n="4" />, Martin Scholz<fnref n="1" />, André Nordhues<fnref n="1" />, Jan Niklas Offenborn<fnref n="1" />, Dimitris Petroutsos⁵,⁶,⁷,⁸, Giovanni Finazzi⁵,⁶,⁷,⁸, Christian Fufezan<fnref n="1" />, Kaiyao Huang<fnref n="4" />, Genji Kurisu²,³ & Michael Hippler<fnref n="1" />

Calcium ($Ca^{2+}$) and redox signalling play important roles in acclimation processes from archaea to eukaryotic organisms. Herein we characterized a unique protein from *Chlamydomonas reinhardtii* that has the competence to integrate $Ca^{2+}$- and redox-related signalling. This protein, designated as calredoxin (CRX), combines four $Ca^{2+}$-binding EF-hands and a thioredoxin (TRX) domain. A crystal structure of CRX, at 1.6 Å resolution, revealed an unusual calmodulin-fold of the $Ca^{2+}$-binding EF-hands, which is functionally linked via an inter-domain communication path with the enzymatically active TRX domain. CRX is chloroplast-localized and interacted with a chloroplast 2-Cys peroxiredoxin (PRX1). $Ca^{2+}$-binding to CRX is critical for its TRX activity and for efficient binding and reduction of PRX1. Thereby, CRX represents a new class of $Ca^{2+}$-dependent 'sensor-responder' proteins. Genetically engineered *Chlamydomonas* strains with strongly diminished amounts of CRX revealed altered photosynthetic electron transfer and were affected in oxidative stress response underpinning a function of CRX in stress acclimation.

<footnote n="1">¹Institute of Plant Biology and Biotechnology, University of Münster, 48143 Münster, Germany. ²Institute for Protein Research, Osaka University, Suita Osaka 565-0871, Japan. ³Core Research for Evolutional Science and Technology (CREST), Japan Science and Technology Agency (JST), Saitama 332-0012, Japan.</footnote>

<footnote n="4">⁴Key Laboratory of Algal Biology, Institute of Hydrobiology, Chinese Academy of Sciences, Wuhan, Hubei 430072, China. ⁵Centre National Recherche Scientifique, Unité Mixte Recherche 5168, Laboratoire Physiologie Cellulaire et Végétale, F-38054 Grenoble, France. ⁶Commissariat à l'Energie Atomique et Energies Alternatives, l'Institut de Recherches en Technologies et Sciences pour le Vivant, F-38054 Grenoble, France. ⁷Université Grenoble 1, F-38041 Grenoble, France. ⁸Institut National Recherche Agronomique, UMR1200, F-38054 Grenoble, France. * These authors contributed equally to this work.</footnote>

<footnote>Correspondence and requests for materials should be addressed to M.H. (email: mhippler@uni-muenster.de).</footnote>
Living organisms respond to dynamic alterations in their environment by matching acclimation strategies to cope with impending stresses. For appropriate cellular responses, cells measure changes in the environment by signal perception and decode signals via signal-transduction pathways to tune their physiology accordingly. In this signal-transduction leitmotif, calcium ($Ca^{2+}$) and redox signalling play important roles in all types of cells from archaea to eukaryotic organisms. For sensing and promotion of $Ca^{2+}$ derived signals, $Ca^{2+}$ sensor proteins are crucial and mostly characterized by the presence of one or more EF-hand $Ca^{2+}$-binding motifs. This motif is highly conserved and EF-hand proteins are encoded in all eukaryotic genomes¹. Notably, ~250 EF-hand proteins are encoded in the Arabidopsis genome, more than in any other organism investigated to date².

In line, $Ca^{2+}$-dependent responses are central in acclimation of plants towards environmental changes. Numerous environmental cues cause fluctuations in the cytosolic-free $Ca^{2+}$ concentration, fundamental in the initiation of the appropriate physiological response of the plant (for review ref. 3). Two types of signalling components decode changes in cellular $Ca^{2+}$ concentration in plants. Type I 'sensor-responder' proteins possess both $Ca^{2+}$-binding and enzymatic 'effector' domains. Type II components, such as calmodulin (CaM), have a $Ca^{2+}$-binding domain but do not exhibit an enzymatic activity and are designated as 'sensor-relay' proteins. CaM are small proteins and possess a pair of $Ca^{2+}$-binding EF-hand motifs⁴. Interestingly, in plants the CaM family is extended by a large number of CaM-like proteins (CMLs), which differ from canonical CaMs by variations in length and/or by their number of EF-hands⁵. The physiological role of CaMs and CMLs in plants has been described thoroughly (for recent reviews see refs 5,6). Although apparent CaM-dependent processes have also been described for chloroplasts (as reviewed in refs 7,8), no chloroplast-localized CaMs and CMLs have been identified so far.

Although, there is increasing evidence that chloroplasts are part of the cellular $Ca^{2+}$ network and contribute to the cytosolic $Ca^{2+}$ signalling, thereby modulating the action of downstream cytosolic 'sensor-responder' and 'sensor-relay' proteins. $Ca^{2+}$ is taken up into the chloroplast in the light but released into the cytosol in the dark (as reviewed in ref. 8). The chloroplast-localized $Ca^{2+}$ sensor protein (CAS) may contribute to cellular $Ca^{2+}$ signalling. CAS function is crucial for stomatal regulation⁹,¹⁰ and chloroplast-mediated activation of immune signalling¹¹ in *A. thaliana*. Moreover, it is involved in photo-acclimation in *C. reinhardtii* by contributing to the control of expression of LHCSR3 (ref. 12). LHCSR3 is crucial for qE (ref. 13), the energy-dependent component of non-photochemical quenching, required for efficient light to heat dissipation. Effective qE requires the acidification of the thylakoid lumen and is therefore functionally connected to photosynthetic cyclic electron flow (CEF), which contributes to the pH-gradient across the thylakoid membrane¹⁴. In microalgae and vascular plants, CEF operates via an NAD(P)H dehydrogenase (NDH)-dependent and/or PROTONGRADIENT REGULATION5 (PGR5)-related pathway¹⁵. The thylakoid protein PGR5-Like 1 (PGRL1) participates in the PGR5-dependent CEF pathway¹⁶⁻¹⁸ and associates together with PGR5 with photosystem I (PSI), thereby facilitating the operation of CEF (ref. 18). Recent data further indicated that PGRL1 might operate as a ferredoxin-plastoquinone reductase¹⁹. Activity of CEF is modulated via the availability of $Ca^{2+}$ (ref. 20). Moreover, there is evidence that CEF is redox controlled²¹ in *C. reinhardtii* and activated by hydrogen peroxide in Arabidopsis²². The redox control could operate via the chloroplast thioredoxin (TRX) system and its redox regulation of the PGR5 and PGRL1 cycle¹⁹,

required for PGRL1 homodimer to monomer conversion. In line, reversible activation/inactivation of CEF has been described to operate with an apparent midpoint potential of -306 mV, consistent with a TRX-mediated redox modulation of a thiol/disulfide couple of PGRL1 (ref. 23).

TRXs are protein oxidoreductases that harbour a redox-active dithiol/disulfide motif in their active site. In its reduced state this motif allows cleavage of a disulfide bond of a target protein via a bimolecular nucleophilic substitution reaction. Oxidized TRXs are reduced by TRX reductases (TRXRs). In contrast to other organisms, plants and cyanobacteria possess two independent TRXRs, a ferredoxin-dependent (FTR) and a NADPH-dependent (NTR) and various types of TRXs²⁴. NTRC, which resides in the chloroplast, is present in cyanobacteria and plants and has a TRX sequence fused to the C-terminal sequence of NTR²⁵. TRXs play a central role in redox regulation and are interlinked with signal-transduction pathways in acclimation responses²⁶. The TRX system is in particular important for chloroplasts and the regulation of photosynthesis²⁴ as evidenced by, for example, its involvement in metabolic regulation and anti-oxidative stress response. Here, the reduction of oxidized 2-Cys peroxiredoxins (PRXs), which catalytically reduce hydrogen peroxide and organic peroxides, by the TRX system, such as ACHT1 (ref. 27) and NTRC²⁸, provides a strong link of the latter to antioxidant metabolism²⁹. This is exemplified by the fact that PRXs also operate as redox sensors³⁰, thereby further strengthening the link between ROS and redox signalling networks.

Herein we provide evidence that $Ca^{2+}$- and redox-dependent regulation is directly linked within one protein, involved in the regulation of photosynthetic redox regulation and oxidative stress defense. This protein is designated as calredoxin (CRX, Cre03.g202950.t1.1) and consists of 4 EF-hands that are functionally connected to a TRX domain. In this work we characterized functional properties of CRX, which is conserved in the green algae lineage³¹ but absent in vascular plants.

## Results

### Recombinant CRX displays $Ca^{2+}$-dependent TRX activity.

Recombinant CRX was isolated and purified via $Ni^{2+}$-NTA affinity chromatography. Binding of divalent cations to recombinant CRX was measured using Microscale Thermophoresis in the presence of calcium and magnesium. For calculation of $Ca^{2+}$ and $Mg^{2+}$ dissociation constants, the ratio between the fluorescence of the fluorescently labelled protein before and after the thermophoretic movement was determined at increasing cation concentrations (Fig. 1a). An increase in $Mg^{2+}$ concentration did not change the fluorescence ratio. On the other hand, an increase in $Ca^{2+}$ concentration led to a clear decrease in the fluorescence ratio, indicating specific binding of $Ca^{2+}$ but not $Mg^{2+}$ to CRX. The $Ca^{2+}$ response data were fitted, resulting in a dissociation constant ($K_d$ value) of 88.2 ± 16.5 nM. To assess the activity of the TRX domain, an assay was employed in which CRX was reduced by NTR and NADPH, and its activity was measured photo-metrically via reduction of 5,5'-dithiobis-(2-nitrobenzoic acid) (DTNB). To determine whether the TRX activity of CRX depends on the presence of calcium, the activity assay was performed at increasing concentrations of $Ca^{2+}$ (Fig. 1b). As shown, the TRX domain of CRX is active and its full activity requires calcium. Fitting the TRX- and $Ca^{2+}$-dependent activity curve by Michaelis-Menten kinetics revealed a half maximum activity at 281.1 ± 153.8 nM free $Ca^{2+}$. Importantly, a site-directed mutant of CRX where both cysteine residues of the TRX active site were altered to Ser (C238S, C241S) exhibited only minor DTNB reduction activity, proving that an active TRX domain is required
for CRX-dependent substrate reduction (Supplementary Fig. 1a). Notably, the overall $Ca^{2+}$ dissociation constant as measured by Microscale Thermophoresis (88 nM) is slightly smaller to the value obtained from the TRX activity assay (281 nM), suggesting that for TRX activity all four EF-hands are occupied by $Ca^{2+}$. Next, we determined the redox midpoint potential of the TRX active domain by measuring the disulfide/dithiol redox state using the monobromobimane fluorescence technique³². The resulting redox titration curve was fitted by the Nernst equation (Fig. 1c). The best fit is represented by a two-electron Nernst curve (see solid line, Fig. 1c) and yields a redox midpoint potential of

$-288.2 \pm 5.3$ mV, which is close to the redox midpoint potential of chloroplast TRX *f* from spinach and pea³².

**CRX is localized in the chloroplast.** ChloroP³³ and PredAlgo³⁴ protein sequence analyses indicated that the nuclear encoded CRX possesses a transit peptide and is accordingly predicted to reside in the chloroplast. To verify this localization, chloroplasts were isolated and assayed for the presence of CRX. Immunoblot analyses revealed that CRX is enriched in chloroplasts versus whole cells and isolated mitochondria (Supplementary Fig. 2). Independently, localization of CRX was assessed by engineering strains that express a *ble-2A-crx-mVenus* construct, a technique utilized previously for subcellular localization studies in *Chlamydomonas*³⁵. Confocal laser-scanning microscopy revealed presence of YFP fluorescence in transformed *Chlamydomonas* cells (Fig. 2a,d). The YFP fluorescence partially overlaps with chlorophyll autofluorescence (Fig. 2b,e) as also visualized in the merged images (Fig. 2c,f). The CRX-YFP fluorescence is nearly identical with recently documented mCherry fluorescence signals stemming from a construct targeted to the chloroplast stroma via the *psad* transit peptide sequence³⁶. Therefore, we conclude that CRX is localized in the chloroplast stroma of *Chlamydomonas*. Chloroplast localization of CRX was further independently confirmed by transient expression of a *crx-mCherry* construct in *Nicotiana benthamiana* epidermal leaf cells (Supplementary Fig. 3a–c). Here, the *Chlamydomonas crx* transit peptide drives chloroplast localization of CRX-mCherry in tobacco, as shown for other transiently expressed *Chlamydomonas* green fluorescent protein constructs having also their endogenous transit peptide sequences²⁰.

**CRX interacts with chloroplast 2-cys peroxiredoxin.** To identify CRX protein-protein interaction partners, we have altered the CRX cysteine residues of the active disulfide to serine and performed thiol-trapping experiments as described before³⁷⁻³⁹. To this end, wild type (WT) and mutated versions of recombinant CRX were immobilized on a CNBr-activated resin and a whole-cell lysate of photoheterotrophically grown *C. reinhardtii* was added to the column to supply potential CRX target proteins to the column. The C241S mutated CRX was most efficient in trapping potential targets due to the missing

**Figure 1 | In vitro characterization of calredoxin.** (a) $K_d$ estimation for $Ca^{2+}$-binding by MST measurements. Fluorescently labelled recombinant WT CRX (CRX) was incubated with defined concentrations of free $Ca^{2+}$ (filled circles) or $Mg^{2+}$ (open circles) and the ratio of detected fluorescence before and after the thermophoretic movement was plotted against the corresponding cation concentration. Data for CRX incubated with $Ca^{2+}$ were fitted according to the law of mass action (black line) and gave a $K_d$ of 88.2 nM $\pm$ 16.5 nM. Each data point represents the mean value of at least three experiments ($\pm$ s.d.). (b) CRX shows $Ca^{2+}$-dependent redox activity. 10 µM recombinant WT CRX (closed circles) was reduced by *E. coli* NTR and NADPH in defined $Ca^{2+}$ concentrations for 10 min at RT. 200 µM DTNB were added as substrate for reduction by CRX and the increase in absorption at 412 nm was recorded to calculate the redox activity (slope 0–80 s after addition of DTNB). Data were normalized on the highest activity measured for each protein purification and fitted by Michaelis–Menten kinetics ($K_d$: 281.1 $\pm$ 153.8 nM). Error bars represent s.d. of three independent measurements. Assay modified after ref. 62. (c) Oxidation-reduction titration of WT CRX. The disulfide/dithiol redox state at each $E_h$ value was monitored using the monobromobimane fluorescence method. The line represents a fit of the data to a two-electron Nernst curve and yielded an $E_h$ of $-288.2 \pm 5.3$ mV. Data were acquired in two independent experiments.

<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>**Figure 2 |** Chloroplast localization of calredoxin in *C. reinhardtii*.</figcaption></figure>

Microscopy images of a transgenic strain expressing ble-2A-calredoxin-mVenus. (a,d) YFP fluorescence (detected with a 525-555 nm filter) and (b,e) chlorophyll fluorescence (AUF, detected with a 690-740 nm filter). Merged images (YFP/AUF) are shown in the bottom row (c,f). The second row (d-f) shows the fluorescence signals in a single cell at higher magnification. Scale bars, 20 µm.

resolving cysteine in the active site of the TRX domain. The potential targets were eluted with 10 mM dithiothreitol (DTT), digested tryptically and analysed by mass spectrometry. The log2 intensities after label-free protein quantification of the C241S sample were plotted either against the intensities of the WT (Fig. 3a) or the C238S (Fig. 3b) sample. Most proteins were detected at equal ratios (black line represents ratio = 1). Two proteins were significantly enriched in the C241S sample: PRX1, Cre06.g257601.t1.2 (open square) and another 2-cys PRX, Cre02.g114600.t1.2 (open circle; see also Supplementary Data 1). In an independent experiment whole-cell lysates stemming from photoautotrophically grown *Chlamydomonas* cells were subjected to the thiol-trapping protocol (Supplementary Fig. 4a-c, Supplementary Data 2). In concurrence, PRX1 was the only candidate significantly higher, abundant in the C241S sample. Notably, PRX1 is a chloroplast-localized 2-cys PRX³⁹. Importantly, trapping experiments in the absence of Ca²⁺ revealed particularly diminished enrichment of PRX1, suggesting that binding of CRX to PRX1 requires Ca²⁺ (Fig. 3c,d). Instead, a new protein, a TRX-like protein with unknown function (TRXL1, Cre03.g157800.t1.1) was higher abundant in C241S-dependent enrichment as compared with the two controls (Fig. 3c,d). An activity assay confirmed that CRX is indeed able to reduce PRX1 (Fig. 3e) and drive detoxification of H₂O₂ via PRX1 *in vitro*. Moreover, the reduction of PRX1 via CRX is dependent on Ca²⁺ as addition of EGTA diminishes the reduction of PRX1 by CRX about 10-fold (Fig. 3e). A titration of the electron transfer between PRX1 and CRX in dependence of Ca²⁺ revealed a half-maximal rate of NADPH oxidation at a concentration of 122.3 ± 64.5 nM free Ca²⁺, a value close to the one found in the Microscale Thermophoresis experiment (Fig. 1a). Notably, the rate of PRX1 reduction via CRX, with a value between 30-44 µmol NADPH min⁻¹ µmol⁻¹ PRX1 (Fig. 3e,f), is almost as efficient as described for NTRC and PRX in Arabidopsis²⁸, thereby supporting a role of CRX in ROS defense as electron donor to PRX1 in chloroplasts of *C. reinhardtii*.

**Crystal structure of CRX.** Recombinant CRX was purified from the soluble supernatant with and without adding calcium ions and both were used for crystallization. Although extensive crystallization trials were carried out, only the sample containing calcium could be crystallized. Analytical X-ray fluorescence spectroscopy of the sample without Ca²⁺ revealed that the purified CRX still had small amounts of Ca²⁺ bound even after EGTA treatment and gel-filtration chromatography, thereby causing structural non-homogeneity and hindering crystallization of CRX without Ca²⁺.

Crystals of CRX were obtained at 20 °C. After optimization of the conditions, plate-shaped crystals were grown within 2–3 days using 1.5 M LiSO₄ and 0.1 M MES-NaOH buffer pH 6 as a precipitant. Crystals directly dipped into liquid nitrogen for the cryo-experiment belonged to the monoclinic space group of P2₁ (Table 1). The maximum resolution of the collected X-ray data set was 1.6 Å. To solve the phase problem, we prepared selenomethionine (SeMet) substituted proteins. The structure was determined using the single-wavelength anomalous dispersion method at 2.8 Å with the SeMet derivative crystal. The phases of the native crystal data were obtained by the molecular replacement method and the structure was refined to 1.6 Å resolution. The crystallographic data and refinement statistics are listed in Table 1(for further details see Supplementary experimental procedures).

The crystal contained two CRX molecules per asymmetric unit and the electron density map was clear enough to refine the model without any non-crystallographic restraint. Molecules A and B contained 313 and 310 amino acid residues, respectively, except for flexible regions, four calcium ions in each CRX molecule and 1,156 water molecules in total. Since the model of molecule A contained more structural information than that of molecule B and those were almost identical (r.m.s.d = 0.71 Å for Cα carbons), we described the structure of monomeric CRX by referring to molecule A unless otherwise mentioned (Fig. 4a).

The CRX molecule consists of two domains connected by a flexible linker region: the CaM domain with four Ca²⁺-binding EF-hand loops and the TRX domain with one disulfide bridge. In the CaM domain, one calcium ion was bound in each EF-hand loop located in the N- and C-subdomains, respectively. There were nine α-helices (from H1 to H9) in the CaM domain, two of which sandwiched each EF-hand motif. An additional structured loop was inserted between the helices H2 and H3 and the last helix, H9, was connected to the flexible linker region. The TRX domain exhibited the typical TRX fold, possessing four α-helices (H10–H13) and one mixed β-sheet. Side chains of Cys238 and Cys241 formed the disulfide bridge between the second β-strand and the H11 helix. We have added 1 mM DTT in the crystallization droplets, nevertheless the resultant structure contained the disulfide bridge, likely due to the oxidation of DTT during crystallization and to the negative redox midpoint potential of CRX.

At the interface of the two domains, extensive inter-domain interactions were found (Fig. 4b,c). Eighteen residues from the CaM domain and 21 from the TRX domain were involved in the direct inter-domain interaction. Among them, Ser123 and Asn125 provided ligands to the third calcium ion (Ca3; Fig. 4b,d) and Ser161 and Gln163 to the fourth (Ca4; Fig. 4b). Additionally, four residues (Thr127, Glu132, Glu147 and Asp157) from the CaM domain and six residues (Glu215, Glu216, Thr236, Lys242, Lys263 and Asn267) from the TRX domain mediated the inter-domain hydrogen bond network through water molecules (partly shown in Fig. 4b,c). The inter-domain networks between the Ca4 calcium ion in the CaM domain and the disulfide bridge in the TRX domain is shown in Fig. 4c. Strikingly, when the TRX domain is expressed without the CaM domain, purified from
<figure><img src="image_8.png" /><figcaption>**Figure 3** | Potential calredoxin interaction partners. (a–d) Results from the CRX affinity chromatography. WT and mutated versions of recombinant CRX were immobilized on a CNBr-activated resin and a whole-cell lysate of heterotrophically grown *C. reinhardtii* was added to the column to supply potential CRX target proteins. The log2 protein intensities after label-free quantification (LFQ) of the C241S sample were plotted either against the intensities of the WT (a,c) or the C238S (b,d) sample. Two proteins were repeatedly (two experiments) significantly more abundant in the C241S sample when $Ca^{2+}$ was present on the column (a,b): PRX1, Cre06.g257601.t1.2 (open square) and another 2-cys peroxiredoxin, Cre02.g114600.t1.2 (open circle). Elimination of $Ca^{2+}$ (c,d) reduced the abundance of these proteins and led to identification of a third potential target protein: TRXL1 (Cre03.g157800.t1.1, open triangle). (e) Interaction of CRX and PRX1 in vitro. 5 µM recombinant CRX was reduced by *E. coli* TRXR and NADPH in the presence of 40 µM $H_2O_2$ at RT. The NADPH absorbance at 340 nm was monitored until a steady decrease was observed. After subsequent addition of 1 µM oxidized recombinant PRX1 (indicated by the arrow) NADPH oxidation was increased in the presence of $Ca^{2+}$ (black filled circle, 44.3 µmol NADPH min$^{-1}$µmol PRX1$^{-1}$) in contrast to CRX without $Ca^{2+}$ (open circle) or without reductase (grey filled circle). CRX C238S, C241S (black filled square) was not able to reduce PRX1. One exemplary measurement is shown. Assay modified after ref. 66. (f) Titration of electron transfer between PRX1 and CRX in dependence of $Ca^{2+}$. Fitting the data to Michaelis–Menten kinetics revealed a half-maximal rate of NADPH oxidation at a concentration of 122.3 ± 64.5 nM free $Ca^{2+}$. Scale bars give s.d. of three measurements.</figcaption></figure>

*Escherichia coli* and functionally analysed, the protein is not active, neither in the DTNB assay nor in the PRX1 reduction assay (Supplementary Fig. 1a). To validate the importance of the inter-domain hydrogen bond network, residues Gln154, Lys242 and Lys263 were modified by site-directed mutagenesis and altered recombinant CRX proteins were expressed and measured for their competence to donate electrons to PRX1 (Fig. 4c, purple labels indicate the residues of structure based mutagenesis, Supplementary Fig. 1b). CRX mutants Gln154Ala and Lys242Leu were less efficient in electron transfer to PRX1 at low $Ca^{2+}$-concentrations as compared with WT CRX but had similar maximal electron transfer rates. Mutant CRX Lys263Ile on the other hand had an almost threefold lower maximal electron transfer rate at saturating $Ca^{2+}$-concentrations. These results indicate that the communication between $Ca^{2+}$-binding and TRX domains indeed tunes the $Ca^{2+}$-dependent enzymatic activity of CRX.

**CRX depletion leads to increased CEF and ROS production.** To analyse the physiological role of CRX, its expression was diminished using an amiRNA strategy. Moreover, screening of an insertional *Chlamydomonas* mutant library identified a CRX

mutant with an insertion in the second intron of the *crx* gene (Supplementary Fig. 5). In the insertional *crx* mutant (*IMcrx*) only very minor CRX expression is detectable (*IMcrx*, Fig. 5a). Transformants expressing the amiRNA construct (knock down (KD) strains) showed a strong decrease in CRX protein amounts in comparison to the empty vector control (Fig. 5b). Here, *amiRNA-crx-23* and *amiRNA-crx-12* displayed a more than fourfold reduction of CRX expression. Notably, expression of CRX in the WT and the vector control strain is more than fourfold induced in high light (HL) and photoautotrophic conditions as compared with low light (LL) and photo-heterotrophic conditions (Fig. 5a,b). This is in line with previous quantitative proteomics results, where CRX expression was found to be induced under autotrophic versus photoheterotrophic conditions³¹ and further confirmed by quantitative mass spectrometric data in this study (Fig. 6a).

To test whether the light-stress response is active in CRX depletion strains, the expression of LHCSR3 was assayed by immunoblotting (Fig. 5a,b). In the *IMcrx* as well as in the *amiRNA-crx* strains, expression of LHCSR3 was induced on shift from LL to HL. After 24 h HL, LHCSR3 is slightly diminished in the *IMcrx* as compared with the control, while LHCSR3 expression in the KD strains and the vector control
<table><caption>**Table 1 | Data collection, phasing and refinement statistics for single-wavelength anomalous dispersion (SeMet) structures**</caption><thead><tr><th></th><th>SeMet</th><th>Calredoxin</th></tr></thead><tbody><tr><td colspan="3"><strong>Data collection</strong></td></tr><tr><td>Space group</td><td>P2<sub>1</sub></td><td>P2<sub>1</sub></td></tr><tr><td>Cell dimensions</td><td></td><td></td></tr><tr><td>a, b, c (Å)</td><td>85.6, 54.8, 89.1</td><td>86.0, 55.2, 89.5</td></tr><tr><td>α, β, γ (°)</td><td>90.0, 101.1, 90.0</td><td>90.0, 101.1, 90.0</td></tr><tr><td>Wavelength (Å)</td><td>0.97500</td><td>0.90000</td></tr><tr><td>Resolution (Å)</td><td>50-2.8 (2.85-2.80)</td><td>50-1.6 (1.63-1.60)</td></tr><tr><td>R<sub>merge</sub> (%)</td><td>9.5 (19)</td><td>6 (55.8)</td></tr><tr><td>l/σ l</td><td>60.2 (20.2)</td><td>34.9 (3.1)</td></tr><tr><td>Completeness (%)</td><td>99.6 (100)</td><td>99.7 (100)</td></tr><tr><td>Redundancy</td><td>7.5 (7.6)</td><td>3.7 (3.6)</td></tr><tr><td colspan="3"><strong>Refinement</strong></td></tr><tr><td>Resolution (Å)</td><td></td><td>50-1.6</td></tr><tr><td>No. reflections</td><td></td><td>103,272</td></tr><tr><td>R<sub>work</sub>/R<sub>free</sub> (%)</td><td></td><td>19.88/23.25</td></tr><tr><td>No. of atoms</td><td></td><td></td></tr><tr><td>Protein</td><td></td><td>4,862</td></tr><tr><td>Ca<sup>2+</sup></td><td></td><td>8</td></tr><tr><td>Water</td><td></td><td>849</td></tr><tr><td colspan="3"><strong>B-factors</strong></td></tr><tr><td>Protein</td><td></td><td>25.038</td></tr><tr><td>Ca<sup>2+</sup></td><td></td><td>24.761</td></tr><tr><td>Water</td><td></td><td>35.824</td></tr><tr><td>r.m.s. deviations</td><td></td><td></td></tr><tr><td>Bond lengths (Å)</td><td></td><td>0.0253</td></tr><tr><td>Bond angles (°)</td><td></td><td>2.2260</td></tr></tbody></table>

Highest resolution shell is shown in parentheses.

strain appears to be similar. Thus the induction of LHCSR3 expression due to an increase in light intensity is operating when CRX amounts are diminished in expression. Measurements of linear photosynthetic electron flow revealed no significant differences between control, IMcrx and KD strains (Fig. 5c,d). However, measurements of CEF clearly showed a significant increase in CEF capacity in the IMcrx and the KD strains in comparison to the respective controls (Fig. 5c,d).

Since PRX1 was identified as a potent interaction partner of CRX1, we investigated whether IMcrx and *amiRNA-crx-23* were affected in ROS defense. Photo-oxidative light stress is known to induce lipid peroxidation in plants<sup>40</sup>. To study lipid peroxidation we employed the thiobarbituric acid assay<sup>41</sup>. Notably, we found that lipid peroxidation is enhanced in both CRX deficient strains (Fig. 5e,f), underpinning a function of CRX in ROS defense and indicating that the functional protein-protein interaction between CRX and PRX1 in detoxification of H<sub>2</sub>O<sub>2</sub> found *in vitro* (Fig. 3e,f) is of physiological importance *in vivo* (Fig. 5e,f).

**Depletion of CRX diminishes HL induction of TRX f.** To this end WT and IMcrx cells were grown either under photoheterotrophic (tris-acetate-phosphate (TAP)) LL conditions (30 µE m<sup>-2</sup> s<sup>-1</sup>) and shifted to autotrophic (high salt medium (HSM)) HL growth conditions for 6 h (HL, 180 µE m<sup>-2</sup> s<sup>-1</sup>) in <sup>15</sup>N or <sup>14</sup>N labelled media. Mixed whole-cell samples were digested with trypsin using FASP<sup>42</sup> technology and peptides were analysed by liquid chromatography-mass spectrometry (LC-MS)/MS. We processed samples from eight conditions including swapping experiments (Fig. 6). The mass spectrometric analyses resulted in the identification of 42,815 distinct peptides (pep ≤ 0.05) and 8,555 proteins permitting quantitation

and calculation of 2,251 proteins with ratios. No CRX peptides could be identified in IMcrx underpinning the strong depletion of CRX in IMcrx (Fig. 6a). These data also revealed that CRX was twofold upregulated after 6 h HL HSM in the WT, in accordance to the immunoblot results (Fig. 5a,b). Quantitative analyses of photosynthetic proteins displayed only slight differences between WT and IMcrx under HL HSM. Herein core proteins of PSI and PSII were slightly diminished in IMcrx, an impact that was not observed under LL TAP growth conditions (Supplementary Fig. 6). Co-regulation analyses of the quantitative data were performed by using pyGCluster<sup>43</sup> including 1,289 proteins present in all clustered conditions. After 250,000 iterations, 146 clusters were found (for a threshold for 0.001) at a minimal cluster length of 4 (Supplementary Fig. 7). These data indicated that overall quantitative responses found via clustering are very similar between the WT and IMcrx. For example, HL responses including proteins involved in carbon concentrating mechanism (cluster 113, Supplementary Fig. 7) were likewise induced. A striking difference in the HL response was observed for TRX f (Cre01.g066552.t1.1), a protein not seen in the clusters. TRX f is a chloroplast TRX with crucial functions in the activation of the Calvin-Benson-Bassham (CBB) cycle enzymes in C. reinhardtii<sup>44</sup>. As for CRX, TRX f was found to be twofold induced under HL in the WT (Fig. 6b). However, the protein was not increased in abundance in IMcrx under HL HSM in comparison to LL TAP control, an expression difference was also seen in its downregulation under IMcrx HSM HL versus WT HSM HL (Fig. 6b). This behaviour is very particular for TRX f, as it was not seen for other quantified chloroplast TRX proteins (Fig. 6c). Notably, PRX1 was only marginally induced under HL and was not differentially expressed between WT and IMcrx.

## Discussion

In this work we characterized a protein designated as CRX, representing a new class of Ca<sup>2+</sup>-dependent 'sensor-responder' proteins. Our structural data revealed that all EF-hands of CRX have the competence to bind Ca<sup>2+</sup> at high affinity (Figs 1, 3 and 4) and are functionally interconnected with the TRX domain that is enzymatically active (Figs 1 and 3e,f).

A closer look at the interface of the two domains in the CRX structure provided deep insight into the inter-domain communication path. There were two hydrogen bond networks between the active site of the TRX domain and the Ca4 binding site. One network involved Helix 7, a part of the fourth EF-hand motif that was linked to Lys242 located underneath the disulfide bridge (Fig. 4c). The other was through the Helix 12 of the TRX domain. Asn269 located in the end of Helix 12 was directly interacting with the calcium ligand, Asp157 and its side chain represents an integral part of a tight hydrogen bond network stretching along Gln154, Leu235, Thr236, Asn267, Ala268 and Asn269 to the amide nitrogen of Cys238. This latter network included van der Waals interactions and hydrogen bonds with and without water molecules. These two inter-domain pathways are likely responsible for the Ca<sup>2+</sup>-dependent TRX activity (Fig. 1b) and the diminished activity towards PRX1 in the absence of Ca<sup>2+</sup> (Fig. 3e,f), which is further underpinned by the fact that CRX mutant versions harbouring the alterations Gln154Ala and Lys242Leu were less efficient in electron transfer to PRX1 at low Ca<sup>2+</sup>-concentrations (Supplementary Fig. 1b). The difference in Ca<sup>2+</sup>-dependent TRX activity could also be interpreted in a way that the active refolding of CRX and stimulation of TRX activity and/or binding of NTR required Ca<sup>2+</sup>. However, it appears that the TRX activity of CRX needs Ca<sup>2+</sup>, as the TRX domain without the CaM domain is functionally inactive (Supplementary Fig. 1a). This is an intriguing aspect as the structure of the TRX
<figure><img src="image_9.png" /><figcaption>**Figure 4** | X-ray structure of calredoxin. (a) Overall structure of CRX. Sequence diagram of CRX is shown at the top. The N- and C-subdomains of the CaM domain are displayed in orange and magenta, respectively. Bound Ca²⁺ ions are represented as green spheres. The TRX domain is highlighted in marine-blue with the disulfide bridge shown as yellow ball-and-stick model. (b) Open-book representation of interactions between the CaM and the TRX domains. Residues involved in direct inter-domain interactions except for water-mediated hydrogen bonds are shown. (c) Inter-domain networks between the Ca4 calcium ion in the CaM domain and the disulfide bridge in the TRX domain. Purple labels indicate the residues of structure based mutagenesis. (d) The stereo image of 2|F<sub>o</sub> | -|F<sub>c</sub>| electron density maps (2σ level) showing amino acid residues around Ca3 (upper) and the disulfide bridge between Cys238 and Cys241 (lower).</figcaption></figure>

domain showed noteworthy similarity to the conventional TRX molecule (r.m.s.d = 1.36 Å for Caα atoms of spinach TRX f), strongly suggesting that the Ca²⁺-dependent activity of CRX is modulated via the inter-domain communication path which is absent in TRX f.

Assuming the TRX's interaction sites are conserved in the TRX domain of CRX, the structures of CRX together with the NTR or FTR were predicted based on the published X-ray structures of their electron transfer complex. In the predicted structures (Supplementary Fig. 8a,b), no overlap was found in the CRX:NTR complex;
<figure><img src="image_10.png" /><figcaption>b</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>c</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>Electron flow with DCMU/HA</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>Contribution of CEF</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>d</figcaption></figure>

<figure><img src="image_15.png" /><figcaption>Electron flow with DCMU/HA</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>Contribution of CEF</figcaption></figure>

<figure><img src="image_17.png" /><figcaption>e</figcaption></figure>

<figure><img src="image_18.png" /><figcaption>f</figcaption></figure>

<figure><img src="image_19.png" /></figure>

**Figure 5** | Calredoxin function in vivo. (a,b) Immunoblot analysis of WT versus IMcrx (a) or empty vector (EV) versus amiRNA-crx-12/23 (KD#12/KD#23) (b) whole-cell extracts. ATPB, LHCSR3 and CRX protein expression were examined in cells grown under photoheterotrophic (TAP) LL (30 µE m⁻²s⁻¹), which were shifted to photoautotrophic (HSM) HL (180 µE m⁻²s⁻¹) growth conditions. Chlorophyll (1.5 µg) were loaded per lane and equate 100%. ATPB was used as loading control. (c,d) Comparison of linear photosynthetic electron flow and CEF in WT versus IMcrx (c) or EV versus KD#12 and KD#23 (d) cells grown under TAP/LL, which were shifted for 6 h either to HSM/LL or HSM/HL growth conditions. LL and HL data (± s.d.) refer to analyses of three biological replicates for WT and IMcrx, with 5 and 11 as well as 6 and 12 measurements, respectively; KD LL Data (± s.d.) refer to analyses of three biological replicates for EV, KD#12 and KD#23 with three measurements each; HL Data (± s.d.) refer to analyses of three (EV), three (KD#12) and seven (KD#23) biological replicates, with three, nine and 10 measurements, respectively. For statistical analysis of indicated data t-test with *P*<0.05, **P**<0.01 and ***P**<0.001 was performed. (e,f) Measurement of lipid peroxidation in calredoxin-deficient strains under HL conditions. WT versus IMcrx (e) and of empty vector (EV) versus KD#23 (f) cells after transition from TAP/LL to HSM/HL growth conditions. Malondialdehyde equivalents were measured from whole cells (4 µg Chl ml⁻¹) before (0 h, grown in TAP at 30 µE m⁻²s⁻¹) and after HL treatment (shifted to HSM for 6, 24, 30 h at 180 µE m⁻²s⁻¹). WT versus IMcrx data refer to four biological replicates with eight measurements each, EV versus KD#23 data refer to analyses of three biological replicates with six measurements each. Error bars represent s.d.'s. Statistical comparison of indicated data was done using t-test with *P*<0.05, **P**<0.01 and ***P**<0.001.

however, in the CRX:FTR, a spatial overlap was found between the FTR and the CaM domain. Based on these in silico modelling, it is suggested that CRX is reduced through an NTR-dependent pathway

in chloroplasts and oxidized by a 2-Cys PRX from a structural point of view. The effective reduction of CRX via the NTR from E. coli is consistent with this interpretation (Figs 1 and 3e).
<figure><img src="image_20.png" /><figcaption>**Figure 6** | TRX *f* is diminished in HL on depletion of CRX. MS-based ¹⁵N metabolic labelling-based quantitation from CRX (Cre03.g202950.t1.1) (a) and TRX *f* (Cre01.g066552.t1.1) (b). Each pie slice represents a quantified peptide. The area of a pie slice is proportional to the log2 ratio of the corresponding peptide between the indicated conditions. Colours represent pyQms quantification score; 0.7 (yellow, false discovery rate (FDR) ≤1%) to 1 (blue, perfect match). (c) Selected knowledge based TRX community. The heat map represents the ratios of the proteins ranges from yellow to green indicating an upregulation and from yellow to red indicating a downregulation. The s.d. is visualized by the size of the box, smaller the box higher the s.d. of the protein ratio (see legend). For quantitative MS analysis, ¹⁴N-/¹⁵N-labelled WT and ¹⁵N-/¹⁴N-labelled IMcrx were mixed based on equal protein amount (total protein amount 100 µg). Eight different conditions were examined as followed: (1) ¹⁴N-labelled WT TAP LL versus ¹⁵N-labelled IM TAP LL, (2) ¹⁵N-labelled WT TAP LL versus ¹⁴N-labelled IM TAP LL, (3) ¹⁴N-labelled WT HSM HL versus ¹⁵N-labelled IM HSM HL, (4) ¹⁵N-labelled WT HSM HL versus ¹⁴N-labelled IM HSM HL, (5) ¹⁴N-labelled WT TAP LL versus ¹⁵N-labelled WT HSM HL, (6) ¹⁴N-labelled IM TAP LL versus ¹⁵N-labelled IM HSM HL, (7) ¹⁴N-labelled WT TAP LL versus ¹⁵N-labelled IM HSM HL and (8) ¹⁴N-labelled IM TAP LL versus ¹⁵N-labelled WT HSM HL.</figcaption></figure>

For the structural analysis of the CRX:PRX complex, we referred to the structure of the TRX domain of ERp46 complexed with the C-terminal peptide of mammalian 2-Cys PRX, PRX4, (TRX:PRX4-peptide)⁴⁵. A mammalian PRX4 existed as a ring-shaped decameric form⁴⁶, while the plant-type PRX forms a noncovalent homodimer⁴⁷. Although it is not clear yet which oligomeric state *Chlamydomonas* PRX1 has, the active site structure is conserved between these PRX structures. From the partial complex structure of TRX:PRX4-peptide, we can build the predicted CRX:PRX-peptide structure (Supplementary Fig. 8c).

The TRX domain of CRX has an open space for PRX binding. Additionally, the decameric PRX4 (PDB ID: 3VWU) can be modelled by fitting its C-terminal peptide to the corresponding peptide in the CRX:PRX-peptide structure (Supplementary Fig. 8d). Further the modelled structure shows a close fit to the open space formed by CRX that implies an additional interaction between CRX and PRX mediated by the CaM-like domain on binding of Ca²⁺. The structural arrangement of the four EF-hands within the Ca²⁺-binding domain of CRX is different with regard to the conventional CaM structure. Referring to the
dumbbell-shaped CaM structure, two α-helices (H4 and H5) connecting the second and the third EF-hand motifs correspond to the central helix that would bend when recognizing the target protein. The position of the target protein in the conventional structure of the CaM:peptide complex is occluded by the additional structured loop between H2 and H3 helices, which suggested no similarity in target recognition by CRX. Clearly, target recognition of CRX in dependence of $Ca^{2+}$ via the four EF-hands needs to be further investigated. However, the second protein-protein interaction module of CRX is its TRX domain, as revealed by the capture of PRX1 in the thiol-trapping experiments (Fig. 3a–d, Supplementary Fig. 4a–c). Remarkably, the thiol-trapping of PRX1 as well as CRX electron transfer towards PRX1 were $Ca^{2+}$-dependent (Fig. 3a–f). Thus, as mentioned above, one important function of the $Ca^{2+}$-dependent structural changes of the four EF-hands is manifested in the modulation of TRX binding capability and activity via the inter-domain communication path.

Thereby, it would translate alterations in chloroplast $Ca^{2+}$ concentration and redox potential in modulation of CRX function, placing the protein at the crossroad of redox and $Ca^{2+}$ signalling in the chloroplast. Notably, CRX is the first CML protein identified in the chloroplast. Its presence and function strongly supports the notion that $Ca^{2+}$ plays an important role in regulating algal photosynthesis but also underpins the view of a general functional interconnection of $Ca^{2+}$ and redox signalling in chloroplast photosynthesis.

Consistent with this hypothesis, our *in vivo* analysis reveals that CEF is induced on reducing or removing CRX activity. This effect on CEF can be rationalized based on the more reducing stromal redox poise and/or enhanced ROS (for example, hydrogen peroxide production (see above)) when CRX activity is diminished. The strong depletion of CRX enhanced lipid peroxidation (Fig. 5e,f), in line with an increased oxidative stress. As described above, CRX is an efficient alternative to the NTRC system<sup>28</sup> to reduce hydrogen peroxide via PRX1 (Fig. 3e,f). Therefore, a strong diminishment of CRX should results in a less efficient reduction of oxidized PRX1, thereby increasing in hydrogen peroxide in the light. Importantly, while CRX as well as PRX1 expression is enhanced at increasing light intensities (Fig. 5a, Fig. 6a and ref. 39), the induction of CRX in HL is much stronger than the increase of PRX1 (Fig. 6c). Besides NTRC as shown in Arabidopsis<sup>28</sup>, TRXx is another chloroplast TRX<sup>48</sup>, is involved in electron transfer to PRX1 in *Chlamydomonas*<sup>49</sup>. TRXx expression is comparable between WT and IMcrx, while NTRC expression in WT is not induced in HL but slightly diminished under HL in IMcrx (Fig. 6c), overall suggesting that in HL, CRX is a significant electron donor to PRX1, probably crucial in minimizing photo-oxidative damage.

The increase in chloroplast $Ca^{2+}$ in the light could be important for activation of the TRX domain of CRX (Figs 1b and 3e). Light and $Ca^{2+}$-dependent activation has also been described for other chloroplast enzymes, for example, some of the CBB cycle (reviewed in ref. 7). In this scenario CRX would provide reducing power for oxidized PRX1 in the light. It is striking that on depletion of CRX the expression of TRX *f* is not induced in HL (Fig. 6b,c). As mentioned, TRX *f* is critical for the activation of CBB cycle enzymes<sup>44</sup> and is the only chloroplast TRX that is glutathionylated<sup>50</sup>, impairing its reduction by FTR leading to compromised light activation of target enzymes<sup>50</sup>. Regulation of TRX *f* via CRX is therefore another and new layer of regulation, which needs to be further explored. Here the diminishment of NTRC in IMcrx is also noteworthy, indicating that the depletion of CRX has a particular impact on the chloroplast TRX system.

Notably, the $Ca^{2+}$-dependent reduction of PRX1 via CRX (Fig. 3f) was found at a $K_d$ value of ~122 nM free $Ca^{2+}$, a value

that corresponds to the resting stromal $Ca^{2+}$-concentration<sup>8</sup>. Thus an increase of stromal $Ca^{2+}$ could indeed further activate CRX and accelerate the rate of PRX1 reduction or reduction of other substrates. Therefore, enhanced ROS in the light due to depletion of CRX is a possible explanation for the induction of CEF in *C. reinhardtii* as described for Arabidopsis<sup>22</sup>. In Arabidopsis, however, hydrogen peroxide induced CEF via the NDH-dependent pathway and not via the PGR5-dependent pathway. It is currently not known whether the NDA2-dependent pathway in *C. reinhardtii*<sup>51</sup> is activated by hydrogen peroxide. In the redox regulation of the PGR5 and PGRL1 cycle, the PGRL1 homodimer to monomer conversion requires reduction of a disulfide bond via the TRX system<sup>19</sup>. Although CRX is depleted and TRX *f* is not induced in HL in IMcrx, CEF is still activated and higher than WT. CEF activation could be accomplished by other chloroplast TRXs such as TRXx, which is not affected in IMcrx (Fig. 6c). On the other hand the diminishment of TRX *f* in HL may impair the activation of CCB. This in turn would feedback on the photosynthetic electron transfer, as stromal ATP and NADPH will increase, and activate CEF (refs 52,53), in line with our findings (Fig. 5c,d). In conclusion, CRX is involved in ROS defense and in controlling TRX *f* expression, thereby modulating the stromal redox poise and concomitantly CEF.

Moreover, the CRX inter-domain communication path could translate alterations in chloroplast $Ca^{2+}$ and redox state into cellular signalling.

## Methods

**Culture conditions.** All *C. reinhardtii* strains were grown photoheterotrophically in TAP medium<sup>54</sup> at 25 °C, 120 r.p.m. shaking and 30 µE m<sup>-2</sup> s<sup>-1</sup>. For HL growth experiments, cells were shifted to photoautotrophic growth conditions in HSM, washed once with HSM, resuspended to a chlorophyll concentration of 4 µg ml<sup>-1</sup> and exposed to HL (180 µE m<sup>-2</sup> s<sup>-1</sup>) for indicated hours.

**Protein expression and purification.** For *in vitro* experiments *E. coli* strain BL21 (DE3; Novagen) was transformed with pET-22b(+) plasmids encoding different versions of CRX (WT; C238S; C241S; C238S, C241S; Q154A; K242L; K263I or the TRX domain only). For protein crystallization, CRX WT and the SeMet derivative were expressed from the same vector but in *E. coli* Rosetta and B834 (DE3) pLysS strain, respectively. Cells were grown at 37 °C until an OD<sub>600</sub> of 0.6 was reached. Expression was induced with 0.5 mM isopropyl-β-D-thiogalactoside at 37 °C for 5 h (*in vitro* experiments) or with 0.1 mM Isopropyl-β-D-thiogalactopyranoside at 20 °C for 18–24 h (crystallization experiments). The His-tagged proteins were purified either natively from *E. coli* lysates (for crystallization, microscale thermophoresis (MST) and PRX1 interaction experiments) or under denaturing conditions from inclusion bodies (DTNB activity assay (except Supplementary Fig. 1), redox titration, thiol-trapping experiment) using $Ni^{2+}$-affinity resin. Eluates were dialyzed against 30 mM MOPS, 100 mM KCl, pH 7.2 and sterile filtered. For protein crystallization after $Ni^{2+}$-affinity purification the protein sample was submitted to ion exchange (Hi-Trap Q) chromatography. The eluates were pooled and further purified by gel filtration (Superdex-200).

**Crystallization and structure solution.** WT CRX was crystallized at 20 °C using the hanging-drop vapour diffusion method. Protein concentration was adjusted to 10 mg ml<sup>-1</sup>. Sample buffer contained 20 mM Tris-HCl pH 8.0, 150 mM NaCl, 1 mM DTT and 150 mM $CaCl_2$. Precipitant solution was 1.5 M $LiSO_4$ dissolved in 0.1 M MES-NaOH buffer pH 6. To obtain the crystals, the drops contained 1 µl of the protein and 1 µl of precipitant solution. SeMet derivative crystals were grown with the same experimental setting. X-ray diffraction images were collected at 100 K using a synchrotron radiation source with 0.90000 Å for native and 0.97500 Å for SeMet. All data sets were collected at the SPring-8 using the BL44XU beamline.

The native and SeMet derivative data sets were processed using the HKL-2000 software. Phase problem was solved by the single-wavelength anomalous dispersion method using Phenix AutoSol in Phenix program suite at 2.8 Å resolution. The SeMet derivative protein had 32 SeMet residues and 27 sites were found. The density modified resultant electron density map was clear enough for interpretation and the asymmetric unit contained two molecules of CRX. The model was build using COOT software along with the NCS-averaged electron density map calculated using DM in the CCP4 package. Crystallographic refinement was performed using RefMAC5. The final structure was validated using wwPDB validation server. The analysis showed that all residues were in the favoured region of Ramachandran plot, with no residues as outliers. Crystallographic data and refinement statistics are shown in Table 1.
**Vector constructions.** RNA isolation from *Chlamydomonas reinhardtii* was performed using TRI reagent solution (Sigma Aldrich) and 300 ng of DNase I treated RNA was used for cDNA synthesis (iScript cDNA Synthesis Kit, Biorad) according to the manufacturer's protocol.

For *in vitro* experiments the *crx* and *pxr1* genes were amplified from cDNA lacking the transit peptide and the stop codon using the primers 5'-ggaattccatatgtcgagcgtcctcc-3' and 5'-cgggatcctgagcggccgcagc-3' for CRX and 5'-ggtattccatatgcttcccacgcag-3' and 5'-cgggatcctgacggcagcagagaagtactcc-3' for PRX1. Each amplicon was digested with NdeI and BamHI (NEB) and cloned into the pET-22b(+) vector (Novagen). Based on the CRX plasmid, site-directed mutagenesis (QuickChange II Site-Directed Mutagenesis Kit, Agilent) was performed (C238S; C241S;C238S, C241S; Q154A; K242L and K263I) according to the manufacturer's protocol. The following primers were used for the mutagenesis reaction: C238S: 5'-gcgcgctcacatggagccggccctgcaaggg-3' and 5'-cccttgcagggccggcctcatgtgacgcgcagc-3'; C241S: 5'-gtgctccggcccccagcaagggcat-3' and 5'-catgcccttgctgggccggcacc-3'; C238S, C241S: 5'-gctcacatggagccggcccccagcaagggcat-3' and 5'-atgccttctgctgggccggctcatgatgagc-3'; Q154A: 5'-gtcgtccttgttcgtacatggccatgaggtaacaacagttc-3' and 5'-gaagctgttgacctcatggccatgtacgaacaggacgac-3'; K242L: 5'-ccgcgcatccggccaggcagggccgcag-3' and 5'-tggccggcctgcctgggcatcagcgg-3'; K263L: 5'-gttggcgttaccgaacagaatcacaacacgatgtggtc-3' and 5'-gaccacatcgttgttgtgattctgttcgtaacgccaaac-3'. Moreover the same plasmid was used as template for the amplification of the TRX domain of CRX with the primers 5'-ggataccatatgctgagcagggcgg-3' and 5'-ccgctcagcaggttggcatcag-3'. The amplicon was digested with NdeI and XhoI (NEB) and cloned into an empty pET-22b(+) vector (Novagen).

To visualize calredoxin expression in *C. reinhardtii* *in vivo* a calredoxin-YFP construct was generated. For this purpose the *crx* gene without the stop codon was amplified using the primers 5'-ctaggctagcatgattgcaattgcac-3' and 5'-ctagtccggacagggccgcagc-3'. Then the amplicon was sequentially digested with NheI and BspI (NEB) and cloned into a pChlamiRNA3-based vector, between the *bleomycin-FMDV* 2A sequence³⁵ and the YFP coding region⁵⁵. Nuclear transformation of *C. reinhardtii* strain CC-851 cw2 mt + was performed using the glass bead protocol (1-2 µg of plasmid DNA per transformation, linearized with *Kpn*I and *Not*I)⁵⁶. Transformants were selected on TAP plates supplemented with zeocin (10 µg ml⁻¹).

Calredoxin-knockdown strains were generated using an amiRNA approach⁵⁷. Calredoxin-specific oligonucleotides were designed using the WMD3 tool (http://wmd3.weig爾world.org/cgi-bin/webapp.cgi). The oligonucleotides 5'-ctagtT CGGATCAGTTTTGTTGTTTCATctgcctgatcggcaccatgggggtggtggtgatcagcgcctaTGAA GAACAAAACTGATCCGAg-3' and 5'-ctagTCGGATCAGTTTTGTTCTTCAtag cgcctgatcaccaccaccccctggtccgatcagcaggaTGAACAACAAAACTGATCCGAA-3' that target the 3'UTR of the *crx* gene were annealed and ligated via the *SpeI* restriction site into the pChlamiRNA2 vector. The resulting plasmid was transformed into the arginine auxotrophic cw15-325 (cwd mt + arg7 nit1 + nit2 +) strain by electroporation according to ref. 58, with some modifications. Plasmid DNA (150-300 ng) was added to 250 µl of a concentrated cell suspension (2 × 10⁸ cells ml⁻¹ in TAP supplemented with 40 mM sucrose). The cell/DNA mix was incubated in an electroporation cuvette for 10 min on ice before electroporation with pulse settings of 1 kV and 25 µF (GENE PULSER II coupled to the PULSE CONTROLLER II, Biorad). Cells were transferred into a 50-ml falcon tube containing 10 ml TAP supplemented with 40 mM sucrose and shaken overnight under LL (10 µE m⁻² s⁻¹). Cells were collected by centrifugation at 3,000g for 5 min and resuspended in 1 ml TAP. Resuspended cells (330 µl) were plated onto 1% (w/v) TAP agar plates.

For the transient expression of the calredoxin-mCherry fusion protein in *N. benthamiana* epidermal leaf cells the *crx* gene without the stop codon was amplified using the primers 5'-ggggtaccatgattgcaattgcactg-3' and 5'-catctcagcagggccgcagc-3'. The amplicon was digested with *Kpn*I and *Xho*I (NEB) and cloned into the binary vector pGPTKII.bar_mcherry (kindly provided by Jörg Kudla).

**Insertional mutagenesis and IMcrx identification.** The IMcrx mutant was identified from an insertion mutant library using a PCR-based screening method⁵⁹. The library was generated using the AphVIII fragment from plasmid pMj013b as the insert, and a rescued mutant of ift46-1 (CC-4375) as the parental strain⁶⁰,⁶¹. The flagella less mutant ift46-1 was rescued by transformation of a ~7.2 kb fragment including the full-length genomic sequence of IFT46 fused to the YFT gene at 3'-end, and the Aph7' gene as the selective marker. The whole library was separated into ~100 super-libraries each containing 1,440 transformants. PCR reactions using super-library DNA as template were performed with a gene specific primer and an insert specific primer LGR06-F (primers are listed in Supplementary Table 1). The super-library that included the IMcrx mutant was determined according to the sequencing data of a positive PCR product (Supplementary Fig. S5. The positive super-library was further separated into 15 small-libraries and each small-library DNA was checked individually by PCR with appropriate *crx* target primer. The single clone of the IMcrx mutant was picked from the small library after identification by colony PCR. Furthermore, insertion of the AphVIII fragment was confirmed with a primer pair binding adjacent to the insertion site (EX2-F, EX3-R, Supplementary Table 1, Supplementary Fig. 5).

**Microscale thermophoresis binding analyses.** Recombinant CRX (100 µl of 40 µM) was labelled with red fluorescent amine-reactive dye according to the

instructions of the MO-L001 Monolith Protein Labeling Kit RED-NHS (amine reactive). A constant concentration of labelled CRX (~115 nM) was mixed with defined concentrations of free Ca²⁺ or Mg²⁺ ions (up to 100 µM) in 30 mM MOPS (pH 7.2), 100 mM KCl, 0.05% Tween-20, 5 mM total EGTA or EDTA respectively. To prepare solutions with Ca²⁺ in the micromolar range we took advantage of available calcium buffer kits guidelines (ThermoScientific). Measurements were performed at 22 °C in standard treated glass capillaries (Monolith NT Capillaries, NanoTemper) and the thermophoresis analysis was performed on a NanoTemper Monolith NT.115 instrument (20% LED; 20% infrared-laser power). The MST data were fitted with the law of mass action using the NanoTemper Analysis software to obtain K_d values for binding between CRX and Ca²⁺.

**Measurement of Ca²⁺-dependent redox activity.** The redox activity of CRX was determined photometrically in a TRXR-dependent assay⁶² (modified). *E. coli* (40 nM) TRXR were used to specifically reduce CRX (10 µM) by 10 min incubation in 30 mM MOPS, 100 mM KCl, pH 7.2 with NADPH (200 µM) at RT. After addition of 200 µM DTNB, the formation of TNB⁻ was followed at 412 nm (Ultsropec 3000, Amersham Biosciences). The reduction rate of DNTB was determined at defined free Ca²⁺ and Mg²⁺ concentrations and over a time course of 0–80 s after addition of DTNB. In Fig. 1b the different activities were normalized on the highest activity measured for each batch of purified CRX, plotted against the free ion concentration and modelled according to Michaelis–Menten kinetics with GraphPad Prism 2.01 software.

**Determination of redox potential.** Oxidation-reduction titration was performed as follows. Briefly, recombinant WT CRX (150 µg ml⁻¹) was incubated for 2 h at room temperature in 100 mM MOPS, 2.5 µM CaCl₂, pH 7.0 that contained a total DTT concentration of 2 mM. Different redox potentials (*E*<sub>h</sub>) were adjusted by mixing appropriate quantities of oxidized and reduced DTT. The CRX disulfide/dithiol redox state at each *E*<sub>h</sub> value was monitored using the monobromobimane fluorescence method: 250 mM monobromobimane (final concentration) was allowed to react with reduced CRX thiols for 20 min at RT. To remove excess dye, proteins were precipitated with 10% TCA in acetone and washed with 1% TCA in acetone. Pellets were resuspended in 100 mM Tris/HCl pH 8, 1% SDS and fluorescence of a 3:50 dilution in Tris/HCl, pH 8 was detected in a Jasco Fluorometer at 470 nm (excitation: 380 nm). Data were fitted with a two-electron Nernst curve using GraphPad Prism 2.01 software to calculate the CRX *E*<sub>m</sub> value.

**Transient expression in Nicotiana benthamiana.** *Agrobacterium tumefaciens* (GV601 strain), containing the calredoxin-mCherry fusion construct (OD<sub>600</sub> 0.5), were coinfilttrated with the p19 strain (OD<sub>600</sub> 0.3) into 6-week-old *N. benthamiana* leaves as described previously⁶³. Microscopic analyses were conducted 3 days after infiltration on leaf discs (Ø 1.2 cm) of the lower epidermis at room temperature with water as imaging medium.

**Microscopy.** Protein fluorescence was detected by confocal laser-scanning microscopy using a Leica TCS-SP5 II setup consisting of an inverted confocal laser-scanning microscope DMI6000 (Leica) equipped with an × 63/1.2 water immersion lens (HCX PL APO lambda blue 63.0 × 1.20 Water UV). The following filters were used: mVenus, excitation 514 nm (Argon laser), scanning 525–550 nm; mCherry, excitation 561 nm (DPSS Laser), scanning 605–638 nm. Chlorophyll auto-fluorescence was detected at 690–740 nm (Supplementary Table 2). Image acquisition was performed using the Leica software (Leica Application Suite – Advanced Fluorescence 2.6.0.7266; Leica Microsystems). Images are shown in RGB mode and brightness and contrast were adjusted using Adobe Photoshop CS3 software.

**Calredoxin affinity chromatography.** A thiol-trapping approach was performed based on refs 37,64. Five milligram of either WT CRX or two single-point mutated versions (C238S and C241S) were covalently linked to HCl-washed CNBr-activated sepharose 4B overnight at 4 °C to produce three separate columns. After successful coupling, remaining active groups were blocked with Tris/HCl, pH 8. Columns were equilibrated overnight in binding buffer (50 mM MOPS, 50 mM NaCl and pH 8) either with 2.5 µM free Ca²⁺ in 5 mM EGTA or with EGTA only. Five milligram of total protein from *C. reinhardtii* were loaded onto each column and incubated for 1 h at 4 °C. After washing with 3 × 5 ml 50 mM MOPS, 500 mM NaCl, pH 8.3, potential CRX targets were eluted with 10 mM DTT in 1 ml steps until absorbance at 280 nm dropped significantly, as determined by NanoDrop measurement. The eluates were concentrated in Vivaspin 4 Columns (MWCO 10 kDa, Sartorius) and subsequently submitted to tryptic digestion (Filter Aided Sample Preparation)⁴² for subsequent mass spectrometric analysis.

**Peptide identification by mass spectrometry.** Tryptically digested peptides were chromatographically separated using an Ultimate 3000 RSLCNano System (Thermo Scientific). The sample was loaded on a trap column (C18 PepMap 100, 300 µM × 5 mm, 5-µm particle size, 100-Å pore size; Thermo Scientific) and desalted for 5 min using 0.05% (v/v) TFA/2% (v/v) acetonitrile in ultrapure water
at a flow rate of 5 µl min⁻¹. Then the trap column was switched in-line with the separation column (Acclaim PepMap100 C18, 75 µm × 15 cm, 2-µm particle size, 100-Å pore size; Thermo Scientific). The mobile phases for peptide elution consisted of 0.1% (v/v) formic acid in ultrapure water (A) and 80% acetonitrile/0.08% formic acid in ultrapure water (B). Peptides were eluted at a flow rate of 300 nl min⁻¹ and employing the following gradient profile: 2.5–35% B over 90 min, 35–60% B over 20 min, 60–99% B over 5 min and 99% B for 10 min. Afterwards the column was re-equilibrated with 97.5% A for 25 min.

The LC system was coupled via a nanospray source to a Q Exactive Plus mass spectrometer (Thermo Scientific) operating in positive ion mode. MS data were acquired in a data-dependent manner, dynamically choosing the 12 most abundant precursor ions from the survey scans (scan range *m/z* 400–1,600, resolution 70,000, AGC target value 1e6, maximum injection time 50 ms) for fragmentation (MS/MS) by higher-energy C-trap dissociation (27% normalized collision energy). Dynamic exclusion was enabled with an exclusion duration of 60 s. AGC target value for MS/MS was 5e4 at 50 ms maximum injection time. The precursor isolation window was 2.0 *m/z* and the resolution for higher-energy C-trap dissociation spectra was set to 17,500. The 'underfill ratio', specifying the minimum percentage of the target ion value to be reached at the maximum fill time, was 1%. Singly charged ions, ions with charge state 8 and above as well as ions with unassigned charge states were excluded from fragmentation. Internal lock mass calibration on *m/z* 445.12003 was enabled.

**Label-free protein quantification.** Peptide identification and protein quantification was performed using MaxQuant 1.5.1.2 (ref. 65) with default parameters for tryptic digestion and Orbitrap detection (maximum peptide mass: 8,000 Da). Spectra were searched against a database containing protein sequences of the *Chlamydomonas* v5.5 gene models (Joint Genome Institute) merged with mitochondrial and chloroplast protein sequences from NCBI databases BK000554.2 and NC_001638.1. The database was further supplemented with commonly observed contaminants (http://www.thegpm.org/crap/). Cysteine carbamidomethylation was set as a fixed modification, the oxidation of methionine, protein N-terminal acetylation, and phosphorylation of serine, threonine and tyrosine were set as variable modifications. A false discovery rate of 1% was applied to peptide and protein identifications. False discovery rate calculation was performed by MaxQuant using a target-decoy approach. The 'match between runs' setting was activated.

Correlation plots were produced using log2 transformed label-free quantification protein intensities.

**PRX1 interaction assay.** The interaction of CRX and PRX1 was measured *in vitro* in a TRXR coupled photometrical assay⁶⁶. Recombinant CRX (5 µM) was reduced by incubation with 1 µM *E. coli* TRXR and 120 µM NADPH in 30 mM MOPS, 100 mM KCl, pH 7.2 and defined free amounts of Ca²⁺ as well as 40 µM H₂O₂. The NADPH consumption was measured at 340 nm until a steady decrease in absorption was detected. Then the reaction was started by addition of 1 µM recombinant PRX1 and recording of the absorption at 340 nm was continued. The rate of NADPH oxidation was calculated from the first 60 s after addition of PRX1.

**Immunoblot analysis.** Immunoblot analyses were performed using equal amounts of chlorophyll (2 µg, Fig. 5) or the equal amounts of total protein (30 µg, Supplementary Fig. S1)³¹. All primary antibodies were used as follows: CRX 1:3000 (generated by Eurogentec using recombinant CRX for immunization of rabbits), ATPB 1:5,000 (Agrisera), CA 1:1000 (Agrisera), LHCSR3 1:2,000⁶⁷. The signals were detected by enhanced chemical luminescence. An exemplary full-blot containing molecular weight markers is shown in Supplementary Fig. 9.

**Spectroscopic measurements.** P700 absorption and electrochromic shift signal measurements were done using an LED pump-probe JTS-10 spectrophotometer (BioLogic). Single turnover measurements were performed using a dye laser emitting at 640 nm, pumped by the second harmonic of a Minilite II Nd:YAG laser (Continuum). *C. reinhardtii* cells were harvested and resuspended to a chlorophyll concentration of 20 µg ml⁻¹ in 20 mM HEPES (pH 7.2) containing 10% Ficoll (w/v) and incubated 10 min in the dark under constant shaking to avoid anaerobiosis.

Electron flow rates were determined as the product of *k*<sub>ox</sub> [P700<sub>red</sub>] as described in ref. 21. Measurements performed in the absence of inhibitors were utilized to calculate linear electron flow, while measurements in the presence of PSII inhibitors 3-(3,4-dichlorophenyl)-1,1-dimethylurea (40 µM) and hydroxylamine (2 mM) were used to determine rates for CEF around PSI. Continuous illumination was provided at 630 nm with an intensity of 130 µE m⁻² s⁻¹. *k*<sub>ox</sub> was determined as the initial change in the 520–546 nm electrochromic shift signal at the onset of actinic illumination normalized to the electrochromic shift signal recorded 140 µs after application of a saturating, single turnover flash, in the presence of the aforementioned PSII inhibitors. To assess [P700<sub>red</sub>], continuous actinic illumination was applied for 10 s, followed by a 30-ms short pulse of strong actinic light to fully oxidize P700. P700 absorption changes were followed at 705 nm, corrected by subtraction of the signal at 740 nm.

**Lipid peroxidation measurement.** Malondialdehyde equivalents were measured as described⁴¹, with minor modifications. Cells were grown under photoheterotrophic (TAP) LL (30 µE m⁻² s⁻¹) conditions until a density of 15 µg Chl ml⁻¹ was reached and then shifted to photoautotrophic (HSM) HL (180 µE m⁻² s⁻¹) growth conditions at a chlorophyll concentration of 4 µg ml⁻¹. At designated time points, two 15 ml aliquots of each culture were taken and incubated with butylated hydroxytoluene (final concentration 0.01% w/v) to prevent further lipid peroxidation. Cells were pelleted at 4,000g for 10 min at 4 °C and resuspended in 1.8 ml of degassed extraction buffer (3.75% TCA (w/v) 0.05 N HCL) with or without 0.1% (w/v) thiobarbituric acid. Samples were incubated in a boiling water bath for 15 min and cooled down to room temperature on ice. Cell extracts were pelleted at 4,000g for 10 min and the absorbance of the supernatant was determined at 440, 532 and 600 nm. Malondialdehyde equivalents were calculated as described ⁴¹ using a molar extinction coefficient for the TBA-malondialdehyde complex of 155,000 nmol⁻¹ ml⁻¹.

**Isolation of chloroplasts and mitochondria.** Chloroplasts and mitochondria were isolated as described in ref. 48.

Protein concentration was determined using the Pierce Protein assay kit (Thermo Scientific, Perbio Science GmbH). Bovine serum albumin was used as standard for construction of the calibration curve.

**Isotopic labelling.** Isotopic ¹⁴N and ¹⁵N labelling was performed with two biological replicates for the WT and the Calredoxin-deficient IMcrx strain. Cells grown for 6 h under photoheterotrophic (TAP) LL (30 µE m⁻² s⁻¹) conditions were compared with cells grown for 6 h under photoautotrophic (HSM) and HL (180 µE m⁻² s⁻¹) conditions. For quantitative MS analysis, ¹⁴N-¹⁵N-labelled WT and ¹⁵N-¹⁴N-labelled IMcrx were mixed based on equal protein amount (total protein amount 100 µg). Eight different conditions were examined as follows: (1) ¹⁴N-labelled WT TAP LL versus ¹⁵N-labelled IM TAP LL, (2) ¹⁵N-labelled WT TAP LL versus ¹⁴N-labelled IM TAP LL, (3) ¹⁴N-labelled WT HSM HL versus ¹⁵N-labelled IM HSM HL, (4) ¹⁵N-labelled WT HSM HL versus ¹⁴N-labelled IM HSM HL, (5) ¹⁴N-labelled WT TAP LL versus ¹⁵N-labelled WT HSM HL, (6) ¹⁴N-labelled IM TAP LL versus ¹⁵N-labelled IM HSM HL, (7) ¹⁴N-labelled WT TAP LL versus ¹⁵N-labelled IM HSM HL and (8) ¹⁴N-labelled IM TAP LL versus ¹⁵N-labelled WT HSM HL. Proteins were tryptically digested overnight on centrifugal ultrafiltration devices (cutoff 30 kDa) according to the FASP method⁴² with an enzyme-to-protein ratio of 1:100.

**LC-MS/MS analyses of WT IMcrx strains.** Liquid chromatography and MS parameters were essentially the same as described above with the following modifications:

Peptide chromatography was carried out on an Acclaim PepMap100 C18 capillary column (75 µm × 50 cm, 2-µm particle size, 100-Å pore size; Thermo Scientific). The following gradient was applied: 2.5–35% B over 195 min, 35–60% B over 40 min, 60–99% B over 5 min, 99% B for 15 min.

MS data acquisition was performed using a data-dependent Top20 method (MS1: scan range *m/z* 375–1600, AGC target value 1e6, maximum injection time 30 ms. MS2 parameters: see above). The precursor isolation window was 1.5 *m/z*. Full scan resolution was set to 70,000 (FWHM). Some samples were additionally analysed with an MS1 resolution of 140,000 (FWHM).

**Identification of peptide spectrum matches and statistical evaluation.** MS2

Spectra were analysed using Ursgal⁶⁸. Post processing was performed using percolator (v. 2.08)⁶⁹ and PSMs were filtered by 5% posterior error probability (PEP). Complete analysis pipeline was executed using the Python framework Ursgal (v 0.3.3)⁶⁸ using default profile 'QExactive +'. Additionally, Methionine oxidation and N-terminal acetylation was included as post translational modification search parameter. Details on parameters can be inspected in the Ursgal log JSONs that have been uploaded with the raw result files. In total, 3,234,459 and 2,687,720 peptide spectrum matches for PEPs ≤5% and ≤1%, respectively, were observed. This translates to 42,815 and 35,926 distinct peptides for PEPs ≤5% and ≤1%, respectively. Detailed listing of all identified peptides can be found in Supplementary Material 1.

**Metabolically labelled peptide quantification.** All peptides (PEP ≤5%) were subjected to quantification using pyQms (refs 31,70) using default parameters. To quantify peptides irrespectively of their identification in a given LC-MS/MS run, retention time alignment was performed as described before³¹,⁷⁰ All quantifications were grouped and processed using piqDB, an internal database structure based on Python, mongoDB and Hadoop. In total, all quantifications resulted in 1,462,589 matched isotope pattern chromatograms (MICs) containing at least one isotope pattern match with a mScore ≥0.7. However, to include a MIC further into the quantitative analysis a total of eight quantifications (MS1 hits) were required with a mScore ≥0.7. Finally, protein quantifications were based on peptide identification with a PEP threshold ≤5% but with a more stringent quantification mScore threshold of 0.8 (for at least one match within the MIC). All circle plots, showing raw quantification results for all peptide (incl. charge) in
all conditions were plotted using piqDB. Detailed listing of all quantified peptides can be found in Supplementary Data 3.

**Unsupervised clustering.** In total, 1,289 proteins were quantified in all six conditions and those were subjected to unsupervised clustering using pyGCluster⁴³. All possible combinations between the distance matrices 'Euclidean' and 'Correlation' and the linkages 'complete', 'average', 'weighted', 'centroid', 'median' and 'ward' were clustered during each of the 250,000 iterations. Finally, all agglomerative hierarchical cluster results were combined using pyGCluster's meta clustering approach resulting in 146 frequent (obCoFreq ≥ 0.1%) clusters forming five communities (Supplementary Fig. 7), namely community 21, 92, 99, 113 and 139.

**Data availability.** Coordinates and structure factors have been deposited in the world-wide Protein Data Bank with accession codes of 5E37. Mass spectrometry data has been uploaded to the ProteomeXchange Consortium with the data set identifier PXD003049 for the affinity trap experiment and PXD003759 for the isotopic labelling experiment. All other data that support the findings of this study are available from the corresponding author on reasonable request.

## References

1. Gifford, J. L., Walsh, M. P. & Vogel, H. J. Structures and metal-ion-binding properties of the Ca2⁺-binding helix-loop-helix EF-hand motifs. *Biochem. J.* **405**, 199-221 (2007).

2. Day, I. S., Reddy, V. S., Shad Ali, G. & Reddy, A. S. Analysis of EF-hand-containing proteins in Arabidopsis. *Genome Biol.* **3**, RESEARCH0056 (2002).

3. White, P. J. & Broadley, M. R. Calcium in plants. *Ann. Bot. Lond.* **92**, 487–511 (2003).

4. McCormack, E. & Braam, J. Calmodulins and related potential calcium sensors of Arabidopsis. *New Phytol.* **159**, 585–598 (2003).

5. Bender, K. W. & Snedden, W. A. Calmodulin-related proteins step out from the shadow of their namesake. *Plant Physiol.* **163**, 486–495 (2013).

6. Bouche, N., Yellin, A., Snedden, W. A. & Fromm, H. Plant-specific calmodulin-binding proteins. *Annu. Rev. Plant Biol.* **56**, 435–466 (2005).

7. Hochmal, A. K., Schulze, S., Trompelt, K. & Hippler, M. Calcium-dependent regulation of photosynthesis. *Biochim. Biophys. Acta* **1847**, 993–1003 (2015).

8. Stael, S. *et al*. Plant organellar calcium signalling: an emerging field. *J. Exp. Bot.* **63**, 1525–1542 (2012).

9. Weinl, S. *et al*. A plastid protein crucial for Ca2⁺-regulated stomatal responses. *New Phytol.* **179**, 675–686 (2008).

10. Nomura, H., Komori, T., Kobori, M., Nakahira, Y. & Shiina, T. Evidence for chloroplast control of external Ca2⁺-induced cytosolic Ca2⁺ transients and stomatal closure. *Plant J.* **53**, 988–998 (2008).

11. Nomura, H. *et al*. Chloroplast-mediated activation of plant immune signalling in Arabidopsis. *Nat. Commun.* **3**, 926 (2012).

12. Petroutsos, D. *et al*. The chloroplast calcium sensor CAS is required for photoacclimation in *Chlamydomonas reinhardtii*. *Plant Cell* **23**, 2950–2963 (2011).

13. Peers, G. *et al*. An ancient light-harvesting protein is critical for the regulation of algal photosynthesis. *Nature* **462**, 518–521 (2009).

14. Busch, A. & Hippler, M. The structure and function of eukaryotic photosystem I. *Biochim. Biophys. Acta* **1807**, 864–877 (2011).

15. Wang, C., Yamamoto, H. & Shikanai, T. Role of cyclic electron transport around photosystem I in regulating proton motive force. *Biochim. Biophys. Acta* **1847**, 931–938 (2015).

16. Petroutsos, D. *et al*. PGRL1 participates in iron-induced remodeling of the photosynthetic apparatus and in energy metabolism in *Chlamydomonas reinhardtii*. *J. Biol. Chem.* **284**, 32770–32781 (2009).

17. Tolleter, D. *et al*. Control of hydrogen photoproduction by the proton gradient generated by cyclic electron flow in *Chlamydomonas reinhardtii*. *Plant Cell* **23**, 2619–2630 (2011).

18. DalCorso, G. *et al*. A complex containing PGRL1 and PGR5 is involved in the switch between linear and cyclic electron flow in Arabidopsis. *Cell* **132**, 273–285 (2008).

19. Hertle, A. P. *et al*. PGRL1 is the elusive ferredoxin-plastoquinone reductase in photosynthetic cyclic electron flow. *Mol. Cell* **49**, 511–523 (2013).

20. Terashima, M. *et al*. Calcium-dependent regulation of cyclic photosynthetic electron transfer by a CAS, ANR1, and PGRL1 complex. *Proc. Natl Acad. Sci. USA* **109**, 17717–17722 (2012).

21. Takahashi, H., Clowez, S., Wollman, F. A., Vallon, O. & Rappaport, F. Cyclic electron flow is redox-controlled but independent of state transition. *Nat. Commun.* **4**, 1954 (2013).

22. Strand, D. D. *et al*. Activation of cyclic electron flow by hydrogen peroxide in vivo. *Proc. Natl Acad. Sci. USA* **112**, 5539–5544 (2015).

23. Strand, D. D., Fisher, N., Davis, G. A. & Kramer, D. M. Redox regulation of the antimycin A sensitive pathway of cyclic electron flow around photosystem I in higher plant thylakoids. *Biochim. Biophys. Acta* **1857**, 1–6 (2016).

24. Nikkanen, L. & Rintamaki, E. Thioredoxin-dependent regulatory networks in chloroplasts under fluctuating light conditions. *Philos. Trans. R. Soc. Lond. B. Biol. Sci.* **369**, 20130224 (2014).

25. Serrato, A. J., Perez-Ruiz, J. M., Spinola, M. C. & Cejudo, F. J. A novel NADPH thioredoxin reductase, localized in the chloroplast, which deficiency causes hypersensitivity to abiotic stress in Arabidopsis thaliana. *J. Biol. Chem.* **279**, 43821–43827 (2004).

26. Collet, J. F. & Messens, J. Structure, function, and mechanism of thioredoxin proteins. *Antioxid. Redox Signal.* **13**, 1205–1216 (2010).

27. Dangoor, I., Peled-Zehavi, H., Wittenberg, G. & Danon, A. A chloroplast light-regulated oxidative sensor for moderate light intensity in Arabidopsis. *Plant Cell* **24**, 1894–1906 (2012).

28. Perez-Ruiz, J. M. *et al*. Rice NTRC is a high-efficiency redox system for chloroplast protection against oxidative damage. *Plant Cell* **18**, 2356–2368 (2006).

29. Cejudo, F. J., Ferrandez, J., Cano, B., Puerto-Galan, L. & Guinea, M. The function of the NADPH thioredoxin reductase C-2-Cys peroxiredoxin system in plastid redox regulation and signalling. *FEBS Lett.* **586**, 2974–2980 (2012).

30. Dietz, K. J. Peroxiredoxins in plants and cyanobacteria. *Antioxid. Redox Signal.* **15**, 1129–1159 (2011).

31. Hohner, R. *et al*. The metabolic status drives acclimation of iron deficiency responses in *Chlamydomonas reinhardtii* as revealed by proteomics based hierarchical clustering and reverse genetics. *Mol. Cell. Proteomics* **12**, 2774–2790 (2013).

32. Hirasawa, M. *et al*. Oxidation-reduction properties of chloroplast thioredoxins, ferredoxin:thioredoxin reductase, and thioredoxin f-regulated enzymes. *Biochemistry* **38**, 5200–5205 (1999).

33. Emanuelsson, O., Nielsen, H. & von Heijne, G. Chlorop, a neural network-based method for predicting chloroplast transit peptides and their cleavage sites. *Protein Sci.* **8**, 978–984 (1999).

34. Tardif, M. *et al*. PredAlgo: a new subcellular localization prediction tool dedicated to green algae. *Mol. Biol. Evol.* **29**, 3625–3639 (2012).

35. Rasala, B. A. *et al*. Robust expression and secretion of Xylanase1 in *Chlamydomonas reinhardtii* by fusion to a selection gene and processing with the FMDV 2A peptide. *PLoS One* **7**, e43349 (2012).

36. Rasala, B. A., Chao, S. S., Pier, M., Barrera, D. J. & Mayfield, S. P. Enhanced genetic tools for engineering multigene traits into green algae. *PLoS One* **9**, e94028 (2014).

37. Balmer, Y. *et al*. Proteomics gives insight into the regulatory function of chloroplast thioredoxins. *Proc. Natl Acad. Sci. USA* **100**, 370–375 (2003).

38. Motohashi, K., Kondoh, A., Stumpp, M. T. & Hisabori, T. Comprehensive survey of proteins targeted by chloroplast thioredoxin. *Proc. Natl Acad. Sci. USA* **98**, 11224–11229 (2001).

39. Goyer, A. *et al*. Isolation and characterization of a thioredoxin-dependent peroxidase from *Chlamydomonas reinhardtii*. *Eur. J. Biochem.* **269**, 272–282 (2002).

40. Vavilin, D. V., Ducruet, J. M., Matorin, D. N., Venediktov, P. S. & Rubin, A. B. Membrane lipid peroxidation, cell viability and Photosystem II activity in the green alga Chlorella pyrenoidosa subjected to various stress conditions. *J. Photochem. Photobiol. B* **42**, 233–239 (1998).

41. Hodges, D. M., DeLong, J. M., Forney, C. F. & Prange, R. K. Improving the thiobarbituric acid-reactive-substances assay for estimating lipid peroxidation in plant tissues containing anthocyanin and other interfering compounds. *Planta* **207**, 604–611 (1999).

42. Wisniewski, J. R., Zougman, A., Nagaraj, N. & Mann, M. Universal sample preparation method for proteome analysis. *Nat. Methods* **6**, 359–362 (2009).

43. Jaeger, D., Barth, J., Niehues, A. & Fufezan, C. pyGCluster, a novel hierarchical clustering approach. *Bioinformatics* **30**, 896–898 (2014).

44. Morisse, S. *et al*. Thioredoxin-dependent redox regulation of chloroplastic phosphoglycerate kinase from *Chlamydomonas reinhardtii*. *J. Biol. Chem.* **289**, 30012–30024 (2014).

45. Kojima, R. *et al*. Radically different thioredoxin domain arrangement of ERp46, an efficient disulfide bond introducer of the mammalian PDI family. *Structure* **22**, 431–443 (2014).

46. Sato, Y. *et al*. Synergistic cooperation of PDI family members in peroxiredoxin 4-driven oxidative protein folding. *Sci. Rep.* **3**, 2456 (2013).

47. Echalier, A. *et al*. Crystal structure and solution NMR dynamics of a D (type II) peroxiredoxin glutaredoxin and thioredoxin dependent: a new insight into the peroxiredoxin oligomerism. *Biochemistry* **44**, 1755–1767 (2005).

48. Terashima, M., Specht, M., Naumann, B. & Hippler, M. Characterizing the anaerobic response of *Chlamydomonas reinhardtii* by quantitative proteomics. *Mol. Cell. Proteomics* **9**, 1514–1532 (2010).

49. Lemaire, S. D., Collin, V., Keryer, E., Quesada, A. & Miginiac-Maslow, M. Characterization of thioredoxin y, a new type of thioredoxin identified in the genome of *Chlamydomonas reinhardtii*. *FEBS Lett.* **543**, 87–92 (2003).

50. Michelet, L. *et al*. Glutathionylation of chloroplast thioredoxin f is a redox signaling mechanism in plants. *Proc. Natl Acad. Sci. USA* **102**, 16478–16483 (2005).
51. Jans, F. et al. A type II NAD(P)H dehydrogenase mediates light-independent plastoquinone reduction in the chloroplast of *Chlamydomonas*. *Proc. Natl Acad. Sci. USA* **105**, 20546–20551 (2008).

52. Johnson, X. & Alric, J. Interaction between starch breakdown, acetate assimilation, and photosynthetic cyclic electron flow in *Chlamydomonas reinhardtii*. *J. Biol. Chem.* **287**, 26445–26452 (2012).

53. Alric, J., Lavergne, J. & Rappaport, F. Redox and ATP control of photosynthetic cyclic electron flow in *Chlamydomonas reinhardtii* (I) aerobic conditions. *Biochim. Biophys. Acta* **1797**, 44–51 (2010).

54. Harris, E. H. *The Chlamydomonas Sourcebook. A Comprehensive Guide To Biology And Laboratory Use* (Academic Press, 1989).

55. Lauersen, K. J., Kruse, O. & Mussgnug, J. H. Targeted expression of nuclear transgenes in *Chlamydomonas reinhardtii* with a versatile, modular vector toolkit. *Appl. Microbiol. Biotechnol.* **99**, 3491–3503 (2015).

56. Kindle, K. L., Schnell, R. A., Fernandez, E. & Lefebvre, P. A. Stable nuclear transformation of *Chlamydomonas* using the *Chlamydomonas* gene for nitrate reductase. *J. Cell Biol.* **109**, 2589–2260 (1989).

57. Molnar, A. et al. Highly specific gene silencing by artificial microRNAs in the unicellular alga *Chlamydomonas reinhardtii*. *Plant J.* **58**, 165–174 (2009).

58. Shimogawara, K., Fujiwara, S., Grossman, A. & Usuda, H. High-efficiency transformation of *Chlamydomonas reinhardtii* by electroporation. *Genetics* **148**, 1821–1828 (1998).

59. Gonzalez-Ballester, D. et al. Reverse genetics in *Chlamydomonas*: a platform for isolating insertional mutants. *Plant Methods* **7**, 24 (2011).

60. Zhang, R. et al. High-throughput genotyping of green algal mutants reveals random distribution of mutagenic insertion sites and endonucleolytic cleavage of transforming DNA. *Plant Cell* **26**, 1398–1409 (2014).

61. Hou, Y. et al. Functional analysis of an individual IFT protein: IFT46 is required for transport of outer dynein arms into flagella. *J. Cell Biol.* **176**, 653–665 (2007).

62. Moore, E. C., Reichard, P. & Thelander, L. Enzymatic synthesis of deoxyribonucleotides. V. purification and properties of thioredoxin reductase from *Escherichia coli*. *B. J. Biol. Chem.* **239**, 3445–3452 (1964).

63. Waadt, R., Schlucking, K., Schroeder, J. I. & Kudla, J. Protein fragment bimolecular fluorescence complementation analyses for the in vivo study of protein-protein interactions and cellular protein complex localizations. *Methods Mol. Biol.* **1062**, 629–658 (2014).

64. Verdoucq, L., Vignols, F., Jacquot, J. P., Chartier, Y. & Meyer, Y. In vivo characterization of a thioredoxin h target protein defines a new peroxiredoxin family. *J. Biol. Chem.* **274**, 19714–19722 (1999).

65. Cox, J. & Mann, M. MaxQuant enables high peptide identification rates, individualized p.p.b.-range mass accuracies and proteome-wide protein quantification. *Nat. Biotechnol.* **26**, 1367–1372 (2008).

66. Nelson, K. J. & Parsonage, D. Measurement of peroxiredoxin activity. *Curr. Protoc. Toxicol. Chapter* **7**, Unit7–10 (2011).

67. Naumann, B. et al. Comparative quantitative proteomics to investigate the remodeling of bioenergetic pathways under iron deficiency in *Chlamydomonas reinhardtii*. *Proteomics* **7**, 3964–3979 (2007).

68. Kremer, L. P., Leufken, J., Oyunchimeg, P., Schulze, S. & Fufezan, C. Ursgal, universal python module combining common bottom-up proteomics tools for large-scale analysis. *J. Proteome Res.* **15**, 788–794 (2016).

69. Kall, L., Canterbury, J. D., Weston, J., Noble, W. S. & MacCoss, M. J. Semi-supervised learning for peptide identification from shotgun proteomics datasets. *Nat. Methods* **4**, 923–925 (2007).

70. Barth, J. et al. The interplay of light and oxygen in the reactive oxygen stress response of *Chlamydomonas reinhardtii* dissected by quantitative mass spectrometry. *Mol. Cell. Proteomics* **13**, 969–989 (2014).

## Acknowledgements

M.H. acknowledges support from the Deutsche Forschungsgemeinschaft (DFG, HI 739/9-1, FOR964)). G.K. acknowledges support from the Funding Program for Next Generation World-Leading Researchers (GS016) from the Cabinet Office of Japan. This work was supported in part by an International Joint Research Promotion Program and International Collaborative Research Program, Osaka University, (G.K. and M.H.). G.F. acknowledges funding from the Marie Curie ITN Accliphot, and the HFSP grants. K.H. acknowledges support from the National Natural Science Foundation of China (31171287). We are grateful to Dr Jörg Kudla for access to the NanoTemper station (J.K. acknowledges support from DFG (FOR964)) and to the staff at beamline BL44XU at SPring-8, Japan for their kind help during the data collection. We thank Alan Wanke for initial help with the project.

## Author contributions

A.K.H., K.Z., R.C., P.G., S.S., G.L., A.N., D.P. and G.F. conducted the experiments; A.K.H., K.Z., R.C., P.G., R.M., H.T., S.S., G.L., M.S., J.N.O., D.P., G.F., C.F., K.H., G.K. and M.H. analysed the data; G.K. and M.H. wrote the manuscript.

## Additional information

**Supplementary Information** accompanies this paper at http://www.nature.com/naturecommunications

**Competing financial interests:** The authors declare no competing financial interests.

**Reprints and permission** information is available online at http://npg.nature.com/reprintsandpermissions.

**How to cite this article:** Hochmal, A.K. et al. Calredoxin represents a novel type of calcium-dependent sensor-responder connected to redox regulation in the chloroplast. *Nat. Commun.* **7**:11847 (doi: 10.1038/ncomms11847 (2016)).

<figure><img src="image_21.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/