Search

The ability to study chemical dynamics on ultrafast timescales has greatly advanced with the introduction of X-ray free electron lasers (XFELs) providing short pulses of intense X-rays tailored to probe atomic structure and electronic configuration. Fully exploiting the full potential of XFELs requires specialized experimental endstations along with the development of techniques and methods to successfully carry out experiments. The liquid jet endstation (LJE) at the Linac Coherent Light Source (LCLS) has been developed to study photochemistry and biochemistry in solution systems using a combination of X-ray solution scattering (XSS), X-ray absorption spectroscopy (XAS), and X-ray emission spectroscopy (XES). The pump-probe setup utilizes an optical laser to excite the sample, which is subsequently probed by a hard X-ray pulse to resolve structural and electronic dynamics at their intrinsic femtosecond timescales. The LJE ensures reliable sample delivery to the X-ray interaction point via various liquid jets, enabling rapid replenishment of thin samples with millimolar concentrations and low sample volumes at the 120 Hz repetition rate of the LCLS beam. This paper provides a detailed description of the LJE design and of the techniques it enables, with an emphasis on the diagnostics required for real-time monitoring of the liquid jet and on the spatiotemporal overlap methods used to optimize the signal. Additionally, various scientific examples are discussed, highlighting the versatility of the LJE.

The advent of X-ray free-electron lasers (XFELs) has revolutionized fundamental science, from atomic to condensed matter physics, from chemistry to biology, giving researchers access to X-rays with unprecedented brightness, coherence and pulse duration. All XFEL facilities built until recently provided X-ray pulses at a relatively low repetition rate, with limited data statistics. Here, results from the first megahertz-repetition-rate X-ray scattering experiments at the Spectroscopy and Coherent Scattering (SCS) instrument of the European XFEL are presented. The experimental capabilities that the SCS instrument offers, resulting from the operation at megahertz repetition rates and the availability of the novel DSSC 2D imaging detector, are illustrated. Time-resolved magnetic X-ray scattering and holographic imaging experiments in solid state samples were chosen as representative, providing an ideal test-bed for operation at megahertz rates. Our results are relevant and applicable to any other non-destructive XFEL experiments in the soft X-ray range.

Oxidation states of individual metal atoms within a metalloprotein can be assigned by examining X-ray absorption edges, which shift to higher energy for progressively more positive valence numbers. Indeed, X-ray crystallography is well suited for such a measurement, owing to its ability to spatially resolve the scattering contributions of individual metal atoms that have distinct electronic environments contributing to protein function. However, as the magnitude of the shift is quite small, about +2 eV per valence state for iron, it has only been possible to measure the effect when performed with monochromated X-ray sources at synchrotron facilities with energy resolutions in the range 2-3 × 10-4 (ΔE/E). This paper tests whether X-ray free-electron laser (XFEL) pulses, which have a broader bandpass (ΔE/E = 3 × 10-3) when used without a monochromator, might also be useful for such studies. The program nanoBragg is used to simulate serial femtosecond crystallography (SFX) diffraction images with sufficient granularity to model the XFEL spectrum, the crystal mosaicity and the wavelength-dependent anomalous scattering factors contributed by two differently charged iron centers in the 110-amino-acid protein, ferredoxin. Bayesian methods are then used to deduce, from the simulated data, the most likely X-ray absorption curves for each metal atom in the protein, which agree well with the curves chosen for the simulation. The data analysis relies critically on the ability to measure the incident spectrum for each pulse, and also on the nanoBragg simulator to predict the size, shape and intensity profile of Bragg spots based on an underlying physical model that includes the absorption curves, which are then modified to produce the best agreement with the simulated data. This inference methodology potentially enables the use of SFX diffraction for the study of metalloenzyme mechanisms and, in general, offers a more detailed approach to Bragg spot data reduction.

