Original Research ARTICLE
Kinetic Modeling of the Arabidopsis Cryptochrome Photocycle: FADHo Accumulation Correlates with Biological Activity
- 1Department of Physics and Astronomy, University of California, Irvine, Irvine, CA, USA
- 2UMR 8256 (B2A), IBPS, University of Paris VI, Paris, France
- 3Department of Physics, Xavier University, Cincinnati, OH, USA
- 4Department of Biology, Xavier University, Cincinnati, OH, USA
Cryptochromes are flavoprotein photoreceptors with multiple signaling roles during plant de-etiolation and development. Arabidopsis cryptochromes (cry1 and cry2) absorb light through an oxidized flavin (FADox) cofactor which undergoes reduction to both FADH° and FADH− redox states. Since the FADH° redox state has been linked to biological activity, it is important to estimate its concentration formed upon illumination in vivo. Here we model the photocycle of isolated cry1 and cry2 proteins with a three-state kinetic model. Our model fits the experimental data for flavin photoconversion in vitro for both cry1 and cry2, providing calculated quantum yields which are significantly lower in cry1 than for cry2. The model was applied to the cryptochrome photocycle in vivo using biological activity in plants as a readout for FADH° concentration. The fit to the in vivo data provided quantum yields for cry1 and cry2 flavin reduction similar to those obtained in vitro, with decreased cry1 quantum yield as compared to cry2. These results validate our assumption that FADH° concentration correlates with biological activity. This is the first reported attempt at kinetic modeling of the cryptochrome photocycle in relation to macroscopic signaling events in vivo, and thereby provides a theoretical framework to the components of the photocycle that are necessary for cryptochrome response to environmental signals.
Plants adapt to their light environment by means of multiple photoreceptors which optimally absorb at different wavelengths of light throughout the visible spectrum. These include specific photoreceptors absorbing in the blue—UV/A such as cryptochromes (Chaves et al., 2011; Wang et al., 2014) and phototropins (Christie et al., 2015), red/far red light absorbing phytochromes (Burgie and Vierstra, 2014; Xu et al., 2015) and UV-B specific receptors (Jenkins, 2014). Light sensitivity is achieved through pigment molecules (chromophores) bound to a protein backbone (apoprotein). The pigments absorb photons at specific wavelengths of light to initiate a primary photochemical reaction. These reactions, in turn, trigger changes within the photoreceptor apoprotein leading to the initiation of biological signaling. Generally, such changes involve conformational change in the protein which allows access to signaling partners and/or modifications such as phosphorylation or ubiquitination (Galvão and Fankhauser, 2015).
In Arabidopsis, two cryptochromes have been shown to mediate significant signaling functions (cry1 and cry2) (Chaves et al., 2011). These proteins are highly conserved within their first 500 amino acid residues, which comprise the N-terminal flavin binding domains. This domain absorbs light and undergoes the primary photochemical reactions involving intra-protein electron and proton transfer to the flavin. By contrast, cry1 and cry2 diverge greatly at their C-terminal domains, which undergo conformational change involved in signaling. Cry1 plays a key role in de-etiolation responses and photomorphogenesis as well as in many aspects of vegetative growth. Cry2 also plays a role in seedling photomorphogenesis, including hypocotyl growth inhibition and cotyledon expansion. However, Cry2 function during de-etiolation is apparent primarily at low blue light intensity and not at high light. This specificity of cry2 for conditions of dim blue light is thought to follow from the fact that cry2, once activated by blue light, undergoes rapid ubiquitination resulting in targeting to the proteosome and degradation (Yu et al., 2007, 2009). In addition to its role in de-etiolation responses, cry2 has been implicated in the Arabidopsis photoperiodic initiation of flowering response wherein long days induce earlier flowering than short days (Valverde et al., 2004). Mechanistically, both cry1 and cry2 have been shown to interact with signaling partners (CIB1, SPA1) (Liu et al., 2008, 2011) in response to illumination, indicative of a light induced conformational change leading to substrate binding. In this respect, cryptochrome functions similarly to other classes of known plant photoreceptors (phytochrome, phototropin, and UVR8 type receptors) which also undergo conformational changes in response to illumination.
Photochemical reactions of cryptochromes are induced by light absorption through the flavin (FAD) chromophore, and have been well characterized (reviewed in Chaves et al., 2011). Briefly, cryptochrome—bound FAD occurs in the oxidized (FADox) state in the dark. Upon illumination, the excited state flavin is reduced via multiple electron and proton transfer events to a mixture of neutral radical (FADH°) and fully reduced (FADH−) flavin redox states. Once formed, the reduced redox state intermediates are relatively stable (on the order of minutes) and undergo reoxidation to the dark (FADox) resting state at rates that are determined by the concentration of molecular oxygen (Müller and Ahmad, 2011). Therefore, the proportion of cryptochrome in any given redox state under constant illumination is determined by the steady state equilibrium reached between the forward (light driven photoreduction to FADH° and FADH−) and reverse (reoxidation to FADox) reactions. A description of this redox cycle is shown in Figure 1.
Figure 1. Reaction scheme of the plant cryptochrome photocycle. In the dark, flavin occurs in the oxidized (blue light absorbing) redox state (FADox). Subsequent to illumination, the flavin is photoreduced to the neutral radical (FADH°) redox form by blue light at a defined rate (k1) determined by the quantum yield. The FADH° can be further reduced to the fully reduced (FADH−) form by blue or green light with at a rate k2. The reduced flavins (FADH° and FADH−) reoxidize to the resting state (FADox) (k1b, k2b) at rates that are independent of light. Equilibrium concentrations of the different redox forms of cryptochrome result from the sum of the four reaction rates (k1, k2, k1b, k2b).
Many lines of evidence have identified the radical (FADH°) redox intermediate is the signaling state (biologically activate) form of cryptochrome. Briefly, the action spectrum for cryptochrome (Ahmad et al., 2002) and in vivo measurements of cryptochrome-bound flavin redox state in living cells indicate that the inactive (light-absorbing) state of the photoreceptor is the oxidized (FADox) redox state (Banerjee et al., 2007; Bouly et al., 2007; Engelhard et al., 2014). Upon illumination, a large conformational change has been shown to occur in vitro and linked to formation of the FADH° redox state (Kondoh et al., 2011). The lifetime of this FADH° redox state correlates with that of the biologically active “lit” state for both cry1 and cry2 in vivo (Herbel et al., 2013). Significantly, both cry1 and cry2 biological activity is diminished under illumination treatments (green light) that deplete the concentration of the FADH° redox state (see Figure 1; Banerjee et al., 2007; Bouly et al., 2007; Herbel et al., 2013). Finally, in each instance where in vivo biological activity has been reported to occur in mutants of cry2 or cry1 in the literature (Li et al., 2011; Gao et al., 2015) the light dependent reduction of the flavin cofactor has also been shown to occur under in vivo conditions (Engelhard et al., 2014; El-Esawi et al., 2015). From these and many other studies (see Chaves et al., 2011), the available evidence is consistent with Arabidopsis cryptochromes being activated by flavin photoreduction and that the neutral radical (FADH°) redox state represents the “lit” state in vivo.
Assigning a given redox form of cryptochrome as the “lit” state raises the challenge of directly linking biological activation to the concentration of this redox form induced upon illumination. This task is further complicated by the fact that, unlike most plant photoreceptors which are converted by light into well-defined “on” and “off” chemical states, cryptochromes in fact exist in three redox states. That is, illumination induces formation of not one but two reduced states (FADH° and FADH−), each with a characteristic (and different) reoxidation rate that occurs continuously throughout the illumination period. Therefore, the concentration of the “activated” redox state results from the equilibrium defined by the sum rate constants (k1, k2, k1b, k2b) that must be considered as a function of the quality (wavelength) and intensity of the light signal.
In this study, we fit a three-state kinetic model of cryptochrome photocycle to absorption spectra of isolated proteins illuminated in vitro. We use the kinetic model to obtain the rate constants (Figure 1) and calculate the quantum yield of the light-driven reactions (FADox to FADH° and FADH° to FADH−). To apply the kinetic model to cryptochrome responses in vivo, it is not possible in living plants to directly measure the concentration of the different redox state intermediates. We therefore consider biological activity as a measure for FADH° concentration formed in vivo. Indeed, the experimental data for cryptochrome responses in vivo at different intensities and wavelengths of light provided an excellent fit for the kinetic model applied to spectra in vitro. This model thereby validates the assumption that FADH° is the signaling state, and provides further insight into many of the known characteristics of cryptochrome responses in vivo.
Materials and Methods
Cryptochrome Protein Samples
Arabidopsis cryptochrome-1 (cry1) and cryptochrome-2 (cry2) proteins were expressed and purified using baculovirus expression constructs in insect cell cultures as previously described (Banerjee et al., 2007; Bouly et al., 2007). Photoreduction experiments were conducted by illuminating protein samples in PBS (Phosphate Buffered Saline: 5 mM NaPO4 pH 7.5 or 8.2, 150 mM NaCl directly in quartz cuvettes, 1 cm path length) at 21°C. Spectra were taken using a Cary 300 UV/Vis spectrophotometer at 21°C. Concentration of cry1 and cry2 protein samples was between 50 and 150 μM for in vitro experiments.
LEDS were from Quadica Developments, Brantford, Ontario. Blue light: Royal-Blue (447.5 nm). Green Light: Lime-Green (567 nm). Spectra are as indicated by the manufacturer.
Seeds from phyAphyB Arabidopsis phytochrome-deficient mutants from ecotype Landsberg erecta (Ahmad and Cashmore, 1997) were sown and germinated as described previously (Herbel et al., 2013). Germinating seedlings were maintained for four days at 22°C in darkness. Etiolated seedlings were then illuminated, harvested into liquid nitrogen, the proteins extracted and then applied to nitrocellulose membranes for Western blotting as described previously (Herbel et al., 2013). Quantitation of the cry2 signal from the Western blots was by ImageJ image analysis software from photographic images of the blots.
The reaction scheme depicted in Figure 1 reports the cryptochrome photocycle studied here. The time evolution of the intermediate states concentrations is described by a set of coupled first-order kinetic equations (Espensen, 1981):
where k1 and k2 are the two forward rate constants, and k1b and k2b are the dark reoxidation rate constants. In Equation (1) square brackets denote the concentrations of the transient states FADox, FADHo, and FADH−. The fit of the kinetic model to absorption spectra of the isolated proteins allows to find the rate constants. To this end we apply the Beer-Lambert law, which relates the concentration of the transient states to the absorbance A (Schmidt, 2005). According to the Beer-Lambert law the absorbance A at a given wavelength λ and at time t is linearly dependent on the concentration of the absorbing species:
where N is the number of different light-absorbing species in the system with concentration ci (M). εi (M−1 cm−1) is the molar extinction coefficient and d (cm) is the thickness of the absorbing medium. In all our experiments we use the same quartz cuvette of path length d = 1 cm.
We record absorption spectra from 400 to 570 nm, by illuminating cry samples with blue (450 ± 10 nm) and green light (560 ± 10 nm). The neutral radical (FADH°) flavin redox state can absorb green as well as blue light, while the FADox absorbs blue (Banerjee et al., 2007; Bouly et al., 2007). The FADH− radical absorbs at wavelengths out of the range considered here (Müller and Ahmad, 2011). Thus, from the Beer-Lambert law we have:
where εox(450) and εH(560) are the molar extinction coefficients, respectively, of FADox and FADHo.
Since cry–bound FAD only occurs in the oxidized (FADox) state in the dark, the initial concentration of cry in the sample before illumination, i.e., at time t = 0, can be found according to A(450, 0) = εox(450)[FADox](0), while the other intermediate states are unpopulated.
Based on the absorption spectra of FADox and FADHo (Liu et al., 2010; Björn, 2015) we estimate that the extinction coefficient of FADHo at 450 nm and 560 nm is the same and equals εH (450) = εH (560) = εox (450)/2.
By normalizing absorbance and concentration to the dark, and considering that εH (450) = εox (450)/2, and εH (560) = εox (450)/2, we obtain from Equation (3) a simplified expression between concentrations and absorbance, which we use in our calculations:
In Equation (4) square brackets and A label, respectively, normalized concentration and normalized absorbance. We numerically solve Equation (1) with a Runge-Kutta method, and obtain the time evolution of the concentration of the transient states corresponding to data acquisition times, with the assumption that at a given time t, the sum of the concentrations of the three states is constant, i.e.,
In Equation (1), while reoxidation rates k1b and k2b are independent of light, the forward rates k1 and k2 are light dependent. At a given wavelength λ the rate constant is given by kλ = σλ Iλ where Iλ is the photon fluence rate (mol m−2 s−1) at wavelength λ and σλ (mol−1 m2) denotes the photoconversion cross-section. σλ is related to the quantum yield ϕλ according to σλ = 2.3 ε(λ) ϕλ (Kendrick and Kronenberg, 1994).
In our model we consider that the efficiency of photoconversion of FADHo upon blue light illumination is negligible (see Section Results). Thus the forward rates k1 and k2 are related to the photon fluence rate of blue light I1 and of green light I2 by, respectively, k1 = σ1 I1 and k2, k2 = σ2I2. I2, where ε1 = εox(450) and ε2 = εH(560).
Two-State Model: Dark Reoxidation k1b, and Quantum Yield ϕ1 under Blue Light
In order to determine the rate constants k1 and k1b, we first consider an abbreviated photocycle consisting of two states FADox and FADHo. This is possible because in atmospheric oxygen and in the presence of mild reductants, isolated Arabidopsis cry1 and cry2 proteins, under blue light illumination, accumulate in primarily two redox forms: FADox and FADH° (Banerjee et al., 2007; Bouly et al., 2007). The two-state kinetic model is described by the following system of differential equations:
We set the initial concentration of cry at time t = 0, i.e., before illumination to [FADox](0) = 1. Furthermore, at a given later time t . Analytical solutions of Equation (5) are straightforward (see Equations S1, S2 in Supplementary Material).
Dark reoxidation rate (k1b)
Since flavin reoxidation from FADH° to FADox can be readily monitored spectroscopically, we first obtain experimental values for the dark reoxidation kinetics from reduced (FADH°) to oxidized (FADox) flavin states for the cry.
We define the reoxidation time as td, and td = 0 is the time in which cry is placed in darkness after being illuminated for a certain time t at fluence rate I1. We record the absorption spectrum after increasing times in darkness td, until complete reoxidation to FADox. When only the dark reoxidation occurs (k1 = 0) the analytical solutions of Equation (2) are:
where cox and co are the initial concentrations of, respectively, FADox and FADHo at the dark time td = 0. We plot the concentrations of FADox and FADHo, obtained from the spectra by applying Equation (4), as a function of the dark recovery time td, and fit the data with, respectively, Equations (6) and (7) to find k1b. For the fitting model we use the Least-Squares algorithm Levenberg-Marquardt provided by Matlab. We calculate the half-life τ1∕2 that, for a first-order reaction, is giving by τ1∕2 = ln(2) .
Forward photoreduction rate (k1). Two-state-based algorithm
For an illumination time much smaller than the reoxidation time k1b = 0 in Equation (2). The analytical solutions for this case can be used to calculate the rate constant k1 from the absorption spectra (see Equations S3, S4 in Supplementary Material). However, for longer illumination times, the reoxidation rate has to be taken into account. For this we implement a simple algorithm based on Equation (5). Such algorithm takes as input either the [FADox] or [FADHo] concentration, the k1b previously found, and outputs the corresponding rate constant k1 according to the two-state kinetic model. The concentration of FADox and FADHo are obtained from the spectra by using Equation (4). In our calculations we use as input both [FADox] and [FADHo] concentration and compare the results to confirm that the two-state model approximates well the blue-light experiments.
Quantum yield of FADox—FADH° conversion
We calculate the quantum yield ϕ1 from the photoconversion cross section σ1 according to σ1 = 2.3 εox(450) ϕ1 (Kendrick and Kronenberg, 1994). To obtain the photoconversion cross section we illuminate cry sample with increasing blue light fluence rates I1 and record absorption spectra. For each photon fluence rate I1 we calculate k1 with the two-state algorithm described above. Plotting I1 vs. k1, and fitting the data with a linear function, k1 = σ1 I1, allows estimation of σ1. From σ1 we calculate ϕ1 by using the experimentally determined extinction coefficient of cry1 and cry2 (see Section Extinction Coefficient). We point out that only one single photon fluence rate value is enough to calculate the quantum yield. However, we prefer to use a dose-response profile to find a linear range of blue light fluence rates that both allows prediction, and confirms the linear correlation between k1 and I1.
Three-State Model: Forward Rate k2 and Quantum Yield ϕ2 under Green Light
Forward photoreaction rate k2. Three-state-based algorithm
To find the rate constant k2 we co-illuminate samples with blue and green light at fluence rates, respectively, I1 and I2, and record the absorption spectrum. From the spectra we obtain the normalized concentration of [FADox] and [FADHo] by using Equation (4). To calculate k2 we implement an algorithm by numerically solving Equation (1). This algorithm takes as input the concentration of either [FADox] or [FADHo], k1, k1b and k2b, and outputs k2, according to the three-state model. We use the k1 and k1b values obtained from the blue light experiments in the present study, and k2b = 0.011 s−1 is provided from previous studies (Müller and Ahmad, 2011). For short illumination times, one can neglect reoxidation rates, and the algorithm in this case resolves Equation (1) with k1b = k2b = 0.
Quantum yield of FADHo—FADH− conversion
We calculate the quantum yield ϕ2 from the photo-conversion cross section σ2 according to σ2 = 2.3 ε2 ϕ2. To this end we perform a series of experiments to obtain a dose-response profile. We illuminate cry samples with same blue light fluence rate I1 and increasing green light fluence rates I2. For each photon fluence rate I2 we calculate the corresponding k2 by using the three-state algorithm. We plot I2 vs. k2 and fit the data with a linear function, k2 = σ2 I2, to estimate the photo-conversion cross section σ2. We then calculate ϕ2 by considering that ε2 = εox(450)/2.
We calculated the extinction coefficient at 450 nm of cry1 and cry2 using absorption of the purified cryptochrome protein at 450 nm together with protein concentration determined by Bradford assay. The extinction coefficient of cry1 resulted εox (450) = 6415.5 M−1 cm−1, and that of cry2 εox (450) = 5094 M−1 cm−1. see Supplementary Figure A. In our calculations we use as units for the extinction coefficient mol−1 m2, thus for cry1 εox(450) = 641.55 mol−1 m2, and for cry2 εox(450) = 509.4 mol−1 m2.
Kinetic Model Applied to In vivo Responses
To apply the kinetic model to in vivo responses we make the assumption that biological activity is directly proportional to the concentration of the FADH° flavin state (Banerjee et al., 2007; Bouly et al., 2007). In the case of cry1, we choose blue-light dependent inhibition of hypocotyl elongation as a “readout” for biological activity, and consider the concentration of FADH° as inversely proportional to the blue-light dependent inhibition of hypocotyl length (L). In the case of cry2 we use light-dependent degradation of cry2 as a readout for biological function, and we consider the concentration of FADH° as inversely proportional to the cry2 protein concentration (C). In both cases, phytochrome-deficient phyAphyB mutants were used for assay of cry-dependent function, in order to avoid the considerable effect of phytochrome, which also absorbs in the blue and green spectral regions and significantly enhances the sensitivity of cryptochrome-dependent signaling pathways (Ahmad and Cashmore, 1997). Therefore, in our studies, only the effects of light on the cry (blue light receptor) are detected as biological activity.
The length L or concentration C are measured as a function of the photon fluence rate, thereby obtaining a dose-dependent biological response profile. We convert this light dose-biological response curve into a light dose-FADHo concentration curve by using the Equation (S5) given in Supplementary Material. In this way we can calculate quantum yields by applying the same method as for in vitro data.
Two-State Model. Quantum Yield of FADox—FADH° Conversion
For each blue light fluence rate I1, we calculate the rate constant k1 by the two-state algorithm, which takes as input the FADHo values and the dark reoxidation k1b provided from previous studies (Herbel et al., 2013). By plotting I1 vs. k1, and fitting the data with expression k1 = σ1 I1, we estimate the photoconversion cross section σ1. We calculate the quantum yield ϕ1 according to σ1 = 2.3 ε1 ϕ1 and by using the experimentally determined extinction coefficient ε1 = εox(450) found for cry in vitro.
Seedlings were co-illuminated with a blue light fluence rate I1, and increasing fluence rates of green light I2. A dose-biological response profile was then converted in dose-FADHo profile. For each photon fluence rate I2 we calculate the rate constant k2 by the three-state-based algorithm, which inputs FADHo, k1 previously obtained with the blue light experiments, and k1b and k2b provided from the literature (Müller and Ahmad, 2011; Herbel et al., 2013). We fit the data with a linear function, k2 = σ2 I2, to estimate σ2 and calculate the quantum yield ϕ2 using the in vitro extinction coefficient ε2.
The goal of this study is to apply a simple kinetic model to the cryptochrome photocycle (Figure 1) that can accurately predict the effects of illumination on redox state interconversion in vitro and relate this model to observations on biological activation in vivo. We first apply the model in vitro, to samples of purified isolated cryptochrome (cry1 and cry2) which were photoreduced under defined illumination conditions. In this way, concentrations of redox state intermediates could be accurately determined and the reaction rates and quantum yields calculated by the model. We next apply the kinetic model to plant cryptochrome responses in vivo to correlate flavin redox state interconversion that could account for biological activity.
Two-State Model for Cryptochrome Photocycle In vitro
For analysis of the cryptochrome photocycle in vitro, samples of purified cry1 and cry2 proteins were photoreduced in vitro and allowed to reoxidize in monochromatic blue light (450 nm). The simpler two-state model (an abbreviated photocycle from FADox to FADH° and back) is valid under these conditions as there is almost no FADH− accumulation (Banerjee et al., 2007; Bouly et al., 2007). We therefore first modeled only the rates k1 and k1b, and used them to obtain quantum yield and half-life of cry under conditions of monochromatic blue light.
Dark Reoxidation Rate (k1b)
A sample of Cry2 (pH = 7.5, with 10 mM Betamercaptoethanol as reducing agent) was illuminated with blue light for 20 s at a fluence rate of I1 = 400 μmol m−2s−1, and then placed in darkness (td = 0). Figure 2A shows the absorption spectrum after increasing times in darkness td, until complete reoxidation to FADox. From the spectra of Figure 2A we obtained the concentrations of FADox and FADHo as a function of the reoxidation time td by using Equation (4). Figure 2B (the red triangles) reports the concentration of FADHo as a function of the dark recovery time td, and the fit (blue curve) of the data with the two-state dark reoxidation model reported in Equation (7). The fit resulted in a reoxidation rate of k1b = 0.003 s−1, or half-life of τ1∕2 = 231 s. We have also fit the increase of FADox with Equation (6) and obtained a similar reoxidation rate (k1b = 0.0036 s−1, τ1∕2 = 192 s with a goodness of the fit R2≈1, confirming that the two-state model well approximates our experiments. The average of the two half-lives is reported in Table 1.
Figure 2. Rate constants and quantum yield for two-state reduction and reoxidation of AtCry2 in PBS pH7.5, 10 mM βME. (A) Isolated purified cry2 protein was illuminated for 20 seconds (s) at 200 μmol m−2 s−1 blue light and placed in darkness (td = 0 s). Normalized absorption spectra are reported at increasing dark intervals. The FADox concentration is obtained from the absorbance at 450 and 560 nm, and FADHo concentration from absorbance at 560 nm according to Equation 4. (B) Normalized concentration of FADHo as a function of the dark reoxidation time (td). The FADHo concentration is obtained from the spectra of panel (A). The red triangles represent the experimental data, and the blue curve is the fit of the experimental data with the two-state dark reoxidation model for FADHo (Equation 7). The calculated reoxidation rate is k1b = 0.003 s−1 (half-life of τ1∕2 = 231 s). The goodness of the fit is large (R2 = 0.996). (C) Isolated purified cry2 protein was illuminated for 20 s at the indicated blue light fluence rates . Normalized absorption spectra are presented. (D) Calculated forward rate constant k1 vs. photon fluence rate I1 (red triangles). For each I1, the rate constant k1 was calculated by using the two-state algorithm, which inputs the concentration of FADHo, obtained from the spectra of panel (C), and k1b from panel (B) and outputs the rate constant k1. The linear fit k1 = σ1 I1 of the data are reported in blue. The photo-conversion cross section resulted σ1 = 1.6 × 10−4 μmol−1 m2. σ1 is related to the quantum yield ϕ1 and the extinction coefficient ε1 according to σ1 = 2.3 ε1 ϕ1. By using the experimentally calculated ϵ1 = 509.4 mol−1 m2 (5094 M−1 cm−1) the quantum yield was ϕ1 = 0.137. For details of calculations see Method Section: Two-State Model: Dark Reoxidation k1b, and Quantum Yield ϕ1 under Blue Light.
Table 1. In vitro and in vivo parameters of the cry1 and cry2 photocycle depicted in Figure 1.
Forward Photoreduction Rate (k1) and Quantum Yield (ϕ1) of FADox—FADH° Conversion
Flavin reduction by light (FADox to FADH°) rate constants, and quantum yields, were estimated from experimental results of photoreduction of isolated cryptochrome proteins by using the two-state kinetic model as described in the Method Section Two-State Model: Dark Reoxidation k1b, and Quantum Yield ϕ1 under Blue Light. In order to derive the forward rate constant k1 as a function of photon fluence rate I1, Figure 2C shows the spectra obtained from cry2 photoreduced at different blue light fluence rates. For each fluence rate I1 we calculated the rate constant k1 by using the two-state-based algorithm. This algorithm takes as input the concentration of FADHo (obtained from Figure 2C according to Equation 4), the dark reoxidation rate k1b previously found (from Figure 2B), and output k1 according to the two-state model. Figure 2D (red triangles) reports the rate constants k1 as a function of the blue light fluence rates I1 (μmol m−2s−1). The blue curve in Figure 2D is the linear fit of the data (k1 = σ1 I1) which allows to estimate the photoconversion cross section σ1. As can be seen (Figure 2D), the goodness of the fit is excellent (R2≈1), and resulted in σ1 = 1.6 × 10−4 μmol−1 m2. We calculated the quantum yield ϕ1 according to σ1 = 2.3 ε1 ϕ1, which resulted ϕ1 = 0.14, using the experimentally derived extinction coefficient of ε1 = 509.4 mol−1 m2 (5094 M−1 cm−1).
By performing the same calculations with FADox as input to the algorithm we obtain similar results, with σ1 = 2.8 × 10−4 μmol−1 m2 (R2 = 0.98), and quantum yield of ϕ1 = 0.24. This similarity confirms that the two-state kinetic model accurately describes the kinetics of cryptochrome flavin reoxidation under blue light derived experimentally. The average of the two quantum yields is reported in Table 1.
Effect of pH on Kinetics (k1 and k1b) of the Cry2 Photocycle
To further test the validity of the two-state modeling approach under blue light, we evaluate the effect of pH change on the kinetics of the cryptochrome photocycle. At pH 8.2, the efficiency of forward electron transfer is reportedly decreased in Arabidopsis cryptochrome-1 (Müller et al., 2014). We accordingly modeled both forward (photoreduction) and back (reoxidation) kinetics of Arabidopsis cry2 at pH 8.2, using the same buffer composition and concentration of reducing agent as at pH 7.5 (see Figure 2).
Firstly, AtCry2 samples were photoreduced and returned to darkness (Figure 3A). Spectra were taken at intervals during the dark reoxidation time (td). Concentration of FADHo, obtained from the spectra of Figure 3A, was then plotted as a function of the reoxidation time td (Figure 3B). The experimental data (red triangles) were fitted with the two-state dark reoxidation model (Equation 7) to find k1b. The resulting reoxidation rate was k1b = 0.004 s−1 (half-life of τ1∕2 = 2.7 min), which is very similar to that found at pH 7.5 (Figures 2A,B). By fitting the concentration of FADox as a function of the dark reoxidation time with Equation (6), we obtained similar results with k1b = 0.003 s−1 (τ1∕2 = 3.2 min, R2≈1). Therefore, the two states approximation accurately models the experimental results in this case as well. Changes of pH do not affect the reoxidation rate of cry2 protein under these experimental conditions.
Figure 3. Rate constants and quantum yield for two-state reduction and reoxidation of AtCry2 in PBS pH 8.2, 10 mM βME. (A) Isolated purified cry2 protein was illuminated for 20 s at 4000 μmol m−2s−1 blue light and placed in darkness (td = 0 s). Normalized absorption spectra are reported at increasing dark intervals. (B) Normalized concentration of FADHo as a function of the dark reoxidation time (td). The FADHo concentration is obtained from the spectra of panel (A). The red triangles represent the experimental data, and the blue curve is the fit of the experimental data with the two-state dark reoxidation model for FADHo (Equation 7). The resulting reoxidation rate was k1b = 0.0043 s−1 (half-life of τ1∕2 = 2.7 min). (C) Isolated purified cry2 protein was illuminated for 10 s at the indicated blue light fluence rates . Normalized absorption spectra are presented. (D) Calculated forward rate constant k1 vs. photon fluence rate I1 (red triangles). For each fluence rate I1 we calculated the rate constant k1 by using the two-state model algorithm, with input of FADHo (from panel C) and k1b from panel (B). The blue curve reports the linear fit of the data (k1 = σ1I1), which gives a photoconversion cross section of σ1 = 4.42 × 10−5 μmol−1 m2, and a quantum yield of ϕ1 = 0.038. For details of calculations see caption of Figure 2.
We next considered the effect of pH on the forward reaction rate (k1). Photoreduction of purified cry2 sample was performed at multiple blue light fluence rates (I1) to obtain decrease at 450 nm and increase at 560 nm indicative of flavin reduction (Figure 3C). From the spectra of Figure 3C we calculated the concentration of FADox and FADHo. For each photon fluence rate I1 we calculated the rate constant k1 by using the concentration of FADHo, and k1b from Figure 3B, as input to the two-state algorithm. Figure 3C reports k1 as function of I1 (red triangles), and the linear fit k1 = σ1I1 (blue curve). The fit gives a photoconversion cross section of σ1 = 4.42 × 10−5μmol−1 m2. The quantum yield was then calculated, as explained in the previous section, and resulted ϕ1 = 0.038. Using the concentration of FADox to find k1, we obtained similar results (σ1 = 5.1 × 10−5 μmol−1 m2, R2 = 0.98, ϕ1 = 0.0435). The quantum yield at pH = 8.2 is therefore one order of magnitude lower than the quantum yield at pH = 7.5.
To summarize the effect of pH change on the cry2 photocycle, our work indicates that the dark reoxidation (k1b) reaction is unaffected by pH change, at least by the change studied here, therefore decrease in quantum yield is due to decreased efficiency in forward light-dependent photoreduction by an order of magnitude. These results are in agreement with prior studies (Müller et al., 2014) validating the accuracy of our model in describing experimental results.
Comparison of Cry1 and Cry2 Photoconversion Efficiency
Cry1 has been described in the literature as playing a principal role at high blue light intensity, whereas cry2 effects, particularly on photomorphogenesis, are more evident at low blue light intensity (Lin et al., 1998). To provide a further test for the relevance of our kinetic model to experimental findings, we compared the calculated reoxidation rate and quantum yield of cry1 and cry2 flavin reduction in vitro, using the two-state model.
Protein samples of cry1 and cry2 were both photoreduced in PBS pH 7.5 with the addition of 5 mM DTT as reductant, since cry1 cannot otherwise be readily reduced under aerobic conditions (Müller and Ahmad, 2011). Dark reoxidation rates for cry1 (Figures 4A,B) and cry2 (Figures 5A,B) were evaluated from spectra taken during the time course of reoxidation by using the two-state model. The reoxidation time td as a function of the FADox concentration is reported in Figures 5B, 6B (red triangles). From the fit (blue curve) the half-life of FADH° to FADox interconversion of cry1 resulted τ1∕2 = 53 s, while that of cry2 was τ1∕2 = 2 min. Similar results are obtained by fitting the decrease in concentration of FADHo with Equation (6) (for cry1 k1b = 0.016 s−1 with R2≈1, and for cry2 k1b = 0.0054 s−1, with R2≈1). Once again, the experimental data showed an excellent fit to the model.
Figure 4. Rate constants and quantum yield for two-state reduction and reoxidation of AtCry1 in 5 mM DTT. (A) Isolated purified cry1 protein was illuminated for 30 s at 200 μmol m−2s−1 blue light and placed in darkness (td = 0 s). Normalized absorption spectra are shown at increasing dark intervals. (B) Normalized concentration of FADox as a function of the dark reoxidation time (td). The FADox concentration is obtained from the spectra of panel (A). The red triangles represent the experimental data, and the blue curve is the fit of the experimental data with Equation (6). The reoxidation rate resulted k1b = 0.013 s−1 (half-life of τ1∕2 = 53 s). (C) Isolated purified cry1 protein was illuminated for 30 s at the indicated blue light fluence rates . Normalized absorption spectra are presented. (D) Calculated forward rate constant k1 vs. photon fluence rate I1 (red triangles). For each I1 of panel (C), the rate constant k1 was calculated by using the two-state model algorithm with concentration of FADox and k1b (from panel B) as input. From the linear fit of the data, k1 = σ1 I1 (shown in blue) the photoconversion cross section resulted in σ1 = 6.4 × 10−5 μmol−1 m2, which gives a quantum yield of ϕ1 = 0.043. For details of calculations see caption of Figure 2.
Figure 5. Rate constants and quantum yield for two-state reduction and reoxidation of AtCry2 in 5 mM DTT. (A) Isolated purified cry2 protein was illuminated for 30 s at 200 μmol m−2s−1 blue light and placed in darkness (td = 0 s). Normalized absorption spectra are reported at increasing dark intervals. (B) Normalized concentration of FADox as a function of the dark reoxidation time (td). The FADox concentration is obtained from the spectra of panel (A). The red triangles represent the experimental data, and the blue curve is the fit of the experimental data with Equation (6). The reoxidation resulted k1b = 0.0058 s−1 (half-life of τ1∕2 = 2 min). (C) Isolated purified cry2 protein was illuminated for 30 s at the indicated blue light fluence rates . Normalized absorption spectra are presented. (D) Calculated forward rate constant k1 vs. photon fluence rate I1 (red triangles). For each fluence rate I1 of panel (C), the rate constant k1 was calculated by using the two-state algorithm, with FADox concentration and k1b from panel (B) as input. From the linear fit of the data k1 = σ1 I1 (shown in blue) the photoconversion cross section resulted σ1 = 2.5 × 10−4 μmol−1 m2, which gives a quantum yield of ϕ1 = 0.213. For details of calculations see caption of Figure 2.
Figure 6. Forward rate constant k2 and quantum yield ϕ2 under green light for AtCry2 in PBS pH7.5, 10 mM βME. (A) Isolated purified cry2 protein was co-illuminate for 20 s at 100 μmol m−2s−1 blue light I1 and an increasing fluence rate of green light I2 as indicated. Normalized absorption spectra are reported. (B) Calculated forward rate constant k2 vs. green light of fluence rate I2 (red triangles). For each photon fluence rate of panel (A) I2, the rate constant k2 was calculated by using the three-state algorithm, with the concentration of [FADHo] and k1 (from Figure 2D) as input. We neglected the reoxidation rates k1b = k2b = 0 (see text). The linear fit k2 = σ2 I2 of the data are reported in blue. σ2 is the photo-conversion cross section, which resulted σ2 = 1.6 × 10−5 μmol−1 m2. From σ2 we calculated a quantum yield of ϕ2 = 0.027, according to σ2 = 2.3 ε2 ϕ2. For details of calculations see Method Section Three-State Model: Forward Rate k2 and Quantum Yield ϕ2 under Green Light.
Next, the light driven FADox to FADH° photoconversion rates were determined experimentally and modeled for cry1 (Figures 4C,D) and cry2 (Figures 5C,D) by using the two-state kinetic model (see Method Section Two-State Model or Cryptochrome Photocycle in vitro). Quantum yield calculated for cry1 was ϕ1 = 0.043 whereas that for cry2 ϕ1 = 0.213.
Similar results were obtained by considering the change in concentration of FADHo to calculate the quantum yield (for cry1 σ1 = 4.8 × 10−5 μmol−1 m2 and ϕ1 = 0.033 with R2 = 0.96, and for cry2 σ1 = 1.76 × 10−4 μmol−1 m2 and ϕ1 = 0.15 with R2 = 0.97). This means that, under comparable illumination and buffer conditions, cry2 has about 10-fold more efficient response to light than does cry1.
In sum, the two-state kinetic model accurately describes the experimental data for both cry1 and cry2 photoreduction. The obtained quantum yield for cry1, which is 10-fold lower as compared to cry2, is furthermore in good agreement with the biological role of cry2 at lower blue light intensities.
Three-State Model of Cryptochrome Photoreduction
Under conditions of monochromatic blue light illumination, it is evident from our above results that a two-state model adequately describes the cryptochrome photocycle. This follows from the fact that only two redox states (FADox and FADH°) accumulate to reasonable proportions in blue light, which is absorbed preferentially by FADox state and because the FADH° state converts relatively inefficiently to the FADH− redox state (see Figure 1 and Bouly et al., 2007; Burney et al., 2012). Therefore, the contribution of the FADH− state to the equilibrium reached by cryptochrome can be neglected and the two-state model can be applied. However, under natural conditions blue light is only one component of the ambient light environment and there is a higher proportion of UV, turquoise, green, and yellow light (all absorbed by the radical FADH°) than of purely blue light in the spectrum. Therefore, the cryptochrome photocycle is more complex than the simpler “on”—“off” two state conversion model for most photoreceptors, and the three-state kinetic model must be considered.
FADH° to FADH− Photoconversion
To complete the modeling of the cryptochrome photocycle, we therefore consider the FADH°—FADH− redox state transition (k2) and its reverse (reoxidation) reaction (k2b). To do so, we take advantage of the spectral properties of the neutral radical (FADH°) flavin redox state, which can absorb green (500—600 nm) as well as blue light. No other redox form of FAD can absorb green light. As can be seen (Figure 1), illumination with green light induces the reduction of FADH° to FADH−(k2 but not k1). Therefore, the rate constant k2 can be experimentally determined by assessing the effect of co-illumination of blue plus green light in comparison to blue light alone. Any change in cry2 absorbance induced as a result of co-illumination with green light must necessarily be due to depletion of the neutral radical form (FADH°) of cry2.
To determine the k2 rate constant for FADH° to FADH− forward light driven photoconversion from the spectral data we have illuminated purified cry2 photoreceptor at a photon fluence rate of blue light I1 = 100 μmol m−2 s−1 for 20 s (Figure 6A). This illumination induces formation of the neutral radical redox state. The samples were then co-illuminated with increasing fluence rates of green light (I2) varying from 0 to 1000 μmol m−2 s−1 (Figure 6A). Decrease in absorbance at 450 nm and 550 nm with increasing concentrations of green light results from FADH° to FADH− redox state transition.
From the spectra of Figure 6A we calculated the concentration of FADHo and FADox by using Equation (4). While the concentration of FADHo decreases with increasing green light fluence rates, FADox remains constant (1% variation with respect to the photoreduction with only blue light), meaning that the reoxidation rates can be neglected. By neglecting reoxidation (k1b = k2b = 0) we calculated, for each fluence rate I2, the rate constant k2 from the three-state algorithm. By taking as input the FADHo values and the rate constant k1 previously found (Figure 2D), this algorithm outputs k2 according to Equation (1). Figure 6B reports (red triangles) I2 as function of k2, and the linear fit k2 = σ2I2 (blue curve). The fit provided a photoconversion cross section of σ2 = 1.6 × 10−5 μmol−1 m2. As can be seen, also in this case the model (blue line) provides an excellent fit with the data (R2 = 0.96). The quantum yield for this reaction was calculated as ϕ2 = 0.027, which is almost an order of magnitude lower than that of the FADox to FADH° interconversion ϕ1 = 0.14 (Figure 2D). To confirm that reoxidation can be neglected, we also calculated k2 with the three-state algorithm by considering the reoxidation rates. For k1b we used the value previously found in the present manuscript (Figure 2B), and for k2d = 0.011 s−1, a value determined from prior publications (Müller and Ahmad, 2011). We obtained the same results, i.e., ϕ2 = 0.027 with σ2 = 1.58 × 10−5 μmol−1 m2, confirming that the depletion of FADHo is indeed due to green light. In sum, under conditions of steady state illumination, the receptor will occupy primarily FADH° and FADox redox forms, with only a minor contribution from the fully reduced (FADH−) redox state (see Supplementary Figure B).
Kinetic Modeling of the Cryptochrome Photocycle In vivo
Ultimately, our goal is to apply the kinetic model, which we used for isolated protein spectra, to predict the equilibrium redox states of cryptochromes adopted in response to illumination in vivo. However, it is not possible to obtain direct measurements of the flavin redox state in living plants. Therefore, to apply the kinetic model to cryptochrome photocycle in vivo, we make the following assumptions. Firstly, we consider that the simpler two-state model (with only FADox and FADH° flavin states) is adequate to describe the cryptochrome photocycle under continuous blue (450 nm) light (Figure 2). This assumption is based on the fact that quantum yield for k2 is 10-fold lower than for k1 (Figure 6) and therefore the concentration of FADH° should always be in large excess to that of FADH− (see Figure 1 for interconversion of redox states). Secondly, we make the assumption that biological activity is directly proportional to the concentration of the FADH° flavin state, as indicated by many studies in the literature (reviewed in Chaves et al., 2011). To convert the “readout” for biological activity to FADHo concentration we used Equation (S5) given in the Supplementary Material. Using these two assumptions, we applied the kinetic model depicted in Figure 1 to in vivo responses for both cry1 and cry2.
Modeling the Cry1 Photocycle under Blue Light In vivo
In the case of cry1, we chose blue-light dependent inhibition of hypocotyl elongation as a “readout” for biological activity (Ahmad et al., 2002). In this phenotype, seedlings are allowed to grow for several days at different fluence rates of blue light I1. The length of the hypocotyl (L) is inversely proportional to the blue light fluence rate, and in this case the cryptochrome receptor mediates shortening of the hypocotyl. We chose this property since it is proportional to the photon fluence rate, and we wished to measure a response as closely tied to primary light absorption characteristics at the photoreceptor as possible. To this end we use mutants deficient in phytochrome (phyAphyB double mutants) since phytochrome also absorbs blue light and has profound effects on hypocotyl growth, likely by acting downstream of the cryptochrome by an independent mechanism (Ahmad and Cashmore, 1997).
For the experimental procedure, Arabidopsis seedlings were grown at different fluence rates of blue light I1. Hypocotyl length (L) was plotted as a function of blue light fluence rate (Figure 7A).
Figure 7. Modeling in vivo Cry1 biological activity in blue light. (A) Double plot of (left, blue vertical axis) of normalized hypocotyl length L, and (right, red vertical axis) [FADHo], as a function of the photon fluence rate I1. The FADHo values were obtained from L by assuming that the hypocotyl length is inversely proportional to FADHo concentration, as explained in Method Section Three-State Model of Cryptochrome Photoreduction. The error bar represents the standard error of the hypocotyl growth measurements for each I1. (B) Rate constant k1 as a function of the photon fluence rate I1 (red triangles). The rate constant k1 was calculated by the two-state model algorithm, with input of the [FADHo] concentration, k1b = 0.0023 s−1 and with output the k1. The blue curve reports the fit of the data, from which a photoconversion cross section of σ1 = 2.08 × 10−5 μmol−1 m2 resulted. We calculate from σ1 (Method Section Two-State Model for Cryptochrome Photocycle In vitro), a quantum yield of ϕ1 = 0.014. The k1b value was taken from the literature (Herbel et al., 2013).
To apply the two-state kinetic model, we converted hypocotyl length (L) to in vivo FADHo concentration, which we double plot as a function of the photon fluence rate in Figure 7A. Figure 7A thus shows our assumption, i.e., that FADHo correlates with the “readout” of biological activity.
For each photon fluence rate I1 in Figure 7A, we calculated the rate constant k1, by using the two-state algorithm with input of the FADHo concentration (Figure 7A) and k1b. The dark reoxidation reactions in vivo k1b was taken from previously obtained experimental values (Herbel et al., 2013). The output value was the k1, as obtained for cry1 in vitro (Figure 4). Figure 7B reports k1 as a function of I1. The blue curve in Figure 7B reports the liner fit k1 = σ1 I1, and resulted in a photoconversion cross section of σ1 = 2.08 × 10−5 μmol−1 m2. We calculated the quantum yield using the extinction coefficient obtained for cry1 in vitro. The quantum yield resulted ϕ1 = 0.014, which is in excellent agreement with values obtained from the in vitro studies (Figure 4).
Modeling the Cry2 Photocycle under Blue Light In vivo
To model the cry2 photocycle in vivo, we have used light-dependent degradation of cry2 as a “readout” for biological function. In the dark, cry2 protein accumulates to high levels in seedlings. However, upon transfer to light, the protein is rapidly degraded within 30 min of the start of illumination (Banerjee et al., 2007; Bouly et al., 2007). This effect relies on conformational change subsequent to light absorption followed by ubiquitination of the receptor, and is therefore directly linked to cry2 primary activation by light and formation of the signaling state. Experiments were performed in phytochrome-deficient mutants (phyAphyB double mutants), to avoid potential effects of phytochrome on cry2 degradation. We therefore consider cry2 protein concentration after illumination as a measure for accumulation of the cry2 signaling state and of FADH° accumulation.
In our experiments, we first irradiated dark-grown seedlings at different blue light fluence rates for 30 min. Cry2 protein expression levels were evaluated on Western blots (see Supplementary Figure C) and the signal quantitated by imaging software ImageJ. In this way the protein concentration values C were plotted as a function of the photon fluence rate (Figure 8A).
Figure 8. Modeling in vivo Cry2 biological activity in blue light. (A) Double plot of (left, blue vertical axis) normalized protein concentration C, and (right, red vertical axis) [FADHo], as a function of photon fluence rate I1. The FADHo values were obtained from C by assuming that FADHo concentration is inversely proportional to C, as explained in Method Section Three-State Model of Cryptochrome Photoreduction. The error bar represents the standard error of the measurements for each I1. (B) Rate constant k1 as a function of the photon fluence rate I1 (red triangles). The rate constant k1 was calculated by the two-state model algorithm, which inputs the [FADHo] concentration from panel (A), k1b = 7.2 × 10−4 s−1, and outputs k1. The blue curve reports the linear fit of the data (k1 = σ1 I1), from which a photoconversion cross section of σ1 = 4.8 × 10−4 μmol−1 m2 resulted. From σ1 we calculated (as explained in Method Section Two-State Model for Cryptochrome Photocycle In vitro), a quantum yield of ϕ1 = 0.41, in close agreement with the calculated quantum yield of cry2 in vitro (Figure 2D and Table 1). The k1b value was taken from the literature (Herbel et al., 2013).
Also in this case, we converted the “readout” of biological activity, i.e., the cry2 protein concentration C, to FADH° concentration, which we double plot in Figure 8A. We used the FADHo concentration to calculate the photoconversion cross section σ1 by plotting I1 vs. the calculated k1 (Figure 8B). The value k1d for the dark reoxidation reaction of cry2 in vivo was already previously obtained experimentally (Herbel et al., 2013). The quantum yield for photoconversion of cry2 was ϕ1 = 0.41, again in agreement with the values obtained from in vitro studies (Figure 2).
Three-State Model for the Cry2 Photocycle In vivo
Finally, we provide a more comprehensive model of the cry2 photocycle in vivo by taking into consideration the third redox state, FADH−, which is induced by green light and which had been successfully modeled in vitro (Figure 6). We again use the cry2 protein degradation assay as a measure for biological activity (see above), but illuminating with green light (560 nm) in order to induce FADH°—FADH− photoconversion. In this way, the effect of the three redox states on cryptochrome photoconversion could be modeled and compared to the biological activity.
For these experiments, all seedlings were illuminated at a sub-saturating fluence rate of blue light (10 μmol m−2s−1) for the duration of the light treatments (30 min). Seedlings were in addition co-illuminated with increasing fluence rates of green light I2. The levels of cry2 protein were analyzed from Western blot images (see Supplementary Figure C for gel image). The cry2 protein concentration was converted in cry2 FADH° concentration and normalized (Figure 9A; increasing concentration of cry2 shows decreased biological activity). For each green light fluence rate I2, we then calculated the forward rate k2, by using the three-state model algorithm. We input the normalized FADHo concentration as a function of the photon fluence rate (from Figure 9A), k1, k1b, and k2b. The k1 value was as determined in this study (Figure 8), while the k1b, k2b values were provided from previous studies (i.e., k1b = 7.2 × 10−4 s−1 (Herbel et al., 2013) and k2b = 0.011s−1 (Müller and Ahmad, 2011). The output k2 was obtained according to Equation (1). Figure 9B reports I2 as a function of k2 (red triangles), and the liner fit of the data, k2 = σ2 I2 (blue curve). From the linear fit, a photoconversion cross section of σ2 = 3 × 10−5 μmol−1 m2 was obtained. A quantum yield of ϕ2 = 0.05 could be calculated by using the extinction coefficient determined for cry2 in vitro (ε2 = 254.7 mol−1 m2, i.e., 2547 M−1 cm−1). Once again, this value in vivo is in close agreement with the quantum yield obtained for cry2 (ϕ2) in vitro (Figure 6B and Table 1).
Figure 9. Three-state Modeling of in vivo Cry2 biological activity. (A) Double plot of (left, green vertical axis) normalized protein concentration C, and (right, red vertical axis) [FADHo], as a function of the green light fluence rate I2. Seedlings were co-illuminated at a fluence rate of blue light I1 = 10 μmol m−2s−1 and increasing fluence rates of green light (I2) as shown in the horizontal axis, for 30 min. The FADHo values were obtained from C by assuming that FADHo concentration is inversely proportional to protein concentration, as explained in the Method Section Three-State Model of Cryptochrome Photoreduction. The error bar represents the standard error of the measurements. (B) Rate constant k2 as a function of I2 (red triangles). The rate constant k2 was calculated by the three-state model algorithm, which input the [FADHo] concentration from panel (A), k2b = 0.011 s−1 (Herbel et al., 2013), k2b = 7.2 × 10−4 s−1 (Müller and Ahmad, 2011), and k1 from Figure 8B. The output is the k2. The blue curve reports the linear fit of the data (k2 = σ2 I2), from which a photoconversion cross section of σ2 = 3 × 10−5 μmol−1 m2resulted. From σ2 we calculated (as explained in Methods Section Two-State Model for Cryptochrome Photocycle In vitro), a quantum yield of ϕ2 = 0.05, in close agreement with the calculated quantum yield of cry2 in vitro (Figure 4D).
This study represents a first attempt to model the cryptochrome photocycle and determine kinetic parameters of relevance to biological function in vivo. The method adopted here seeks only to model interconversion of the redox states of cryptochromes (FADox, FADH°, and FADH−). It is a much simplified method as compared to global analysis, which has been traditionally used primarily to identify absorbing species from spectra and find reoxidation rates (Müller and Ahmad, 2011). The two-state model studied here is a simple integrable model and the method used here is less computationally expensive than global analysis, and has a similar level of accuracy. We relate absorbance and concentration of two absorbing species (FADox and FADH°) by using the Beer-Lambert law, which allows a more accurate fit of the two-state kinetic model to spectra than have been obtained in prior studies (Burney et al., 2009). Most importantly, our method can be readily adapted to estimate quantum yields in vivo, which is not the case for global analysis.
Quantum yields have been traditionally obtained at a given photon fluence rate by using different evaluation techniques that were at the same time correlated with receptor photoconversion. Pr/Pfr photoconversion has been followed in phytochromes and related to function in this way, for example (Hermann et al., 1985; Kelly and Lagarias, 1985; Mancinelli, 1988). Quantum yield of Arabidopsis cry1 has been calculated in this way for the efficiency of primary electron transfer reactions (Giovani et al., 2003; Müller et al., 2014), which were not, however, related to biological function.
Here, by contrast, we exploit the linear correlation between photon fluence rate and forward rate constant to find the photoconversion cross sections, which we then use to calculate the quantum yields. This method allows us to experimentally find a linear range of intensities where prediction of the two or three states intermediate concentrations can be made.
We have estimated the quantum yield of both cry1 and cry2 in vitro from spectra of flavin photoreduction at a range of light intensities. Increasing the buffer pH from 7.5 to 8.2 resulted in a 10-fold decrease in quantum yield of flavin reduction in cry2 derived from the kinetic model, consistent with prior studies showing reduced amplitude of primary electron transfer from flavin under such conditions (Müller et al., 2014). The rate of reoxidation of the FADH° redox state for both cry1 and cry2 was similar under all tested conditions (Figures 2–5), consistent with the dependence of this parameter on the concentration of molecular oxygen (Müller and Ahmad, 2011). In sum, our simple kinetic model accurately describes the redox state transitions of cryptochromes in vitro in a manner consistent with effects of known modulating factors.
To model the kinetics of the cryptochrome photocycle in vivo, we have made the assumption that biological activity can be used as a measure of FADH° concentration. This follows from numerous studies on the cryptochrome photocycle that have correlated FADH° formation and flavin reduction with biological signaling and activity in vivo (Chaves et al., 2011). These include recent studies showing decrease in blue light-dependent biological activity in cry2 (protein degradation assay) in mutants that impair light-dependent radical formation in vitro (Li et al., 2011—see also Engelhard et al., 2014). Furthermore, a dose-response curve of cry1 biological activity (hypocotyl growth inhibition) showed decrease in light sensitivity of several orders of magnitude in cry1 photoreduction mutants as compared to the appropriate (expressing similar concentrations of wild type cry1) control seedlings (see supplement in Gao et al., 2015). For the rate of reoxidation from FADH° to FADox in our model, we used values derived from in vivo studies of the cry1 and cry2 flavin state lifetimes (Herbel et al., 2013). Fitting the kinetic model to biological data, we were able to calculate quantum yields for the biological response in vivo which were somewhat similar to the corresponding values calculated from our spectral analysis in vitro (see Table 1—compare for example ϕ1 values measured for cry1 and cry2 of 0.038 and 0.19 in vitro at pH7.5 to ϕ1 values obtained in vivo of 0.014 and 0.4, respectively).
This similarity (within 2.5-fold) was not a required result for our model to have validity, as many external variables may affect photoreceptor responsivity in vivo. For example, the rate constants of the redox state interconversion events (k1, k2, k1b, and k2b) could vary greatly in vivo through tuning by the cellular environment (Engelhard et al., 2014). Such variability due to cellular environment was also pointed out for the case of the phytochrome (Mancinelli, 1988) where direct measurements of the active state (Pfr) concentration is also not possible in vivo (but see Rausenberger et al., 2010). On a more general note, predicted values of the state of photoreceptors calculated from in vitro photochemical parameters and the spectral photon flux distribution may vary from the actual ones found in vivo as a result of light scattering in whole tissues and shielding by other plant pigments (Mancinelli, 1988). A further problem is the phenomenon of signal amplification through secondary pathways. Indeed, it is possible to derive quantum efficiencies for photoreceptor light sensing in vivo that are >1, if the biological readout selected for analysis is subject to significant signal amplification. We address these issues by evaluating phenotypes in etiolated seedlings with a minimum of cell layers and accessory pigments, and also selecting a genetic background (phyAphyB double mutant) providing a minimum of signal amplification.
Despite these potential limitations, the fact that quantum yield for biological activity matches so closely with the calculated values for cryptochrome flavin reduction in vitro supports our assumption that flavin redox state determines biological activity. This is particularly striking in the case of k2 for cry2 (response to bichromatic green/blue illumination), for which no other explanation can be reasonably given. Furthermore, qualitative effects in vivo such as the relative efficiencies in vivo (cry1 as compared to cry2—see below; response to blue as compared to biochromatic blue plus green light) could be clearly determined by this approach.
The quantum yield of cry2 photoconversion (see results summarized in Table 1) is within the range of other photoreceptors such as phytochrome A, which has a photoconversion quantum yield of the order of 0.14 (Gensch et al., 1996) or of the LOV2 domain in phototropins of 0.2 (Kasahara et al., 2002). Cry1 is within the range of quantum yield calculated for LOV1 of 0.026 (Kasahara et al., 2002), which responds at higher blue light intensities. In terms of quantum yields of sensory receptors, it should be considered that cry1 in particular regulates growth processes that occur in full sunlight during de-etiolation and vegetative growth of plants. Therefore, it need not have high quantum yield such as phototropins and phytochrome A, which are specialized for responses at extremely low light intensities (Smith, 1995). To the contrary, too high a photon sensitivity would eliminate the ability to respond to higher intensities of light, as the receptor would reach saturation too quickly. Along these lines the quantum yield of cry2 is around 10-fold higher than cry1, consistent with published effects of cry2 at lower light intensities than for cry1 (Lin et al., 1998).
One of the unique characteristic of the cryptochrome photocycle is the fact that it exists in three states rather than solely as a two states system. The in vitro obtained quantum yield for cry2 photoconversion were 0.188 (FADox to FADH°; Figure 2) and 0.027 (FADH° to FADH−; Figure 6). This means that under continuous illumination in monochromatic blue light, essentially only the radical (FADH°) redox form should accumulate in response to illumination. This is in distinction to photolyases which undergo full reduction under continuous illumination (Burney et al., 2012). However, under natural conditions of full sunlight, there is a significant contribution from other wavelengths of light, including green light. In particular, in the case of shading under plant canopies, the ratio of green to blue light can be quite high as green light passes through shading leaves whereas blue light is absorbed. Under these conditions, the biological activity of cryptochrome could be considerably modulated by the second (FADH° to FADH−) redox state transition, as has indeed been suggested in prior studies (Bouly et al., 2007).
Variations of the three-state photocycle appear to be of general relevance to cryptochrome activation. For instance, in the case of algal cryptochromes (Beel et al., 2012) it appears that the FADH° redox state is the (dark) resting state and the FADH− redox state is the signaling state. Therefore, because the FADH° redox state also absorbs green and red light, algal cryptochromes are activated by UV, blue, green, and red light rather than just blue light such as for plant cryptochromes. Conversely, the drosophila cry appears to adopt just a two-state (FADox to FAD°−) photocycle, with FAD°− as the signaling state (Berndt et al., 2007). Intriguingly, the avian cry1a, which is thought to be implicated in sensing of the geomagnetic field (Ritz et al., 2000), also appears to adopt a three-state photocycle wherein the FADH− but not the FADH° redox state is the signaling state (Nießner et al., 2013). The kinetic model described in this work should therefore be applicable to all of these cryptochrome photocycles.
In sum, the plant cryptochrome photocycle provides the basis for a rapid and versatile response to the light environment. The receptor can respond within a matter of minutes to changes in the light intensity, even at very high light, by a simple shift in the equilibrium concentration of the FADH° redox state. Because cry1 and cry2 have overlapping functions but the quantum yield of cry2 is 10-fold higher, the combined action of cry1 and cry2 is capable of sensing and responding to light intensity in an almost linear fashion over two orders of magnitude. In this way, the cryptochromes appear to divide the task of responding to the light environment along the lines of what is seen for phytochromes, where the most abundant but photolabile phytochrome (phyA) is specialized for response to dim light whereas stable but less abundant phytochromes (phy B,C,D,E) respond to light at higher intensities (Smith, 1995). Finally, the cryptochromes also have the unique feature of a three-state photocycle, which provides for differential responsivity depending on the wavelength composition. Our model should therefore be useful in unraveling the details of the cryptochrome photocycle in any biological system.
MA: Designed and performed experiments, wrote paper. JL: Performed experiments. MP: Performed modeling, wrote paper. TR: Assisted modeling. JW: Performed experiments.
Funding was from AFOSR (FA9550-14-0-0409).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
We are indebted to Alain d'Harlingue and Pierre-Etienne Bouchet for help with plant growth and technical assistance and to Austin Glascoe for help with recombinant protein production.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/article/10.3389/fpls.2016.00888
Ahmad, M., and Cashmore, A. R. (1997). The blue light receptor cryptochrome 1 shows functional dependence on phytochrome A or phytochrome B in Arabidopsis thaliana. Plant J. 11, 421–427. doi: 10.1046/j.1365-313X.1997.11030421.x
Ahmad, M., Grancher, N., Heil, M., Black, R. C., Giovani, B., Galland, P., et al. (2002). Action spectrum for cryptochromedependent hypocotyl growth inhibition in Arabidopsis. Plant Physiol. 129, 774–785. doi: 10.1104/pp.010969
Banerjee, R., Schleicher, E., Meier, S., Muñoz Viana, R., Pokorny, R., Ahmad, M., et al. (2007). The signaling state of Arabidopsis cryptochrome 2 contains flavin semiquinone. J. Biol. Chem. 282, 14916–14922. doi: 10.1074/jbc.M700616200
Beel, B., Prager, K., Spexard, M., Sasso, S., Weiss, D., Müller, N., et al. (2012). A flavin binding cryptochrome photoreceptor responds to both blue and red light in Chlamydomonas reinhardtii. Plant Cell 24, 2992–3008. doi: 10.1105/tpc.112.098947
Berndt, A., Kottke, T., Breitkreuz, H., Dvorsky, R., Hennig, S., Alexander, M., et al. (2007). A novel photoreaction mechanism for the circadian blue light photoreceptor Drosophila cryptochrome. J. Biol. Chem. 282, 13011–13021. doi: 10.1074/jbc.M608872200
Bouly, J.-P., Schleicher, E., Dionisio-Sese, M., Vandenbussche, F., Van Der Straeten, D., Bakrim, N., et al. (2007). Cryptochrome bluelight photoreceptors are activated through interconversion of flavinredox states. J. Biol. Chem. 282, 9383–9391. doi: 10.1074/jbc.M609842200
Burney, S., Hoang, N., Caruso, M., Dudkin, E. A., Ahmad, M., and Bouly, J. P. (2009). Conformational change induced by ATP binding correlates with enhanced biological function of Arabidopsis cryptochrome. FEBS Lett. 583, 1427–1433. doi: 10.1016/j.febslet.2009.03.040
Burney, S., Wenzel, R., Kottke, T., Roussel, T., Hoang, N., Bouly, J. P., et al. (2012). Single amino acid substitution reveals latent photolyase activity in Arabidopsis cry1. Angew. Chem. Int. Ed. Engl. 51, 9356–9360. doi: 10.1002/anie.201203476
Chaves, I., Pokorny, R., Byrdin, M., Hoang, N., Ritz, T., Brettel, K., et al. (2011). The cryptochromes: blue light photoreceptors in plants and animals. Annu. Rev. Plant Biol. 62, 335–364. doi: 10.1146/annurev-arplant-042110-103759
El-Esawi, M., Glascoe, A., Engle, D., Ritz, T., Link, J., and Ahmad, M. (2015). Cellular metabolites modulate in vivo signaling of Arabidopsis cryptochrome-1. Plant Signal Behav. 10:e1063758. doi: 10.1080/15592324.2015.1063758
Engelhard, C., Wang, X., Robles, D., Moldt, J., Essen, L. O., Batschauer, A., et al. (2014). Cellular metabolites enhance light sensitivity through alternate electron transfer pathways in Arabidopsis cryptochrome. Plant Cell 26, 4519–4531. doi: 10.1105/tpc.114.129809
Gao, J., Wang, X., Zhang, M., Bian, M., Deng, W., Zuo, Z., et al. (2015). Trp triad-dependent rapid photoreduction is not required for the function of Arabidopsis CRY1. Proc. Natl. Acad. Sci. U.S.A. 112, 9135–9140. doi: 10.1073/pnas.1504404112
Gensch, T., Churio, S. M., Braslavsky, S. E., and Schaffner, K. (1996). Primary quantum yield and volume change of phytochrome-a phototransformation determined by laser-induced optoacoustic spectroscopy. Photochem. Photobiol. 63, 719–725. doi: 10.1111/j.1751-1097.1996.tb09621.x
Hermann, G., Appenroth, K. J., Müller, E., and Augsten, H. (1985). Reexamination of the partial quantum yields for the phytochrome photo conversion by adequate photokinetic methods. Biochem. Physiol. Pflanzen 180, 533–542. doi: 10.1016/S0015-3796(85)80101-3
Kasahara, M., Swartz, T. E., Olney, M. A., Onodera, A., Mochizuki, N., Fukuzawa, H., et al. (2002). Photochemical properties of the flavin mononucleotide-binding domains of the phototropins from Arabidopsis, Rice, and Chlamydomonas reinhardtii. Plant Physiol. 129, 762–773. doi: 10.1104/pp.002410
Kondoh, M., Shiraishi, C., Müller, P., Ahmad, M., Hitomi, K., Getzoff, E. D., et al. (2011). Light-induced conformational changes in full-length Arabidopsis thaliana cryptochrome. J. Mol. Biol. 413, 128–137. doi: 10.1016/j.jmb.2011.08.031
Li, X., Wang, Q., Yu, X., Liu, H., Yang, H., Zhao, C., et al. (2011). Arabidopsis cryptochrome 2 (CRY2) functions by the photoactivation mechanism distinct from the tryptophan (trp) triad-dependent photoreduction. Proc. Natl. Acad. Sci. U.S.A. 108, 20844–20849. doi: 10.1073/pnas.1114579108
Lin, C., Yang, H., Guo, H., Mockler, T., Chen, J., and Cashmore, A. R. (1998). Enhancement of blue-light sensitivity of Arabidopsis seedlings by a blue light receptor cryptochrome 2. Proc. Natl. Acad. Sci. U.S.A. 95, 2686–2890. doi: 10.1073/pnas.95.5.2686
Liu, B., Zuo, Z., Liu, H., Liu, X., and Lin, C. (2011). Arabidopsis cryptochrome 1 interacts with SPA1 to suppress COP1 activity in response to blue light. Genes Dev. 25, 1029–1034. doi: 10.1101/gad.2025011
Liu, H., Yu, X., Li, K., Klejnot, J., Yang, H., Lisiero, D., et al. (2008). Photoexcited CRY2 interacts with CIB1 to regulate transcription and floral initiation in Arabidopsis. Science 322, 1535–1539. doi: 10.1126/science.1163927
Mancinelli, A. L. (1988). Some thoughts about the use of predicted values of the state of phytochrome in plant photomorphogenesis research. Plant Cell Environ. 11, 429–439. doi: 10.1111/j.1365-3040.1988.tb01780.x
Müller, P., and Ahmad, M. (2011). Light-activated cryptochrome reacts with molecular oxygen to form a flavin-superoxide radical pair consistent with magnetoreception. J. Biol. Chem. 286, 21033–21040. doi: 10.1074/jbc.M111.228940
Müller, P., Bouly, J. P., Balland, K., Hitomi, V., Getzoff, E., Ritz, T., et al. (2014). ATP binding turns plant cryptochrome into an efficient natural photoswitch. Sci. Rep. 4:5175. doi: 10.1038/srep05175
Nießner, C., Denzau, S., Stapput, K., Ahmad, M., Peichl, L., Wiltschko, W., et al. (2013). Activated Cryptochrome 1a concurs with magnetic orientation in birds. J. R. Soc. Interface 10:20130638. doi: 10.1098/rsif.2013.0638
Rausenberger, J., Hussong, A., Kircher, S., Kirchenbauer, D., Timmer, J., Nagy, F., et al. (2010). An integrative model for phytochrome B mediated photomorphogenesis: from protein dynamics to physiology. PLoS ONE 5:e10721. doi: 10.1371/annotation/4563eaf4-e45b-4d9e-ab06-5f1794bf11e3
Valverde, F., Mouradov, A., Soppe, W., Ravenscroft, D., Samach, A., and Coupland, G. (2004). Photoreceptor regulation of CONSTANS protein in photoperiodic flowering. Science 303, 1003–1006. doi: 10.1126/science.1091761
Yu, X., Klejnot, J., Zhao, X., Shalitin, D., Maymon, M., Yang, H., et al. (2007). Arabidopsis cryptochrome 2 completes its posttranslational life cycle in the nucleus. Plant Cell 19, 3146–3156. doi: 10.1105/tpc.107.053017
Yu, X., Sayegh, R., Maymon, M., Warpeha, K., Klejnot, J., Yang, H., et al. (2009). Formation of nuclear bodies of Arabidopsis CRY2 in response to blue light is associated with its blue light-dependent degradation. Plant Cell 21, 118–130. doi: 10.1105/tpc.108.061663
Keywords: cryptochrome, flavoprotein, kinetic modeling, signaling, photoreduction
Citation: Procopio M, Link J, Engle D, Witczak J, Ritz T and Ahmad M (2016) Kinetic Modeling of the Arabidopsis Cryptochrome Photocycle: FADHo Accumulation Correlates with Biological Activity. Front. Plant Sci. 7:888. doi: 10.3389/fpls.2016.00888
Received: 25 March 2016; Accepted: 06 June 2016;
Published: 28 June 2016.
Edited by:Enamul Huq, The University of Texas at Austin, USA
Copyright © 2016 Procopio, Link, Engle, Witczak, Ritz and Ahmad. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Margaret Ahmad, firstname.lastname@example.org