Autoindexing with outlier rejection and identification of superimposed lattices

Rebamipide and Derivatives are Potent, Selective Inhibitors of Histidine Phosphatase Activity of the Suppressor of T Cell Receptor Signaling Proteins

The suppressor of T cell receptor signaling (Sts) proteins are negative regulators of immune signaling. Genetic inactivation of these proteins leads to significant resistance to infection. From a 590,000 compound high-throughput screen, we identified the 2-(1H)-quinolinone derivative, rebamipide, as a putative inhibitor of Sts phosphatase activity. Rebamipide, and a small library of derivatives, are competitive, selective inhibitors of Sts-1 with IC50 values from low to submicromolar. SAR analysis indicates that the quinolinone, the acid, and the amide moieties are all essential for activity. A crystal structure confirmed the SAR and reveals key interactions between this class of compound and the protein. Although rebamipide has poor cell permeability, we demonstrated that a liposomal preparation can inactivate the phosphatase activity of Sts-1 in cells. These studies demonstrate that Sts-1 enzyme activity can be pharmacologically inactivated and provide foundational tools and insights for the development of immune-enhancing therapies that target the Sts proteins.

A reference-based multi-lattice indexing method integrating prior information correction and iterative refinement in protein crystallography

A new multi-lattice indexing method based on the principle of whole-pattern matching given cell dimensions and space-group symmetry is presented for macromolecular crystallography. The proposed method, termed the multi-crystal data processing suite (MCDPS), features a local correction for prior information accompanied by iterative refinement of experimental parameters, both of which are numerically and experimentally demonstrated to be critical for accurately identifying multiple crystal lattices. Further analysis of data reduction and structure determination with conventional single-crystal programs reveals that the processed multi-lattice data sets are comparable in quality to typical single-crystal ones in terms of crystallographic metrics. Importantly, it is confirmed that careful exclusion of overlapping reflections prior to scaling is necessary to guarantee an accurate data reduction result. The potential for multi-lattice indexing in solving the general macroscopic twinning problem is also explored.

MXCuBE2: the dawn of MXCuBE Collaboration

MXCuBE2 is the second-generation evolution of the MXCuBE beamline control software, initially developed and used at ESRF - the European Synchrotron. MXCuBE2 extends, in an intuitive graphical user interface (GUI), the functionalities and data collection methods available to users while keeping all previously available features and allowing for the straightforward incorporation of ongoing and future developments. MXCuBE2 introduces an extended abstraction layer that allows easy interfacing of any kind of macromolecular crystallography (MX) hardware component, whether this is a diffractometer, sample changer, detector or optical element. MXCuBE2 also works in strong synergy with the ISPyB Laboratory Information Management System, accessing the list of samples available for a particular experimental session and associating, either from instructions contained in ISPyB or from user input via the MXCuBE2 GUI, different data collection types to them. The development of MXCuBE2 forms the core of a fruitful collaboration which brings together several European synchrotrons and a software development factory and, as such, defines a new paradigm for the development of beamline control platforms for the European MX user community.

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.

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.

Protein crystallization and initial neutron diffraction studies of the photosystem II subunit PsbO

The PsbO protein of photosystem II stabilizes the active-site manganese cluster and is thought to act as a proton antenna. To enable neutron diffraction studies, crystals of the β-barrel core of PsbO were grown in capillaries. The crystals were optimized by screening additives in a counter-diffusion setup in which the protein and reservoir solutions were separated by a 1% agarose plug. Crystals were cross-linked with glutaraldehyde. Initial neutron diffraction data were collected from a 0.25 mm3 crystal at room temperature using the MaNDi single-crystal diffractometer at the Spallation Neutron Source, Oak Ridge National Laboratory.

FELIX: an algorithm for indexing multiple crystallites in X-ray free-electron laser snapshot diffraction images