Oxidation states of individual metal atoms within a metalloprotein can be assigned by examining X-ray absorption edges, which shift to higher energy for progressively more positive valence numbers. Indeed, X-ray crystallography is well suited for such a measurement, owing to its ability to spatially resolve the scattering contributions of individual metal atoms that have distinct electronic environments contributing to protein function. However, as the magnitude of the shift is quite small, about +2 eV per valence state for iron, it has only been possible to measure the effect when performed with monochromated X-ray sources at synchrotron facilities with energy resolutions in the range 2-3 × 10-4 (ΔE/E). This paper tests whether X-ray free-electron laser (XFEL) pulses, which have a broader bandpass (ΔE/E = 3 × 10-3) when used without a monochromator, might also be useful for such studies. The program nanoBragg is used to simulate serial femtosecond crystallography (SFX) diffraction images with sufficient granularity to model the XFEL spectrum, the crystal mosaicity and the wavelength-dependent anomalous scattering factors contributed by two differently charged iron centers in the 110-amino-acid protein, ferredoxin. Bayesian methods are then used to deduce, from the simulated data, the most likely X-ray absorption curves for each metal atom in the protein, which agree well with the curves chosen for the simulation. The data analysis relies critically on the ability to measure the incident spectrum for each pulse, and also on the nanoBragg simulator to predict the size, shape and intensity profile of Bragg spots based on an underlying physical model that includes the absorption curves, which are then modified to produce the best agreement with the simulated data. This inference methodology potentially enables the use of SFX diffraction for the study of metalloenzyme mechanisms and, in general, offers a more detailed approach to Bragg spot data reduction.

Oxidation states of individual metal atoms within a metalloprotein can be assigned by examining X-ray absorption edges, which shift to higher energy for progressively more positive valence numbers. Indeed, X-ray crystallography is well suited for such a measurement, owing to its ability to spatially resolve the scattering contributions of individual metal atoms that have distinct electronic environments contributing to protein function. However, as the magnitude of the shift is quite small, about +2 eV per valence state for iron, it has only been possible to measure the effect when performed with monochromated X-ray sources at synchrotron facilities with energy resolutions in the range 2-3 × 10-4 (ΔE/E). This paper tests whether X-ray free-electron laser (XFEL) pulses, which have a broader bandpass (ΔE/E = 3 × 10-3) when used without a monochromator, might also be useful for such studies. The program nanoBragg is used to simulate serial femtosecond crystallography (SFX) diffraction images with sufficient granularity to model the XFEL spectrum, the crystal mosaicity and the wavelength-dependent anomalous scattering factors contributed by two differently charged iron centers in the 110-amino-acid protein, ferredoxin. Bayesian methods are then used to deduce, from the simulated data, the most likely X-ray absorption curves for each metal atom in the protein, which agree well with the curves chosen for the simulation. The data analysis relies critically on the ability to measure the incident spectrum for each pulse, and also on the nanoBragg simulator to predict the size, shape and intensity profile of Bragg spots based on an underlying physical model that includes the absorption curves, which are then modified to produce the best agreement with the simulated data. This inference methodology potentially enables the use of SFX diffraction for the study of metalloenzyme mechanisms and, in general, offers a more detailed approach to Bragg spot data reduction.

How changes in enzyme structure and dynamics facilitate passage along the reaction coordinate is a fundamental unanswered question. Here, we use time-resolved mix-and-inject serial crystallography (MISC) at an X-ray free electron laser (XFEL), ambient-temperature X-ray crystallography, computer simulations, and enzyme kinetics to characterize how covalent catalysis modulates isocyanide hydratase (ICH) conformational dynamics throughout its catalytic cycle. We visualize this previously hypothetical reaction mechanism, directly observing formation of a thioimidate covalent intermediate in ICH microcrystals during catalysis. ICH exhibits a concerted helical displacement upon active-site cysteine modification that is gated by changes in hydrogen bond strength between the cysteine thiolate and the backbone amide of the highly strained Ile152 residue. These catalysis-activated motions permit water entry into the ICH active site for intermediate hydrolysis. Mutations at a Gly residue (Gly150) that modulate helical mobility reduce ICH catalytic turnover and alter its pre-steady-state kinetic behavior, establishing that helical mobility is important for ICH catalytic efficiency. These results demonstrate that MISC can capture otherwise elusive aspects of enzyme mechanism and dynamics in microcrystalline samples, resolving long-standing questions about the connection between nonequilibrium protein motions and enzyme catalysis.

