Influence of pump laser fluence on ultrafast myoglobin structural dynamics

High-intensity femtosecond pulses from an X-ray free-electron laser enable pump-probe experiments for the investigation of electronic and nuclear changes during light-induced reactions. On timescales ranging from femtoseconds to milliseconds and for a variety of biological systems, time-resolved serial femtosecond crystallography (TR-SFX) has provided detailed structural data for light-induced isomerization, breakage or formation of chemical bonds and electron transfer1,2. However, all ultrafast TR-SFX studies to date have employed such high pump laser energies that nominally several photons were absorbed per chromophore3-17. As multiphoton absorption may force the protein response into non-physiological pathways, it is of great concern18,19 whether this experimental approach20 allows valid conclusions to be drawn vis-à-vis biologically relevant single-photon-induced reactions18,19. Here we describe ultrafast pump-probe SFX experiments on the photodissociation of carboxymyoglobin, showing that different pump laser fluences yield markedly different results. In particular, the dynamics of structural changes and observed indicators of the mechanistically important coherent oscillations of the Fe-CO bond distance (predicted by recent quantum wavepacket dynamics21) are seen to depend strongly on pump laser energy, in line with quantum chemical analysis. Our results confirm both the feasibility and necessity of performing ultrafast TR-SFX pump-probe experiments in the linear photoexcitation regime. We consider this to be a starting point for reassessing both the design and the interpretation of ultrafast TR-SFX pump-probe experiments20 such that mechanistically relevant insight emerges.

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.

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.

Large language models generate functional protein sequences across diverse families

Deep-learning language models have shown promise in various biotechnological applications, including protein design and engineering. Here we describe ProGen, a language model that can generate protein sequences with a predictable function across large protein families, akin to generating grammatically and semantically correct natural language sentences on diverse topics. The model was trained on 280 million protein sequences from >19,000 families and is augmented with control tags specifying protein properties. ProGen can be further fine-tuned to curated sequences and tags to improve controllable generation performance of proteins from families with sufficient homologous samples. Artificial proteins fine-tuned to five distinct lysozyme families showed similar catalytic efficiencies as natural lysozymes, with sequence identity to natural proteins as low as 31.4%. ProGen is readily adapted to diverse protein families, as we demonstrate with chorismate mutase and malate dehydrogenase.

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.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and x-ray diffraction at room temperature

The nonstructural protein 3 (NSP3) macrodomain of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (Mac1) removes adenosine diphosphate (ADP) ribosylation posttranslational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the coronavirus disease 2019 pandemic. Here, we determined neutron and x-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose-bound states. We characterize extensive solvation in the Mac1 active site and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase the potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a reevaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

The mechanisms of catalysis and ligand binding for the SARS-CoV-2 NSP3 macrodomain from neutron and X-ray diffraction at room temperature

The NSP3 macrodomain of SARS CoV 2 (Mac1) removes ADP-ribosylation post-translational modifications, playing a key role in the immune evasion capabilities of the virus responsible for the COVID-19 pandemic. Here, we determined neutron and X-ray crystal structures of the SARS-CoV-2 NSP3 macrodomain using multiple crystal forms, temperatures, and pHs, across the apo and ADP-ribose-bound states. We characterize extensive solvation in the Mac1 active site, and visualize how water networks reorganize upon binding of ADP-ribose and non-native ligands, inspiring strategies for displacing waters to increase potency of Mac1 inhibitors. Determining the precise orientations of active site water molecules and the protonation states of key catalytic site residues by neutron crystallography suggests a catalytic mechanism for coronavirus macrodomains distinct from the substrate-assisted mechanism proposed for human MacroD2. These data provoke a re-evaluation of macrodomain catalytic mechanisms and will guide the optimization of Mac1 inhibitors.

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.

Hydroxyl radical mediated damage of proteins in low oxygen solution investigated using X-ray footprinting mass spectrometry

In the method of X-ray footprinting mass spectrometry (XFMS), proteins at micromolar concentration in solution are irradiated with a broadband X-ray source, and the resulting hydroxyl radical modifications are characterized using liquid chromatography mass spectrometry to determine sites of solvent accessibility. These data are used to infer structural changes in proteins upon interaction with other proteins, folding, or ligand binding. XFMS is typically performed under aerobic conditions; dissolved molecular oxygen in solution is necessary in many, if not all, the hydroxyl radical modifications that are generally reported. In this study we investigated the result of X-ray induced modifications to three different proteins under aerobic versus low oxygen conditions, and correlated the extent of damage with dose calculations. We observed a concentration-dependent protecting effect at higher protein concentration for a given X-ray dose. For the typical doses used in XFMS experiments there was minimal X-ray induced aggregation and fragmentation, but for higher doses we observed formation of covalent higher molecular weight oligomers, as well as fragmentation, which was affected by the amount of dissolved oxygen in solution. The higher molecular weight products in the form of dimers, trimers, and tetramers were present in all sample preparations, and, upon X-ray irradiation, these oligomers became non-reducible as seen in SDS-PAGE. The results provide an important contribution to the large body of X-ray radiation damage literature in structural biology research, and will specifically help inform the future planning of XFMS, and well as X-ray crystallography and small-angle X-ray scattering experiments.

