Search

The employment of directed acyclic graphs to advance the tracking, control and appraisal of crystallographic phasing strategies is discussed., Crystallographic phasing strategies increasingly require the exploration and ranking of many hypotheses about the number, types and positions of atoms, molecules and/or molecular fragments in the unit cell, each with only a small chance of being correct. Accelerating this move has been improvements in phasing methods, which are now able to extract phase information from the placement of very small fragments of structure, from weak experimental phasing signal or from combinations of molecular replacement and experimental phasing information. Describing phasing in terms of a directed acyclic graph allows graph-management software to track and manage the path to structure solution. The crystallographic software supporting the graph data structure must be strictly modular so that nodes in the graph are efficiently generated by the encapsulated functionality. To this end, the development of new software, Phasertng, which uses directed acyclic graphs natively for input/output, has been initiated. In Phasertng, the codebase of Phaser has been rebuilt, with an emphasis on modularity, on scripting, on speed and on continuing algorithm development. As a first application of phasertng, its advantages are demonstrated in the context of phasertng.xtricorder, a tool to analyse and triage merged data in preparation for molecular replacement or experimental phasing. The description of the phasing strategy with directed acyclic graphs is a generalization that extends beyond the functionality of Phasertng, as it can incorporate results from bioinformatics and other crystallographic tools, and will facilitate multifaceted search strategies, dynamic ranking of alternative search pathways and the exploitation of machine learning to further improve phasing strategies.

The prerequisites for density modification of maps from electron cryomicroscopy are examined and a procedure for incorporating model-based information is presented., Density modification uses expectations about features of a map such as a flat solvent and expected distributions of density in the region of the macromolecule to improve individual Fourier terms representing the map. This process transfers information from one part of a map to another and can improve the accuracy of a map. Here, the assumptions behind density modification for maps from electron cryomicroscopy are examined and a procedure is presented that allows the incorporation of model-based information. Density modification works best in cases where unfiltered, unmasked maps with clear boundaries between the macromolecule and solvent are visible, and where there is substantial noise in the map, both in the region of the macromolecule and the solvent. It also is most effective if the characteristics of the map are relatively constant within regions of the macromolecule and the solvent. Model-based information can be used to improve density modification, but model bias can in principle occur. Here, model bias is reduced by using ensemble models that allow an estimation of model uncertainty. A test of model bias is presented that suggests that even if the expected density in a region of a map is specified incorrectly by using an incorrect model, the incorrect expectations do not strongly affect the final map.

AMPylation is an inactivating modification that matches the activity of the major endoplasmic reticulum (ER) chaperone BiP to the burden of unfolded proteins. A single ER-localised Fic protein, FICD (HYPE), catalyses both AMPylation and deAMPylation of BiP. However, the basis for the switch in FICD’s activity is unknown. We report on the transition of FICD from a dimeric enzyme, that deAMPylates BiP, to a monomer with potent AMPylation activity. Mutations in the dimer interface or in residues tracing an inhibitory relay from the dimer interface to the enzyme’s active site favour BiP AMPylation in vitro and in cells. Mechanistically, monomerisation relieves a repressive effect allosterically-propagated from the dimer interface to the inhibitory Glu234, thereby permitting AMPylation-competent binding of MgATP. Whereas, a reciprocal signal propagated from the nucleotide binding site, provides a mechanism for coupling the oligomeric-state and enzymatic activity of FICD to the energy status of the ER.Impact StatementUnique amongst known chaperones, the endoplasmic reticulum (ER)-localized Hsp70, BiP, is subject to transient inactivation under conditions of low ER stress by reversible, covalent modification – AMPylation. The enzyme responsible for this modification, FICD, is in fact a bifunctional enzyme with a single active site capable of both AMPylation and deAMPylation. Here we elucidate, by biochemical, biophysical and structural means, the mechanism by which this enzyme is able to switch enzymatic modality: by regulation of its oligomeric state. The oligomeric state-dependent reciprocal regulation of FICD activity is, in turn, sensitive to the ATP/ADP ratio. This allosteric pathway potentially facilitates the sensing of unfolded protein load in the ER and permits the transduction of this signal into a post-translational buffering of ER chaperone activity.