A nonlinear least-squares method for refining a parametric expression describing the estimated errors of reflection intensities in serial crystallographic (SX) data is presented. This approach, which is similar to that used in the rotation method of crystallographic data collection at synchrotrons, propagates error estimates from photon-counting statistics to the merged data. Here, it is demonstrated that the application of this approach to SX data provides better SAD phasing ability, enabling the autobuilding of a protein structure that had previously failed to be built. Estimating the error in the merged reflection intensities requires the understanding and propagation of all of the sources of error arising from the measurements. One type of error, which is well understood, is the counting error introduced when the detector counts X-ray photons. Thus, if other types of random errors (such as readout noise) as well as uncertainties in systematic corrections (such as from X-ray attenuation) are completely understood, they can be propagated along with the counting error, as appropriate. In practice, most software packages propagate as much error as they know how to model and then include error-adjustment terms that scale the error estimates until they explain the variance among the measurements. If this is performed carefully, then during SAD phasing likelihood-based approaches can make optimal use of these error estimates, increasing the chance of a successful structure solution. In serial crystallography, SAD phasing has remained challenging, with the few examples of de novo protein structure solution each requiring many thousands of diffraction patterns. Here, the effects of different methods of treating the error estimates are estimated and it is shown that using a parametric approach that includes terms proportional to the known experimental uncertainty, the reflection intensity and the squared reflection intensity to improve the error estimates can allow SAD phasing even from weak zinc anomalous signal.

A nonlinear least-squares method for refining a parametric expression describing the estimated errors of reflection intensities in serial crystallographic (SX) data is presented. This approach, which is similar to that used in the rotation method of crystallographic data collection at synchrotrons, propagates error estimates from photon-counting statistics to the merged data. Here, it is demonstrated that the application of this approach to SX data provides better SAD phasing ability, enabling the autobuilding of a protein structure that had previously failed to be built. Estimating the error in the merged reflection intensities requires the understanding and propagation of all of the sources of error arising from the measurements. One type of error, which is well understood, is the counting error introduced when the detector counts X-ray photons. Thus, if other types of random errors (such as readout noise) as well as uncertainties in systematic corrections (such as from X-ray attenuation) are completely understood, they can be propagated along with the counting error, as appropriate. In practice, most software packages propagate as much error as they know how to model and then include error-adjustment terms that scale the error estimates until they explain the variance among the measurements. If this is performed carefully, then during SAD phasing likelihood-based approaches can make optimal use of these error estimates, increasing the chance of a successful structure solution. In serial crystallography, SAD phasing has remained challenging, with the few examples of de novo protein structure solution each requiring many thousands of diffraction patterns. Here, the effects of different methods of treating the error estimates are estimated and it is shown that using a parametric approach that includes terms proportional to the known experimental uncertainty, the reflection intensity and the squared reflection intensity to improve the error estimates can allow SAD phasing even from weak zinc anomalous signal.

The DIALS diffraction-modeling software package has been applied to serial crystallography data. Diffraction modeling is an exercise in determining the experimental parameters, such as incident beam wavelength, crystal unit cell and orientation, and detector geometry, that are most consistent with the observed positions of Bragg spots. These parameters can be refined by nonlinear least-squares fitting. In previous work, it has been challenging to refine both the positions of the sensors (metrology) on multipanel imaging detectors such as the CSPAD and the orientations of all of the crystals studied. Since the optimal models for metrology and crystal orientation are interdependent, alternate cycles of panel refinement and crystal refinement have been required. To simplify the process, a sparse linear algebra technique for solving the normal equations was implemented, allowing the detector panels to be refined simultaneously against the diffraction from thousands of crystals with excellent computational performance. Separately, it is shown how to refine the metrology of a second CSPAD detector, positioned at a distance of 2.5 m from the crystal, used for recording low-angle reflections. With the ability to jointly refine the detector position against the ensemble of all crystals used for structure determination, it is shown that ensemble refinement greatly reduces the apparent nonisomorphism that is often observed in the unit-cell distributions from still-shot serial crystallography. In addition, it is shown that batching the images by timestamp and re-refining the detector position can realistically model small, time-dependent variations in detector position relative to the sample, and thereby improve the integrated structure-factor intensity signal and heavy-atom anomalous peak heights.

