Program VUE: analysing distributions of cryo-EM projections using uniform spherical grids.

Three-dimensional cryo electron microscopy reconstructions are obtained by extracting information from a large number of projections of the object. These projections correspond to different 'views' or 'orientations', i.e. directions in which these projections show the reconstructed object. Uneven distribution of these views and the presence of dominating preferred orientations may distort the reconstructed spatial images. This work describes the program VUE (views on uniform grids for cryo electron microscopy), designed to study such distributions. Its algorithms, based on uniform virtual grids on a sphere, allow an easy calculation and accurate quantitative analysis of the frequency distribution of the views. The key computational element is the Lambert azimuthal equal-area projection of a spherical uniform grid onto a disc. This projection keeps the surface area constant and represents the frequency distribution with no visual bias. Since it has multiple tunable parameters, the program is easily adaptable to individual needs, and to the features of a particular project or of the figure to be produced. It can help identify problems related to an uneven distribution of views. Optionally, it can modify the list of projections, distributing the views more uniformly. The program can also be used as a teaching tool.

Accounting for nonuniformity of bulk-solvent: A mosaic model.

A flat mask-based model is almost universally used in macromolecular crystallography to account for disordered (bulk) solvent. This model assumes any voxel of the crystal unit cell that is not occupied by the atomic model is occupied by the solvent. The properties of this solvent are assumed to be exactly the same across the whole volume of the unit cell. While this is a reasonable approximation in practice, there are a number of scenarios where this model becomes suboptimal. In this work, we enumerate several of these scenarios and describe a new generalized approach to modeling the bulk-solvent which we refer to as mosaic bulk-solvent model. The mosaic bulk-solvent model allows nonuniform features of the solvent in the crystal to be accounted for in a computationally efficient way. It is implemented in the computational crystallography toolbox and the Phenix software.

Universal parameters of bulk-solvent masks.

The bulk solvent is a major component of biomacromolecular crystals that contributes significantly to the observed diffraction intensities. Accurate modelling of the bulk solvent has been recognized as important for many crystallographic calculations. Owing to its simplicity and modelling power, the flat (mask-based) bulk-solvent model is used by most modern crystallographic software packages to account for disordered solvent. In this model, the bulk-solvent contribution is defined by a binary mask and a scale (scattering) function. The mask is calculated on a regular grid using the atomic model coordinates and their chemical types. The grid step and two radii, solvent and shrinkage, are the three parameters that govern the mask calculation. They are highly correlated and their choice is a compromise between the computer time needed to calculate the mask and the accuracy of the mask. It is demonstrated here that this choice can be optimized using a unique value of 0.6 Šfor the grid step irrespective of the data resolution, and the radii values adjusted correspondingly. The improved values were tested on a large sample of Protein Data Bank entries derived from X-ray diffraction data and are now used in the computational crystallography toolbox (CCTBX) and in Phenix as the default choice.

Local heterogeneity analysis of crystallographic and cryo-EM maps using shell-approximation.

In X-ray crystallography and cryo-EM, experimental maps can be heterogeneous, showing different level of details in different regions. In this work we interpret heterogeneity in terms of two parameters, assigned individually for each atom, combining the conventional atomic displacement parameter with the resolution of the atomic image in the map. We propose a local real-space procedure to estimate the values of these heterogeneity parameters, assuming that a fragment of the density map and atomic positions are given. The procedure is based on an analytic representation of the atomic image, as a function of the inhomogeneity parameters and atomic coordinates. In this article, we report the results of the tests both with maps simulated and those derived from experimental data. For simulated maps containing regions with different resolutions, the method determines the local map resolution around the atomic centers and the values of the displacement parameter with reasonable accuracy. For experimental maps, obtained as a Fourier synthesis of a given global resolution, estimated values of the local resolution are close to the global one, and the values of the estimated displacement parameters are close to the respective values of the closest atoms in the refined model. Shown successful applications of the proposed method to experimental crystallographic and cryo-EM maps can be seen as a practical proof of method.

Efficient structure-factor modeling for crystals with multiple components.

