Changes in an enzyme ensemble during catalysis observed by high-resolution XFEL crystallography

Enzymes populate ensembles of structures necessary for catalysis that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography at an x-ray free electron laser to observe catalysis in a designed mutant isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations, and formation of the thioimidate intermediate selects for catalytically competent substates. The influence of cysteine ionization on the ICH ensemble is validated by determining structures of the enzyme at multiple pH values. Large molecular dynamics simulations in crystallo and time-resolved electron density maps show that Asp17 ionizes during catalysis and causes conformational changes that propagate across the dimer, permitting water to enter the active site for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics.

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.

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.

Going around the Kok cycle of the water oxidation reaction with femtosecond X-ray crystallography

The water oxidation reaction in photosystem II (PS II) produces most of the molecular oxygen in the atmosphere, which sustains life on Earth, and in this process releases four electrons and four protons that drive the downstream process of CO2 fixation in the photosynthetic apparatus. The catalytic center of PS II is an oxygen-bridged Mn4Ca complex (Mn4CaO5) which is progressively oxidized upon the absorption of light by the chlorophyll of the PS II reaction center, and the accumulation of four oxidative equivalents in the catalytic center results in the oxidation of two waters to dioxygen in the last step. The recent emergence of X-ray free-electron lasers (XFELs) with intense femtosecond X-ray pulses has opened up opportunities to visualize this reaction in PS II as it proceeds through the catalytic cycle. In this review, we summarize our recent studies of the catalytic reaction in PS II by following the structural changes along the reaction pathway via room-temperature X-ray crystallography using XFELs. The evolution of the electron density changes at the Mn complex reveals notable structural changes, including the insertion of OX from a new water molecule, which disappears on completion of the reaction, implicating it in the O-O bond formation reaction. We were also able to follow the structural dynamics of the protein coordinating with the catalytic complex and of channels within the protein that are important for substrate and product transport, revealing well orchestrated conformational changes in response to the electronic changes at the Mn4Ca cluster.

Mapping protein dynamics at high spatial resolution with temperature-jump X-ray crystallography

Understanding and controlling protein motion at atomic resolution is a hallmark challenge for structural biologists and protein engineers because conformational dynamics are essential for complex functions such as enzyme catalysis and allosteric regulation. Time-resolved crystallography offers a window into protein motions, yet without a universal perturbation to initiate conformational changes the method has been limited in scope. Here we couple a solvent-based temperature jump with time-resolved crystallography to visualize structural motions in lysozyme, a dynamic enzyme. We observed widespread atomic vibrations on the nanosecond timescale, which evolve on the submillisecond timescale into localized structural fluctuations that are coupled to the active site. An orthogonal perturbation to the enzyme, inhibitor binding, altered these dynamics by blocking key motions that allow energy to dissipate from vibrations into functional movements linked to the catalytic cycle. Because temperature jump is a universal method for perturbing molecular motion, the method demonstrated here is broadly applicable for studying protein dynamics.

Structure of a ribonucleotide reductase R2 protein radical

Aerobic ribonucleotide reductases (RNRs) initiate synthesis of DNA building blocks by generating a free radical within the R2 subunit; the radical is subsequently shuttled to the catalytic R1 subunit through proton-coupled electron transfer (PCET). We present a high-resolution room temperature structure of the class Ie R2 protein radical captured by x-ray free electron laser serial femtosecond crystallography. The structure reveals conformational reorganization to shield the radical and connect it to the translocation path, with structural changes propagating to the surface where the protein interacts with the catalytic R1 subunit. Restructuring of the hydrogen bond network, including a notably short O-O interaction of 2.41 angstroms, likely tunes and gates the radical during PCET. These structural results help explain radical handling and mobilization in RNR and have general implications for radical transfer in proteins.

Approximating lattice similarity

A method is proposed for choosing unit cells for a group of crystals so that they all appear as nearly similar as possible to a selected cell. Related unit cells with varying cell parameters or indexed with different lattice centering can be accommodated.

Changes in an Enzyme Ensemble During Catalysis Observed by High Resolution XFEL Crystallography