A novel algorithm for indexing multiple crystals in snapshot X-ray diffraction images, especially suited for serial crystallography data, is presented. The algorithm, FELIX, utilizes a generalized parametrization of the Rodrigues-Frank space, in which all crystal systems can be represented without singularities. The new algorithm is shown to be capable of indexing more than ten crystals per image in simulations of cubic, tetragonal and monoclinic crystal diffraction patterns. It is also used to index an experimental serial crystallography dataset from lysozyme microcrystals. The increased number of indexed crystals is shown to result in a better signal-to-noise ratio, and fewer images are needed to achieve the same data quality as when indexing one crystal per image. The relative orientations between the multiple crystals indexed in an image show a slight tendency of the lysozme microcrystals to adhere on ([Formula: see text]10) facets.

Diffraction-geometry refinement in the DIALS framework

A comprehensive description of the methods used within the DIALS framework for diffraction-geometry refinement using predicted reflection centroids is given. Examples of the advanced features of the software are provided.

Rapid data collection and modern computing resources provide the opportunity to revisit the task of optimizing the model of diffraction geometry prior to integration. A comprehensive description is given of new software that builds upon established methods by performing a single global refinement procedure, utilizing a smoothly varying model of the crystal lattice where appropriate. This global refinement technique extends to multiple data sets, providing useful constraints to handle the problem of correlated parameters, particularly for small wedges of data. Examples of advanced uses of the software are given and the design is explained in detail, with particular emphasis on the flexibility and extensibility it entails.

Acoustic Injectors for Drop-On-Demand Serial Femtosecond Crystallography

X-ray free-electron lasers (XFELs) provide very intense X-ray pulses suitable for macromolecular crystallography. Each X-ray pulse typically lasts for tens of femtoseconds and the interval between pulses is many orders of magnitude longer. Here we describe two novel acoustic injection systems that use focused sound waves to eject picoliter to nanoliter crystal-containing droplets out of microplates and into the X-ray pulse from which diffraction data are collected. The on-demand droplet delivery is synchronized to the XFEL pulse scheme, resulting in X-ray pulses intersecting up to 88% of the droplets. We tested several types of samples in a range of crystallization conditions, wherein the overall crystal hit ratio (e.g., fraction of images with observable diffraction patterns) is a function of the microcrystal slurry concentration. We report crystal structures from lysozyme, thermolysin, and stachydrine demethylase (Stc2). Additional samples were screened to demonstrate that these methods can be applied to rare samples.

Assessment of microcrystal quality by transmission electron microscopy for efficient serial femtosecond crystallography

Serial femtosecond crystallography (SFX) employing high-intensity X-ray free-electron laser (XFEL) sources has enabled structural studies on microcrystalline protein samples at non-cryogenic temperatures. However, the identification and optimization of conditions that produce well diffracting microcrystals remains an experimental challenge. Here, we report parallel SFX and transmission electron microscopy (TEM) experiments using fragmented microcrystals of wild type (WT) homoprotocatechuate 2,3-dioxygenase (HPCD) and an active site variant (H200Q). Despite identical crystallization conditions and morphology, as well as similar crystal size and density, the indexing efficiency of the diffraction data collected using the H200Q variant sample was over 7-fold higher compared to the diffraction results obtained using the WT sample. TEM analysis revealed an abundance of protein aggregates, crystal conglomerates and a smaller population of highly ordered lattices in the WT sample as compared to the H200Q variant sample. While not reported herein, the 1.75 Å resolution structure of the H200Q variant was determined from ∼16 min of beam time, demonstrating the utility of TEM analysis in evaluating sample monodispersity and lattice quality, parameters critical to the efficiency of SFX experiments.

New methods for indexing multi-lattice diffraction data

A new indexing method is presented which is capable of indexing multiple crystal lattices from narrow wedges of diffraction data. The method takes advantage of a simplification of Fourier transform-based methods that is applicable when the unit-cell dimensions are known a priori. The efficacy of this method is demonstrated with both semi-synthetic multi-lattice data and real multi-lattice data recorded from crystals of ∼1 µm in size, where it is shown that up to six lattices can be successfully indexed and subsequently integrated from a 1° wedge of data. Analysis is presented which shows that improvements in data-quality indicators can be obtained through accurate identification and rejection of overlapping reflections prior to scaling.

Exploiting fast detectors to enter a new dimension in room-temperature crystallography

