The primary structural photoresponse of phytochrome proteins captured by a femtosecond X-ray laser

  1. Elin Claesson
  2. Weixiao Yuan Wahlgren
  3. Heikki Takala
  4. Suraj Pandey
  5. Leticia Castillon
  6. Valentyna Kuznetsova
  7. Léocadie Henry
  8. Matthijs Panman
  9. Melissa Carrillo
  10. Joachim Kübel
  11. Rahul Nanekar
  12. Linnéa Isaksson
  13. Amke Nimmrich
  14. Andrea Cellini
  15. Dmitry Morozov
  16. Michał Maj
  17. Moona Kurttila
  18. Robert Bosman
  19. Eriko Nango
  20. Rie Tanaka
  21. Tomoyuki Tanaka
  22. Luo Fangjia
  23. So Iwata
  24. Shigeki Owada
  25. Keith Moffat
  26. Gerrit Groenhof
  27. Emina A Stojković
  28. Janne A Ihalainen
  29. Marius Schmidt  Is a corresponding author
  30. Sebastian Westenhoff  Is a corresponding author
  1. Department of Chemistry and Molecular Biology, University of Gothenburg, Sweden
  2. Department of Biological and Environmental Science, Nanoscience Center, University of Jyvaskyla, Finland
  3. Department of Anatomy, Faculty of Medicine, University of Helsinki, Finland
  4. Physics Department, University of Wisconsin-Milwaukee, United States
  5. Department of Biology, Northeastern Illinois University, United States
  6. Department of Chemistry, Nanoscience Center, University of Jyvaskyla, Finland
  7. Department of Cell Biology, Graduate School of Medicine, Kyoto University, Japan
  8. RIKEN SPring-8 Center, Japan
  9. Japan Synchrotron Radiation Research Institute, Japan
  10. Department of Biochemistry and Molecular Biology and Institute for Biophysical Dynamics, University of Chicago, United States

Abstract

Phytochrome proteins control the growth, reproduction, and photosynthesis of plants, fungi, and bacteria. Light is detected by a bilin cofactor, but it remains elusive how this leads to activation of the protein through structural changes. We present serial femtosecond X-ray crystallographic data of the chromophore-binding domains of a bacterial phytochrome at delay times of 1 ps and 10 ps after photoexcitation. The data reveal a twist of the D-ring, which leads to partial detachment of the chromophore from the protein. Unexpectedly, the conserved so-called pyrrole water is photodissociated from the chromophore, concomitant with movement of the A-ring and a key signaling aspartate. The changes are wired together by ultrafast backbone and water movements around the chromophore, channeling them into signal transduction towards the output domains. We suggest that the observed collective changes are important for the phytochrome photoresponse, explaining the earliest steps of how plants, fungi and bacteria sense red light.

eLife digest

Plants adapt to the availability of light throughout their lives because it regulates so many aspects of their growth and reproduction. To detect the level of light, plant cells use proteins called phytochromes, which are also found in some bacteria and fungi. Phytochrome proteins change shape when they are exposed to red light, and this change alters the behaviour of the cell. The red light is absorbed by a molecule known as chromophore, which is connected to a region of the phytochrome called the PHY-tongue. This region undergoes one of the key structural changes that occur when the phytochrome protein absorbs light, turning from a flat sheet into a helix.

Claesson, Wahlgren, Takala et al. studied the structure of a bacterial phytochrome protein almost immediately after shining a very brief flash of red light using a laser. The experiments revealed that the structure of the protein begins to change within a trillionth of a second: specifically, the chromophore twists, which disrupts its attachment to the protein, freeing the protein to change shape. Claesson, Wahlgren, Takala et al. note that this structure is likely a very short-lived intermediate state, which however triggers more changes in the overall shape change of the protein.

One feature of the rearrangement is the disappearance of a particular water molecule. This molecule can be found at the core of many different phytochrome structures and interacts with several parts of the chromophore and the phytochrome protein. It is unclear why the water molecule is lost, but given how quickly this happens after the red light is applied it is likely that this disappearance is an integral part of the reshaping process.

Together these events disrupt the interactions between the chromophore and the PHY-tongue, enabling the PHY-tongue to change shape and alter the structure of the phytochrome protein. Understanding and controlling this process could allow scientists to alter growth patterns in plants, such as crops or weeds.

Introduction

Phytochrome photosensor proteins are crucial for the optimal development of all vegetation on Earth (Butler et al., 1959; Gan et al., 2014; Quail et al., 1995). Prototypical phytochromes can exist in two photochemical states with differential cellular signaling activity, called red light-absorbing (Pr) and far-red light-absorbing (Pfr) state (Figure 1—figure supplement 1). As a result, phytochromes can distinguish two colors of light, providing plants, fungi, and bacteria with primitive two-color vision. Light is detected by a bilin chromophore, which is covalently linked to the photosensory core of the protein (Wagner et al., 2005), comprising of PAS (Per/Arndt/Sim), GAF (cGMP phosphodiesterase/adenyl cyclase/FhlA) and PHY (phytochrome-specific) domains. Two propionate side chains additionally anchor the chromophore non-covalently to the protein (Figure 1b). The signaling sites of the phytochrome are found in its C- and N-terminal output domains, which vary between species. Important for the signaling is a stretch of amino acids in the PHY domain, called the PHY-tongue, which changes from a β-sheet in Pr into an α-helix in Pfr state (Essen et al., 2008; Yang et al., 2008; Takala et al., 2014Sanchez et al., 2019) The chromophore connects to the PHY-tongue via a strictly conserved aspartatic acid, which is expected to play a crucial role in signal transduction.

Figure 1 with 5 supplements see all
Photoinduced observed difference electron density features are focused on the chromophore binding pocket.

(a) The observed difference electron density map at 1 ps is displayed together with the DrBphPdark structure. Red and green electron density peaks, contoured at 4.5 σ, denote negative and positive densities, respectively. Monomer A is colored blue and monomer B is in aqua. (b) Schematic illustration of the biliverdin chromophore. The hydrogen-bonding networks between the propionate groups and the protein are marked with dashed lines. In DrBphP, the chromophore is covalently linked to a cysteine residue in the PAS domain (solid line).

Key to phytochrome function is the primary photoresponse on picosecond time scales. Here, light signals are translated into conformational changes. The changes arise in the electronically excited bilin, but must then be transduced to the surrounding protein residues. This prepares the protein for a formation of the first intermediate (Lumi-R for prototypical phytochromes), in which isomerization of the D-ring has likely occurred (Rüdiger et al., 1983; Dasgupta et al., 2009; Yang et al., 2012; Rockwell et al., 2009; Ihalainen et al., 2018). The mechanism that leads to the first intermediate is currently not well understood, because crystallographic observations of phytochromes directly after photoexcitation have not been available.

Results