A recent publication by Seng et al. in this journal reports the crystallographic structure of refolded, full-length SMN protein and two disease-relevant derivatives thereof. Here, we would like to suggest that at least two of the structures reported in that study are incorrect. We present evidence that one of the associated crystallographic datasets is derived from a crystal of the bacterial Sm-like protein Hfq and that a second dataset is derived from a crystal of the bacterial Gab protein. Both proteins are frequent contaminants of bacterially overexpressed proteins which might have been co-purified during metal affinity chromatography. A third structure presented in the Seng et al. paper cannot be examined further because neither the atomic coordinates, nor the diffraction intensities were made publicly available. The Tudor domain protein SMN has been shown to be a component of the SMN complex, which mediates the assembly of RNA-protein complexes of uridine-rich small nuclear ribonucleoproteins (UsnRNPs). Importantly, this activity is reduced in SMA patients, raising the possibility that the aetiology of SMA is linked to RNA metabolism. Structural studies on diverse components of the SMN complex, including fragments of SMN itself have contributed greatly to our understanding of the cellular UsnRNP assembly machinery. Yet full-length SMN has so far evaded structural elucidation. The Seng et al. study claimed to have closed this gap, but based on the results presented here, the only conclusion that can be drawn is that the Seng et al. study is largely invalid and should be retracted from the literature. published

The Z mutation (E342K) of α1-antitrypsin (α1-AT), carried by 4% of Northern Europeans, predisposes to early onset of emphysema due to decreased functional α1-AT in the lung and to liver cirrhosis due to accumulation of polymers in hepatocytes. However, it remains unclear why the Z mutation causes intracellular polymerization of nascent Z α1-AT and why 15% of the expressed Z α1-AT is secreted into circulation as functional, but polymerogenic, monomers. Here, we solve the crystal structure of the Z-monomer and have engineered replacements to assess the conformational role of residue Glu-342 in α1-AT. The results reveal that Z α1-AT has a labile strand 5 of the central β-sheet A (s5A) with a consequent equilibrium between a native inhibitory conformation, as in its crystal structure here, and an aberrant conformation with s5A only partially incorporated into the central β-sheet. This aberrant conformation, induced by the loss of interactions from the Glu-342 side chain, explains why Z α1-AT is prone to polymerization and readily binds to a 6-mer peptide, and it supports that annealing of s5A into the central β-sheet is a crucial step in the serpins' metastable conformational formation. The demonstration that the aberrant conformation can be rectified through stabilization of the labile s5A by binding of a small molecule opens a potential therapeutic approach for Z α1-AT deficiency.

Descriptions are given of the maximum-likelihood gyre method implemented in Phaser for optimizing the orientation and relative position of rigid-body fragments of a model after the orientation of the model has been identified, but before the model has been positioned in the unit cell, and also the related gimble method for the refinement of rigid-body fragments of the model after positioning. Gyre refinement helps to lower the root-mean-square atomic displacements between model and target molecular-replacement solutions for the test case of antibody Fab(26-10) and improves structure solution with ARCIMBOLDO_SHREDDER., This research was supported by the Wellcome Trust (Principal Research Fellowship to RJR, grant 082961/Z/07/Z) and by grant BB/L006014/1 from the BBSRC, UK. The research was facilitated by Wellcome Trust Strategic Award 100140 to the Cambridge Institute for Medical Research. MS and CM received financial support from CCP4 for a sabbatical in the group of RJR. CM is grateful to MINECO for her BES-2015- 071397 scholarship associated with the Structural Biology Maria de Maeztu Unit of Excellence. IU was supported by grants BIO2015-64216-P, BIO2013-49604-EXP and MDM2014-0435-01 (the Spanish Ministry of Economy and Competitiveness) and Generalitat de Catalunya (2014SGR997).

The Editors of Acta Cryst. D have received a Letter to the Editor from Alex Wlodawer and Zbigniew Dauter (Wlodawer & Dauter, 2017[Wlodawer, A. & Dauter, Z. (2017). Acta Cryst. D73, 379-380.]) concerning the use of the term atomic resolution. As this is an important topic, a number of representatives from the structural biology community were asked to respond to the Letter and their comments are given here.

