Chemical crystallography by serial femtosecond X-ray diffraction.

Inorganic-organic hybrid materials represent a large share of newly reported structures, owing to their simple synthetic routes and customizable properties1. This proliferation has led to a characterization bottleneck: many hybrid materials are obligate microcrystals with low symmetry and severe radiation sensitivity, interfering with the standard techniques of single-crystal X-ray diffraction2,3 and electron microdiffraction4-11. Here we demonstrate small-molecule serial femtosecond X-ray crystallography (smSFX) for the determination of material crystal structures from microcrystals. We subjected microcrystalline suspensions to X-ray free-electron laser radiation12,13 and obtained thousands of randomly oriented diffraction patterns. We determined unit cells by aggregating spot-finding results into high-resolution powder diffractograms. After indexing the sparse serial patterns by a graph theory approach14, the resulting datasets can be solved and refined using standard tools for single-crystal diffraction data15-17. We describe the ab initio structure solutions of mithrene (AgSePh)18-20, thiorene (AgSPh) and tethrene (AgTePh), of which the latter two were previously unknown structures. In thiorene, we identify a geometric change in the silver-silver bonding network that is linked to its divergent optoelectronic properties20. We demonstrate that smSFX can be applied as a general technique for structure determination of beam-sensitive microcrystalline materials at near-ambient temperature and pressure.

In Situ Structural Observation of a Substrate- and Peroxide-Bound High-Spin Ferric-Hydroperoxo Intermediate in the P450 Enzyme CYP121.

The P450 enzyme CYP121 from Mycobacterium tuberculosis catalyzes a carbon-carbon (C-C) bond coupling cyclization of the dityrosine substrate containing a diketopiperazine ring, cyclo(l-tyrosine-l-tyrosine) (cYY). An unusual high-spin (S = 5/2) ferric intermediate maximizes its population in less than 5 ms in the rapid freeze-quenching study of CYP121 during the shunt reaction with peracetic acid or hydrogen peroxide in acetic acid solution. We show that this intermediate can also be observed in the crystalline state by EPR spectroscopy. By developing an on-demand-rapid-mixing method for time-resolved serial femtosecond crystallography with X-ray free-electron laser (tr-SFX-XFEL) technology covering the millisecond time domain and without freezing, we structurally monitored the reaction in situ at room temperature. After a 200 ms peracetic acid reaction with the cocrystallized enzyme-substrate microcrystal slurry, a ferric-hydroperoxo intermediate is observed, and its structure is determined at 1.85 Å resolution. The structure shows a hydroperoxyl ligand between the heme and the native substrate, cYY. The oxygen atoms of the hydroperoxo are 2.5 and 3.2 Å from the iron ion. The end-on binding ligand adopts a near-side-on geometry and is weakly associated with the iron ion, causing the unusual high-spin state. This compound 0 intermediate, spectroscopically and structurally observed during the catalytic shunt pathway, reveals a unique binding mode that deviates from the end-on compound 0 intermediates in other heme enzymes. The hydroperoxyl ligand is only 2.9 Å from the bound cYY, suggesting an active oxidant role of the intermediate for direct substrate oxidation in the nonhydroxylation C-C bond coupling chemistry.

Untangling the sequence of events during the S2 → S3 transition in photosystem II and implications for the water oxidation mechanism.

In oxygenic photosynthesis, light-driven oxidation of water to molecular oxygen is carried out by the oxygen-evolving complex (OEC) in photosystem II (PS II). Recently, we reported the room-temperature structures of PS II in the four (semi)stable S-states, S1, S2, S3, and S0, showing that a water molecule is inserted during the S2 → S3 transition, as a new bridging O(H)-ligand between Mn1 and Ca. To understand the sequence of events leading to the formation of this last stable intermediate state before O2 formation, we recorded diffraction and Mn X-ray emission spectroscopy (XES) data at several time points during the S2 → S3 transition. At the electron acceptor site, changes due to the two-electron redox chemistry at the quinones, QA and QB, are observed. At the donor site, tyrosine YZ and His190 H-bonded to it move by 50 µs after the second flash, and Glu189 moves away from Ca. This is followed by Mn1 and Mn4 moving apart, and the insertion of OX(H) at the open coordination site of Mn1. This water, possibly a ligand of Ca, could be supplied via a "water wheel"-like arrangement of five waters next to the OEC that is connected by a large channel to the bulk solvent. XES spectra show that Mn oxidation (τ of ∼350 µs) during the S2 → S3 transition mirrors the appearance of OX electron density. This indicates that the oxidation state change and the insertion of water as a bridging atom between Mn1 and Ca are highly correlated.