To address this gap of knowledge, we recorded time-resolved serial femtosecond X-ray crystallographic (SFX) data of the PAS-GAF domains of the phytochrome from Deinococcus radiodurans (DrBphPCBD) at 1 ps and 10 ps after femtosecond optical excitation. The experiments were performed in Japan, using the SPring-8 Angstrom Compact Free Electron Laser (SACLA) tuned to 7 KeV (Tono et al., 2015). For homogeneous excitation of the crystals, we photoexcited micrometer-sized crystals in a grease jet with a photon density of 1.7 mJ/mm2 (1/e2 measure, see Materials and methods) into the flank of the absorption peak at 640 nm (Figure 1—figure supplement 1). Taking into account the significant light scattering in the grease-buffer mixture (Figure 1—figure supplement 1), we estimate that the average number of photons per chromophore is 0.5–1 (see Materials and methods). We recorded the SFX data at 1 ps for several excitation fluences (Figure 2). Lowering the excitation density tenfold from 1.7 mJ/mm2 photons to 0.2 mJ/mm2 resulted in a joint reduction of all difference signals. Critical signals, like the twist of the D-ring and the photodissociation of the pyrrole water from the chromophore sustained when lowering the excitation densities, indicating that the signal arises predominately from one-photon excitation. The refined structure in dark (DrBphPdark), 2.07 Å resolution (Table 1), was very similar to our previous dark structure solved by SFX (5K5B, RMSD 0.646 Å and 0.610 Å for monomers A and B) (Edlund et al., 2016), but the present crystals contained two monomers in the asymmetric unit (Figure 1a, Figure 1—figure supplement 2, Figure 1—figure supplement 3). The refined 1 ps structure was solved to 2.21 Å (Table 1).

SFX data as a function of excitation fluence.

The DrBphPdarkstructure (green) is shown together with the observed difference electron density, contoured at 3.5 σ, at 1 ps collected using (a) 0.2 mJ/mm2, (b) 0.4 mJ/mm2, (c) 1.3 mJ/mm2, and (d) 1.7 mJ/mm2. All spot sizes were computed assuming Gaussian line shapes with the (1/e2) convention. The data shown in panel A-C were collected at SACLA in May 2019, whereas the data shown in panel D was collected in October 2018. The same experimental setup was used in both occasions. The laser energy of the experiment in 2018 can be found in the Materials and methods section. The energies for the experiment in 2019 were 16 µJ, 42 µJ, and 106 µJ (panels A-C, respectively). During the experiment in 2019, the femtosecond laser beam was misaligned by 50 µm distance from the interaction spot between X-rays and jet in the direction of flow. The laser intensities were corrected for this displacement assuming a Gaussian line shape. The excitation fluence is similar to previous femtosecond time-resolved SFX experiments (Nogly et al., 2018; Pande et al., 2016; Barends et al., 2015; Coquelle et al., 2018); however, we found high scattering in the grease/buffer mixture (Figure 1—figure supplement 1). Since the crystallographic signals were reduced when lowering the excitation fluence and disappeared completely when reaching 1/10 of the maximum value, we conclude that the excitation fluence that actually reaches the crystals in the grease matrix is much lower than the incoming photon fluence and that the photoexcitation is in the single-photon regime.

Table 1
Crystallographic table.
Darkone ps10 ps
PDB code6T3L6T3U
Data collection
Temperature (K)293293293
Space GroupP212121P212121P212121
Cell dimensions (a, b, c)
a, b, c (Å)54.98 116.69 117.8654.98 116.69 117.8654.98 116.69 117.86
α, β, γ ()90.0 90.0 90.090.0 90.0 90.090.0 90.0 90.0
Data resolution overall (Å)‡45.77–2.0741.46–2.2145.77–2.14
(2.10–2.07)(2.25–2.21)(2.17–2.14)
Rsplit (%)5.79 (120.05)10.59 (114.64)5.70 (121.86)
SNR (I/σ(I))9.21 (0.83)6.10 (0.88)10.11 (0.99)
CC(1/2)0.99 (0.33)0.98 (0.38)0.99 (0.344)
Completeness (%)‡100 (100)100 (100)100 (100)
Multiplicity‡461.35 (65.9)106.11 (34.3)347.36 (62.1)
Number of hits14907442853159997
Number of indexed patterns707262115070335
Indexing rate(%)&47.4449.3543.96
Number of total reflections24017763531017917823530
Number of unique reflections520603931643279
Refinement
Resolution (Å)45.82–2.0736.94–2.21
(2.12–2.07)(2.27–2.21)
Rwork / Rfree0.162/0.1910.230/0.256
(0.317/0.346)(0.411/0.443)
Number of atoms51235135
Average B factor (Å2)76.4478.63
R.m.s deviations
Bond lengths (Å)0.0070.006
Bond angles ()1.2511.152
  1. † Rsplit = 1/2Σhkl|Ieven-Iodd|1/2Σhkl|Ieven+Iodd| .

    ‡ Highest resolution shell is shown in parentheses. 

  2. § Ratio of the number of indexed images to the total number of hits.

From the time-resolved data, we calculated Fourier difference electron density maps (|Fo|light-|Fo|dark), which report on the change of structure due to optical excitation (see Materials and methods). Briefly, the diffraction data for light and dark were scaled to each other and subtracted, assuming preservation of the phases (see Materials and methods for details). The map at 1 ps indicates many significant changes in difference electron density (Figure 1a) above the background level of 3.0 standard deviations (σ) (Figure 1—figure supplement 4). The changes cluster around the chromophore, with the strongest negative densities for the pyrrole water (monomer A: −8.2σ, B: −9.4σ, Table 2). The map at 10 ps contains similar significant features, but at weaker overall intensity (pyrrole water A: −5.0σ, B:−6.5σ) (Figure 1—figure supplement 4 and Figure 1—figure supplement 5). We ascribe this to a lower population of the activated state at 10 ps compared to 1 ps. Monomer A has a lower signal strength than monomer B, but provided a clearer difference map around the chromophore. We refined a structural model (DrBphP1ps) using extrapolated structure factors (Figure 3—figure supplement 1; Pande et al., 2016). The refinement of the structure against the 1 ps data was successful using a photoexcitation density of 8%. However, we aborted our attempts to refine a structural model against the 10 ps data, as the model would have become unreliable due to an even lower photoactivation yield. We focus our discussion on monomer A and the 1 ps time point, although all conclusions are supported by monomer B and the features observed in the difference maps at 10 ps (Table 2, Figure 1—figure supplement 5).