The DIALS diffraction-modeling software package has been applied to serial crystallography data. Diffraction modeling is an exercise in determining the experimental parameters, such as incident beam wavelength, crystal unit cell and orientation, and detector geometry, that are most consistent with the observed positions of Bragg spots. These parameters can be refined by nonlinear least-squares fitting. In previous work, it has been challenging to refine both the positions of the sensors (metrology) on multipanel imaging detectors such as the CSPAD and the orientations of all of the crystals studied. Since the optimal models for metrology and crystal orientation are interdependent, alternate cycles of panel refinement and crystal refinement have been required. To simplify the process, a sparse linear algebra technique for solving the normal equations was implemented, allowing the detector panels to be refined simultaneously against the diffraction from thousands of crystals with excellent computational performance. Separately, it is shown how to refine the metrology of a second CSPAD detector, positioned at a distance of 2.5 m from the crystal, used for recording low-angle reflections. With the ability to jointly refine the detector position against the ensemble of all crystals used for structure determination, it is shown that ensemble refinement greatly reduces the apparent nonisomorphism that is often observed in the unit-cell distributions from still-shot serial crystallography. In addition, it is shown that batching the images by timestamp and re-refining the detector position can realistically model small, time-dependent variations in detector position relative to the sample, and thereby improve the integrated structure-factor intensity signal and heavy-atom anomalous peak heights.

The M2 proton channel of influenza A is a drug target that is essential for the reproduction of the flu virus. It is also a model system for the study of selective, unidirectional proton transport across a membrane. Ordered water molecules arranged in "wires" inside the channel pore have been proposed to play a role in both the conduction of protons to the four gating His37 residues and the stabilization of multiple positive charges within the channel. To visualize the solvent in the pore of the channel at room temperature while minimizing the effects of radiation damage, data were collected to a resolution of 1.4 Å using an X-ray free-electron laser (XFEL) at three different pH conditions: pH 5.5, pH 6.5, and pH 8.0. Data were collected on the Inwardopen state, which is an intermediate that accumulates at high protonation of the His37 tetrad. At pH 5.5, a continuous hydrogen-bonded network of water molecules spans the vertical length of the channel, consistent with a Grotthuss mechanism model for proton transport to the His37 tetrad. This ordered solvent at pH 5.5 could act to stabilize the positive charges that build up on the gating His37 tetrad during the proton conduction cycle. The number of ordered pore waters decreases at pH 6.5 and 8.0, where the Inwardopen state is less stable. These studies provide a graphical view of the response of water to a change in charge within a restricted channel environment.

A major goal for X-ray free-electron laser (XFEL) based science is to elucidate structures of biological molecules without the need for crystals. Filament systems may provide some of the first single macromolecular structures elucidated by XFEL radiation, since they contain one-dimensional translational symmetry and thereby occupy the diffraction intensity region between the extremes of crystals and single molecules. Here, we demonstrate flow alignment of as few as 100 filaments (Escherichia coli pili, F-actin, and amyloid fibrils), which when intersected by femtosecond X-ray pulses result in diffraction patterns similar to those obtained from classical fiber diffraction studies. We also determine that F-actin can be flow-aligned to a disorientation of approximately 5 degrees. Using this XFEL-based technique, we determine that gelsolin amyloids are comprised of stacked β-strands running perpendicular to the filament axis, and that a range of order from fibrillar to crystalline is discernable for individual α-synuclein amyloids.