Observations of phase changes in monoolein during high viscous injection.

Serial crystallography of membrane proteins often employs high-viscosity injectors (HVIs) to deliver micrometre-sized crystals to the X-ray beam. Typically, the carrier medium is a lipidic cubic phase (LCP) media, which can also be used to nucleate and grow the crystals. However, despite the fact that the LCP is widely used with HVIs, the potential impact of the injection process on the LCP structure has not been reported and hence is not yet well understood. The self-assembled structure of the LCP can be affected by pressure, dehydration and temperature changes, all of which occur during continuous flow injection. These changes to the LCP structure may in turn impact the results of X-ray diffraction measurements from membrane protein crystals. To investigate the influence of HVIs on the structure of the LCP we conducted a study of the phase changes in monoolein/water and monoolein/buffer mixtures during continuous flow injection, at both atmospheric pressure and under vacuum. The reservoir pressure in the HVI was tracked to determine if there is any correlation with the phase behaviour of the LCP. The results indicated that, even though the reservoir pressure underwent (at times) significant variation, this did not appear to correlate with observed phase changes in the sample stream or correspond to shifts in the LCP lattice parameter. During vacuum injection, there was a three-way coexistence of the gyroid cubic phase, diamond cubic phase and lamellar phase. During injection at atmospheric pressure, the coexistence of a cubic phase and lamellar phase in the monoolein/water mixtures was also observed. The degree to which the lamellar phase is formed was found to be strongly dependent on the co-flowing gas conditions used to stabilize the LCP stream. A combination of laboratory-based optical polarization microscopy and simulation studies was used to investigate these observations.

Anisotropic X-Ray Scattering of Transiently Oriented Water.

We study the structural dynamics of liquid water by time-resolved anisotropic x-ray scattering under the optical Kerr effect condition. In this way, we can separate the anisotropic scattering decay of 160 fs from the delayed temperature increase of ∼0.1 K occurring at 1 ps and quantify transient changes in the O-O pair distribution function. Polarizable molecular dynamics simulations reproduce well the experiment, indicating transient alignment of molecules along the electric field, which shortens the nearest-neighbor distances. In addition, analysis of the simulated water local structure provides evidence that two hypothesized fluctuating water configurations exhibit different polarizability.

Deep residual networks for crystallography trained on synthetic data.

The use of artificial intelligence to process diffraction images is challenged by the need to assemble large and precisely designed training data sets. To address this, a codebase called Resonet was developed for synthesizing diffraction data and training residual neural networks on these data. Here, two per-pattern capabilities of Resonet are demonstrated: (i) interpretation of crystal resolution and (ii) identification of overlapping lattices. Resonet was tested across a compilation of diffraction images from synchrotron experiments and X-ray free-electron laser experiments. Crucially, these models readily execute on graphics processing units and can thus significantly outperform conventional algorithms. While Resonet is currently utilized to provide real-time feedback for macromolecular crystallography users at the Stanford Synchrotron Radiation Lightsource, its simple Python-based interface makes it easy to embed in other processing frameworks. This work highlights the utility of physics-based simulation for training deep neural networks and lays the groundwork for the development of additional models to enhance diffraction collection and analysis.

Interpreting macromolecular diffraction through simulation.