Table 2
Difference electron density features listed for certain atoms.
ObjectLabel1ps1ps10ps10ps
ABAB
Pyrrole Water
Pyrrole Water (-)XV−8.2σ−9.4σ−5.0σ−6.5σ
Pyrrole Water (+) Alt. 1XVI4.8σ4.2σ3.0σ3.4σ
Pyrrole Water (+) Alt. 2XVII4.4σ5.8σ3.0σ3.3σ
D-ring
D-ring N-H/C = O (-)I−4.2σ−7.3σ−6.4σ−5.1σ
D-ring Methyl (-)II−4.3σ−3.5σ−2.9σ−4.9σ
D-ring Vinyl (-)III−4.0σ−4.8σ−3.5σ−3.4σ
D-ring N-H/C = O (+)IV4.6σ6.4σ4.6σ4.1σ
D-ring Methyl (+)V3.5σ3.4σ--
D-ring Vinyl (+)VI4.1σ4.6σ3.4σ2.6σ
C-ring
C-propionate (-)IX−6.4σ−5.8σ−5.9σ−4.7σ
C-ring (-)VII−5.5σ−4.4σ−5.1σ−3.2σ
C-propionate (+)X6.9σ6.4σ7.4σ4.9σ
C-ring (+)VIII4.9σ5.2σ3.2σ5.4σ
A-ring
A-ring C = O (-)XVIII−4.3σ−4.7σ−3.6σ−2.9σ
A-ring C = O (+)XIX5.2σ5.0σ5.7σ4.7σ
His260
His260 Sidechain(-)XI−6.7σ−5.2σ−6.3σ−4.1σ
His260 Sidechain(+)XII4.3σ4.5σ3.5σ3.2σ
Tyr263
Tyr263 Sidechain(-)XIII−6.9σ−7.0σ−4.7σ−6.2σ
Tyr263 Sidechain(+)XIV4.8σ4.9σ3.0σ4.3σ
Asp207
Asp207 Sidechain (-)XX−6.2σ−6.7σ−5.6σ−3.4σ
Asp207 Sidechain(+)XXI5.7σ5.0σ4.2σ-
Asp207 Backbone (-)XXII−6.2σ−5.4v−5.1σ−4.5σ
Asp207 Backbone (+)XIII4.6σ4.4σ4.5σ3.5σ
Tyr176
Tyr176 Sidechain(-)−4.9σ−3.9σ−5.1σ−3.8σ
Phe203
Phe203 Sidechain(-)−4.1σ−4.6σ−4.6σ-

First, we inspect the D-ring region at 1 ps (Figure 3). We observe strong negative difference density features on the atoms of the D-ring (marked I, II, III), correlating with density gains at both faces of the ring (IV, V, VI). These features strongly indicate that the D-ring twists. The positive feature IV homes the N-H and C = O groups, whereas V and VI indicate densities for the methyl and vinyl groups in the twisted ring (Figure 3c). Excellent agreement was obtained between the observed difference map and the difference map calculated from DrBphP1ps (Fc1ps-Fcdark), when the D-ring twists (C14-C15-C16-ND) from around 20° in the dark to 60° monomer A) and 90° monomer B) at 1 ps (Figure 3c and e). Although the twisting movement is clearly indicated by the difference map, we judge the precision of the angle to be low and approximately ±25°.

Figure 3 with 4 supplements see all
Observed and calculated difference electron densities reveal a twist of the D-ring and significant protein rearrangements.

The observed difference electron density with the refined DrBphPdark(blue) and DrBphP1ps(beige) structures, shown for (a) the B-, C-, and D-ring surroundings, (b) the strictly conserved His260 and Tyr263, and (c) the D-ring. The calculated difference electron density shown for (d) His260 and Tyr263 and (e) the D-ring. The D-ring twists counter-clockwise when viewed along C15-C16 bond toward the C-ring. The observed difference electron density is contoured at 3.3 σ. And the calculated difference electron density is contoured at 3.5 and 5.0 σ for panel d and e, respectively. Monomer A is shown in this figure.

Concomitant with the twist of the D-ring, the C-ring translates by approximately 0.69 Å as indicated by the correlated negative (VII) and positive (VIII) electron density (Figure 3a). Furthermore, the C-ring propionate chain detaches from its conserved anchoring residues Ser272 and Ser274 (IX and X, Figure 3a). The strictly conserved His260 retracts from its position (XI and XII) and Tyr263 moves upward at 1 ps (XIII and XIV, Figure 3b). The water network connecting the C-ring propionate, the D-ring C = O, and His290 rearranges accordingly (Figure 3a). The excellent agreement between calculated and observed difference maps confirms these observations (Figure 3d and e, Figure 3—figure supplement 2). We conclude that the twist of the D-ring causes detachment of the C-ring propionate from the protein scaffold by dislocation of the C-ring, facilitated by the associated hydrogen bonding network.

Turning our attention to the B-ring, we find that the B-ring propionate breaks its salt bridge to Arg254 (Figure 3—figure supplement 3). However, this is not caused by movements of the chromophore backbone, as we observe little change on the B-ring itself. Instead, we find that a water bridge between the B- and C-ring propionates is broken as indicated by negative difference electron densities on the waters (Figure 3—figure supplement 3). Additionally, the highly conserved helix from Ser257 to Val269, moves away from the chromophore by an average of 0.36 Å in monomer A and 0.62 Å in monomer B (distances relative to the pyrrole water, Figure 3—figure supplement 4). The changes of the D-ring are transduced to Ser257 via the side chains of His260 and Tyr263, and as a result, the hydrogen bond of Ser257 to the B-ring propionate group breaks. The amino acids in the stretch are over 50% conserved (Figure 3—figure supplement 4), suggesting that it has evolved to transfer an ultrafast signal. We conclude that relaxation of the protein is necessary for the detachment of the B-ring propionate from the protein scaffold.

Next to the changes around the D-ring, the maps reveal strong difference electron density on the A-ring (XVIII and XIX), Asp207 (XX to XXIII) (Figure 4a) and the pyrrole water (XV) (Figure 4b). When interpreted and modelled as downward movement of the A-ring and Asp207 and photodissociation of the pyrrole water from the chromophore, excellent agreement between calculated and observed difference electron density is obtained (Figure 4—figure supplement 1). The A-ring is covalently attached to the protein backbone in phytochromes (Song et al., 2014), which renders complete isomerization impossible, but is sufficiently flexible to accommodate the proposed changes. The pyrrole water may either move to feature (XVI), or occupy an anisotropic, worm-shaped feature which extends from the A-ring to the D-ring (XVII) (Figure 4b).

Figure 4 with 1 supplement see all
Photodissociation of the pyrrole water, displacement of the A-ring and its effect on the proteins scaffold.

(a). The observed difference electron density displayed with the DrBphPdark (blue) and DrBphP1ps (beige) structures around the A-ring, Asp207 and pyrrole water (PW). The structural model was inconclusive as to whether the A-ring twists around the double bond between the B- and A-ring, or whether it tilts downward hinged on the connection between C- and B-ring. (b). The regions of the pyrrole water (PW) and the area between the pyrrole rings show negative and positive densities, respectively. The observed difference electron density is contoured at 3.3σ. (c). Density displayed for the backbone below the A-ring, including side chains of the strictly conserved Ile208 and Tyr176 as well as the surrounding water network. Monomer A is shown in this figure.