Enzymes populate ensembles of structures with intrinsically different catalytic proficiencies that are difficult to experimentally characterize. We use time-resolved mix-and-inject serial crystallography (MISC) at an X-ray free electron laser (XFEL) to observe catalysis in a designed mutant (G150T) isocyanide hydratase (ICH) enzyme that enhances sampling of important minor conformations. The active site exists in a mixture of conformations and formation of the thioimidate catalytic intermediate selects for catalytically competent substates. A prior proposal for active site cysteine charge-coupled conformational changes in ICH is validated by determining structures of the enzyme over a range of pH values. A combination of large molecular dynamics simulations of the enzyme in crystallo and time-resolved electron density maps shows that ionization of the general acid Asp17 during catalysis causes additional conformational changes that propagate across the dimer interface, connecting the two active sites. These ionization-linked changes in the ICH conformational ensemble permit water to enter the active site in a location that is poised for intermediate hydrolysis. ICH exhibits a tight coupling between ionization of active site residues and catalysis-activated protein motions, exemplifying a mechanism of electrostatic control of enzyme dynamics.

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.

XFEL Microcrystallography of Self-Assembling Silver n-Alkanethiolates

New synthetic hybrid materials and their increasing complexity have placed growing demands on crystal growth for single-crystal X-ray diffraction analysis. Unfortunately, not all chemical systems are conducive to the isolation of single crystals for traditional characterization. Here, small-molecule serial femtosecond crystallography (smSFX) at atomic resolution (0.833 Å) is employed to characterize microcrystalline silver n-alkanethiolates with various alkyl chain lengths at X-ray free electron laser facilities, resolving long-standing controversies regarding the atomic connectivity and odd-even effects of layer stacking. smSFX provides high-quality crystal structures directly from the powder of the true unknowns, a capability that is particularly useful for systems having notoriously small or defective crystals. We present crystal structures of silver n-butanethiolate (C4), silver n-hexanethiolate (C6), and silver n-nonanethiolate (C9). We show that an odd-even effect originates from the orientation of the terminal methyl group and its role in packing efficiency. We also propose a secondary odd-even effect involving multiple mosaic blocks in the crystals containing even-numbered chains, identified by selected-area electron diffraction measurements. We conclude with a discussion of the merits of the synthetic preparation for the preparation of microdiffraction specimens and compare the long-range order in these crystals to that of self-assembled monolayers.

Water Networks in Photosystem II Using Crystalline Molecular Dynamics Simulations and Room-Temperature XFEL Serial Crystallography

Structural dynamics of water and its hydrogen-bonding networks play an important role in enzyme function via the transport of protons, ions, and substrates. To gain insights into these mechanisms in the water oxidation reaction in Photosystem II (PS II), we have performed crystalline molecular dynamics (MD) simulations of the dark-stable S1 state. Our MD model consists of a full unit cell with 8 PS II monomers in explicit solvent (861 894 atoms), enabling us to compute the simulated crystalline electron density and to compare it directly with the experimental density from serial femtosecond X-ray crystallography under physiological temperature collected at X-ray free electron lasers (XFELs). The MD density reproduced the experimental density and water positions with high fidelity. The detailed dynamics in the simulations provided insights into the mobility of water molecules in the channels beyond what can be interpreted from experimental B-factors and electron densities alone. In particular, the simulations revealed fast, coordinated exchange of waters at sites where the density is strong, and water transport across the bottleneck region of the channels where the density is weak. By computing MD hydrogen and oxygen maps separately, we developed a novel Map-based Acceptor-Donor Identification (MADI) technique that yields information which helps to infer hydrogen-bond directionality and strength. The MADI analysis revealed a series of hydrogen-bond wires emanating from the Mn cluster through the Cl1 and O4 channels; such wires might provide pathways for proton transfer during the reaction cycle of PS II. Our simulations provide an atomistic picture of the dynamics of water and hydrogen-bonding networks in PS II, with implications for the specific role of each channel in the water oxidation reaction.

Self-Supervised Deep Learning for Model Correction in the Computational Crystallography Toolbox

The Computational Crystallography Toolbox (cctbx) is open-source software that allows for processing of crystallographic data, including from serial femtosecond crystallography (SFX), for macromolecular structure determination. We aim to use the modules in cctbx to determine the oxidation state of individual metal atoms in a macromolecule. Changes in oxidation state are reflected in small shifts of the atom's X-ray absorption edge. These energy shifts can be extracted from the diffraction images recorded in serial femtosecond crystallography, given knowledge of a forward physics model. However, as the diffraction changes only slightly due to the absorption edge shift, inaccuracies in the forward physics model make it extremely challenging to observe the oxidation state. In this work, we describe the potential impact of using self-supervised deep learning to correct the scientific model in cctbx and provide uncertainty quantification. We provide code for forward model simulation and data analysis, built from cctbx modules, at https://github.com/gigantocypris/SPREAD, which can be integrated with machine learning. We describe open questions in algorithm development to help spur advances through dialog between crystallographers and machine learning researchers. New methods could help elucidate charge transfer processes in many reactions, including key events in photosynthesis.