Diffraction intensities from a crystallographic experiment include contributions from the entire unit cell of the crystal: the macromolecule, the solvent around it and eventually other compounds. These contributions cannot typically be well described by an atomic model alone, i.e. using point scatterers. Indeed, entities such as disordered (bulk) solvent, semi-ordered solvent (e.g. lipid belts in membrane proteins, ligands, ion channels) and disordered polymer loops require other types of modeling than a collection of individual atoms. This results in the model structure factors containing multiple contributions. Most macromolecular applications assume two-component structure factors: one component arising from the atomic model and the second one describing the bulk solvent. A more accurate and detailed modeling of the disordered regions of the crystal will naturally require more than two components in the structure factors, which presents algorithmic and computational challenges. Here an efficient solution of this problem is proposed. All algorithms described in this work have been implemented in the computational crystallography toolbox (CCTBX) and are also available within Phenix software. These algorithms are rather general and do not use any assumptions about molecule type or size nor about those of its components.

Direct calculation of cryo-EM and crystallographic model maps for real-space refinement.

This work addresses the problem of the calculation of limited-resolution maps from an atomic model in cryo-electron microscopy and in X-ray and neutron crystallography, including cases where the resolution varies from one molecular region to another. Such maps are necessary in real-space refinement for comparison with the experimental maps. For an appropriate numeric comparison, the calculated maps should reproduce not only the structural features contained in the experimental maps but also the principal map distortions. These model maps can be obtained with no use of Fourier transforms but, similar to density distributions, as a sum of individual atomic contributions. Such contributions, referred to as atomic density images, are atomic densities morphed to reflect distortions of the experimental map, in particular the loss of resolution. They are described by functions composed of a central peak surrounded by Fourier ripples. For practical calculations, atomic images should be cut at some distance. It is shown that to reach a reasonable accuracy such a distance should be significantly larger than the distance customarily applied when calculating density distributions. This is a consequence of the slow rate with which the amplitude of the Fourier ripples decreases. Such a large distance means that at least a few ripples should be included in calculations in order to obtain a map that is sufficiently accurate. Oscillating functions describing these atomic contributions depend, for a given atomic type, on the resolution and on the atomic displacement parameter values. To express both the central peak and the Fourier ripples of the atomic images, these functions are represented by the sums of especially designed terms, each concentrated in a spherical shell and depending analytically on the atomic parameters. In this work, the strength of the dependence of the accuracy of resulting map on the accuracy of the atomic displacement parameters and on the truncation distance, i.e. the number of ripples included in atomic density images, is analyzed. This analysis is completed by practical aspects of the calculation of maps of inhomogeneous resolution. Tests show that the calculation of limited-resolution maps from an atomic model as a sum of atomic contributions requires a large truncation radius extending beyond the central peak of an atomic image and the first Fourier ripples. The article discusses the practical details of such calculations expressing atomic contributions as analytic functions of the atomic coordinates, the atomic displacement parameters and the local resolution.

Analytic modeling of inhomogeneous-resolution maps in cryo-electron microscopy and crystallography.

Refinement of macromolecular atomic models versus experimental maps in crystallography and cryo-electron microscopy is a critical step in structure solution. For an appropriate comparison, model maps should mimic the imperfections in the experimental maps, mainly atomic disorder and limited resolution, which are often inhomogeneous over the molecular region. In the suggested method, these model maps are calculated as the sum of atomic contributions expressed through a specifically designed function describing a solitary spherical wave. Thanks to this function, atomic contributions are analytically expressed through their atomic displacement parameter and local resolution, a value now associated with each atom. Such a full analytic dependence of inhomogeneous-resolution map values on model parameters permits the refinement of all of these parameters together.

Toward Real Real-Space Refinement of Atomic Models.

High-quality atomic models providing structural information are the results of their refinement versus diffraction data (reciprocal-space refinement), or versus experimental or experimentally based maps (real-space refinement). A proper real-space refinement can be achieved by comparing such a map with a map calculated from the atomic model. Similar to density distributions, the maps of a limited and even inhomogeneous resolution can also be calculated as sums of terms, known as atomic images, which are three-dimensional peaky functions surrounded by Fourier ripples. These atomic images and, consequently, the maps for the respective models, can be expressed analytically as functions of coordinates, atomic displacement parameters, and the local resolution. This work discusses the practical feasibility of such calculation for the real-space refinement of macromolecular atomic models.