Furthermore, correlated negative and positive electron density features are observed on backbone atoms of the highly conserved stretch from Pro201 to His209. These difference electron density features indicate that the residues move away from the centre of the chromophore by an average of 0.54 Å and 0.57 Å in monomers A and B, respectively (Figure 3—figure supplement 4b). The stretch includes Asp207 and it is located between the A-ring and the PHY domain, which makes it plausible that the changes in the chromophore cause this protein rearrangement. The changes are complemented by significant rearrangements of a stretch of waters and a conserved Tyr176 (Figure 4c).

Discussion

The structure of DrBphPCBD 1 ps after photoexcition reveals changes of the biliverdin chromophore and the surrounding residues. We find a twist of the D-ring, displacement of the C-ring, and associated changes of the water network which connects the D-ring, the C-ring propionate, and His290. Further, we identify a disruption of the salt bridge between the B-ring propionate and the Arg254, and significant changes around the A-ring, Asp207 and the pyrrole water. The changes are retained at 10 ps, even though at a lower population (Figure 1—figure supplement 5). The extensive and coordinated structural changes in the binding pocket (Figure 3) manifest a liberation of the chromophore from the protein scaffold, which we propose to be necessary for the conformational rearrangements to occur in the downstream photoconversion to Pfr.

Infrared spectral data indicate significant reorganization of the chromophore and several amino acids including the PHY-tongue region as early as in Lumi-R state, which is the first known ground state intermediate in the photoconversion from Pr to Pfr (Ihalainen et al., 2018; van Thor et al., 2007). However, the structure of the bilin and the binding pocket in Lumi-R is not known, because structural information is missing. Since the quantum yield of reaching the Lumi-R state is low (on the order of 10%), spectroscopic investigation of the mechanism is difficult, and it is currently not fully established how the Lumi-R state is reached. Crystallographic data does not report on whether the chromophore is electronically excited or not, hence we cannot determine whether the structure that we observe is in a relaxed excited state or in a ground state. The time delay (1 ps) supports that our structure presents an intermediate enroute to the Lumi-R state.

The D-ring of the bilin chromophore isomerizes around the C15-C16 bond from (Z) in Pr to (E) in Pfr (Rüdiger et al., 1983; Burgie et al., 2016; Takala et al., 2014; Yang et al., 2008; Essen et al., 2008). Circular dichroism spectroscopy and solid-state NMR spectroscopy have indicated that the position of the D-ring inverts from an ’α’-facial (Pr) to a ’β’-facial (Pfr) position in cyanobacterial and plant phytochromes, whereas it stays ’α’-facial in bacterial phytochromes (Rockwell et al., 2009; Song et al., 2011; Song et al., 2018). Based on anticipated steric clashes with the C-ring methyl group, it has been proposed that the D-ring rotates counter-clockwise in plant and cyanobacterial phytochromes, but clockwise in bacterial phytochromes (Rockwell et al., 2009). Moreover, spectroscopy has shown that the D-ring of the bilin chromophore is already isomerized in the Lumi-R state (van Thor et al., 2007; Heyne et al., 2002; Yang et al., 2012).

Seemingly contradictory, we now observe that the D-ring is rotated counter-clockwise by tens of degrees for the bacterial DrBphP at 1 ps time delay (Figure 3c and e). The conformation is strongly supported by the difference map. It contains two positive peaks (V and VI in Figure 3c), which indicate the new position of the vinyl and methyl group of the D-ring. We tested models in which the D-ring was rotated in a clockwise direction, but the agreement with the experimental difference map decreased. Thus, it may be that the D-ring indeed rotates counter-clockwise in bacterial phytochromes, similar to plant and cyanobacterial phytochromes. For complete isomerization, this would mean that the C-ring moves out of the way during the rotation. We observe significant movements of the C-ring, which may be an indication for that such a mechanism is possible. Raising a note of caution, we cannot fully exclude that the truncation of our phytochrome construct or the crystal packing influences the direction of rotation. NMR studies have reported conformational heterogeneity in the chromophore binding pocket of phytochromes in solution (Song et al., 2011; Lim et al., 2018; Song et al., 2018; Gustavsson et al., 2020). Crystallization could select one of the conformations, which may have a preferred rotation in the counter-clockwise direction. More experiments are needed to clarify this question.

It is interesting to compare the structural changes at 1 ps time delay to the changes observed in the conversion between Pr and Pfr (Burgie et al., 2016; Takala et al., 2014; Stojković et al., 2014). Major changes include a flipped D-ring, changes in conserved residues of the chromophore-binding pocket, for example Tyr176, His201 and Phe203, and refolding of the PHY tongue. The PHY tongue is not included in our construct, but Tyr176 and Phe203 are associated with difference electron density features in our maps (Table 2). However, the movements are much smaller at 1 ps compared to the Pr-to-Pfr transition. This is not unexpected, given the short time delays, but it shows that the residues are tightly coupled to the chromophore. Interestingly, the Pr and Pfr structures also reveal a sliding movement of the entire chromophore (Yang et al., 2011; Burgie et al., 2016; Takala et al., 2014). This requires that the propionic groups have to break their bonds to the protein scaffold. Our data indicate that this is part of the primary photoresponse.

The photodissociation of the pyrrole water from the chromophore is a surprising finding. The pyrrole water is ubiquitously found in phytochrome structures (Essen et al., 2008; Yang et al., 2008; Otero et al., 2016; Burgie et al., 2016; Wagner et al., 2005; Burgie et al., 2014; Schmidt et al., 2018; Yang et al., 2011). Our fluence dependent SFX data show that the negative density on the pyrrole water is the last signal to disappear when lowering the photon excitation densities 10-fold (Figure 2). This makes us confident that the photodissociation reaction is not caused by multi-photon effects. The removal of the water requires significant energy, because the hydrogen bonds to the A-, B-, and C-rings of the chromophore and the backbone C = O group of Asp207 have to be broken. We do not think that the twist of the D-ring causes this through direct steric interactions, because there is no contact between the pyrrole water and the D-ring. Rather, it may be triggered by an excited state charge redistribution between the pyrrole water and the chromophore, for example by ultrafast proton or electron transfer (Toh et al., 2010). Such charge re-distributions are typically facilitated by changes in geometry (Nosenko et al., 2008) and may therefore be caused indirectly by structural changes of the A-, C-, or D-rings, but this requires further investigation.

Conformational changes of the A-ring, Asp207 and the pyrrole water have not been considered to occur on picosecond time scales. The strictly conserved Asp207 is a key residue for signal transduction because it connects the chromophore to the PHY-tongue in Pr and Pfr (Essen et al., 2008; Yang et al., 2008; Takala et al., 2014). Its displacement suggests, together with the relocation of the residue stretch surrounding it, that disruption of the GAF-PHY interface may occur as early as 1 ps after photoexcitation (Figure 3—figure supplement 4b). With a hydrogen bond to the pyrrole water and in tight steric contact with the A-ring, Asp207 thereby acts as an extended arm of the chromophore. We propose that the photodissociation of the pyrrole water from the bilin and the change of the A-ring are integral parts of ultrafast phytochrome signaling toward the PHY domain.