Structural evidence for intermediates during O2 formation in photosystem II

In natural photosynthesis, the light-driven splitting of water into electrons, protons and molecular oxygen forms the first step of the solar-to-chemical energy conversion process. The reaction takes place in photosystem II, where the Mn4CaO5 cluster first stores four oxidizing equivalents, the S0 to S4 intermediate states in the Kok cycle, sequentially generated by photochemical charge separations in the reaction center and then catalyzes the O-O bond formation chemistry1-3. Here, we report room temperature snapshots by serial femtosecond X-ray crystallography to provide structural insights into the final reaction step of Kok's photosynthetic water oxidation cycle, the S3→[S4]→S0 transition where O2 is formed and Kok's water oxidation clock is reset. Our data reveal a complex sequence of events, which occur over micro- to milliseconds, comprising changes at the Mn4CaO5 cluster, its ligands and water pathways as well as controlled proton release through the hydrogen-bonding network of the Cl1 channel. Importantly, the extra O atom Ox, which was introduced as a bridging ligand between Ca and Mn1 during the S2→S3 transition4-6, disappears or relocates in parallel with Yz reduction starting at approximately 700 μs after the third flash. The onset of O2 evolution, as indicated by the shortening of the Mn1-Mn4 distance, occurs at around 1,200 μs, signifying the presence of a reduced intermediate, possibly a bound peroxide.

Capturing the sequence of events during the water oxidation reaction in photosynthesis using XFELs

Ever since the discovery that Mn was required for oxygen evolution in plants by Pirson in 1937 and the period-four oscillation in flash-induced oxygen evolution by Joliot and Kok in the 1970s, understanding of this process has advanced enormously using state-of-the-art methods. The most recent in this series of innovative techniques was the introduction of X-ray free-electron lasers (XFELs) a decade ago, which led to another quantum leap in the understanding in this field, by enabling operando X-ray structural and X-ray spectroscopy studies at room temperature. This review summarizes the current understanding of the structure of Photosystem II (PS II) and its catalytic centre, the Mn4 CaO5 complex, in the intermediate Si (i = 0-4)-states of the Kok cycle, obtained using XFELs.

Redox-controlled reorganization and flavin strain within the ribonucleotide reductase R2b-NrdI complex monitored by serial femtosecond crystallography

Redox reactions are central to biochemistry and are both controlled by and induce protein structural changes. Here, we describe structural rearrangements and crosstalk within the Bacillus cereus ribonucleotide reductase R2b-NrdI complex, a di-metal carboxylate-flavoprotein system, as part of the mechanism generating the essential catalytic free radical of the enzyme. Femtosecond crystallography at an X-ray free electron laser was utilized to obtain structures at room temperature in defined redox states without suffering photoreduction. Together with density functional theory calculations, we show that the flavin is under steric strain in the R2b-NrdI protein complex, likely tuning its redox properties to promote superoxide generation. Moreover, a binding site in close vicinity to the expected flavin O₂ interaction site is observed to be controlled by the redox state of the flavin and linked to the channel proposed to funnel the produced superoxide species from NrdI to the di-manganese site in protein R2b. These specific features are coupled to further structural changes around the R2b-NrdI interaction surface. The mechanistic implications for the control of reactive oxygen species and radical generation in protein R2b are discussed.

De novo determination of mosquitocidal Cry11Aa and Cry11Ba structures from naturally-occurring nanocrystals

Cry11Aa and Cry11Ba are the two most potent toxins produced by mosquitocidal Bacillus thuringiensis subsp. israelensis and jegathesan, respectively. The toxins naturally crystallize within the host; however, the crystals are too small for structure determination at synchrotron sources. Therefore, we applied serial femtosecond crystallography at X-ray free electron lasers to in vivo-grown nanocrystals of these toxins. The structure of Cry11Aa was determined de novo using the single-wavelength anomalous dispersion method, which in turn enabled the determination of the Cry11Ba structure by molecular replacement. The two structures reveal a new pattern for in vivo crystallization of Cry toxins, whereby each of their three domains packs with a symmetrically identical domain, and a cleavable crystal packing motif is located within the protoxin rather than at the termini. The diversity of in vivo crystallization patterns suggests explanations for their varied levels of toxicity and rational approaches to improve these toxins for mosquito control.

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.

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 S₂ to S₃ 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.