Structures of two novel crystal forms of Aspergillus oryzae alpha amylase (taka-amylase)

The structures of Aspergillus oryzae α-amylase were determined in a tetragonal crystal, having one molecule as asymmetric unit, and a monoclinic crystal with two molecules as asymmetric unit. Both crystal forms were obtained from trace contaminants of an old commercial lipase preparation. Structures were determined and refined to 1.65 Å and 1.43 Å resolution respectively. The latter crystal has a non-crystallographic (NCS) twofold axis within the asymmetric unit. Glycosylation at Asn197 is evident, and in the tetragonal crystal can be seen to include three, partially disordered sugar residues following the initial N-acetyl glucosamine (NAG). Superposition of the tetragonal crystal model on the α-amylases from Bacillus subtilis (PDB:1BAG), pig pancreas (PDB:3L2L), and barley (PDB:1AMY), show a high degree of coincidence, particularly for the (β/α)₈-barrel domains, and especially within the active site. Using this structural agreement between amylases, we extrapolated the binding model of a six residue, limit dextrin found in pig pancreas α-amylase to the A. oryzae enzyme model, which predicts substrate interacting amino acid residues.

Fragment binding to the Nsp3 macrodomain of SARS-CoV-2 identified through crystallographic screening and computational docking

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) macrodomain within the nonstructural protein 3 counteracts host-mediated antiviral adenosine diphosphate-ribosylation signaling. This enzyme is a promising antiviral target because catalytic mutations render viruses nonpathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of 2533 diverse fragments resulted in 214 unique macrodomain-binders. An additional 60 molecules were selected from docking more than 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several fragment hits were confirmed by solution binding using three biophysical techniques (differential scanning fluorimetry, homogeneous time-resolved fluorescence, and isothermal titration calorimetry). The 234 fragment structures explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

Fragment Binding to the Nsp3 Macrodomain of SARS-CoV-2 Identified Through Crystallographic Screening and Computational Docking

The SARS-CoV-2 macrodomain (Mac1) within the non-structural protein 3 (Nsp3) counteracts host-mediated antiviral ADP-ribosylation signalling. This enzyme is a promising antiviral target because catalytic mutations render viruses non-pathogenic. Here, we report a massive crystallographic screening and computational docking effort, identifying new chemical matter primarily targeting the active site of the macrodomain. Crystallographic screening of diverse fragment libraries resulted in 214 unique macrodomain-binding fragments, out of 2,683 screened. An additional 60 molecules were selected from docking over 20 million fragments, of which 20 were crystallographically confirmed. X-ray data collection to ultra-high resolution and at physiological temperature enabled assessment of the conformational heterogeneity around the active site. Several crystallographic and docking fragment hits were validated for solution binding using three biophysical techniques (DSF, HTRF, ITC). Overall, the 234 fragment structures presented explore a wide range of chemotypes and provide starting points for development of potent SARS-CoV-2 macrodomain inhibitors.

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.

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.

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.

Covalent Modification and Regulation of the Nuclear Receptor Nurr1 by a Dopamine Metabolite

Nurr1, a nuclear receptor essential for the development, maintenance, and survival of midbrain dopaminergic neurons, is a potential therapeutic target for Parkinson's disease, a neurological disorder characterized by the degeneration of these same neurons. Efforts to identify Nurr1 agonists have been hampered by the recognition that it lacks several classic regulatory elements of nuclear receptor function, including the canonical ligand-binding pocket. Here we report that the dopamine metabolite 5,6-dihydroxyindole (DHI) binds directly to and modulates the activity of Nurr1. Using biophysical assays and X-ray crystallography, we show that DHI binds to the ligand-binding domain within a non-canonical pocket, forming a covalent adduct with Cys566. In cultured cells and zebrafish, DHI stimulates Nurr1 activity, including the transcription of target genes underlying dopamine homeostasis. These findings suggest avenues for developing synthetic Nurr1 ligands to ameliorate the symptoms and progression of Parkinson's disease.

Challenge data set for macromolecular multi-microcrystallography