We demonstrate that within 1 ps, the D-ring twists, that the chromophore is liberated from the protein (Figure 5a) and that movements of the pyrrole water, the A-ring and Asp207 lead to signaling directed toward the PHY-tongue (Figure 5b). When mapped on the structure of the complete photosensory core module (Takala et al., 2014), both changes work together to destabilize the Arg466:Asp207 salt bridge. Tyr263 moves up, caused by the twist of the D-ring, and Asp207 moves down, caused by changes of the A-ring, retracting both residues from the salt bridge.

Two correlated photochemical events guide the primary photorespone of phytochrome proteins.

(a). The structures (DrBphPdark, blue and DrBphP1ps, beige) indicate that rotation of the D-ring initiates breakage of non-covalent bonds of the propionates to the protein scaffold. Even the C- and A-rings are displaced significantly and the pyrrole water is dislocated from its original location at 1 ps (shade). (b). The same structures are overlayed with the complete photosensory core in Pr state (PDB ID 4O0P, pink) (Takala et al., 2014). The scissor-like separation of Asp207 and Tyr263 could result in breakage of the hydrogen bonds to Arg466 of the conserved PRXSF motif located in the PHY-tongue region.

Our data reveal a highly collective primary photoresponse for phytochromes. This is consistent with the fact that most point mutations of conserved residues alter, but do not inhibit, photoconversion (Wagner et al., 2008). The ultrafast structural changes are more extensive than in bacteriorhodopsin, photoactive yellow proteins, and in a fluorescent protein (Pande et al., 2016; Nogly et al., 2018; Coquelle et al., 2018). While previously observed ultrafast backbone movements have been interpreted as ’protein quakes’ for myoglobin and bacteriorhodopsin (Barends et al., 2015; Nogly et al., 2018; Toh et al., 2010), the present backbone motion in the phytochrome binding pocket are much more directed (Figure 3—figure supplement 4). The changes occur in highly conserved regions of the protein and are part of the collective signaling response of the entire binding pocket.

Phytochromes have to be able to stabilize the bilin and to direct its photoisomerization from two photochemical ground states, Pr and Pfr. These differ both structurally and electronically, which precludes a single reaction trajectory for isomerization in the two directions. With this in mind, the observed primary photoresponse is reasonable. The structural signal is highly delocalized already at 1 ps, causing near-simultaneous liberation of the chromophore and initial signal transduction. We propose that these reaction trajectories stabilize each other, navigating the protein into a productive reaction path. The multidimensional reaction trajectory is consistent with the low quantum yields for photoconversion (Lamparter et al., 1997), which are characteristic for the phytochrome superfamily. Whereas the twisting motion of the D-ring has been the working model for phytochrome activation and is now confirmed, the photodissociation of the pyrrole water is highly surprising. We propose that both chemical events work together and enable phytochrome proteins to translate light information into structural signals, guiding the growth and development of plants, fungi, and bacteria on Earth.

Materials and methods

Protein purification and crystallization

Request a detailed protocol

The His6-tagged PAS-GAF domain from D. radiodurans (aa 1–321) in vector pET21b(+) (Wagner et al., 2005) was expressed and purified as previously described (Lehtivuori et al., 2013; Takala et al., 2014). The recombinant protein was expressed in Escherichia coli strain BL21(DE3), either with or without Ho1 to yield holo- or apoprotein, respectively. Cells were lysed with Emulsiflex and cleared by centrifugation (20,000 rpm, 30 min, +4°C). Full biliverdin incorporation was ensured by adding 8 mg of biliverdin hydrochloride (Frontier Scientific) per litre of cell culture to the cell lysate, followed by overnight incubation on ice. The protein was then purified at room temperature with HisTrap HP column (GE Healtcare) in 30 mM Tris, 50 mM NaCl and 5 mM imidazole (pH 8) and eluted with increasing imidazole concentration (gradient elution over 5–500 mM). Size-exclusion chromatography was then conducted with a HiLoad 26/600 Superdex 200 pg column (GE Healthcare) in buffer (30 mM Tris pH 8.0). Finally, the protein was concentrated to 30–50 mg/mL and flash-frozen in liquid nitrogen.

Crystals were set up under green safe light and grown in dark. Batch crystallization was performed as described (Edlund et al., 2016). 50 µL of purified protein (25–30 mg/mL) was added to 450 µL of reservoir solution (60 mM Sodium acetate pH 4.95, 3.3% PEG 400, 1 mM DTT and 30% 2-methyl-2,4-pentanediol) and immediately mixed. Initial microcrystals were grown on a tipping table at 4 °C for 48 hr. Once the microcrystals were formed, additional protein was added to increase crystal size. The microcrystals were first pelleted by brief centrifugation and 400 µL of supernatant was removed. 200 µL of diluted protein (14 mg/mL in 30 mM Tris pH 8.0) was then added to the microcrystals along with 200 µL of fresh reservoir solution. After 48 hr incubation on a tipping table at room temperature, crystals of diffraction quality (20–70 µm long needles) were formed (Figure 1—figure supplement 2 and Table 1).

Transient absorption experiment of microcrystals

Request a detailed protocol

Transient absorption experiments were performed on a home build setup based on a Ti:sapphire femtosecond laser system (1 kHz, 800 nm). The main beam was split into pump and probe beams. The pump beam was sent through the home build noncollinear optical parametric amplifier to produce excitation pulses at 640 nm central wavelength. The probe beam was focused on a 2 mm sapphire plate to generate broadband (400–760 nm) white light which was split by 50/50 beamsplitter to reference and probe beams. The mutual polarization of the pump and probe beams was set to the magic angle (54.7°) by Berek compensator. The probe beam was focused on a sample cuvette that was continuously translated in vertical axis to prevent sample degradation. The microcrystals were washed with crystallization buffer five time in order to remove the solubilized proteins. 2.5 µL of microcrystals including a small amount of crystallization buffer were placed between two CaF2 windows without a spacer. The OD of the sample was about 0.6. Time-resolved absorption changes were measured by detecting probe and reference beams dispersed on the double-diode array; the time delay between pump and probe pulses was set by a computer controlled delay line placed in the probe beam path. All measurements were carried out in room temperature.

Light scattering of the grease jet

Request a detailed protocol

In order to estimate the light intensity of the optical laser in the grease jet, the optical transmission of the grease, grease mixed with crystallization buffer, microcrystals in the grease matrix and pure microcrystals were measured with a transmittance diode-array UV-Vis spectrometer (Cary 8454, Agilent Technologies) (Figure 1—figure supplement 1c). 2.5 µl each of sample was placed between two CaF2 windows with a 50 µm Teflon spacer and squeezed together, resembling the characteristics of the jet during the XFEL experiments. The raw spectra of microcrystals (Figure 1—figure supplement 1c) in the buffer has been measured between two CaF2 windows without spacer to minimize the absorption loss, the pathlength was estimated to be ≤ 50µm.