Room temperature XFEL crystallography reveals asymmetry in the vicinity of the two phylloquinones in photosystem I

Photosystem I (PS I) has a symmetric structure with two highly similar branches of pigments at the center that are involved in electron transfer, but shows very different efficiency along the two branches. We have determined the structure of cyanobacterial PS I at room temperature (RT) using femtosecond X-ray pulses from an X-ray free electron laser (XFEL) that shows a clear expansion of the entire protein complex in the direction of the membrane plane, when compared to previous cryogenic structures. This trend was observed by complementary datasets taken at multiple XFEL beamlines. In the RT structure of PS I, we also observe conformational differences between the two branches in the reaction center around the secondary electron acceptors A_1A and A_1B. The π-stacked Phe residues are rotated with a more parallel orientation in the A-branch and an almost perpendicular confirmation in the B-branch, and the symmetry breaking PsaB-Trp673 is tilted and further away from A_1A. These changes increase the asymmetry between the branches and may provide insights into the preferential directionality of electron transfer.

X-ray free-electron laser studies reveal correlated motion during isopenicillin N synthase catalysis

Isopenicillin N synthase (IPNS) catalyzes the unique reaction of l-δ-(α-aminoadipoyl)-l-cysteinyl-d-valine (ACV) with dioxygen giving isopenicillin N (IPN), the precursor of all natural penicillins and cephalosporins. X-ray free-electron laser studies including time-resolved crystallography and emission spectroscopy reveal how reaction of IPNS:Fe(II):ACV with dioxygen to yield an Fe(III) superoxide causes differences in active site volume and unexpected conformational changes that propagate to structurally remote regions. Combined with solution studies, the results reveal the importance of protein dynamics in regulating intermediate conformations during conversion of ACV to IPN. The results have implications for catalysis by multiple IPNS-related oxygenases, including those involved in the human hypoxic response, and highlight the power of serial femtosecond crystallography to provide insight into long-range enzyme dynamics during reactions presently impossible for nonprotein catalysts.

An on-demand, drop-on-drop method for studying enzyme catalysis by serial crystallography

Serial femtosecond crystallography has opened up many new opportunities in structural biology. In recent years, several approaches employing light-inducible systems have emerged to enable time-resolved experiments that reveal protein dynamics at high atomic and temporal resolutions. However, very few enzymes are light-dependent, whereas macromolecules requiring ligand diffusion into an active site are ubiquitous. In this work we present a drop-on-drop sample delivery system that enables the study of enzyme-catalyzed reactions in microcrystal slurries. The system delivers ligand solutions in bursts of multiple picoliter-sized drops on top of a larger crystal-containing drop inducing turbulent mixing and transports the mixture to the X-ray interaction region with temporal resolution. We demonstrate mixing using fluorescent dyes, numerical simulations and time-resolved serial femtosecond crystallography, which show rapid ligand diffusion through microdroplets. The drop-on-drop method has the potential to be widely applicable to serial crystallography studies, particularly of enzyme reactions with small molecule substrates.

Challenges in solving structures from radiation-damaged tomograms of protein nanocrystals assessed by simulation

Structure-determination methods are needed to resolve the atomic details that underlie protein function. X-ray crystallography has provided most of our knowledge of protein structure, but is constrained by the need for large, well ordered crystals and the loss of phase information. The rapidly developing methods of serial femtosecond crystallography, micro-electron diffraction and single-particle reconstruction circumvent the first of these limitations by enabling data collection from nanocrystals or purified proteins. However, the first two methods also suffer from the phase problem, while many proteins fall below the molecular-weight threshold required for single-particle reconstruction. Cryo-electron tomography of protein nanocrystals has the potential to overcome these obstacles of mainstream structure-determination methods. Here, a data-processing scheme is presented that combines routines from X-ray crystallography and new algorithms that have been developed to solve structures from tomograms of nanocrystals. This pipeline handles image-processing challenges specific to tomographic sampling of periodic specimens and is validated using simulated crystals. The tolerance of this workflow to the effects of radiation damage is also assessed. The simulations indicate a trade-off between a wider tilt range to facilitate merging data from multiple tomograms and a smaller tilt increment to improve phase accuracy. Since phase errors, but not merging errors, can be overcome with additional data sets, these results recommend distributing the dose over a wide angular range rather than using a finer sampling interval to solve the protein structure.

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.