A synthetic data set demonstrating a particularly challenging case of indexing ambiguity in the context of radiation damage was generated. This set shall serve as a standard benchmark and reference point for the ongoing development of new methods and new approaches to robust structure solution when single-crystal methods are insufficient. Of the 100 short wedges of data, only the first 36 are currently necessary to solve the structure by `cheating', or using the correct reference structure as a guide. The total wall-clock time and number of crystals required to solve the structure without cheating is proposed as a metric for the efficacy and efficiency of a given multi-crystal automation pipeline.

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.

Guide-bound structures of an RNA-targeting A-cleaving CRISPR-Cas13a enzyme

CRISPR adaptive immune systems protect bacteria from infections by deploying CRISPR RNA (crRNA)-guided enzymes to recognize and cut foreign nucleic acids. Type VI-A CRISPR-Cas systems include the Cas13a enzyme, an RNA-activated RNase capable of crRNA processing and single-stranded RNA degradation upon target-transcript binding. Here we present the 2.0-Å resolution crystal structure of a crRNA-bound Lachnospiraceae bacterium Cas13a (LbaCas13a), representing a recently discovered Cas13a enzyme subtype. This structure and accompanying biochemical experiments define the Cas13a catalytic residues that are directly responsible for crRNA maturation. In addition, the orientation of the foreign-derived target-RNA-specifying sequence in the protein interior explains the conformational gating of Cas13a nuclease activation. These results describe how Cas13a enzymes generate functional crRNAs and how catalytic activity is blocked before target-RNA recognition, with implications for both bacterial immunity and diagnostic applications.

Se-SAD serial femtosecond crystallography datasets from selenobiotinyl-streptavidin

We provide a detailed description of selenobiotinyl-streptavidin (Se-B SA) co-crystal datasets recorded using the Coherent X-ray Imaging (CXI) instrument at the Linac Coherent Light Source (LCLS) for selenium single-wavelength anomalous diffraction (Se-SAD) structure determination. Se-B SA was chosen as the model system for its high affinity between biotin and streptavidin where the sulfur atom in the biotin molecule (C₁₀H₁₆N₂O₃S) is substituted with selenium. The dataset was collected at three different transmissions (100, 50, and 10%) using a serial sample chamber setup which allows for two sample chambers, a front chamber and a back chamber, to operate simultaneously. Diffraction patterns from Se-B SA were recorded to a resolution of 1.9 Å. The dataset is publicly available through the Coherent X-ray Imaging Data Bank (CXIDB) and also on LCLS compute nodes as a resource for research and algorithm development.

Advances in X-ray free electron laser (XFEL) diffraction data processing applied to the crystal structure of the synaptotagmin-1 / SNARE complex

X-ray free electron lasers (XFELs) reduce the effects of radiation damage on macromolecular diffraction data and thereby extend the limiting resolution. Previously, we adapted classical post-refinement techniques to XFEL diffraction data to produce accurate diffraction data sets from a limited number of diffraction images (Uervirojnangkoorn et al., 2015), and went on to use these techniques to obtain a complete data set from crystals of the synaptotagmin-1 / SNARE complex and to determine the structure at 3.5 Å resolution (Zhou et al., 2015). Here, we describe new advances in our methods and present a reprocessed XFEL data set of the synaptotagmin-1 / SNARE complex. The reprocessing produced small improvements in electron density maps and the refined atomic model. The maps also contained more information than those of a lower resolution (4.1 Å) synchrotron data set. Processing a set of simulated XFEL diffraction images revealed that our methods yield accurate data and atomic models.

DOI:http://dx.doi.org/10.7554/eLife.18740.001

Mapping the conformational landscape of a dynamic enzyme by multitemperature and XFEL crystallography

Determining the interconverting conformations of dynamic proteins in atomic detail is a major challenge for structural biology. Conformational heterogeneity in the active site of the dynamic enzyme cyclophilin A (CypA) has been previously linked to its catalytic function, but the extent to which the different conformations of these residues are correlated is unclear. Here we compare the conformational ensembles of CypA by multitemperature synchrotron crystallography and fixed-target X-ray free-electron laser (XFEL) crystallography. The diffraction-before-destruction nature of XFEL experiments provides a radiation-damage-free view of the functionally important alternative conformations of CypA, confirming earlier synchrotron-based results. We monitored the temperature dependences of these alternative conformations with eight synchrotron datasets spanning 100-310 K. Multiconformer models show that many alternative conformations in CypA are populated only at 240 K and above, yet others remain populated or become populated at 180 K and below. These results point to a complex evolution of conformational heterogeneity between 180-–240 K that involves both thermal deactivation and solvent-driven arrest of protein motions in the crystal. The lack of a single shared conformational response to temperature within the dynamic active-site network provides evidence for a conformation shuffling model, in which exchange between rotamer states of a large aromatic ring in the middle of the network shifts the conformational ensemble for the other residues in the network. Together, our multitemperature analyses and XFEL data motivate a new generation of temperature- and time-resolved experiments to structurally characterize the dynamic underpinnings of protein function.

DOI:http://dx.doi.org/10.7554/eLife.07574.001

Proteins are the workhorses of the cell. The shape that a protein molecule adopts enables it to carry out its role. However, a protein’s shape, or 'conformation', is not static. Instead, a protein can shift between different conformations. This is particularly true for enzymes – the proteins that catalyze chemical reactions. The region of an enzyme where the chemical reaction happens, known as the active site, often has to change its conformation to allow catalysis to proceed. Changes in temperature can also make a protein shift between alternative conformations. Understanding how a protein shifts between conformations gives insight into how it works.

A common method for studying protein conformation is X-ray crystallography. This technique uses a beam of X-rays to figure out where the atoms of the protein are inside a crystal made of millions of copies of that protein. At room temperature or biological temperature, X-rays can rapidly damage the protein. Because of this, most crystal structures are determined at very low temperatures to minimize damage. But cooling to low temperatures changes the conformations that the protein adopts, and usually causes fewer conformations to be present.

Keedy, Kenner, Warkentin, Woldeyes et al. have used X-ray crystallography from a very low temperature (-173°C or 100 K) to above room temperature (up to 27°C or 300 K) to explore the alternative conformations of an enzyme called cyclophilin A. These alternative conformations include those that have previously been linked to this enzyme’s activity. Starting at a low temperature, parts of the enzyme were seen to shift from having a single conformation to many conformations above a threshold temperature. Unexpectedly, different parts of the enzyme have different threshold temperatures, suggesting that there isn’t a single transition across the whole protein. Instead, it appears the way a protein’s conformation changes in response to temperature is more complex than was previously realized. This result suggests that conformations in different parts of a protein are coupled to each other in complex ways.

Keedy, Kenner, Warkentin, Woldeyes et al. then performed X-ray crystallography at room temperature using an X-ray free-electron laser (XFEL). This technique can capture the protein’s structure before radiation damage occurs, and confirmed that the alternative conformations observed were not affected by radiation damage.

The combination of X-ray crystallography at multiple temperatures, new analysis methods for identifying and measuring alternative conformations, and XFEL crystallography should help future studies to characterize conformational changes in other proteins.

DOI:http://dx.doi.org/10.7554/eLife.07574.002

Macromolecular X-ray structure determination using weak, single-wavelength anomalous data

We describe a likelihood-based method for determining the substructure of anomalously scattering atoms in macromolecular crystals that allows successful structure determination by single-wavelength anomalous diffraction (SAD) X-ray analysis with weak anomalous signal. With the use of partial models and electron density maps in searches for anomalously scattering atoms, testing of alternative values of parameters and parallelized automated model-building, this method has the potential to extend the applicability of the SAD method in challenging cases.

Goniometer-based femtosecond crystallography with X-ray free electron lasers

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.

The R-factor gap in macromolecular crystallography: an untapped potential for insights on accurate structures

In macromolecular crystallography, the agreement between observed and predicted structure factors (R_cryst and R_free) is seldom better than 20%. This is much larger than the estimate of experimental error (R_merge). The difference between R_cryst and R_merge is the R-factor gap. There is no such gap in small-molecule crystallography, for which calculated structure factors are generally considered more accurate than the experimental measurements. Perhaps the true noise level of macromolecular data is higher than expected? Or is the gap caused by inaccurate phases that trap refined models in local minima? By generating simulated diffraction patterns using the program MLFSOM, and including every conceivable source of experimental error, we show that neither is the case. Processing our simulated data yielded values that were indistinguishable from those of real data for all crystallographic statistics except the final R_cryst and R_free. These values decreased to 3.8% and 5.5% for simulated data, suggesting that the reason for high R-factors in macromolecular crystallography is neither experimental error nor phase bias, but rather an underlying inadequacy in the models used to explain our observations. The present inability to accurately represent the entire macromolecule with both its flexibility and its protein-solvent interface may be improved by synergies between small-angle X-ray scattering, computational chemistry and crystallography. The exciting implication of our finding is that macromolecular data contain substantial hidden and untapped potential to resolve ambiguities in the true nature of the nanoscale, a task that the second century of crystallography promises to fulfill.

Crystallographic data processing for free-electron laser sources

A processing pipeline for diffraction data acquired using the `serial crystallography' methodology with a free-electron laser source is described with reference to the crystallographic analysis suite CrystFEL and the pre-processing program Cheetah. A detailed analysis of the nature and impact of indexing ambiguities is presented. Simulations of the Monte Carlo integration scheme, which accounts for the partially recorded nature of the diffraction intensities, are presented and show that the integration of partial reflections could be made to converge more quickly if the bandwidth of the X-rays were to be increased by a small amount or if a slight convergence angle were introduced into the incident beam.