Implementation of the riding hydrogen model in CCTBX to support the next generation of X-ray and neutron joint refinement in Phenix.

A fundamental prerequisite for implementing new procedures of atomic model refinement against neutron diffraction data is the efficient handling of hydrogen atoms. The riding hydrogen model, which constrains hydrogen atom parameters to those of the non-hydrogen atoms, is a plausible parameterization for refinements. This work describes the implementation of the riding hydrogen model in the Computational Crystallography Toolbox and in Phenix. Riding hydrogen atoms can be found in several different configurations that are characterized by specific geometries. For each configuration, the hydrogen atom parameterization and the expressions for the gradients of refinement target function with respect to non-hydrogen parameters are described.

Including crystallographic symmetry in quantum-based refinement: Q|R#2.

Three-dimensional structure models refined using low-resolution data from crystallographic or electron cryo-microscopy experiments can benefit from high-quality restraints derived from quantum-chemical methods. However, nonperiodic atom-centered quantum-chemistry codes do not inherently account for nearest-neighbor interactions of crystallographic symmetry-related copies in a satisfactory way. Here, these nearest-neighbor effects have been included in the model by expanding to a super-cell and then truncating the super-cell to only include residues from neighboring cells that are interacting with the asymmetric unit. In this way, the fragmentation approach can adequately and efficiently include nearest-neighbor effects. It has previously been shown that a moderately sized X-ray structure can be treated using quantum methods if a fragmentation approach is applied. In this study, a target protein (PDB entry 4gif) was partitioned into a number of large fragments. The use of large fragments (typically hundreds of atoms) is tractable when a GPU-based package such as TeraChem is employed or cheaper (semi-empirical) methods are used. The QM calculations were run at the HF-D3/6-31G level. The models refined using a recently developed semi-empirical method (GFN2-xTB) were compared and contrasted. To validate the refinement procedure for a non-P1 structure, a standard set of crystallographic metrics were used. The robustness of the implementation is shown by refining 13 additional protein models across multiple space groups and a summary of the refinement metrics is presented.

Macromolecular structure determination using X-rays, neutrons and electrons: recent developments in Phenix.

Diffraction (X-ray, neutron and electron) and electron cryo-microscopy are powerful methods to determine three-dimensional macromolecular structures, which are required to understand biological processes and to develop new therapeutics against diseases. The overall structure-solution workflow is similar for these techniques, but nuances exist because the properties of the reduced experimental data are different. Software tools for structure determination should therefore be tailored for each method. Phenix is a comprehensive software package for macromolecular structure determination that handles data from any of these techniques. Tasks performed with Phenix include data-quality assessment, map improvement, model building, the validation/rebuilding/refinement cycle and deposition. Each tool caters to the type of experimental data. The design of Phenix emphasizes the automation of procedures, where possible, to minimize repetitive and time-consuming manual tasks, while default parameters are chosen to encourage best practice. A graphical user interface provides access to many command-line features of Phenix and streamlines the transition between programs, project tracking and re-running of previous tasks.

Mask-based approach to phasing of single-particle diffraction data. II. Likelihood-based selection criteria.