High-Resolution XFEL Structure of the Soluble Methane Monooxygenase Hydroxylase Complex with its Regulatory Component at Ambient Temperature in Two Oxidation States

Soluble methane monooxygenase (sMMO) is a multicomponent metalloenzyme that catalyzes the conversion of methane to methanol at ambient temperature using a nonheme, oxygen-bridged dinuclear iron cluster in the active site. Structural changes in the hydroxylase component (sMMOH) containing the diiron cluster caused by complex formation with a regulatory component (MMOB) and by iron reduction are important for the regulation of O₂ activation and substrate hydroxylation. Structural studies of metalloenzymes using traditional synchrotron-based X-ray crystallography are often complicated by partial X-ray-induced photoreduction of the metal center, thereby obviating determination of the structure of the enzyme in pure oxidation states. Here, microcrystals of the sMMOH:MMOB complex from Methylosinus trichosporium OB3b were serially exposed to X-ray free electron laser (XFEL) pulses, where the ≤35 fs duration of exposure of an individual crystal yields diffraction data before photoreduction-induced structural changes can manifest. Merging diffraction patterns obtained from thousands of crystals generates radiation damage-free, 1.95 Å resolution crystal structures for the fully oxidized and fully reduced states of the sMMOH:MMOB complex for the first time. The results provide new insight into the manner by which the diiron cluster and the active site environment are reorganized by the regulatory protein component in order to enhance the steps of oxygen activation and methane oxidation. This study also emphasizes the value of XFEL and serial femtosecond crystallography (SFX) methods for investigating the structures of metalloenzymes with radiation sensitive metal active sites.

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.

Artificial Iron Proteins: Modeling the Active Sites in Non-Heme Dioxygenases

An important class of non-heme dioxygenases contains a conserved Fe binding site that consists of a 2-His-1-carboxylate facial triad. Results from structural biology show that, in the resting state, these proteins are six-coordinate with aqua ligands occupying the remaining three coordination sites. We have utilized biotin-streptavidin (Sav) technology to design new artificial Fe proteins (ArMs) that have many of the same structural features found within active sites of these non-heme dioxygenases. An Sav variant was isolated that contains the S112E mutation, which installed a carboxylate side chain in the appropriate position to bind to a synthetic FeII complex confined within Sav. Structural studies using X-ray diffraction (XRD) methods revealed a facial triad binding site that is composed of two N donors from the biotinylated ligand and the monodentate coordination of the carboxylate from S112E. Two aqua ligands complete the primary coordination sphere of the FeII center with both involved in hydrogen bond networks within Sav. The corresponding FeIII protein was also prepared and structurally characterized to show a six-coordinate complex with two exogenous acetato ligands. The FeIII protein was further shown to bind an exogenous azido ligand through replacement of one acetato ligand. Spectroscopic studies of the ArMs in solution support the results found by XRD.

Comparing serial X-ray crystallography and microcrystal electron diffraction (MicroED) as methods for routine structure determination from small macromolecular crystals

Innovative new crystallographic methods are facilitating structural studies from ever smaller crystals of biological macromolecules. In particular, serial X-ray crystallography and microcrystal electron diffraction (MicroED) have emerged as useful methods for obtaining structural information from crystals on the nanometre to micrometre scale. Despite the utility of these methods, their implementation can often be difficult, as they present many challenges that are not encountered in traditional macromolecular crystallography experiments. Here, XFEL serial crystallography experiments and MicroED experiments using batch-grown microcrystals of the enzyme cyclophilin A are described. The results provide a roadmap for researchers hoping to design macromolecular microcrystallography experiments, and they highlight the strengths and weaknesses of the two methods. Specifically, we focus on how the different physical conditions imposed by the sample-preparation and delivery methods required for each type of experiment affect the crystal structure of the enzyme.

Towards the spatial resolution of metalloprotein charge states by detailed modeling of XFEL crystallographic diffraction

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.