A major goal for X-ray free-electron laser (XFEL) based science is to elucidate structures of biological molecules without the need for crystals. Filament systems may provide some of the first single macromolecular structures elucidated by XFEL radiation, since they contain one-dimensional translational symmetry and thereby occupy the diffraction intensity region between the extremes of crystals and single molecules. Here, we demonstrate flow alignment of as few as 100 filaments (Escherichia coli pili, F-actin, and amyloid fibrils), which when intersected by femtosecond X-ray pulses result in diffraction patterns similar to those obtained from classical fiber diffraction studies. We also determine that F-actin can be flow-aligned to a disorientation of approximately 5 degrees. Using this XFEL-based technique, we determine that gelsolin amyloids are comprised of stacked β-strands running perpendicular to the filament axis, and that a range of order from fibrillar to crystalline is discernable for individual α-synuclein amyloids.

The M2 proton channel of influenza A is a drug target that is essential for the reproduction of the flu virus. It is also a model system for the study of selective, unidirectional proton transport across a membrane. Ordered water molecules arranged in "wires" inside the channel pore have been proposed to play a role in both the conduction of protons to the four gating His37 residues and the stabilization of multiple positive charges within the channel. To visualize the solvent in the pore of the channel at room temperature while minimizing the effects of radiation damage, data were collected to a resolution of 1.4 Å using an X-ray free-electron laser (XFEL) at three different pH conditions: pH 5.5, pH 6.5, and pH 8.0. Data were collected on the Inwardopen state, which is an intermediate that accumulates at high protonation of the His37 tetrad. At pH 5.5, a continuous hydrogen-bonded network of water molecules spans the vertical length of the channel, consistent with a Grotthuss mechanism model for proton transport to the His37 tetrad. This ordered solvent at pH 5.5 could act to stabilize the positive charges that build up on the gating His37 tetrad during the proton conduction cycle. The number of ordered pore waters decreases at pH 6.5 and 8.0, where the Inwardopen state is less stable. These studies provide a graphical view of the response of water to a change in charge within a restricted channel environment.

The reaction of peroxides with peroxidases oxidizes the heme iron from Fe(III) to Fe(IV)=O and a porphyrin or aromatic side chain to a cationic radical. X-ray-generated hydrated electrons rapidly reduce Fe(IV), thereby requiring very short exposures using many crystals, and, even then, some reduction cannot be avoided. The new generation of X-ray free electron lasers capable of generating intense X-rays on the tenths of femtosecond time scale enables structure determination with no reduction or X-ray damage. Here, we report the 1.5-Å crystal structure of cytochrome c peroxidase (CCP) compound I (CmpI) using data obtained with the Stanford Linear Coherent Light Source (LCLS). This structure is consistent with previous structures. Of particular importance is the active site water structure that can mediate the proton transfer reactions required for both CmpI formation and reduction of Fe(IV)=O to Fe(III)-OH. The structures indicate that a water molecule is ideally positioned to shuttle protons between an iron-linked oxygen and the active site catalytic His. We therefore have carried out both computational and kinetic studies to probe the reduction of Fe(IV)=O. Kinetic solvent isotope experiments show that the transfer of a single proton is critical in the peroxidase rate-limiting step, which is very likely the proton-coupled reduction of Fe(IV)=O to Fe(III)-OH. We also find that the pKa of the catalytic His substantially increases in CmpI, indicating that this active site His is the source of the proton required in the reduction of Fe(IV)=O to Fe(IV)-OH.

New opportunities for studying (sub)microcrystalline materials with small unit cells, both organic and inorganic, will open up when the X-ray free electron laser (XFEL) presently being constructed in Switzerland (SwissFEL) comes online in 2017. Our synchrotron-based experiments mimicking the 4%-energy-bandpass mode of the SwissFEL beam show that it will be possible to record a diffraction pattern of up to 10 randomly oriented crystals in a single snapshot, to index the resulting reflections, and to extract their intensities reliably. The crystals are destroyed with each XFEL pulse, but by combining snapshots from several sets of crystals, a complete set of data can be assembled, and crystal structures of materials that are difficult to analyze otherwise will become accessible. Even with a single shot, at least a partial analysis of the crystal structure will be possible, and with 10-50 femtosecond pulses, this offers tantalizing possibilities for time-resolved studies.