Implementation and performance of SIBYLS: a dual endstation small-angle X-ray scattering and macromolecular crystallography beamline at the Advanced Light Source

The SIBYLS beamline (12.3.1) of the Advanced Light Source at Lawrence Berkeley National Laboratory, supported by the US Department of Energy and the National Institutes of Health, is optimized for both small-angle X-ray scattering (SAXS) and macromolecular crystallography (MX), making it unique among the world's mostly SAXS or MX dedicated beamlines. Since SIBYLS was commissioned, assessments of the limitations and advantages of a combined SAXS and MX beamline have suggested new strategies for integration and optimal data collection methods and have led to additional hardware and software enhancements. Features described include a dual mode monochromator [containing both Si(111) crystals and Mo/B(4)C multilayer elements], rapid beamline optics conversion between SAXS and MX modes, active beam stabilization, sample-loading robotics, and mail-in and remote data collection. These features allow users to gain valuable insights from both dynamic solution scattering and high-resolution atomic diffraction experiments performed at a single synchrotron beamline. Key practical issues considered for data collection and analysis include radiation damage, structural ensembles, alternative conformers and flexibility. SIBYLS develops and applies efficient combined MX and SAXS methods that deliver high-impact results by providing robust cost-effective routes to connect structures to biology and by performing experiments that aid beamline designs for next generation light sources.