A new type of mask-selection criterion is suggested for mask-based phasing. In this phasing approach, a large number of connected molecular masks are randomly generated. Structure-factor phases corresponding to a trial mask are accepted as an admissible solution of the phase problem if the mask satisfies some specified selection rules that are key to success. The admissible phase sets are aligned and averaged to give a preliminary solution of the phase problem. The new selection rule is based on the likelihood of the generated mask. It is defined as the probability of reproducing the observed structure-factor magnitudes by placing atoms randomly into the mask. While the result of the direct comparison of mask structure-factor magnitudes with observed ones using a correlation coefficient is highly dominated by a few very strong low-resolution reflections, a new method gives higher weight to relatively weak high-resolution reflections that allows them to be phased accurately. This mask-based phasing procedure with likelihood-based selection has been applied to simulated single-particle diffraction data of the photosystem II monomer. The phase set obtained resulted in a 16 Šresolution Fourier synthesis (more than 4000 reflections) with 98% correlation with the exact phase set and 69% correlation for about 2000 reflections in the highest resolution shell (20-16 Å). This work also addresses another essential problem of phasing methods, namely adequate estimation of the resolution achieved. A model-trapping analysis of the phase sets obtained by the mask-based phasing procedure suggests that the widely used `50% shell correlation' criterion may be too optimistic in some cases.

New tools for the analysis and validation of cryo-EM maps and atomic models.

Recent advances in the field of electron cryomicroscopy (cryo-EM) have resulted in a rapidly increasing number of atomic models of biomacromolecules that have been solved using this technique and deposited in the Protein Data Bank and the Electron Microscopy Data Bank. Similar to macromolecular crystallography, validation tools for these models and maps are required. While some of these validation tools may be borrowed from crystallography, new methods specifically designed for cryo-EM validation are required. Here, new computational methods and tools implemented in PHENIX are discussed, including d99 to estimate resolution, phenix.auto_sharpen to improve maps and phenix.mtriage to analyze cryo-EM maps. It is suggested that cryo-EM half-maps and masks should be deposited to facilitate the evaluation and validation of cryo-EM-derived atomic models and maps. The application of these tools to deposited cryo-EM atomic models and maps is also presented.

From deep TLS validation to ensembles of atomic models built from elemental motions. II. Analysis of TLS refinement results by explicit interpretation.

TLS modelling was developed by Schomaker and Trueblood to describe atomic displacement parameters through concerted (rigid-body) harmonic motions of an atomic group [Schomaker & Trueblood (1968), Acta Cryst. B24, 63-76]. The results of a TLS refinement are T, L and S matrices that provide individual anisotropic atomic displacement parameters (ADPs) for all atoms belonging to the group. These ADPs can be calculated analytically using a formula that relates the elements of the TLS matrices to atomic parameters. Alternatively, ADPs can be obtained numerically from the parameters of concerted atomic motions corresponding to the TLS matrices. Both procedures are expected to produce the same ADP values and therefore can be used to assess the results of TLS refinement. Here, the implementation of this approach in PHENIX is described and several illustrations, including the use of all models from the PDB that have been subjected to TLS refinement, are provided.

Real-space refinement in PHENIX for cryo-EM and crystallography.

This article describes the implementation of real-space refinement in the phenix.real_space_refine program from the PHENIX suite. The use of a simplified refinement target function enables very fast calculation, which in turn makes it possible to identify optimal data-restraint weights as part of routine refinements with little runtime cost. Refinement of atomic models against low-resolution data benefits from the inclusion of as much additional information as is available. In addition to standard restraints on covalent geometry, phenix.real_space_refine makes use of extra information such as secondary-structure and rotamer-specific restraints, as well as restraints or constraints on internal molecular symmetry. The re-refinement of 385 cryo-EM-derived models available in the Protein Data Bank at resolutions of 6 Šor better shows significant improvement of the models and of the fit of these models to the target maps.

Aminoglycoside interactions and impacts on the eukaryotic ribosome.

Aminoglycosides are chemically diverse, broad-spectrum antibiotics that target functional centers within the bacterial ribosome to impact all four principle stages (initiation, elongation, termination, and recycling) of the translation mechanism. The propensity of aminoglycosides to induce miscoding errors that suppress the termination of protein synthesis supports their potential as therapeutic interventions in human diseases associated with premature termination codons (PTCs). However, the sites of interaction of aminoglycosides with the eukaryotic ribosome and their modes of action in eukaryotic translation remain largely unexplored. Here, we use the combination of X-ray crystallography and single-molecule FRET analysis to reveal the interactions of distinct classes of aminoglycosides with the 80S eukaryotic ribosome. Crystal structures of the 80S ribosome in complex with paromomycin, geneticin (G418), gentamicin, and TC007, solved at 3.3- to 3.7-Å resolution, reveal multiple aminoglycoside-binding sites within the large and small subunits, wherein the 6'-hydroxyl substituent in ring I serves as a key determinant of binding to the canonical eukaryotic ribosomal decoding center. Multivalent binding interactions with the human ribosome are also evidenced through their capacity to affect large-scale conformational dynamics within the pretranslocation complex that contribute to multiple aspects of the translation mechanism. The distinct impacts of the aminoglycosides examined suggest that their chemical composition and distinct modes of interaction with the ribosome influence PTC read-through efficiency. These findings provide structural and functional insights into aminoglycoside-induced impacts on the eukaryotic ribosome and implicate pleiotropic mechanisms of action beyond decoding.

The integrative role of cryo electron microscopy in molecular and cellular structural biology.

After gradually moving away from preparation methods prone to artefacts such as plastic embedding and negative staining for cell sections and single particles, the field of cryo electron microscopy (cryo-EM) is now heading off at unprecedented speed towards high-resolution analysis of biological objects of various sizes. This 'revolution in resolution' is happening largely thanks to new developments of new-generation cameras used for recording the images in the cryo electron microscope which have much increased sensitivity being based on complementary metal oxide semiconductor devices. Combined with advanced image processing and 3D reconstruction, the cryo-EM analysis of nucleoprotein complexes can provide unprecedented insights at molecular and atomic levels and address regulatory mechanisms in the cell. These advances reinforce the integrative role of cryo-EM in synergy with other methods such as X-ray crystallography, fluorescence imaging or focussed-ion beam milling as exemplified here by some recent studies from our laboratory on ribosomes, viruses, chromatin and nuclear receptors. Such multi-scale and multi-resolution approaches allow integrating molecular and cellular levels when applied to purified or in situ macromolecular complexes, thus illustrating the trend of the field towards cellular structural biology.

From deep TLS validation to ensembles of atomic models built from elemental motions. Addenda and corrigendum.

Researcher feedback has indicated that in Urzhumtsev et al. [(2015) Acta Cryst. D71, 1668-1683] clarification of key parts of the algorithm for interpretation of TLS matrices in terms of elemental atomic motions and corresponding ensembles of atomic models is required. Also, it has been brought to the attention of the authors that the incorrect PDB code was reported for one of test models. These issues are addressed in this article.

Mask-based approach to phasing of single-particle diffraction data.

A Monte Carlo-type approach for low- and medium-resolution phasing of single-particle diffraction data is suggested. Firstly, the single-particle phase problem is substituted with the phase problem for an imaginary crystal. A unit cell of this crystal contains a single isolated particle surrounded by a large volume of bulk solvent. The developed phasing procedure then generates a large number of connected and finite molecular masks, calculates their Fourier coefficients, selects the sets with magnitudes that are highly correlated with the experimental values and finally aligns the selected phase sets and calculates the averaged phase values. A test with the known structure of monomeric photosystem II resulted in phases that have 97% correlation with the exact phases in the full 25 Å resolution shell (1054 structure factors) and correlations of 99, 94, 81 and 79% for the resolution shells ∞-60, 60-40, 40-30 and 30-25 Å, respectively. The same procedure may be used for crystallographic ab initio phasing.

Helical Symmetry of Nucleic Acids: Obstacle or Help in Structure Solution?

Crystallographic molecular replacement method is the key tool to define an atomic structure of nucleic acids. Frequently nucleic acids are packed forming continuous helices in the crystal. This arrangement of individual molecules in "infinite" pseudo helical structures in crystal may be the reason why the molecular replacement fails to find a unique position of the search atomic model as the method requires. The Patterson function, calculated as a Fourier series with diffraction intensities, has auxiliary peaks for such a molecular packing. Those near the origin peak indicate the orientation of the helices. The coordinates of other peaks are related to the molecular position and the rotation angle between two such "infinite" helices. Thus, the peak analysis allows getting molecular position even without a search model. An intelligent selecting and averaging of the phase sets corresponding to multiple probable positions of the search model again result in a unique solution but in the form of a Fourier synthesis and not a model. This synthesis can be used then to build an atomic model as it is the case for usual phasing methods.