A genetic algorithm has been developed to optimize the phases of the strongest reflections in SIR/SAD data. This is shown to facilitate density modification and model building in several test cases., Experimental phasing of diffraction data from macromolecular crystals involves deriving phase probability distributions. These distributions are often bimodal, making their weighted average, the centroid phase, improbable, so that electron-density maps computed using centroid phases are often non-­interpretable. Density modification brings in information about the characteristics of electron density in protein crystals. In successful cases, this allows a choice between the modes in the phase probability distributions, and the maps can cross the borderline between non-interpretable and interpretable. Based on the suggestions by Vekhter [Vekhter (2005 ▶), Acta Cryst. D61, 899–902], the impact of identifying optimized phases for a small number of strong reflections prior to the density-modification process was investigated while using the centroid phase as a starting point for the remaining reflections. A genetic algorithm was developed that optimizes the quality of such phases using the skewness of the density map as a target function. Phases optimized in this way are then used in density modification. In most of the tests, the resulting maps were of higher quality than maps generated from the original centroid phases. In one of the test cases, the new method sufficiently improved a marginal set of experimental SAD phases to enable successful map interpretation. A computer program, SISA, has been developed to apply this method for phase improvement in macromolecular crystallo­graphy.

The statistical effects of translational noncrystallographic symmetry can be characterized by maximizing parameters describing the noncrystallographic symmetry in a likelihood function, thereby unmasking the competing statistical effects of twinning., In the case of translational noncrystallographic symmetry (tNCS), two or more copies of a component in the asymmetric unit of the crystal are present in a similar orientation. This causes systematic modulations of the reflection intensities in the diffraction pattern, leading to problems with structure determination and refinement methods that assume, either implicitly or explicitly, that the distribution of intensities is a function only of resolution. To characterize the statistical effects of tNCS accurately, it is necessary to determine the translation relating the copies, any small rotational differences in their orientations, and the size of random coordinate differences caused by conformational differences. An algorithm to estimate these parameters and refine their values against a likelihood function is presented, and it is shown that by accounting for the statistical effects of tNCS it is possible to unmask the competing statistical effects of twinning and tNCS and to more robustly assess the crystal for the presence of twinning.

Crystallographic studies of ligands bound to biological macromolecules (proteins and nucleic acids) represent\ud an important source of information concerning drug-target interactions, providing atomic level insights\ud into the physical chemistry of complex formation between macromolecules and ligands. Of the\ud more than 115,000 entries extant in the Protein Data Bank (PDB) archive, ~75% include at least one non-polymeric\ud ligand. Ligand geometrical and stereochemical quality, the suitability of ligand models for in silico drug\ud discovery and design, and the goodness-of-fit of ligand models to electron-density maps vary widely across\ud the archive. We describe the proceedings and conclusions from the first Worldwide PDB/Cambridge Crystallographic\ud Data Center/Drug Design Data Resource (wwPDB/CCDC/D3R) Ligand Validation Workshop\ud held at the Research Collaboratory for Structural Bioinformatics at Rutgers University on July 30–31, 2015.\ud Experts in protein crystallography from academe and industry came together with non-profit and for-profit\ud software providers for crystallography and with experts in computational chemistry and data archiving to\ud discuss and make recommendations on best practices, as framed by a series of questions central to structural\ud studies of macromolecule-ligand complexes. What data concerning bound ligands should be archived\ud in the PDB? How should the ligands be best represented? How should structural models of macromoleculeligand\ud complexes be validated? What supplementary information should accompany publications of structural\ud studies of biological macromolecules? Consensus recommendations on best practices developed in\ud response to each of these questions are provided, together with some details regarding implementation.\ud Important issues addressed but not resolved at the workshop are also enumerated.

The foundations and current features of a widely used graphical user interface for macromolecular crystallography are described., A new Python-based graphical user interface for the PHENIX suite of crystallography software is described. This interface unifies the command-line programs and their graphical displays, simplifying the development of new interfaces and avoiding duplication of function. With careful design, graphical interfaces can be displayed automatically, instead of being manually constructed. The resulting package is easily maintained and extended as new programs are added or modified.