Converting three-space matrices to equivalent six-space matrices for Delone scalars in S6

Given a matrix for transforming vectors in the three-space of unit-cell edge vectors, the corresponding matrix to transform vectors in the six-space of Delone scalars is derived.

The transformations from the primitive cells of the centered Bravais lattices to the corresponding centered cells have conventionally been listed as three-by-three matrices that transform three-space lattice vectors. Using those three-by-three matrices when working in the six-dimensional space of lattices represented as Selling scalars as used in Delone (Delaunay) reduction, one could transform to the three-space representation, apply the three-by-three matrices and then back-transform to the six-space representation, but it is much simpler to have the equivalent six-by-six matrices and apply them directly. The general form of the transformation from the three-space matrix to the corresponding matrix operating on Selling scalars (expressed in space S⁶) is derived, and the particular S⁶matrices for the centered Delone types are listed. (Note: in his later publications, Boris Delaunay used the Russian version of his surname, Delone.)

Best practices for high data-rate macromolecular crystallography (HDRMX)

In macromolecular crystallography, higher flux, smaller beams, and faster detectors open the door to experiments with very large numbers of very small samples that can reveal polymorphs and dynamics but require re-engineering of approaches to the clustering of images both at synchrotrons and XFELs (X-ray free electron lasers). The need for the management of orders of magnitude more images and limitations of file systems favor a transition from simple one-file-per-image systems such as CBF to image container systems such as HDF5. This further increases the load on computers and networks and requires a re-examination of the presentation of metadata. In this paper, we discuss three important components of this problem-improved approaches to the clustering of images to better support experiments on polymorphs and dynamics, recent and upcoming changes in metadata for Eiger images, and software to rapidly validate images in the revised Eiger format.

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.

A space for lattice representation and clustering

Algorithms for defining the difference between two lattices are described. They are based on the work of Selling and Delone (Delaunay).

Algorithms for quantifying the differences between two lattices are used for Bravais lattice determination, database lookup for unit cells to select candidates for molecular replacement, and recently for clustering to group together images from serial crystallography. It is particularly desirable for the differences between lattices to be computed as a perturbation-stable metric, i.e. as distances that satisfy the triangle inequality, so that standard tree-based nearest-neighbor algorithms can be used, and for which small changes in the lattices involved produce small changes in the distances computed. A perturbation-stable metric space related to the reduction algorithm of Selling and to the Bravais lattice determination methods of Delone is described. Two ways of representing the space, as six-dimensional real vectors or equivalently as three-dimensional complex vectors, are presented and applications of these metrics are discussed. (Note: in his later publications, Boris Delaunay used the Russian version of his surname, Delone.)

Selling reduction versus Niggli reduction for crystallographic lattices

The unit-cell reduction described by Selling and used by Delone (whose early publications were under the spelling Delaunay) is explained in a simple form. The transformations needed to implement the reduction are listed. The simplicity of this reduction contrasts with the complexity of Niggli reduction.

Structures of the intermediates of Kok's photosynthetic water oxidation clock

Inspired by the period-four oscillation in flash-induced oxygen evolution of photosystem II discovered by Joliot in 1969, Kok performed additional experiments and proposed a five-state kinetic model for photosynthetic oxygen evolution, known as Kok's S-state clock or cycle^1,2. The model comprises four (meta)stable intermediates (S₀, S₁, S₂ and S₃) and one transient S₄ state, which precedes dioxygen formation occurring in a concerted reaction from two water-derived oxygens bound at an oxo-bridged tetra manganese calcium (Mn₄CaO₅) cluster in the oxygen-evolving complex^3-7. This reaction is coupled to the two-step reduction and protonation of the mobile plastoquinone Q_B at the acceptor side of PSII. Here, using serial femtosecond X-ray crystallography and simultaneous X-ray emission spectroscopy with multi-flash visible laser excitation at room temperature, we visualize all (meta)stable states of Kok's cycle as high-resolution structures (2.04-2.08 Å). In addition, we report structures of two transient states at 150 and 400 µs, revealing notable structural changes including the binding of one additional 'water', Ox, during the S₂→S₃ state transition. Our results suggest that one water ligand to calcium (W3) is directly involved in substrate delivery. The binding of the additional oxygen Ox in the S₃ state between Ca and Mn1 supports O-O bond formation mechanisms involving O5 as one substrate, where Ox is either the other substrate oxygen or is perfectly positioned to refill the O5 position during O₂ release. Thus, our results exclude peroxo-bond formation in the S₃ state, and the nucleophilic attack of W3 onto W2 is unlikely.