A departure from a linear or an exponential intensity decay in the diffracting power of protein crystals as a function of absorbed dose is reported. The observation of a lag phase raises the possibility of collecting significantly more data from crystals held at room temperature before an intolerable intensity decay is reached. A simple model accounting for the form of the intensity decay is reintroduced and is applied for the first time to high frame-rate room-temperature data collection.

Accurate macromolecular structures using minimal measurements from X-ray free-electron lasers

X-ray free-electron laser (XFEL) sources enable the use of crystallography to solve three-dimensional macromolecular structures under native conditions and without radiation damage. Results to date, however, have been limited by the challenge of deriving accurate Bragg intensities from a heterogeneous population of microcrystals, while at the same time modeling the X-ray spectrum and detector geometry. Here we present a computational approach designed to extract meaningful high-resolution signals from fewer diffraction measurements.

New Python-based methods for data processing

Current pixel-array detectors produce diffraction images at extreme data rates (of up to 2 TB h(-1)) that make severe demands on computational resources. New multiprocessing frameworks are required to achieve rapid data analysis, as it is important to be able to inspect the data quickly in order to guide the experiment in real time. By utilizing readily available web-serving tools that interact with the Python scripting language, it was possible to implement a high-throughput Bragg-spot analyzer (cctbx.spotfinder) that is presently in use at numerous synchrotron-radiation beamlines. Similarly, Python interoperability enabled the production of a new data-reduction package (cctbx.xfel) for serial femtosecond crystallography experiments at the Linac Coherent Light Source (LCLS). Future data-reduction efforts will need to focus on specialized problems such as the treatment of diffraction spots on interleaved lattices arising from multi-crystal specimens. In these challenging cases, accurate modeling of close-lying Bragg spots could benefit from the high-performance computing capabilities of graphics-processing units.

Autoindexing diffraction images with iMosflm

The principles of one-dimensional FFT-based autoindexing of diffraction images are described together with practical issues that may arise. A procedure for indexing multiple lattices as implemented in iMosflm is presented.

An overview of autoindexing diffraction images based on one-­dimensional fast Fourier transforms is presented. The implementation of the algorithm in the Mosflm/iMosflm program suite is described with a discussion of practical issues that may arise and ways of assessing the success or failure of the procedure. Recent developments allow indexing of images that show multiple lattices, and several examples demonstrate the success of this approach in real cases.

Outrunning free radicals in room-temperature macromolecular crystallography

A significant increase in the lifetime of room-temperature macromolecular crystals is reported through the use of a high-brilliance X-ray beam, reduced exposure times and a fast-readout detector. This is attributed to the ability to collect diffraction data before hydroxyl radicals can propagate through the crystal, fatally disrupting the lattice. Hydroxyl radicals are shown to be trapped in amorphous solutions at 100 K. The trend in crystal lifetime was observed in crystals of a soluble protein (immunoglobulin γ Fc receptor IIIa), a virus (bovine enterovirus serotype 2) and a membrane protein (human A(2A) adenosine G-protein coupled receptor). The observation of a similar effect in all three systems provides clear evidence for a common optimal strategy for room-temperature data collection and will inform the design of future synchrotron beamlines and detectors for macromolecular crystallography.

X-ray-excited optical luminescence of protein crystals: a new tool for studying radiation damage during diffraction data collection

During X-ray irradiation protein crystals radiate energy in the form of small amounts of visible light. This is known as X-ray-excited optical luminescence (XEOL). The XEOL of several proteins and their constituent amino acids has been characterized using the microspectrophotometers at the Swiss Light Source and Diamond Light Source. XEOL arises primarily from aromatic amino acids, but the effects of local environment and quenching within a crystal mean that the XEOL spectrum of a crystal is not the simple sum of the spectra of its constituent parts. Upon repeated exposure to X-rays XEOL spectra decay non-uniformly, suggesting that XEOL is sensitive to site-specific radiation damage. However, rates of XEOL decay were found not to correlate to decays in diffracting power, making XEOL of limited use as a metric for radiation damage to protein crystals.

Crystal structure of perakine reductase, founding member of a novel aldo-keto reductase (AKR) subfamily that undergoes unique conformational changes during NADPH binding