Fluence calculations

Request a detailed protocol

The optical laser parameters used for the experiment were as follows: wavelength was 640 nm, the laser-spot dimensions at the focus was 100 × 80 μm2 FWHM (170 × 136 μm2 at 1/e2 intensity), the pulse energy was 40 µJ, the nominal pulse duration was 70 fs (not confirmed at the sample position), and the repetition rate was 30 Hz. The energy of a photon at 640 nm is 3.1×10-19J. Using the photon density of the laser at 1/e2 convention of 1.7 mJ/mm2, we obtain a photon fluence of 5.48×1015photonsmm-2. The extinction coefficient of biliverdin in the phytochrome at 640 nm (ϵ640) is 27.7×103M-1cm-1 and the cross section is then σ640=ln(10)ϵ6401000/NA=1.06×10-14mm2molecule-1, where NA is Avogadro’s number. Multiplying the photon fluence with the cross section yields 58 photons per molecule. Light scattering in the carrier matrix decreases the effective fluence of photons that interact with the crystals. Our absorption spectra of the grease (Figure 1—figure supplement 1) indicate that the grease is transparent when untreated, but attenuates the light intensity by 2 orders of magnitude in the visible region when mixed with crystallization buffer or crystals (pathlength 50 µm). This indicates that almost every photon is scattered. Therefore, even when neglecting the scattering of the jet surface, crystals will be exposed to a photon fluence that is significantly reduced. A reduction of the photon fluence by 2 orders of magnitude is a realistic assumption as we used grease jets with a diameter of 75 µm or 100 µm. Another factor that contributes to the reduction of the number of photons per chromophore and non-homogeneous illumination of the microcrystals is the orientation of the crystals and the high chromophore density in them. The first few chromophores in the light path will shade the remaining chromophores in the needle-shaped crystals. Since the X-rays probe every molecule in their path with approximately same likelihood, the average photon fluence per probed chromophore is reduced. Assuming that the effective fluence inside the grease jet is reduced by a factor of 100, we estimate an average number of 0.5–1 photons per chromophore. This is consistent with the photoexcitation yield of 8% and with our experimental finding that the difference signal vanishes under the noise signal when reducing the photon fluence by a factor of 10 (Figure 2).

SFX data acquisition

Request a detailed protocol

Serial femtosecond crystallographic data were collected at SPring-8 Angstrom Compact Free Electron Laser (SACLA) in two beamtimes in October 2018 and May 2019. The microcrystals were pelleted by brief centrifugation and the crystal pellet was mixed with 180 µL of grease. The grease/crystal mixture was loaded into a 4 mm sample reservoir for data acquisition. The sample was delivered to the X-ray beam at a flow rate of 2.5 µL/min or 4.2 µL/min for 75 µm and 100 µm diameter nozzles, respectively. The time resolution of the experiment was limited by the jitter of the XFEL of 100 fs r.m.s. The experimental settings were nominally the same for the 1 ps and 10 ps delay times and all data were recorded during 7 hr of beamtime. We also recorded data at 3 ps delay time, but these generated electron density maps of poor quality due to an unknown reason and were therefore not analyzed further.

Data processing

Request a detailed protocol

The background of the detector was estimated by averaging the first 150 dark images in each run and then subtracted from each diffraction pattern. Diffraction images with Bragg spots (the ‘hits’) were found by a version of Cheetah adapted for SACLA (Nakane et al., 2016; Barty et al., 2014). These hits were indexed by the program CrystFEL (version 0.6.3) (White et al., 2012). Indexing was performed using Dirax and Mosflm (Duisenberg, 1992; Battye et al., 2011). Spot finding in each diffraction image was done with the peakfinder8 algorithm using the parameters (min SNR = 4.5, threshold = 100, minimum pixel counts = 3). The indexed patterns were merged and scaled using partialator in CrystFEL and hkl files were produced. The figure of merits (Table 1) were calculated by using compare_hkl and check_hkl in CrystFEL. The histograms of the unit cell parameters are presented in Figure 3. All diffraction images have been deposited to CXIDB under ID 121.

Refinement of dark structure

Request a detailed protocol

The initial phases were solved by molecular replacement with Phaser (McCoy et al., 2007) and the PAS-GAF crystal structure (PDB ID 5K5B) (Edlund et al., 2016) as a search model. The structure was refined with REFMAC version 5.8.0135 (Murshudov et al., 2011) with a weight factor for the geometry restraints of 0.05, accompanied by model building steps with Coot 0.8.2 (Emsley et al., 2010). The final structure (DrBphPdark) had Rwork/Rfree of 0.161/0.192 and no Ramachandran outliers (Table 1). The coordinates and structure factors have been deposited in the Protein Data Bank under the accession code 6T3L.

Computation of difference electron density maps

Request a detailed protocol

The difference structure factors (ΔF) are computed from the measured structure factor amplitudes in dark and for preset delay times between laser and X-ray pulses as |ΔFo|=w(|Folight|-|Fodark|) and with phases taken from the dark structural model (DrBphPdark). |Fodark| and |Folight| were brought to the absolute scale by first scaling |Fodark| to |Fcdark| and then scaling |Folight| to |Fodark| using the CCP4 program Scaleit (Winn et al., 2011). Difference Fourier density maps were calculated with a low resolution scaling cut-off at 18 Å . A weighting factor (w) was determined for each reflection to reduce the influence of outliers (Ren et al., 2001). From the weighted ΔF, a difference electron density map (Δρ) is calculated using the program ‘fft’ from the ccp4 suite of programs (Winn et al., 2011). Since ΔF are on the absolute scale, Δρ is on half the absolute scale as a result of the difference Fourier approximation (Henderson and Moffat, 1971; Pandey et al., 2020).

Structure refinement of light structure

Request a detailed protocol

Extrapolated structural factors were assembled from amplitudes computed as |Fe|=|Fcdark|+α*|ΔFo|. The Fcdark denotes the calculated structure factors of the refined dark structure (DrBphPdark). The phases were taken from DrBphPdark is inversely related to the population of the photoinduced state by (100/α)*2 (Pandey et al., 2020; Henderson and Moffat, 1971). We estimated α based on Fe map features in the chromophore-binding pocket. Too high values for α lead to physically unrealistic negative electron density. We converged to α=25, which corresponds to 8% photoexcitation yield.