The point-spread function of fiber-coupled area detectors

The point-spread function (PSF) of a fiber-optic taper-coupled CCD area detector was measured over five decades of intensity using a 20 µm X-ray beam and ~2000-fold averaging. The `tails' of the PSF clearly revealed that it is neither Gaussian nor Lorentzian, but instead resembles the solid angle subtended by a pixel at a point source of light held a small distance (~27 µm) above the pixel plane. This converges to an inverse cube law far from the beam impact point. Further analysis revealed that the tails are dominated by the fiber-optic taper, with negligible contribution from the phosphor, suggesting that the PSF of all fiber-coupled CCD-type detectors is best described as a Moffat function.

Time-resolved protein nanocrystallography using an X-ray free-electron laser

We demonstrate the use of an X-ray free electron laser synchronized with an optical pump laser to obtain X-ray diffraction snapshots from the photoactivated states of large membrane protein complexes in the form of nanocrystals flowing in a liquid jet. Light-induced changes of Photosystem I-Ferredoxin co-crystals were observed at time delays of 5 to 10 µs after excitation. The result correlates with the microsecond kinetics of electron transfer from Photosystem I to ferredoxin. The undocking process that follows the electron transfer leads to large rearrangements in the crystals that will terminally lead to the disintegration of the crystals. We describe the experimental setup and obtain the first time-resolved femtosecond serial X-ray crystallography results from an irreversible photo-chemical reaction at the Linac Coherent Light Source. This technique opens the door to time-resolved structural studies of reaction dynamics in biological systems.

Radiation damage in protein serial femtosecond crystallography using an x-ray free-electron laser

X-ray free-electron lasers deliver intense femtosecond pulses that promise to yield high resolution diffraction data of nanocrystals before the destruction of the sample by radiation damage. Diffraction intensities of lysozyme nanocrystals collected at the Linac Coherent Light Source using 2 keV photons were used for structure determination by molecular replacement and analyzed for radiation damage as a function of pulse length and fluence. Signatures of radiation damage are observed for pulses as short as 70 fs. Parametric scaling used in conventional crystallography does not account for the observed effects.

Structure-factor analysis of femtosecond microdiffraction patterns from protein nanocrystals

A complete set of structure factors has been extracted from hundreds of thousands of femtosecond single-shot X-ray microdiffraction patterns taken from randomly oriented nanocrystals. The method of Monte Carlo integration over crystallite size and orientation was applied to experimental data from Photosystem I nanocrystals. This arrives at structure factors from many partial reflections without prior knowledge of the particle-size distribution. The data were collected at the Linac Coherent Light Source (the first hard-X-ray laser user facility), to which was fitted a hydrated protein nanocrystal injector jet, according to the method of serial crystallography. The data are single 'still' diffraction snapshots, each from a different nanocrystal with sizes ranging between 100 nm and 2 µm, so the angular width of Bragg peaks was dominated by crystal-size effects. These results were compared with single-crystal data recorded from large crystals of Photosystem I at the Advanced Light Source and the quality of the data was found to be similar. The implications for improving the efficiency of data collection by allowing the use of very small crystals, for radiation-damage reduction and for time-resolved diffraction studies at room temperature are discussed.

Automated electron-density sampling reveals widespread conformational polymorphism in proteins

Although proteins populate large structural ensembles, X-ray diffraction data are traditionally interpreted using a single model. To search for evidence of alternate conformers, we developed a program, Ringer, which systematically samples electron density around the dihedral angles of protein side chains. In a diverse set of 402 structures, Ringer identified weak, nonrandom electron-density features that suggest of the presence of hidden, lowly populated conformations for >18% of uniquely modeled residues. Although these peaks occur at electron-density levels traditionally regarded as noise, statistically significant (P < 10(-5)) enrichment of peaks at successive rotameric chi angles validates the assignment of these features as unmodeled conformations. Weak electron density corresponding to alternate rotamers also was detected in an accurate electron density map free of model bias. Ringer analysis of the high-resolution structures of free and peptide-bound calmodulin identified shifts in ensembles and connected the alternate conformations to ligand recognition. These results show that the signal in high-resolution electron density maps extends below the traditional 1 sigma cutoff, and crystalline proteins are more polymorphic than current crystallographic models. Ringer provides an objective, systematic method to identify previously undiscovered alternate conformations that can mediate protein folding and function.

The minimum crystal size needed for a complete diffraction data set

In this work, classic intensity formulae were united with an empirical spot-fading model in order to calculate the diameter of a spherical crystal that will scatter the required number of photons per spot at a desired resolution over the radiation-damage-limited lifetime. The influences of molecular weight, solvent content, Wilson B factor, X-ray wavelength and attenuation on scattering power and dose were all included. Taking the net photon count in a spot as the only source of noise, a complete data set with a signal-to-noise ratio of 2 at 2 A resolution was predicted to be attainable from a perfect lysozyme crystal sphere 1.2 microm in diameter and two different models of photoelectron escape reduced this to 0.5 or 0.34 microm. These represent 15-fold to 700-fold less scattering power than the smallest experimentally determined crystal size to date, but the gap was shown to be consistent with the background scattering level of the relevant experiment. These results suggest that reduction of background photons and diffraction spot size on the detector are the principal paths to improving crystallographic data quality beyond current limits.

A beginner's guide to radiation damage

Many advances in the understanding of radiation damage to protein crystals, particularly at cryogenic temperatures, have been made in recent years, but with this comes an expanding literature, and, to the new breed of protein crystallographer who is not really interested in X-ray physics or radiation chemistry but just wants to solve a biologically relevant structure, the technical nature and breadth of this literature can be daunting. The purpose of this paper is to serve as a rough guide to radiation damage issues, and to provide references to the more exacting and detailed work. No attempt has been made to report precise numbers (a factor of two is considered satisfactory), and, since there are aspects of radiation damage that are demonstrably unpredictable, the 'worst case scenario' as well as the 'average crystal' are discussed in terms of the practicalities of data collection.

Determination of X-ray flux using silicon pin diodes

Accurate measurement of photon flux from an X-ray source, a parameter required to calculate the dose absorbed by the sample, is not yet routinely available at macromolecular crystallography beamlines. The development of a model for determining the photon flux incident on pin diodes is described here, and has been tested on the macromolecular crystallography beamlines at both the Swiss Light Source, Villigen, Switzerland, and the Advanced Light Source, Berkeley, USA, at energies between 4 and 18 keV. These experiments have shown that a simple model based on energy deposition in silicon is sufficient for determining the flux incident on high-quality silicon pin diodes. The derivation and validation of this model is presented, and a web-based tool for the use of the macromolecular crystallography and wider synchrotron community is introduced.

Molecular functions of the histone acetyltransferase chaperone complex Rtt109-Vps75

Histone acetylation and nucleosome remodeling regulate DNA damage repair, replication and transcription. Rtt109, a recently discovered histone acetyltransferase (HAT) from Saccharomyces cerevisiae, functions with the histone chaperone Asf1 to acetylate lysine K56 on histone H3 (H3K56), a modification associated with newly synthesized histones. In vitro analysis of Rtt109 revealed that Vps75, a Nap1 family histone chaperone, could also stimulate Rtt109-dependent acetylation of H3K56. However, the molecular function of the Rtt109-Vps75 complex remains elusive. Here we have probed the molecular functions of Vps75 and the Rtt109-Vps75 complex through biochemical, structural and genetic means. We find that Vps75 stimulates the k_cat of histone acetylation by ∼100-fold relative to Rtt109 alone and enhances acetylation of K9 in the H3 histone tail. Consistent with the In vitro evidence, cells lacking Vps75 showed a substantial reduction (60%) in H3K9 acetylation during S phase. X-ray structural, biochemical and genetic analyses of Vps75 indicate a unique, structurally dynamic Nap1-like fold that suggests a potential mechanism of Vps75-dependent activation of Rtt109. Together, these data provide evidence for a multifunctional HAT-chaperone complex that acetylates histone H3 and deposits H3-H4 onto DNA, linking histone modification and nucleosome assembly.

Structure of a designed, right-handed coiled-coil tetramer containing all biological amino acids

The previous design of an unprecedented family of two-, three-, and four-helical, right-handed coiled coils utilized nonbiological amino acids to efficiently pack spaces in the oligomer cores. Here we show that a stable, right-handed parallel tetrameric coiled coil, called RH4B, can be designed entirely using biological amino acids. The X-ray crystal structure of RH4B was determined to 1.1 Angstrom resolution using a designed metal binding site to coordinate a single Yb(2+) ion per 33-amino acid polypeptide chain. The resulting experimental phases were particularly accurate, and the experimental electron density map provided an especially clear, unbiased view of the molecule. The RH4B structure closely matched the design, with equivalent core rotamers and an overall root-mean-square deviation for the N-terminal repeat of the tetramer of 0.24 Angstrom. The clarity and resolution of the electron density map, however, revealed alternate rotamers and structural differences between the three sequence repeats in the molecule. These results suggest that the RH4B structure populates an unanticipated variety of structures.

A steric block in translation caused by the antibiotic spectinomycin

The widely used antibiotic spectinomycin inhibits bacterial protein synthesis by blocking translocation of messenger RNA and transfer RNAs on the ribosome. Here, we show that in crystals of the Escherichia coli 70S ribosome spectinomycin binding traps a distinct swiveling state of the head domain of the small ribosomal subunit. Spectinomycin also alters the rate and completeness of reverse translocation in vitro. These structural and biochemical data indicate that in solution spectinomycin sterically blocks swiveling of the head domain of the small ribosomal subunit and thereby disrupts the translocation cycle.

Structural basis for aminoglycoside inhibition of bacterial ribosome recycling

Aminoglycosides are widely used antibiotics that cause messenger RNA decoding errors, block mRNA and transfer RNA translocation, and inhibit ribosome recycling. Ribosome recycling follows the termination of protein synthesis and is aided by ribosome recycling factor (RRF) in bacteria. The molecular mechanism by which aminoglycosides inhibit ribosome recycling is unknown. Here we show in X-ray crystal structures of the Escherichia coli 70S ribosome that RRF binding causes RNA helix H69 of the large ribosomal subunit, which is crucial for subunit association, to swing away from the subunit interface. Aminoglycosides bind to H69 and completely restore the contacts between ribosomal subunits that are disrupted by RRF. These results provide a structural explanation for aminoglycoside inhibition of ribosome recycling.

Structural insight into KCNQ (Kv7) channel assembly and channelopathy

Kv7.x (KCNQ) voltage-gated potassium channels form the cardiac and auditory I(Ks) current and the neuronal M-current. The five Kv7 subtypes have distinct assembly preferences encoded by a C-terminal cytoplasmic assembly domain, the A-domain Tail. Here, we present the high-resolution structure of the Kv7.4 A-domain Tail together with biochemical experiments that show that the domain is a self-assembling, parallel, four-stranded coiled coil. Structural analysis and biochemical studies indicate conservation of the coiled coil in all Kv7 subtypes and that a limited set of interactions encode assembly specificity determinants. Kv7 mutations have prominent roles in arrhythmias, deafness, and epilepsy. The structure together with biochemical data indicate that A-domain Tail arrhythmia mutations cluster on the solvent-accessible surface of the subunit interface at a likely site of action for modulatory proteins. Together, the data provide a framework for understanding Kv7 assembly specificity and the molecular basis of a distinct set of Kv7 channelopathies.

Basis for a ubiquitin-like protein thioester switch toggling E1-E2 affinity

Ubiquitin-like proteins (UBLs) are conjugated by dynamic E1-E2-E3 enzyme cascades. E1 enzymes activate UBLs by catalysing UBL carboxy-terminal adenylation, forming a covalent E1 throught UBL thioester intermediate, and generating a thioester-linked E2 throught UBL product, which must be released for subsequent reactions. Here we report the structural analysis of a trapped UBL activation complex for the human NEDD8 pathway, containing NEDD8's heterodimeric E1 (APPBP1-UBA3), two NEDD8s (one thioester-linked to E1, one noncovalently associated for adenylation), a catalytically inactive E2 (Ubc12), and MgATP. The results suggest that a thioester switch toggles E1-E2 affinities. Two E2 binding sites depend on NEDD8 being thioester-linked to E1. One is unmasked by a striking E1 conformational change. The other comes directly from the thioester-bound NEDD8. After NEDD8 transfer to E2, reversion to an alternate E1 conformation would facilitate release of the E2 throught NEDD8 thioester product. Thus, transferring the UBL's thioester linkage between successive conjugation enzymes can induce conformational changes and alter interaction networks to drive consecutive steps in UBL cascades.

XANES measurements of the rate of radiation damage to selenomethionine side chains

The radiation-induced disordering of selenomethionine (SeMet) side chains represents a significant impediment to protein structure solution. Not only does the increased B-factor of these sites result in a serious drop in phasing power, but some sites decay much faster than others in the same unit cell. These radio-labile SeMet side chains decay faster than high-order diffraction spots with dose, making it difficult to detect this kind of damage by inspection of the diffraction pattern. The selenium X-ray absorbance near-edge spectrum (XANES) from samples containing SeMet was found to change significantly after application of X-ray doses of 10-100 MGy. Most notably, the sharp ;white line' feature near the canonical Se edge disappears. The change was attributed to breakage of the Cgamma-Se bond in SeMet. This spectral change was used as a probe to measure the decay rate of SeMet with X-ray dose in cryo-cooled samples. Two protein crystal types and 15 solutions containing free SeMet amino acid were examined. The damage rate was influenced by the chemical and physical condition of the sample, and the half-decaying dose for the selenium XANES signal ranged from 5 to 43 MGy. These decay rates were 34- to 3.8-fold higher than the rate at which the Se atoms interacted directly with X-ray photons, so the damage mechanism must be a secondary effect. Samples that cooled to a more crystalline state generally decayed faster than samples that cooled to an amorphous solid. The single exception was a protein crystal where a nanocrystalline cryoprotectant had a protective effect. Lowering the pH, especially with ascorbic or nitric acids, had a protective effect, and SeMet lifetime increased monotonically with decreasing sample temperature (down to 93 K). The SeMet lifetime in one protein crystal was the same as that of the free amino acid, and the longest SeMet lifetime measured was found in the other protein crystal type. This protection was found to arise from the folded structure of the protein molecule. A mechanism to explain observed decay rates involving the damaging species following the electric field lines around protein molecules is proposed.

Functional and structural characterization of a thiol peroxidase from Mycobacterium tuberculosis

A thiol peroxidase (Tpx) from Mycobacterium tuberculosis was functionally analyzed. The enzyme shows NADPH-linked peroxidase activity using a thioredoxin-thioredoxin reductase system as electron donor, and anti-oxidant activity in a thiol-dependent metal-catalyzed oxidation system. It reduces H2O2, t-butyl hydroperoxide, and cumene hydroperoxide, and is inhibited by sulfhydryl reagents. Mutational studies revealed that the peroxidatic (Cys60) and resolving (Cys93) cysteine residues are critical amino acids for catalytic activity. The X-ray structure determined to a resolution of 1.75 A shows a thioredoxin fold similar to that of other peroxiredoxin family members. Superposition with structural homologues in oxidized and reduced forms indicates that the M. tuberculosis Tpx is a member of the atypical two-Cys peroxiredoxin family. In addition, the short distance that separates the Calpha atoms of Cys60 and Cys93 and the location of these cysteine residues in unstructured regions may indicate that the M. tuberculosis enzyme is oxidized, though the side-chain of Cys60 is poorly visible. It is solely in the reduced Streptococcus pneumoniae Tpx structure that both residues are part of two distinct helical segments. The M. tuberculosis Tpx is dimeric both in solution and in the crystal structure. Amino acid residues from both monomers delineate the active site pocket.

Structures of the bacterial ribosome at 3.5 A resolution

We describe two structures of the intact bacterial ribosome from Escherichia coli determined to a resolution of 3.5 angstroms by x-ray crystallography. These structures provide a detailed view of the interface between the small and large ribosomal subunits and the conformation of the peptidyl transferase center in the context of the intact ribosome. Differences between the two ribosomes reveal a high degree of flexibility between the head and the rest of the small subunit. Swiveling of the head of the small subunit observed in the present structures, coupled to the ratchet-like motion of the two subunits observed previously, suggests a mechanism for the final movements of messenger RNA (mRNA) and transfer RNAs (tRNAs) during translocation.

Structural and mechanistic studies of VPS4 proteins

VPS4 ATPases function in multivesicular body formation and in HIV-1 budding. Here, we report the crystal structure of monomeric apo human VPS4B/SKD1 (hVPS4B), which is composed of five distinct elements: a poorly ordered N-terminal MIT domain that binds ESCRT-III substrates, large (mixed alpha/beta) and small (alpha) AAA ATPase domains that closely resemble analogous domains in the p97 D1 ATPase cassette, a three-stranded antiparallel beta domain inserted within the small ATPase domain, and a novel C-terminal helix. Apo hVPS4B and yeast Vps4p (yVps4p) proteins dimerized in solution, and assembled into larger complexes (10-12 subunits) upon ATP binding. Human and yeast adaptor proteins (LIP5 and yVta1p, respectively) bound the beta domains of the fully assembled hVPS4B and yVps4p proteins. We therefore propose that Vps4 proteins cycle between soluble, inactive low molecular weight complexes and active, membrane-associated double-ring structures that bind ATP and coassemble with LIP5/Vta1. Finally, HIV-1 budding was inhibited by mutations in a loop that projects into the center of the modeled hVPS4B rings, suggesting that hVPS4B may release the assembled ESCRT machinery by pulling ESCRT-III substrates up into the central pore.