This chapter discusses the use of diffraction simulators to improve experimental outcomes in macromolecular crystallography, in particular for future experiments aimed at diffuse scattering. Consequential decisions for upcoming data collection include the selection of either a synchrotron or free electron laser X-ray source, rotation geometry or serial crystallography, and fiber-coupled area detector technology vs. pixel-array detectors. The hope is that simulators will provide insights to make these choices with greater confidence. Simulation software, especially those packages focused on physics-based calculation of the diffraction, can help to predict the location, size, shape, and profile of Bragg spots and diffuse patterns in terms of an underlying physical model, including assumptions about the crystal's mosaic structure, and therefore can point to potential issues with data analysis in the early planning stages. Also, once the data are collected, simulation may offer a pathway to improve the measurement of diffraction, especially with weak data, and might help to treat problematic cases such as overlapping patterns.

Real-Time XFEL Data Analysis at SLAC and NERSC: a Trial Run of Nascent Exascale Experimental Data Analysis.

X-ray scattering experiments using Free Electron Lasers (XFELs) are a powerful tool to determine the molecular structure and function of unknown samples (such as COVID-19 viral proteins). XFEL experiments are a challenge to computing in two ways: i) due to the high cost of running XFELs, a fast turnaround time from data acquisition to data analysis is essential to make informed decisions on experimental protocols; ii) data collection rates are growing exponentially, requiring new scalable algorithms. Here we report our experiences analyzing data from two experiments at the Linac Coherent Light Source (LCLS) during September 2020. Raw data were analyzed on NERSC's Cori XC40 system, using the Superfacility paradigm: our workflow automatically moves raw data between LCLS and NERSC, where it is analyzed using the software package CCTBX. We achieved real time data analysis with a turnaround time from data acquisition to full molecular reconstruction in as little as 10 min -- sufficient time for the experiment's operators to make informed decisions. By hosting the data analysis on Cori, and by automating LCLS-NERSC interoperability, we achieved a data analysis rate which matches the data acquisition rate. Completing data analysis with 10 mins is a first for XFEL experiments and an important milestone if we are to keep up with data collection trends.

Structural dynamics in the water and proton channels of photosystem II during the S2 to S3 transition.

Light-driven oxidation of water to molecular oxygen is catalyzed by the oxygen-evolving complex (OEC) in Photosystem II (PS II). This multi-electron, multi-proton catalysis requires the transport of two water molecules to and four protons from the OEC. A high-resolution 1.89 Å structure obtained by averaging all the S states and refining the data of various time points during the S2 to S3 transition has provided better visualization of the potential pathways for substrate water insertion and proton release. Our results indicate that the O1 channel is the likely water intake pathway, and the Cl1 channel is the likely proton release pathway based on the structural rearrangements of water molecules and amino acid side chains along these channels. In particular in the Cl1 channel, we suggest that residue D1-E65 serves as a gate for proton transport by minimizing the back reaction. The results show that the water oxidation reaction at the OEC is well coordinated with the amino acid side chains and the H-bonding network over the entire length of the channels, which is essential in shuttling substrate waters and protons.

Beyond integration: modeling every pixel to obtain better structure factors from stills.

Most crystallographic data processing methods use pixel integration. In serial femtosecond crystallography (SFX), the intricate interaction between the reciprocal lattice point and the Ewald sphere is integrated out by averaging symmetrically equivalent observations recorded across a large number (104-106) of exposures. Although sufficient for generating biological insights, this approach converges slowly, and using it to accurately measure anomalous differences has proved difficult. This report presents a novel approach for increasing the accuracy of structure factors obtained from SFX data. A physical model describing all observed pixels is defined to a degree of complexity such that it can decouple the various contributions to the pixel intensities. Model dependencies include lattice orientation, unit-cell dimensions, mosaic structure, incident photon spectra and structure factor amplitudes. Maximum likelihood estimation is used to optimize all model parameters. The application of prior knowledge that structure factor amplitudes are positive quantities is included in the form of a reparameterization. The method is tested using a synthesized SFX dataset of ytterbium(III) lysozyme, where each X-ray laser pulse energy is centered at 9034 eV. This energy is 100 eV above the Yb3+ L-III absorption edge, so the anomalous difference signal is stable at 10 electrons despite the inherent energy jitter of each femtosecond X-ray laser pulse. This work demonstrates that this approach allows the determination of anomalous structure factors with very high accuracy while requiring an order-of-magnitude fewer shots than conventional integration-based methods would require to achieve similar results.

