Search

Summary: The upcoming exascale computing systems Frontier and Aurora will draw much of their computing power from GPU accelerators. The hardware for these systems will be provided by AMD and Intel, respectively, each supporting their own GPU programming model. The challenge for applications that harness one of these exascale systems will be to avoid lock‐in and to preserve performance portability. We report here on our results of using Kokkos to accelerate a real‐world application on NERSC's Perlmutter Phase 1 (using NVIDIA A100 accelerators) and Crusher, the testbed system for OLCF's Frontier (using AMD MI250X). By porting to Kokkos, we successfully ran the same X‐ray tracing code on both systems and achieved speed‐ups between 13 % and 66 % compared to the original CUDA code. These results are a highly encouraging demonstration of using Kokkos to accelerate production science code. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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. Synthetic macromolecular crystallography diffraction‐image data were generated to demonstrate the challenges of combining data from multiple crystals with indexing ambiguity in the context of heavy radiation damage. The nature of the problems encountered using contemporary data‐processing programs is summarized. [ABSTRACT FROM AUTHOR]

In macromolecular crystallography, the agreement between observed and predicted structure factors (Rcryst and R free ) is seldom better than 20%. This is much larger than the estimate of experimental error (R merge ). The difference between Rcryst 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 Rcryst 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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/B4C 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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

The article presents a study of the structure factor of femtosecond microdiffraction forms form protein nanocrystals. A set of structure factors was extracted from multiple femtosecond single-shot X-ray microdiffraction patterns obtained from nanocrystals. Monte Carlo integration was employed to experimental data from Photosystem I nanocrystals at the Linac Coherent Light Source. Results showed that crystal-size effects dominated the angular width of Bragg peaks.

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. [ABSTRACT FROM AUTHOR]

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. [ABSTRACT FROM AUTHOR]

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 Å resolution using a designed metal binding site to coordinate a single Yb2+ 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 Å. 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. [ABSTRACT FROM AUTHOR]

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 Cγ—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. [ABSTRACT FROM AUTHOR]

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 α/β) and small (α) AAA ATPase domains that closely resemble analogous domains in the p97 D1 ATPase cassette, a three-stranded antiparallel β 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 β 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. [ABSTRACT FROM AUTHOR]

At the Advanced Light Source, three protein crystallography beamlines have been built that use as a source one of the three 6 T single-pole superconducting bending magnets (superbends) that were recently installed in the ring. The use of such single-pole super-conducting bend magnets enables the development of a hard X-ray program on a relatively low-energy 1.9 GeV ring without taking up insertion-device straight sections. The source is of relatively low power but, owing to the small electron beam emittance, it has high brightness. X-ray optics are required to preserve the brightness and to match the illumination requirements for protein crystallography. This was achieved by means of a collimating premirror bent to a plane parabola, a double-crystal monochromator followed by a toroidal mirror that focuses in the horizontal direction with a 2:1 demagnification. This optical arrangement partially balances aberrations from the collimating and toroidal mirrors such that a tight focused spot size is achieved. The optical properties of the beamline are an excellent match to those required by the small protein crystals that are typically measured. The design and performance of these new beamlines are described. [ABSTRACT FROM AUTHOR]

Class I major histocompatibility complex (MHC) molecules bind peptides derived from degraded proteins for display to T cells of the immune system. Peptides bind to MHC proteins with varying affinities, depending upon their sequence and length. We demonstrate that the thermal stability of the MHC-peptide complex depends directly on peptide binding affinity. We use this correlation to develop a convenient method to determine peptide dissociation constants by measuring MHC-peptide complex stability using thermal denaturation profiles monitored by circular dichroism. [ABSTRACT FROM AUTHOR]

A response by Robin L. Owen and colleagues to a letter to the editor about their article regarding the determination of x-ray flux with the use of silicon diodes in the 2009 issue is presented.

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.