SAD data can be used in Phaser to solve novel structures, supplement molecular-replacement phase information or identify anomalous scatterers from a final refined model., Phaser is a program that implements likelihood-based methods to solve macromolecular crystal structures, currently by molecular replacement or single-wavelength anomalous diffraction (SAD). SAD phasing is based on a likelihood target derived from the joint probability distribution of observed and calculated pairs of Friedel-related structure factors. This target combines information from the total structure factor (primarily non-anomalous scattering) and the difference between the Friedel mates (anomalous scattering). Phasing starts from a substructure, which is usually but not necessarily a set of anomalous scatterers. The substructure can also be a protein model, such as one obtained by molecular replacement. Additional atoms are found using a log-likelihood gradient map, which shows the sites where the addition of scattering from a particular atom type would improve the likelihood score. An automated completion algorithm adds new sites, choosing optionally among different atom types, adds anisotropic B-factor parameters if appropriate and deletes atoms that refine to low occupancy. Log-likelihood gradient maps can also identify which atoms in a refined protein structure are anomalous scatterers, such as metal or halide ions. These maps are more sensitive than conventional model-phased anomalous difference Fouriers and the iterative completion algorithm is able to find a significantly larger number of convincing sites.

Ten measures of experimental electron-density-map quality are examined and the skewness of electron density is found to be the best indicator of actual map quality. A Bayesian approach to estimating map quality is developed and used in the PHENIX AutoSol wizard to make decisions during automated structure solution., Estimates of the quality of experimental maps are important in many stages of structure determination of macromolecules. Map quality is defined here as the correlation between a map and the corresponding map obtained using phases from the final refined model. Here, ten different measures of experimental map quality were examined using a set of 1359 maps calculated by re-analysis of 246 solved MAD, SAD and MIR data sets. A simple Bayesian approach to estimation of map quality from one or more measures is presented. It was found that a Bayesian estimator based on the skewness of the density values in an electron-density map is the most accurate of the ten individual Bayesian estimators of map quality examined, with a correlation between estimated and actual map quality of 0.90. A combination of the skewness of electron density with the local correlation of r.m.s. density gives a further improvement in estimating map quality, with an overall correlation coefficient of 0.92. The PHENIX AutoSol wizard carries out automated structure solution based on any combination of SAD, MAD, SIR or MIR data sets. The wizard is based on tools from the PHENIX package and uses the Bayesian estimates of map quality described here to choose the highest quality solutions after experimental phasing.

The MBT repeat has been recently identified as a key domain capable of methyl-lysine histone recognition. Functional work has pointed to a role for MBT domain-containing proteins in transcriptional repression of developmental control genes such as Hox genes. In this study, L3MBTL2, a human homolog of Drosophila Sfmbt critical for Hox gene silencing, is demonstrated to preferentially recognize lower methylation states of several histone-derived peptides through its fourth MBT repeat. High-resolution crystallographic analysis of the four MBT repeats of this protein reveals its unique asymmetric rhomboid architecture, as well as binding mechanism, which preclude the interaction of the first three MBT repeats with methylated peptides. Structural elucidation of an L3MBTL2-H4K20me1 complex and comparison with other MBT-histone peptide complexes also suggests that an absence of distinct surface contours surrounding the methyl-lysine-binding pocket may underlie the lack of sequence specificity observed for members of this protein family.

The highly automated PHENIX AutoBuild wizard is described. The procedure can be applied equally well to phases derived from isomorphous/anomalous and molecular-replacement methods., The PHENIX AutoBuild wizard is a highly automated tool for iterative model building, structure refinement and density modification using RESOLVE model building, RESOLVE statistical density modification and phenix.refine structure refinement. Recent advances in the AutoBuild wizard and phenix.refine include automated detection and application of NCS from models as they are built, extensive model-completion algorithms and automated solvent-molecule picking. Model-completion algorithms in the AutoBuild wizard include loop building, crossovers between chains in different models of a structure and side-chain optimization. The AutoBuild wizard has been applied to a set of 48 structures at resolutions ranging from 1.1 to 3.2 Å, resulting in a mean R factor of 0.24 and a mean free R factor of 0.29. The R factor of the final model is dependent on the quality of the starting electron density and is relatively independent of resolution.