Perakine reductase (PR) catalyzes the NADPH-dependent reduction of the aldehyde perakine to yield the alcohol raucaffrinoline in the biosynthetic pathway of ajmaline in Rauvolfia, a key step in indole alkaloid biosynthesis. Sequence alignment shows that PR is the founder of the new AKR13D subfamily and is designated AKR13D1. The x-ray structure of methylated His(6)-PR was solved to 2.31 Å. However, the active site of PR was blocked by the connected parts of the neighbor symmetric molecule in the crystal. To break the interactions and obtain the enzyme-ligand complexes, the A213W mutant was generated. The atomic structure of His(6)-PR-A213W complex with NADPH was determined at 1.77 Å. Overall, PR folds in an unusual α(8)/β(6) barrel that has not been observed in any other AKR protein to date. NADPH binds in an extended pocket, but the nicotinamide riboside moiety is disordered. Upon NADPH binding, dramatic conformational changes and movements were observed: two additional β-strands in the C terminus become ordered to form one α-helix, and a movement of up to 24 Å occurs. This conformational change creates a large space that allows the binding of substrates of variable size for PR and enhances the enzyme activity; as a result cooperative kinetics are observed as NADPH is varied. As the founding member of the new AKR13D subfamily, PR also provides a structural template and model of cofactor binding for the AKR13 family.

Ferritin structure from Mycobacterium tuberculosis: comparative study with homologues identifies extended C-terminus involved in ferroxidase activity

Ferritins are recognized as key players in the iron storage and detoxification processes. Iron acquisition in the case of pathogenic bacteria has long been established as an important virulence mechanism. Here, we report a 3.0 Å crystal structure of a ferritin, annotated as Bacterioferritin B (BfrB), from Mycobacterium tuberculosis (Mtb), the causative agent of tuberculosis that continues to be one of the world's deadliest diseases. Similar to the other members of ferritin family, the Mtb BfrB subunit exhibits the characteristic fold of a four-helical bundle that possesses the ferroxidase catalytic centre. We compare the structure of Mtb BfrB with representatives of the ferritin family belonging to the archaea, eubacteria and eukarya. Unlike most other ferritins, Mtb BfrB has an extended C-terminus. To dissect the role of this extended C-terminus, truncated Mtb BfrB was purified and biochemical studies implicate this region in ferroxidase activity and iron release in addition to providing stability to the protein. Functionally important regions in a protein of known 3D-structure can be determined by estimating the degree of conservation of the amino-acid sites with its close homologues. Based on the comparative studies, we identify the slowly evolving conserved sites as well as the rapidly evolving variable sites and analyze their role in relation to structure and function of Mtb BfrB. Further, electrostatic computations demonstrate that although the electrostatic environment of catalytic residues is preserved within the family, extensive variability is exhibited by residues defining the channels and pores, in all likelihood keeping up with the diverse functions executed by these ferritins in varied environments.

Data processing and analysis with the autoPROC toolbox

A typical diffraction experiment will generate many images and data sets from different crystals in a very short time. This creates a challenge for the high-throughput operation of modern synchrotron beamlines as well as for the subsequent data processing. Novice users in particular may feel overwhelmed by the tables, plots and numbers that the different data-processing programs and software packages present to them. Here, some of the more common problems that a user has to deal with when processing a set of images that will finally make up a processed data set are shown, concentrating on difficulties that may often show up during the first steps along the path of turning the experiment (i.e. data collection) into a model (i.e. interpreted electron density). Difficulties such as unexpected crystal forms, issues in crystal handling and suboptimal choices of data-collection strategies can often be dealt with, or at least diagnosed, by analysing specific data characteristics during processing. In the end, one wants to distinguish problems over which one has no immediate control once the experiment is finished from problems that can be remedied a posteriori. A new software package, autoPROC, is also presented that combines third-party processing programs with new tools and an automated workflow script that is intended to provide users with both guidance and insight into the offline processing of data affected by the difficulties mentioned above, with particular emphasis on the automated treatment of multi-sweep data sets collected on multi-axis goniostats.