The primary structural photoresponse of phytochrome proteins captured by a femtosecond X-ray laser
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., 2014; Sanchez 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.
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 ( 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).
From the time-resolved data, we calculated Fourier difference electron density maps (), 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).
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 (), 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°.
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).
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.
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 protocolThe -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 protocolTransient 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 protocolIn 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 protocolThe 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 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 . Using the photon density of the laser at convention of 1.7 mJ/mm2, we obtain a photon fluence of . The extinction coefficient of biliverdin in the phytochrome at 640 nm () is and the cross section is then , where 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 protocolSerial 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 protocolThe 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 protocolThe 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 protocolThe difference structure factors () are computed from the measured structure factor amplitudes in dark and for preset delay times between laser and X-ray pulses as and with phases taken from the dark structural model (DrBphPdark). and were brought to the absolute scale by first scaling to and then scaling to 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 () was determined for each reflection to reduce the influence of outliers (Ren et al., 2001). From the weighted , a difference electron density map () is calculated using the program ‘fft’ from the ccp4 suite of programs (Winn et al., 2011). Since 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 protocolExtrapolated structural factors were assembled from amplitudes computed as . The 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 (Pandey et al., 2020; Henderson and Moffat, 1971). We estimated based on map features in the chromophore-binding pocket. Too high values for lead to physically unrealistic negative electron density. We converged to , which corresponds to 8% photoexcitation yield.
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 (). The computation of this difference-difference maps require scaling of the maps to each other. To do so, the highest and lowest intensities of were scaled to the corresponding maximum and minimum of and the observed 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 and 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 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.
-
CXIDBThe primary structural photoresponse of phytochrome proteins captured by a femtosecond X-ray laser.https://doi.org/10.11577/1607859
-
RCSB Protein Data BankID 6T3L. PAS-GAF fragment from Deinococcus radiodurans phytochrome in dark state.
-
RCSB Protein Data BankID 6T3U. PAS-GAF fragment from Deinococcus radiodurans phytochrome 1ps after photoexcitation.
References
-
PHENIX: a comprehensive Python-based system for macromolecular structure solutionActa Crystallographica Section D Biological Crystallography 66:213–221.https://doi.org/10.1107/S0907444909052925
-
Cheetah: software for high-throughput reduction and analysis of serial femtosecond X-ray diffraction dataJournal of Applied Crystallography 47:1118–1131.https://doi.org/10.1107/S1600576714007626
-
iMOSFLM: a new graphical interface for diffraction-image processing with MOSFLMActa Crystallographica. Section D, Biological Crystallography 67:271–281.https://doi.org/10.1107/S0907444910048675
-
Crystallographic and electron microscopic analyses of a bacterial phytochrome reveal local and global rearrangements during photoconversionJournal of Biological Chemistry 289:24573–24587.https://doi.org/10.1074/jbc.M114.571661
-
Indexing in single-crystal diffractometry with an obstinate list of reflectionsJournal of Applied Crystallography 25:92–96.https://doi.org/10.1107/S0021889891010634
-
Features and development of cootActa Crystallographica. Section D, Biological Crystallography 66:486–501.https://doi.org/10.1107/S0907444910007493
-
Modulation of structural heterogeneity controls phytochrome photoswitchingBiophysical Journal 118:415–421.https://doi.org/10.1016/j.bpj.2019.11.025
-
The difference fourier technique in protein crystallography: errors and their treatmentActa Crystallographica Section B Structural Crystallography and Crystal Chemistry 27:1414–1420.https://doi.org/10.1107/S0567740871004060
-
Chromophore-Protein interplay during the phytochrome photocycle revealed by Step-Scan FTIR spectroscopyJournal of the American Chemical Society 140:12396–12404.https://doi.org/10.1021/jacs.8b04659
-
Fluorescence properties of the chromophore-binding domain of bacteriophytochrome from Deinococcus radioduransThe Journal of Physical Chemistry B 117:11049–11057.https://doi.org/10.1021/jp312061b
-
Phaser crystallographic softwareJournal of Applied Crystallography 40:658–674.https://doi.org/10.1107/S0021889807021206
-
REFMAC5 for the refinement of macromolecular crystal structuresActa Crystallographica Section D Biological Crystallography 67:355–367.https://doi.org/10.1107/S0907444911001314
-
Data processing pipeline for serial femtosecond crystallography at SACLAJournal of Applied Crystallography 49:1035–1041.https://doi.org/10.1107/S1600576716005720
-
Proton transfer with a twist? femtosecond dynamics of 7-(2-pyridyl)indole in condensed phase and in supersonic jetsAngewandte Chemie International Edition 47:6037–6040.https://doi.org/10.1002/anie.200801350
-
The D-ring, not the A-ring, rotates in Synechococcus OS-B′ PhytochromeJournal of Biological Chemistry 289:2552–2562.https://doi.org/10.1074/jbc.M113.520031
-
FTIR spectroscopy revealing Light-Dependent refolding of the conserved tongue region of bacteriophytochromeThe Journal of Physical Chemistry Letters 5:2512–2515.https://doi.org/10.1021/jz501189t
-
Formation of the early photoproduct lumi-R of cyanobacterial phytochrome cph1 observed by ultrafast mid-infrared spectroscopyJournal of the American Chemical Society 129:126–132.https://doi.org/10.1021/ja0660709
-
CrystFEL : a software suite for snapshot serial crystallographyJournal of Applied Crystallography 45:335–341.https://doi.org/10.1107/S0021889812002312
-
Overview of the CCP4 suite and current developmentsActa Crystallographica. Section D, Biological Crystallography 67:235–242.https://doi.org/10.1107/S0907444910045749
-
Real-time tracking of phytochrome's orientational changes during Pr photoisomerizationJournal of the American Chemical Society 134:1408–1411.https://doi.org/10.1021/ja209413d
Article and author information
Author details
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,817
- views
-
- 447
- downloads
-
- 87
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Structural Biology and Molecular Biophysics
Cyclic nucleotide-binding domain (CNBD) ion channels play crucial roles in cellular-signaling and excitability and are regulated by the direct binding of cyclic adenosine- or guanosine-monophosphate (cAMP, cGMP). However, the precise allosteric mechanism governing channel activation upon ligand binding, particularly the energetic changes within domains, remains poorly understood. The prokaryotic CNBD channel SthK offers a valuable model for investigating this allosteric mechanism. In this study, we investigated the conformational dynamics and energetics of the SthK C-terminal region using a combination of steady-state and time-resolved transition metal ion Förster resonance energy transfer (tmFRET) experiments. We engineered donor-acceptor pairs at specific sites within a SthK C-terminal fragment by incorporating a fluorescent noncanonical amino acid donor and metal ion acceptors. Measuring tmFRET with fluorescence lifetimes, we determined intramolecular distance distributions in the absence and presence of cAMP or cGMP. The probability distributions between conformational states without and with ligand were used to calculate the changes in free energy (ΔG) and differences in free energy change (ΔΔG) in the context of a simple four-state model. Our findings reveal that cAMP binding produces large structural changes, with a very favorable ΔΔG. In contrast to cAMP, cGMP behaved as a partial agonist and only weakly promoted the active state. Furthermore, we assessed the impact of protein oligomerization and ionic strength on the structure and energetics of the conformational states. This study demonstrates the effectiveness of time-resolved tmFRET in determining the conformational states and the ligand-dependent energetics of the SthK C-terminal region.
-
- Structural Biology and Molecular Biophysics
Ciliary rootlets are striated bundles of filaments that connect the base of cilia to internal cellular structures. Rootlets are critical for the sensory and motile functions of cilia. However, the mechanisms underlying these functions remain unknown, in part due to a lack of structural information of rootlet organization. In this study, we obtain 3D reconstructions of membrane-associated and purified rootlets from mouse retina using cryo-electron tomography. We show that flexible protrusions on the rootlet surface, which emanate from the cross-striations, connect to intracellular membranes. In purified rootlets, the striations were classified into amorphous (A)-bands, associated with accumulations on the rootlet surface, and discrete (D)-bands corresponding to punctate lines of density that run through the rootlet. These striations connect a flexible network of longitudinal filaments. Subtomogram averaging suggests the filaments consist of two intertwined coiled coils. The rootlet’s filamentous architecture, with frequent membrane-connecting cross-striations, lends itself well for anchoring large membranes in the cell.