A foreign-DNA–destroying machine Bacteria have an adaptive immune system, called CRISPR, that identifies invading viruses through their DNA or RNA sequences and cuts them up (see the Perspective by Zhang and Sontheimer). Jackson et al. and Mulepati et al. have determined the structure of the large protein complex, called Cascade, that targets the invading nucleic acids and does the cutting. The seahorse-shaped structure reveals how the 11 subcomponents of Cascade assemble into the final protein complex. The structure also shows how Cascade presents the short CRISPR-derived RNAs so that they can bind and target foreign DNA. Science , this issue p. 1473 and p. 1479 ; see also p. 1452

The crystal structure of the region spanning residues 95–146 of the rotavirus nonstructural protein NSP4 from the asymptomatic human strain ST3 was determined at a resolution of 2.5 Å. Severe diffraction anisotropy, rotational pseudosymmetry and twinning complicated the refinement of this structure. A systematic explanation confirming the crystal pathologies and describing how the structure was successfully refined is given in this report.

An overview of the CCP4 software suite for macromolecular crystallography is given., The CCP4 (Collaborative Computational Project, Number 4) software suite is a collection of programs and associated data and software libraries which can be used for macromolecular structure determination by X-ray crystallography. The suite is designed to be flexible, allowing users a number of methods of achieving their aims. The programs are from a wide variety of sources but are connected by a common infrastructure provided by standard file formats, data objects and graphical interfaces. Structure solution by macromolecular crystallo­graphy is becoming increasingly automated and the CCP4 suite includes several automation pipelines. After giving a brief description of the evolution of CCP4 over the last 30 years, an overview of the current suite is given. While detailed descriptions are given in the accompanying articles, here it is shown how the individual programs contribute to a complete software package.

The three-dimensional structure of glycosomal glyceraldehyde-3-phosphate dehydrogenase [D-glyceraldehyde-3-phosphate:NAD+ oxidoreductase (phosphorylating), EC 1.12.1.12] from the sleeping-sickness parasite Trypanosoma brucei was solved by molecular replacement at 3.2-angstrom resolution with an x-ray data set collected by the Laue method. For data collection, three crystals were exposed to the polychromatic synchrotron x-ray beam for a total of 20.5 sec. The structure was solved by using the Bacillus stearothermophilus enzyme model [Skarzynski, T., Moody, P. C. E. & Wonacott, A. J. (1987) J. Mol. Biol. 193, 171-187] with a partial data set which was 37% complete. The crystals contain six subunits per asymmetric unit, which allowed us to overcome the absence of >60% of the reflections by 6-fold density averaging. After molecular dynamics refinement, the current molecular model has an R factor of 17.6%. Comparing the structure of the trypanosome enzyme with that of the homologous human muscle enzyme, which was determined at 2.4-angstrom resolution, reveals important structural differences in the NAD binding region. These are of great interest for the design of specific inhibitors of the parasite enzyme.

Structures of biomolecular systems are increasingly computed by integrative modeling that relies on varied types of experimental data and theoretical information. We describe here the proceedings and conclusions from the first wwPDB Hybrid/Integrative Methods Task Force Workshop held at the European Bioinformatics Institute in Hinxton, UK, on October 6 and 7, 2014. At the workshop, experts in various experimental fields of structural biology, experts in integrative modeling and visualization, and experts in data archiving addressed a series of questions central to the future of structural biology. How should integrative models be represented? How should the data and integrative models be validated? What data should be archived? How should the data and models be archived? What information should accompany the publication of integrative models? The workshop was supported by funding to PDBe by Wellcome Trust 088944; RCSB PDB by NSF DBI 1338415; PDBj by JST-NBDC; BMRB by NLM P41 LM05799; EMDataBank by NIH GM079429; and tax-deductible donations made to the wwPDB Foundation in support of wwPDB outreach activities.