Fe represents the pure structure factor of the photo-activated state (Figure 3—figure supplement 1). Refinement of a structural model was then performed in real and reciprocal space, using Coot (Emsley et al., 2010) and Phenix (Adams et al., 2010). The equilibrium values for the restraints used in the refinement of the biliverdin chromophore were taken from a minimal energy biliverdin ground state geometry that was obtained at the B3LYP/6–31G* level of density functional theory. Torsional restraints for the excited state geometry with the twisted D-ring were obtained at the SA(5)- CASSCF(12,12)/cc-pVDZ level of ab initio theory. We removed the torsion restrains for the C/D-ring (C14-C15-C16-C17; C14-C15-C16-ND; C13-C14-C15-C16; NC-C14-C15-C16) and for the A/B-ring (C3-C4-C5-C6; NA-C4-C5-C6; C4-C5-C6-C7; C4-C5-C6-NB) during refinement. The overall aim of the refinement was to maximize the agreement between the observed and calculated difference maps. To evaluate the agreement, we subtracted the calculated from the observed difference electron density (Δρo-Δρc). The computation of this difference-difference maps require scaling of the maps to each other. To do so, the highest and lowest intensities of Δρo were scaled to the corresponding maximum and minimum of Δρc and the observed Δρo were interpolated linearly according to this scaling. The resulting difference-difference electron density map was used to identify sites, which required further optimization in subsequent refinement steps. Calculation of Pearson Correlation Coefficient (PCC) values between the Δρo and Δρc were applied to guide refinement of specific regions, such as the D-ring and the whole chromophore region. To do so, the correlation was determined based on electron density within a sphere with a radius of 3.5 Å or 10 Å centred on the D-ring or pyrrole water, respectively. As a final step in the refinement procedure, we refined the models with REFMAC version 5.8.0135 (Murshudov et al., 2011) with high geometry restraints (weight factor 0.005). This was done against phased extrapolated structure factors, using the phases of the refined light and dark structure for computation of phased ΔF as described (Pande et al., 2016). The structures did not change much, although the R factors dropped in this last step of refinement to Rwork/Rfree of 0.230/0.256 (Table 1). The coordinates and structure factors have been deposited in the Protein Data Bank under the accession code 6T3U.

Data availability

Crystallography data have been submitted to protein data bank (PDB) dark:ID: D_1292104678 and PDB ID: 6T3L 1ps:ID: D_1292104679 and PDB ID: 6T3U. Raw diffraction images are deposited under ID 121 at CXIDB.

The following data sets were generated
    1. Claesson E
    2. Takala H
    3. Wahlgren W
    4. Pandey S
    5. Schmidt M
    6. Westenhoff S
    (2020) RCSB Protein Data Bank
    ID 6T3L. PAS-GAF fragment from Deinococcus radiodurans phytochrome in dark state.
    1. Claesson E
    2. Takala H
    3. Wahlgren W
    4. Pandey S
    5. Schmidt M
    6. Westenhoff S
    (2020) RCSB Protein Data Bank
    ID 6T3U. PAS-GAF fragment from Deinococcus radiodurans phytochrome 1ps after photoexcitation.

References

    1. Emsley P
    2. Lohkamp B
    3. Scott WG
    4. Cowtan K
    (2010) Features and development of coot
    Acta Crystallographica. Section D, Biological Crystallography 66:486–501.
    https://doi.org/10.1107/S0907444910007493

Article and author information

Author details

  1. Elin Claesson

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Planned experiments, Purified the protein and prepared the microcrystals, Acquired the data at SACLA, Computed difference electron density maps and refined the structural models, Wrote the paper with input from all authors
    Contributed equally with
    Weixiao Yuan Wahlgren and Heikki Takala
    Competing interests
    No competing interests declared
  2. Weixiao Yuan Wahlgren

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Purified the protein and prepared the microcrystals, Acquired the data at SACLA, Computed difference electron density maps and refined the structural models
    Contributed equally with
    Elin Claesson and Heikki Takala
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-0413-165X
  3. Heikki Takala

    1. Department of Biological and Environmental Science, Nanoscience Center, University of Jyvaskyla, Jyvaskyla, Finland
    2. Department of Anatomy, Faculty of Medicine, University of Helsinki, Helsinki, Finland
    Contribution
    Purified the protein and prepared the microcrystals, Acquired the data at SACLA, Computed difference electron density maps and refined the structural models
    Contributed equally with
    Elin Claesson and Weixiao Yuan Wahlgren
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-2518-8583
  4. Suraj Pandey

    Physics Department, University of Wisconsin-Milwaukee, Milwaukee, United States
    Contribution
    Acquired the data at SACLA and processed the crystallography data
    Competing interests
    No competing interests declared
  5. Leticia Castillon

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Purified the protein and prepared the microcrystals, Acquired the data at SACLA
    Competing interests
    No competing interests declared
  6. Valentyna Kuznetsova

    Department of Biological and Environmental Science, Nanoscience Center, University of Jyvaskyla, Jyvaskyla, Finland
    Contribution
    Acquired the data at SACLA, Measured the UV-Vis spectroscopic data
    Competing interests
    No competing interests declared
  7. Léocadie Henry

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  8. Matthijs Panman

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0003-3853-123X
  9. Melissa Carrillo

    Department of Biology, Northeastern Illinois University, Chicago, United States
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  10. Joachim Kübel

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  11. Rahul Nanekar

    Department of Biological and Environmental Science, Nanoscience Center, University of Jyvaskyla, Jyvaskyla, Finland
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  12. Linnéa Isaksson

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  13. Amke Nimmrich

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  14. Andrea Cellini

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  15. Dmitry Morozov

    Department of Chemistry, Nanoscience Center, University of Jyvaskyla, Jyvaskyla, Finland
    Contribution
    Provided the geometry restrains biliverdin in the excited state
    Competing interests
    No competing interests declared
  16. Michał Maj

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  17. Moona Kurttila

    Department of Biological and Environmental Science, Nanoscience Center, University of Jyvaskyla, Jyvaskyla, Finland
    Contribution
    Measured the UV-Vis spectroscopic data
    Competing interests
    No competing interests declared
  18. Robert Bosman

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  19. Eriko Nango

    1. Department of Cell Biology, Graduate School of Medicine, Kyoto University, Kyoto, Japan
    2. RIKEN SPring-8 Center, Hyogo, Japan
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  20. Rie Tanaka

    1. Department of Cell Biology, Graduate School of Medicine, Kyoto University, Kyoto, Japan
    2. RIKEN SPring-8 Center, Hyogo, Japan
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  21. Tomoyuki Tanaka

    1. Department of Cell Biology, Graduate School of Medicine, Kyoto University, Kyoto, Japan
    2. RIKEN SPring-8 Center, Hyogo, Japan
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  22. Luo Fangjia

    1. Department of Cell Biology, Graduate School of Medicine, Kyoto University, Kyoto, Japan
    2. RIKEN SPring-8 Center, Hyogo, Japan
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  23. So Iwata

    1. Department of Cell Biology, Graduate School of Medicine, Kyoto University, Kyoto, Japan
    2. RIKEN SPring-8 Center, Hyogo, Japan
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  24. Shigeki Owada

    1. RIKEN SPring-8 Center, Hyogo, Japan
    2. Japan Synchrotron Radiation Research Institute, Hyogo, Japan
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  25. Keith Moffat

    Department of Biochemistry and Molecular Biology and Institute for Biophysical Dynamics, University of Chicago, Chicago, United States
    Contribution
    Conceived the experiments
    Competing interests
    No competing interests declared
  26. Gerrit Groenhof

    Department of Chemistry, Nanoscience Center, University of Jyvaskyla, Jyvaskyla, Finland
    Contribution
    Provided the geometry restrains biliverdin in the excited state
    Competing interests
    No competing interests declared
  27. Emina A Stojković

    Department of Biology, Northeastern Illinois University, Chicago, United States
    Contribution
    Acquired the data at SACLA
    Competing interests
    No competing interests declared
  28. Janne A Ihalainen

    Department of Biological and Environmental Science, Nanoscience Center, University of Jyvaskyla, Jyvaskyla, Finland
    Contribution
    Conceived the experiments, Acquired the data at SACLA
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-8741-1587
  29. Marius Schmidt

    Physics Department, University of Wisconsin-Milwaukee, Milwaukee, United States
    Contribution
    Conceived the experiments, Acquired the data at SACLA, Computed difference electron density maps and refined the structural models, Wrote the paper with input from all authors
    For correspondence
    smarius@uwm.edu
    Competing interests
    No competing interests declared
  30. Sebastian Westenhoff

    Department of Chemistry and Molecular Biology, University of Gothenburg, Gothenburg, Sweden
    Contribution
    Conceived the experiments, Planned experiments, Acquired the data at SACLA, Computed difference electron density maps and refined the structural models, Wrote the paper with input from all authors
    For correspondence
    sebastian.westenhoff.2@gu.se
    Competing interests
    No competing interests declared
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-6961-8015

