EMRinger: side chain-directed model and map validation for 3D cryo-electron microscopy.

Advances in high-resolution cryo-electron microscopy (cryo-EM) require the development of validation metrics to independently assess map quality and model geometry. We report EMRinger, a tool that assesses the precise fitting of an atomic model into the map during refinement and shows how radiation damage alters scattering from negatively charged amino acids. EMRinger (https://github.com/fraser-lab/EMRinger) will be useful for monitoring progress in resolving and modeling high-resolution features in cryo-EM.

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.

Structural and Biophysical Characterization of the Mycobacterium tuberculosis Protein Rv0577, a Protein Associated with Neutral Red Staining of Virulent Tuberculosis Strains and Homologue of the Streptomyces coelicolor Protein KbpA.

Mycobacterium tuberculosis protein Rv0577 is a prominent antigen in tuberculosis patients, the component responsible for neutral red staining of virulent strains of M. tuberculosis, a putative component in a methylglyoxal detoxification pathway, and an agonist of toll-like receptor 2. It also has an amino acid sequence that is 36% identical to that of Streptomyces coelicolor AfsK-binding protein A (KbpA), a component in the complex secondary metabolite pathways in the Streptomyces genus. To gain insight into the biological function of Rv0577 and the family of KpbA kinase regulators, the crystal structure for Rv0577 was determined to a resolution of 1.75 Å, binding properties with neutral red and deoxyadenosine were surveyed, backbone dynamics were measured, and thermal stability was assayed by circular dichroism spectroscopy. The protein is composed of four approximate repeats with a ßαßßß topology arranged radially in consecutive pairs to form two continuous eight-strand ß-sheets capped on both ends with an α-helix. The two ß-sheets intersect in the center at roughly a right angle and form two asymmetric deep "saddles" that may serve to bind ligands. Nuclear magnetic resonance chemical shift perturbation experiments show that neutral red and deoxyadenosine bind to Rv0577. Binding to deoxyadenosine is weaker with an estimated dissociation constants of 4.1 ± 0.3 mM for saddle 1. Heteronuclear steady-state {1H}-15N nuclear Overhauser effect, T1, and T2 values were generally uniform throughout the sequence with only a few modest pockets of differences. Circular dichroism spectroscopy characterization of the thermal stability of Rv0577 indicated irreversible unfolding upon heating with an estimated melting temperature of 56 °C.

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.

Improved low-resolution crystallographic refinement with Phenix and Rosetta.

Refinement of macromolecular structures against low-resolution crystallographic data is limited by the ability of current methods to converge on a structure with realistic geometry. We developed a low-resolution crystallographic refinement method that combines the Rosetta sampling methodology and energy function with reciprocal-space X-ray refinement in Phenix. On a set of difficult low-resolution cases, the method yielded improved model geometry and lower free R factors than alternate refinement methods.

Using support vector machines to improve elemental ion identification in macromolecular crystal structures.

In the process of macromolecular model building, crystallographers must examine electron density for isolated atoms and differentiate sites containing structured solvent molecules from those containing elemental ions. This task requires specific knowledge of metal-binding chemistry and scattering properties and is prone to error. A method has previously been described to identify ions based on manually chosen criteria for a number of elements. Here, the use of support vector machines (SVMs) to automatically classify isolated atoms as either solvent or one of various ions is described. Two data sets of protein crystal structures, one containing manually curated structures deposited with anomalous diffraction data and another with automatically filtered, high-resolution structures, were constructed. On the manually curated data set, an SVM classifier was able to distinguish calcium from manganese, zinc, iron and nickel, as well as all five of these ions from water molecules, with a high degree of accuracy. Additionally, SVMs trained on the automatically curated set of high-resolution structures were able to successfully classify most common elemental ions in an independent validation test set. This method is readily extensible to other elemental ions and can also be used in conjunction with previous methods based on a priori expectations of the chemical environment and X-ray scattering.

Indexing amyloid peptide diffraction from serial femtosecond crystallography: new algorithms for sparse patterns.

Still diffraction patterns from peptide nanocrystals with small unit cells are challenging to index using conventional methods owing to the limited number of spots and the lack of crystal orientation information for individual images. New indexing algorithms have been developed as part of the Computational Crystallography Toolbox (cctbx) to overcome these challenges. Accurate unit-cell information derived from an aggregate data set from thousands of diffraction patterns can be used to determine a crystal orientation matrix for individual images with as few as five reflections. These algorithms are potentially applicable not only to amyloid peptides but also to any set of diffraction patterns with sparse properties, such as low-resolution virus structures or high-throughput screening of still images captured by raster-scanning at synchrotron sources. As a proof of concept for this technique, successful integration of X-ray free-electron laser (XFEL) data to 2.5 Šresolution for the amyloid segment GNNQQNY from the Sup35 yeast prion is presented.

Energy-dispersive X-ray emission spectroscopy using an X-ray free-electron laser in a shot-by-shot mode.

The ultrabright femtosecond X-ray pulses provided by X-ray free-electron lasers open capabilities for studying the structure and dynamics of a wide variety of systems beyond what is possible with synchrotron sources. Recently, this "probe-before-destroy" approach has been demonstrated for atomic structure determination by serial X-ray diffraction of microcrystals. There has been the question whether a similar approach can be extended to probe the local electronic structure by X-ray spectroscopy. To address this, we have carried out femtosecond X-ray emission spectroscopy (XES) at the Linac Coherent Light Source using redox-active Mn complexes. XES probes the charge and spin states as well as the ligand environment, critical for understanding the functional role of redox-active metal sites. Kß(1,3) XES spectra of Mn(II) and Mn(2)(III,IV) complexes at room temperature were collected using a wavelength dispersive spectrometer and femtosecond X-ray pulses with an individual dose of up to >100 MGy. The spectra were found in agreement with undamaged spectra collected at low dose using synchrotron radiation. Our results demonstrate that the intact electronic structure of redox active transition metal compounds in different oxidation states can be characterized with this shot-by-shot method. This opens the door for studying the chemical dynamics of metal catalytic sites by following reactions under functional conditions. The technique can be combined with X-ray diffraction to simultaneously obtain the geometric structure of the overall protein and the local chemistry of active metal sites and is expected to prove valuable for understanding the mechanism of important metalloproteins, such as photosystem II.

Room temperature femtosecond X-ray diffraction of photosystem II microcrystals.

Most of the dioxygen on earth is generated by the oxidation of water by photosystem II (PS II) using light from the sun. This light-driven, four-photon reaction is catalyzed by the Mn(4)CaO(5) cluster located at the lumenal side of PS II. Various X-ray studies have been carried out at cryogenic temperatures to understand the intermediate steps involved in the water oxidation mechanism. However, the necessity for collecting data at room temperature, especially for studying the transient steps during the O-O bond formation, requires the development of new methodologies. In this paper we report room temperature X-ray diffraction data of PS II microcrystals obtained using ultrashort (< 50 fs) 9 keV X-ray pulses from a hard X-ray free electron laser, namely the Linac Coherent Light Source. The results presented here demonstrate that the "probe before destroy" approach using an X-ray free electron laser works even for the highly-sensitive Mn(4)CaO(5) cluster in PS II at room temperature. We show that these data are comparable to those obtained in synchrotron radiation studies as seen by the similarities in the overall structure of the helices, the protein subunits and the location of the various cofactors. This work is, therefore, an important step toward future studies for resolving the structure of the Mn(4)CaO(5) cluster without any damage at room temperature, and of the reaction intermediates of PS II during O-O bond formation.

Improved crystal orientation and physical properties from single-shot XFEL stills.

X-ray diffraction patterns from still crystals are inherently difficult to process because the crystal orientation is not uniquely determined by measuring the Bragg spot positions. Only one of the three rotational degrees of freedom is directly coupled to spot positions; the other two rotations move Bragg spots in and out of the reflecting condition but do not change the direction of the diffracted rays. This hinders the ability to recover accurate structure factors from experiments that are dependent on single-shot exposures, such as femtosecond diffract-and-destroy protocols at X-ray free-electron lasers (XFELs). Here, additional methods are introduced to optimally model the diffraction. The best orientation is obtained by requiring, for the brightest observed spots, that each reciprocal-lattice point be placed into the exact reflecting condition implied by Bragg's law with a minimal rotation. This approach reduces the experimental uncertainties in noisy XFEL data, improving the crystallographic R factors and sharpening anomalous differences that are near the level of the noise.

Flexible torsion-angle noncrystallographic symmetry restraints for improved macromolecular structure refinement.

One of the great challenges in refining macromolecular crystal structures is a low data-to-parameter ratio. Historically, knowledge from chemistry has been used to help to improve this ratio. When a macromolecule crystallizes with more than one copy in the asymmetric unit, the noncrystallographic symmetry relationships can be exploited to provide additional restraints when refining the working model. However, although globally similar, NCS-related chains often have local differences. To allow for local differences between NCS-related molecules, flexible torsion-based NCS restraints have been introduced, coupled with intelligent rotamer handling for protein chains, and are available in phenix.refine for refinement of models at all resolutions.

Automated identification of elemental ions in macromolecular crystal structures.

Many macromolecular model-building and refinement programs can automatically place solvent atoms in electron density at moderate-to-high resolution. This process frequently builds water molecules in place of elemental ions, the identification of which must be performed manually. The solvent-picking algorithms in phenix.refine have been extended to build common ions based on an analysis of the chemical environment as well as physical properties such as occupancy, B factor and anomalous scattering. The method is most effective for heavier elements such as calcium and zinc, for which a majority of sites can be placed with few false positives in a diverse test set of structures. At atomic resolution, it is observed that it can also be possible to identify tightly bound sodium and magnesium ions. A number of challenges that contribute to the difficulty of completely automating the process of structure completion are discussed.

Ligand placement based on prior structures: the guided ligand-replacement method.

The process of iterative structure-based drug design involves the X-ray crystal structure determination of upwards of 100 ligands with the same general scaffold (i.e. chemotype) complexed with very similar, if not identical, protein targets. In conjunction with insights from computational models and assays, this collection of crystal structures is analyzed to improve potency, to achieve better selectivity and to reduce liabilities such as absorption, distribution, metabolism, excretion and toxicology. Current methods for modeling ligands into electron-density maps typically do not utilize information on how similar ligands bound in related structures. Even if the electron density is of sufficient quality and resolution to allow de novo placement, the process can take considerable time as the size, complexity and torsional degrees of freedom of the ligands increase. A new module, Guided Ligand Replacement (GLR), was developed in Phenix to increase the ease and success rate of ligand placement when prior protein-ligand complexes are available. At the heart of GLR is an algorithm based on graph theory that associates atoms in the target ligand with analogous atoms in the reference ligand. Based on this correspondence, a set of coordinates is generated for the target ligand. GLR is especially useful in two situations: (i) modeling a series of large, flexible, complicated or macrocyclic ligands in successive structures and (ii) modeling ligands as part of a refinement pipeline that can automatically select a reference structure. Even in those cases for which no reference structure is available, if there are multiple copies of the bound ligand per asymmetric unit GLR offers an efficient way to complete the model after the first ligand has been placed. In all of these applications, GLR leverages prior knowledge from earlier structures to facilitate ligand placement in the current structure.

Automating crystallographic structure solution and refinement of protein-ligand complexes.

High-throughput drug-discovery and mechanistic studies often require the determination of multiple related crystal structures that only differ in the bound ligands, point mutations in the protein sequence and minor conformational changes. If performed manually, solution and refinement requires extensive repetition of the same tasks for each structure. To accelerate this process and minimize manual effort, a pipeline encompassing all stages of ligand building and refinement, starting from integrated and scaled diffraction intensities, has been implemented in Phenix. The resulting system is able to successfully solve and refine large collections of structures in parallel without extensive user intervention prior to the final stages of model completion and validation.

Accessing protein conformational ensembles using room-temperature X-ray crystallography.

Modern protein crystal structures are based nearly exclusively on X-ray data collected at cryogenic temperatures (generally 100 K). The cooling process is thought to introduce little bias in the functional interpretation of structural results, because cryogenic temperatures minimally perturb the overall protein backbone fold. In contrast, here we show that flash cooling biases previously hidden structural ensembles in protein crystals. By analyzing available data for 30 different proteins using new computational tools for electron-density sampling, model refinement, and molecular packing analysis, we found that crystal cryocooling remodels the conformational distributions of more than 35% of side chains and eliminates packing defects necessary for functional motions. In the signaling switch protein, H-Ras, an allosteric network consistent with fluctuations detected in solution by NMR was uncovered in the room-temperature, but not the cryogenic, electron-density maps. These results expose a bias in structural databases toward smaller, overpacked, and unrealistically unique models. Monitoring room-temperature conformational ensembles by X-ray crystallography can reveal motions crucial for catalysis, ligand binding, and allosteric regulation.

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.

Phaser.MRage: automated molecular replacement.

Phaser.MRage is a molecular-replacement automation framework that implements a full model-generation workflow and provides several layers of model exploration to the user. It is designed to handle a large number of models and can distribute calculations efficiently onto parallel hardware. In addition, phaser.MRage can identify correct solutions and use this information to accelerate the search. Firstly, it can quickly score all alternative models of a component once a correct solution has been found. Secondly, it can perform extensive analysis of identified solutions to find protein assemblies and can employ assembled models for subsequent searches. Thirdly, it is able to use a priori assembly information (derived from, for example, homologues) to speculatively place and score molecules, thereby customizing the search procedure to a certain class of protein molecule (for example, antibodies) and incorporating additional biological information into molecular replacement.

phenix.mr_rosetta: molecular replacement and model rebuilding with Phenix and Rosetta.

The combination of algorithms from the structure-modeling field with those of crystallographic structure determination can broaden the range of templates that are useful for structure determination by the method of molecular replacement. Automated tools in phenix.mr_rosetta simplify the application of these combined approaches by integrating Phenix crystallographic algorithms and Rosetta structure-modeling algorithms and by systematically generating and evaluating models with a combination of these methods. The phenix.mr_rosetta algorithms can be used to automatically determine challenging structures. The approaches used in phenix.mr_rosetta are described along with examples that show roles that structure-modeling can play in molecular replacement.

Towards automated crystallographic structure refinement with phenix.refine.

phenix.refine is a program within the PHENIX package that supports crystallographic structure refinement against experimental data with a wide range of upper resolution limits using a large repertoire of model parameterizations. It has several automation features and is also highly flexible. Several hundred parameters enable extensive customizations for complex use cases. Multiple user-defined refinement strategies can be applied to specific parts of the model in a single refinement run. An intuitive graphical user interface is available to guide novice users and to assist advanced users in managing refinement projects. X-ray or neutron diffraction data can be used separately or jointly in refinement. phenix.refine is tightly integrated into the PHENIX suite, where it serves as a critical component in automated model building, final structure refinement, structure validation and deposition to the wwPDB. This paper presents an overview of the major phenix.refine features, with extensive literature references for readers interested in more detailed discussions of the methods.