Several sets of amino acid surface areas and transfer free energies were used to derive a total of nine sets of atomic solvation parameters (ASPs). We tested the accuracy of each of these sets of parameters in predicting the experimentally determined transfer free energies of the amino acid derivatives from which the parameters were derived. In all cases, the calculated and experimental values correlated well. We then chose three parameter sets and examined the effect of adding an energetic correction for desolvation based on these three parameter sets to the simple potential function used in our multiple start Monte Carlo docking method. A variety of protein-protein interactions and docking results were examined. In the docking simulations studied, the desolvation correction was only applied during the final energy calculation of each simulation. For most of the docking results we analyzed, the use of an octanol-water-based ASP set marginally improved the energetic ranking of the low-energy dockings, whereas the other ASP sets we tested disturbed the ranking of the low-energy dockings in many of the same systems. We also examined the correlation between the experimental free energies of association and our calculated interaction energies for a series of proteinase-inhibitor complexes. Again, the octanol-water-based ASP set was compatible with our standard potential function, whereas ASP sets derived from other solvent systems were not.

Corticosteroids are transported in the blood by a serpin, corticosteroid-binding globulin (CBG), and their normally equilibrated release can be further triggered by the cleavage of the reactive loop of CBG. We report here the crystal structures of cleaved human CBG (cCBG) at 1.8-A resolution and its complex with cortisol at 2.3-A resolution. As expected, on cleavage, CBG undergoes the irreversible S-to-R serpin transition, with the cleaved reactive loops being fully incorporated into the central β-sheet. A connecting loop of helix D, which is in a helix-like conformation in native CBG, unwinds and grossly perturbs the hormone binding site following β-sheet expansion in the cCBG structure but shifts away from the binding site by more than 8 A following the binding of cortisol. Unexpectedly, on cortisol binding, the hormone binding site of cCBG adopts a configuration almost identical with that of the native conformer. We conclude that CBG has adapted an allosteric mechanism of the serpins to allow equilibrated release of the hormones by a flip–flop movement of the intact reactive loop into and out of the β-sheet. The change in the hormone binding affinity results from a change in the flexibility or plasticity of the connecting loop, which modulates the configuration of the binding site.

Crystals of glyceraldehyde phosphate dehydrogenase from the glycosome of Trypanosoma brucei brucei have been grown, and a partial data set has been collected using synchrotron radiation. The crystals diffract initially to 2.3 A resolution. The space group is P2(1)2(1)2, with cell dimensions a = 135 A, b = 255 A, c = 115 A, so there are probably at least two tetramers in the asymmetric unit.

The crystal structure of lipoamide dehydrogenase from Azotobacter vinelandii has been determined by a combination of molecular replacement and isomorphous replacement techniques yielding eventually a good-quality 2.8 Å electron density map. Initially, the structure determination was attempted by molecular replacement procedures alone using a model of human glutathione reductase, which has 26% sequence identity with this bacterial dehydrogenase. The rotation function yielded the correct orientation of the model structure both when the glutathione reductase dimer and monomer were used as starting model. The translation function could not be solved, however. Consequently, data for two heavy-atom derivatives were collected using the Hamburg synchotron facilities. The derivatives had several sites in common, which was presumably a major reason why the electron density map obtained by isomorphous information alone was of poor quality. Application of solvent flattening procedures cleaned up the map considerably, however, showing clearly the outline of the lipoamide dehydrogenase dimer, which has a molecular weight of 100,000. Application of the “phased translation function”, which combines the phase information of both isomorphous and molecular replacement, led to an unambiguous determination of the position of the model structure in the lipoamide dehydrogenase unit cell. The non-crystallographic 2-fold axis of the dimer was optimized by several cycles of constrained-restrained least-squares refinement and subsequently used for phase improvement by 2-fold density averaging. After ten cycles at 3.5 Å, the resolution was gradually extended to 2.8 Å in another 140 cycles. The 2.8 Å electron density distribution obtained in this manner was of much improved quality and allowed building of an atomic model of A. vinelandii lipoamide dehydrogenase. It appears that in the orthorhombic crystals used each dimer is involved in contacts with eight surrounding dimers, leaving unexplained why the crystals are rather fragile. Contacts between subunits within one dimer, which are quite extensive, can be divided into two regions separated by a cavity. In one of the contact regions, the level of sequence identity with glutathione reductase is very low but it is quite high in the other. The folding of the polypeptide chain in each subunit is quite similar to that of glutathione reductase, as is the extended conformation of the co-enzyme FAD. The structure of A. vinelandii lipoamide dehydrogenase solved forms a starting point for investigating details of the catalytic mechanism as well as studying the interactions of this enzyme with its partners in the pyruvate dehydrogenase multi-enzyme complex.