Free-electron laser data for multiple-particle fluctuation scattering analysis

Fluctuation X-ray scattering (FXS) is an emerging experimental technique in which solution scattering data are collected using X-ray exposures below rotational diffusion times, resulting in angularly anisotropic X-ray snapshots that provide several orders of magnitude more information than traditional solution scattering data. Such experiments can be performed using the ultrashort X-ray pulses provided by a free-electron laser source, allowing one to collect a large number of diffraction patterns in a relatively short time. Here, we describe a test data set for FXS, obtained at the Linac Coherent Light Source, consisting of close to 100 000 multi-particle diffraction patterns originating from approximately 50 to 200 Paramecium Bursaria Chlorella virus particles per snapshot. In addition to the raw data, a selection of high-quality pre-processed diffraction patterns and a reference SAXS profile are provided.

Towards characterization of photo-excited electron transfer and catalysis in natural and artificial systems using XFELs

The ultra-bright femtosecond X-ray pulses provided by X-ray Free Electron Lasers (XFELs) open capabilities for studying the structure and dynamics of a wide variety of biological and inorganic systems beyond what is possible at synchrotron sources. Although the structure and chemistry at the catalytic sites have been studied intensively in both biological and inorganic systems, a full understanding of the atomic-scale chemistry requires new approaches beyond the steady state X-ray crystallography and X-ray spectroscopy at cryogenic temperatures. Following the dynamic changes in the geometric and electronic structure at ambient conditions, while overcoming X-ray damage to the redox active catalytic center, is key for deriving reaction mechanisms. Such studies become possible by using the intense and ultra-short femtosecond X-ray pulses from an XFEL, where sample is probed before it is damaged. We have developed methodology for simultaneously collecting X-ray diffraction data and X-ray emission spectra, using an energy dispersive spectrometer, at ambient conditions, and used this approach to study the room temperature structure and intermediate states of the photosynthetic water oxidizing metallo-protein, photosystem II. Moreover, we have also used this setup to simultaneously collect the X-ray emission spectra from multiple metals to follow the ultrafast dynamics of light-induced charge transfer between multiple metal sites. A Mn-Ti containing system was studied at an XFEL to demonstrate the efficacy and potential of this method.

Improving signal strength in serial crystallography with DIALS geometry refinement

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.

A convolutional neural network-based screening tool for X-ray serial crystallography

A new tool is introduced for screening macromolecular X-ray crystallography diffraction images produced at an X-ray free-electron laser light source. Based on a data-driven deep learning approach, the proposed tool executes a convolutional neural network to detect Bragg spots. Automatic image processing algorithms described can enable the classification of large data sets, acquired under realistic conditions consisting of noisy data with experimental artifacts. Outcomes are compared for different data regimes, including samples from multiple instruments and differing amounts of training data for neural network optimization.

DIALS: implementation and evaluation of a new integration package

The DIALS project is a collaboration between Diamond Light Source, Lawrence Berkeley National Laboratory and CCP4 to develop a new software suite for the analysis of crystallographic X-ray diffraction data, initially encompassing spot finding, indexing, refinement and integration. The design, core algorithms and structure of the software are introduced, alongside results from the analysis of data from biological and chemical crystallography experiments.

High-speed fixed-target serial virus crystallography

We report a method for serial X-ray crystallography at X-ray free-electron lasers (XFELs), which allows for full use of the current 120-Hz repetition rate of the Linear Coherent Light Source (LCLS). Using a micropatterned silicon chip in combination with the high-speed Roadrunner goniometer for sample delivery, we were able to determine the crystal structures of the picornavirus bovine enterovirus 2 (BEV2) and the cytoplasmic polyhedrosis virus type 18 polyhedrin, with total data collection times of less than 14 and 10 min, respectively. Our method requires only micrograms of sample and should therefore broaden the applicability of serial femtosecond crystallography to challenging projects for which only limited sample amounts are available. By synchronizing the sample exchange to the XFEL repetition rate, our method allows for most efficient use of the limited beam time available at XFELs and should enable a substantial increase in sample throughput at these facilities.

Flow-aligned, single-shot fiber diffraction using a femtosecond X-ray free-electron laser

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.