Funding

European Research Council (279944)

  • Sebastian Westenhoff

Academy of Finland (285461)

  • Heikki Takala

Academy of Finland (296135)

  • Janne A Ihalainen

Jane and Aatos Erkko Foundation

  • Janne A Ihalainen

National Science Foundation ((NSF)-MCB-RUI 141336)

  • Emina A Stojković

National Science Foundation (NSF-MCB-EAGER 183951)

  • Emina A Stojković

Horizon 2020 (H2020-INFRAEDI-02-2018-82383)

  • Gerrit Groenhof

Horizon 2020 (H2020-EINFRA-2015-1-67572)

  • Gerrit Groenhof

National Science Foundation (NSF-1231306)

  • Marius Schmidt

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Acknowledgements

The experiments at SACLA were performed at BL3 with the approval of the Japan Synchrotron Radiation Research Institute (JASRI) (Proposal No. 2018A8055 and 2019A8007). SW acknowledges the European Research Council for support (grant number: 279944). This work was supported by Academy of Finland grants 285461 and 296135 (HT and JAI, respectively) and Jane and Aatos Erkko foundation (JAI). This research is partially supported by Platform Project for Supporting Drug Discovery and Life Science Research (Basis for Supporting Innovative Drug Discovery and Life Science Research (BINDS)) from Japan Agency for Medical Research and Development (AMED). We thank Dr. Takanori Nakane for assistance with data processing during the beamtime and Heli Lehtivuori for the spectroscopic measurements with the microcrystals. This work was supported by NSF Science and Technology Centers grant NSF-1231306 (‘Biology with X-ray Lasers’) (MS), the National Science Foundation (NSF)-MCB-RUI 1413360 and NSF-MCB-EAGER 1839513 Research Grants to EAS. This work has been done as part of the BioExcel CoE (www.bioexcel.eu), a project funded by the European Union contracts H2020-INFRAEDI-02-2018-823830 and H2020-EINFRA-2015-1-675728 (GG).

Copyright

© 2020, Claesson et al.

This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 3,821
    views
  • 448
    downloads
  • 87
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Elin Claesson
  2. Weixiao Yuan Wahlgren
  3. Heikki Takala
  4. Suraj Pandey
  5. Leticia Castillon
  6. Valentyna Kuznetsova
  7. Léocadie Henry
  8. Matthijs Panman
  9. Melissa Carrillo
  10. Joachim Kübel
  11. Rahul Nanekar
  12. Linnéa Isaksson
  13. Amke Nimmrich
  14. Andrea Cellini
  15. Dmitry Morozov
  16. Michał Maj
  17. Moona Kurttila
  18. Robert Bosman
  19. Eriko Nango
  20. Rie Tanaka
  21. Tomoyuki Tanaka
  22. Luo Fangjia
  23. So Iwata
  24. Shigeki Owada
  25. Keith Moffat
  26. Gerrit Groenhof
  27. Emina A Stojković
  28. Janne A Ihalainen
  29. Marius Schmidt
  30. Sebastian Westenhoff
(2020)
The primary structural photoresponse of phytochrome proteins captured by a femtosecond X-ray laser
eLife 9:e53514.
https://doi.org/10.7554/eLife.53514

Share this article

https://doi.org/10.7554/eLife.53514

Further reading

    1. Structural Biology and Molecular Biophysics
    Johannes Elferich, Lingli Kong ... Nikolaus Grigorieff
    Research Advance

    Images taken by transmission electron microscopes are usually affected by lens aberrations and image defocus, among other factors. These distortions can be modeled in reciprocal space using the contrast transfer function (CTF). Accurate estimation and correction of the CTF is essential for restoring the high-resolution signal in cryogenic electron microscopy (cryoEM). Previously, we described the implementation of algorithms for this task in the cisTEM software package (Grant et al., 2018). Here we show that taking sample characteristics, such as thickness and tilt, into account can improve CTF estimation. This is particularly important when imaging cellular samples, where measurement of sample thickness and geometry derived from accurate modeling of the Thon ring pattern helps judging the quality of the sample. This improved CTF estimation has been implemented in CTFFIND5, a new version of the cisTEM program CTFFIND. We evaluated the accuracy of these estimates using images of tilted aquaporin crystals and eukaryotic cells thinned by focused ion beam milling. We estimate that with micrographs of sufficient quality CTFFIND5 can measure sample tilt with an accuracy of 3° and sample thickness with an accuracy of 5 nm.

    1. Structural Biology and Molecular Biophysics
    Mart GF Last, Leoni Abendstein ... Thomas H Sharp
    Tools and Resources

    Segmentation is a critical data processing step in many applications of cryo-electron tomography. Downstream analyses, such as subtomogram averaging, are often based on segmentation results, and are thus critically dependent on the availability of open-source software for accurate as well as high-throughput tomogram segmentation. There is a need for more user-friendly, flexible, and comprehensive segmentation software that offers an insightful overview of all steps involved in preparing automated segmentations. Here, we present Ais: a dedicated tomogram segmentation package that is geared towards both high performance and accessibility, available on GitHub. In this report, we demonstrate two common processing steps that can be greatly accelerated with Ais: particle picking for subtomogram averaging, and generating many-feature segmentations of cellular architecture based on in situ tomography data. Featuring comprehensive annotation, segmentation, and rendering functionality, as well as an open repository for trained models at aiscryoet.org, we hope that Ais will help accelerate research and dissemination of data involving cryoET.