New opportunities for studying (sub)microcrystalline materials with small unit cells, both organic and inorganic, will open up when the X-ray free electron laser (XFEL) presently being constructed in Switzerland (SwissFEL) comes online in 2017. Our synchrotron-based experiments mimicking the 4%-energy-bandpass mode of the SwissFEL beam show that it will be possible to record a diffraction pattern of up to 10 randomly oriented crystals in a single snapshot, to index the resulting reflections, and to extract their intensities reliably. The crystals are destroyed with each XFEL pulse, but by combining snapshots from several sets of crystals, a complete set of data can be assembled, and crystal structures of materials that are difficult to analyze otherwise will become accessible. Even with a single shot, at least a partial analysis of the crystal structure will be possible, and with 10-50 femtosecond pulses, this offers tantalizing possibilities for time-resolved studies.

Still diffraction patterns from peptide nanocrystals with small unit cells are challenging to index using conventional methods owing to the limited number of spots and the lack of crystal orientation information for individual images. New indexing algorithms have been developed as part of the Computational Crystallography Toolbox (cctbx) to overcome these challenges. Accurate unit-cell information derived from an aggregate data set from thousands of diffraction patterns can be used to determine a crystal orientation matrix for individual images with as few as five reflections. These algorithms are potentially applicable not only to amyloid peptides but also to any set of diffraction patterns with sparse properties, such as low-resolution virus structures or high-throughput screening of still images captured by raster-scanning at synchrotron sources. As a proof of concept for this technique, successful integration of X-ray free-electron laser (XFEL) data to 2.5 Å resolution for the amyloid segment GNNQQNY from the Sup35 yeast prion is presented.

Still diffraction patterns from peptide nanocrystals with small unit cells are challenging to index using conventional methods owing to the limited number of spots and the lack of crystal orientation information for individual images. New indexing algorithms have been developed as part of the Computational Crystallography Toolbox (cctbx) to overcome these challenges. Accurate unit-cell information derived from an aggregate data set from thousands of diffraction patterns can be used to determine a crystal orientation matrix for individual images with as few as five reflections. These algorithms are potentially applicable not only to amyloid peptides but also to any set of diffraction patterns with sparse properties, such as low-resolution virus structures or high-throughput screening of still images captured by raster-scanning at synchrotron sources. As a proof of concept for this technique, successful integration of X-ray free-electron laser (XFEL) data to 2.5 Å resolution for the amyloid segment GNNQQNY from the Sup35 yeast prion is presented.

The emerging method of femtosecond crystallography (FX) may extend the diffraction resolution accessible from small radiation-sensitive crystals and provides a means to determine catalytically accurate structures of acutely radiation-sensitive metalloenzymes. Automated goniometer-based instrumentation developed for use at the Linac Coherent Light Source enabled efficient and flexible FX experiments to be performed on a variety of sample types. In the case of rod-shaped Cpl hydrogenase crystals, only five crystals and about 30 min of beam time were used to obtain the 125 still diffraction patterns used to produce a 1.6-Å resolution electron density map. For smaller crystals, high-density grids were used to increase sample throughput; 930 myoglobin crystals mounted at random orientation inside 32 grids were exposed, demonstrating the utility of this approach. Screening results from cryocooled crystals of β2-adrenoreceptor and an RNA polymerase II complex indicate the potential to extend the diffraction resolution obtainable from very radiation-sensitive samples beyond that possible with undulator-based synchrotron sources.

It has long been known that toxins produced by Bacillus thuringiensis (Bt) are stored in the bacterial cells in crystalline form. Here we describe the structure determination of the Cry3A toxin found naturally crystallized within Bt cells. When whole Bt cells were streamed into an X-ray free-electron laser beam we found that scattering from other cell components did not obscure diffraction from the crystals. The resolution limits of the best diffraction images collected from cells were the same as from isolated crystals. The integrity of the cells at the moment of diffraction is unclear; however, given the short time (∼ 5 µs) between exiting the injector to intersecting with the X-ray beam, our result is a 2.9-Å-resolution structure of a crystalline protein as it exists in a living cell. The study suggests that authentic in vivo diffraction studies can produce atomic-level structural information.