Segmented flow generator for serial crystallography at the European X-ray free electron laser.

Serial femtosecond crystallography (SFX) with X-ray free electron lasers (XFELs) allows structure determination of membrane proteins and time-resolved crystallography. Common liquid sample delivery continuously jets the protein crystal suspension into the path of the XFEL, wasting a vast amount of sample due to the pulsed nature of all current XFEL sources. The European XFEL (EuXFEL) delivers femtosecond (fs) X-ray pulses in trains spaced 100 ms apart whereas pulses within trains are currently separated by 889 ns. Therefore, continuous sample delivery via fast jets wastes >99% of sample. Here, we introduce a microfluidic device delivering crystal laden droplets segmented with an immiscible oil reducing sample waste and demonstrate droplet injection at the EuXFEL compatible with high pressure liquid delivery of an SFX experiment. While achieving ~60% reduction in sample waste, we determine the structure of the enzyme 3-deoxy-D-manno-octulosonate-8-phosphate synthase from microcrystals delivered in droplets revealing distinct structural features not previously reported.

SAD phasing of XFEL data depends critically on the error model.

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.

High-viscosity injector-based pink-beam serial crystallography of microcrystals at a synchrotron radiation source.

Since the first successful serial crystallography (SX) experiment at a synchrotron radiation source, the popularity of this approach has continued to grow showing that third-generation synchrotrons can be viable alternatives to scarce X-ray free-electron laser sources. Synchrotron radiation flux may be increased ∼100 times by a moderate increase in the bandwidth ('pink beam' conditions) at some cost to data analysis complexity. Here, we report the first high-viscosity injector-based pink-beam SX experiments. The structures of proteinase K (PK) and A2A adenosine receptor (A2AAR) were determined to resolutions of 1.8 and 4.2 Šusing 4 and 24 consecutive 100 ps X-ray pulse exposures, respectively. Strong PK data were processed using existing Laue approaches, while weaker A2AAR data required an alternative data-processing strategy. This demonstration of the feasibility presents new opportunities for time-resolved experiments with microcrystals to study structural changes in real time at pink-beam synchrotron beamlines worldwide.

Angular correlations of photons from solution diffraction at a free-electron laser encode molecular structure.

During X-ray exposure of a molecular solution, photons scattered from the same molecule are correlated. If molecular motion is insignificant during exposure, then differences in momentum transfer between correlated photons are direct measurements of the molecular structure. In conventional small- and wide-angle solution scattering, photon correlations are ignored. This report presents advances in a new biomolecular structural analysis technique, correlated X-ray scattering (CXS), which uses angular intensity correlations to recover hidden structural details from molecules in solution. Due to its intense rapid pulses, an X-ray free electron laser (XFEL) is an excellent tool for CXS experiments. A protocol is outlined for analysis of a CXS data set comprising a total of half a million X-ray exposures of solutions of small gold nanoparticles recorded at the Spring-8 Ångström Compact XFEL facility (SACLA). From the scattered intensities and their correlations, two populations of nanoparticle domains within the solution are distinguished: small twinned, and large probably non-twinned domains. It is shown analytically how, in a solution measurement, twinning information is only accessible via intensity correlations, demonstrating how CXS reveals atomic-level information from a disordered solution of like molecules.

Observation of correlated X-ray scattering at atomic resolution.

Tools to study disordered systems with local structural order, such as proteins in solution, remain limited. Such understanding is essential for e.g. rational drug design. Correlated X-ray scattering (CXS) has recently attracted new interest as a way to leverage next-generation light sources to study such disordered matter. The CXS experiment measures angular correlations of the intensity caused by the scattering of X-rays from an ensemble of identical particles, with disordered orientation and position. Averaging over 15 496 snapshot images obtained by exposing a sample of silver nanoparticles in solution to a micro-focused synchrotron radiation beam, we report on experimental efforts to obtain CXS signal from an ensemble in three dimensions. A correlation function was measured at wide angles corresponding to atomic resolution that matches theoretical predictions. These preliminary results suggest that other CXS experiments on disordered ensembles--such as proteins in solution--may be feasible in the future.