Search

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 interpretation of cryo-EM maps often includes the docking of known or predicted structures of the components, which is particularly useful when the map resolution is worse than 4 Å. Although it can be effective to search the entire map to find the best placement of a component, the process can be slow when the maps are large. However, frequently there is a well-founded hypothesis about where particular components are located. In such cases, a local search using a map subvolume will be much faster because the search volume is smaller, and more sensitive because optimizing the search volume for the rotation-search step enhances the signal to noise. A Fourier-space likelihood-based local search approach, based on the previously published em_placement software, has been implemented in the new emplace_local program. Tests confirm that the local search approach enhances the speed and sensitivity of the computations. An interactive graphical interface in the ChimeraX molecular-graphics program provides a convenient way to set up and evaluate docking calculations, particularly in defining the part of the map into which the components should be placed., (open access.)

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein-nucleic acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: Escherichia coli beta-galactosidase with inhibitor, SARS-CoV-2 virus RNA-dependent RNA polymerase with covalently bound nucleotide analog and SARS-CoV-2 virus ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. The quality of submitted ligand models and surrounding atoms were analyzed by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics and contact scores. A composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution., (© 2024. The Author(s), under exclusive licence to Springer Nature America, Inc.)

A defining pathological feature of most neurodegenerative diseases is the assembly of proteins into amyloid that form disease-specific structures 1 . In Alzheimer's disease, this is characterized by the deposition of β-amyloid and tau with disease-specific conformations. The in situ structure of amyloid in the human brain is unknown. Here, using cryo-fluorescence microscopy-targeted cryo-sectioning, cryo-focused ion beam-scanning electron microscopy lift-out and cryo-electron tomography, we determined in-tissue architectures of β-amyloid and tau pathology in a postmortem Alzheimer's disease donor brain. β-amyloid plaques contained a mixture of fibrils, some of which were branched, and protofilaments, arranged in parallel arrays and lattice-like structures. Extracellular vesicles and cuboidal particles defined the non-amyloid constituents of β-amyloid plaques. By contrast, tau inclusions formed parallel clusters of unbranched filaments. Subtomogram averaging a cluster of 136 tau filaments in a single tomogram revealed the polypeptide backbone conformation and filament polarity orientation of paired helical filaments within tissue. Filaments within most clusters were similar to each other, but were different between clusters, showing amyloid heterogeneity that is spatially organized by subcellular location. The in situ structural approaches outlined here for human donor tissues have applications to a broad range of neurodegenerative diseases., (© 2024. The Author(s).)

Five new Co-editors are appointed to the Editorial Board of Acta Cryst. D - Structural Biology.

Artificial intelligence-based protein structure prediction methods such as AlphaFold have revolutionized structural biology. The accuracies of these predictions vary, however, and they do not take into account ligands, covalent modifications or other environmental factors. Here, we evaluate how well AlphaFold predictions can be expected to describe the structure of a protein by comparing predictions directly with experimental crystallographic maps. In many cases, AlphaFold predictions matched experimental maps remarkably closely. In other cases, even very high-confidence predictions differed from experimental maps on a global scale through distortion and domain orientation, and on a local scale in backbone and side-chain conformation. We suggest considering AlphaFold predictions as exceptionally useful hypotheses. We further suggest that it is important to consider the confidence in prediction when interpreting AlphaFold predictions and to carry out experimental structure determination to verify structural details, particularly those that involve interactions not included in the prediction., (© 2023. The Author(s).)

Michael James is remembered., (open access.)

Haloalkane dehalogenase (HLD) enzymes employ an S N 2 nucleophilic substitution mechanism to erase halogen substituents in diverse organohalogen compounds. Subfamily I and II HLDs are well-characterized enzymes, but the mode and purpose of multimerization of subfamily III HLDs are unknown. Here we probe the structural organization of DhmeA, a subfamily III HLD-like enzyme from the archaeon Haloferax mediterranei, by combining cryo-electron microscopy (cryo-EM) and x-ray crystallography. We show that full-length wild-type DhmeA forms diverse quaternary structures, ranging from small oligomers to large supramolecular ring-like assemblies of various sizes and symmetries. We optimized sample preparation steps, enabling three-dimensional reconstructions of an oligomeric species by single-particle cryo-EM. Moreover, we engineered a crystallizable mutant (DhmeA ΔGG ) that provided diffraction-quality crystals. The 3.3 Å crystal structure reveals that DhmeA ΔGG forms a ring-like 20-mer structure with outer and inner diameter of ~200 and ~80 Å, respectively. An enzyme homodimer represents a basic repeating building unit of the crystallographic ring. Three assembly interfaces (dimerization, tetramerization, and multimerization) were identified to form the supramolecular ring that displays a negatively charged exterior, while its interior part harboring catalytic sites is positively charged. Localization and exposure of catalytic machineries suggest a possible processing of large negatively charged macromolecular substrates., (© 2023 The Authors. Protein Science published by Wiley Periodicals LLC on behalf of The Protein Society.)

This editorial acknowledges the transformative impact of new machine-learning methods, such as the use of AlphaFold, but also makes the case for the continuing need for experimental structural biology.

The Collaborative Computational Project No. 4 (CCP4) is a UK-led international collective with a mission to develop, test, distribute and promote software for macromolecular crystallography. The CCP4 suite is a multiplatform collection of programs brought together by familiar execution routines, a set of common libraries and graphical interfaces. The CCP4 suite has experienced several considerable changes since its last reference article, involving new infrastructure, original programs and graphical interfaces. This article, which is intended as a general literature citation for the use of the CCP4 software suite in structure determination, will guide the reader through such transformations, offering a general overview of the new features and outlining future developments. As such, it aims to highlight the individual programs that comprise the suite and to provide the latest references to them for perusal by crystallographers around the world., (open access.)

Optimized docking of models into cryo-EM maps requires exploiting an understanding of the signal expected in the data to minimize the calculation time while maintaining sufficient signal. The likelihood-based rotation function used in crystallography can be employed to establish plausible orientations in a docking search. A phased likelihood translation function yields scores for the placement and rigid-body refinement of oriented models. Optimized strategies for choices of the resolution of data from the cryo-EM maps to use in the calculations and the size of search volumes are based on expected log-likelihood-gain scores computed in advance of the search calculation. Tests demonstrate that the new procedure is fast, robust and effective at placing models into even challenging cryo-EM maps., (open access.)

Fast, reliable docking of models into cryo-EM maps requires understanding of the errors in the maps and the models. Likelihood-based approaches to errors have proven to be powerful and adaptable in experimental structural biology, finding applications in both crystallography and cryo-EM. Indeed, previous crystallographic work on the errors in structural models is directly applicable to likelihood targets in cryo-EM. Likelihood targets in Fourier space are derived here to characterize, based on the comparison of half-maps, the direction- and resolution-dependent variation in the strength of both signal and noise in the data. Because the signal depends on local features, the signal and noise are analysed in local regions of the cryo-EM reconstruction. The likelihood analysis extends to prediction of the signal that will be achieved in any docking calculation for a model of specified quality and completeness. A related calculation generalizes a previous measure of the information gained by making the cryo-EM reconstruction., (open access.)

Cas12a is a programmable nuclease for adaptive immunity against invading nucleic acids in CRISPR-Cas systems. Here, we report the crystal structures of apo Cas12a from Lachnospiraceae bacterium MA2020 (Lb2) and the Lb2Cas12a+crRNA complex, as well as the cryo-EM structure and functional studies of the Lb2Cas12a+crRNA+DNA complex. We demonstrate that apo Lb2Cas12a assumes a unique, elongated conformation, whereas the Lb2Cas12a+crRNA binary complex exhibits a compact conformation that subsequently rearranges to a semi-open conformation in the Lb2Cas12a+crRNA+DNA ternary complex. Notably, in solution, apo Lb2Cas12a is dynamic and can exist in both elongated and compact forms. Residues from Met493 to Leu523 of the WED domain undergo major conformational changes to facilitate the required structural rearrangements. The REC lobe of Lb2Cas12a rotates 103° concomitant with rearrangement of the hinge region close to the WED and RuvC II domains to position the RNA-DNA duplex near the catalytic site. Our findings provide insight into crRNA recognition and the mechanism of target DNA cleavage., Competing Interests: The authors have declared that no competing interests exist., (Copyright: © 2023 Jianwei et al. This is an open access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.)

Experimental structure determination can be accelerated with artificial intelligence (AI)-based structure-prediction methods such as AlphaFold. Here, an automatic procedure requiring only sequence information and crystallographic data is presented that uses AlphaFold predictions to produce an electron-density map and a structural model. Iterating through cycles of structure prediction is a key element of this procedure: a predicted model rebuilt in one cycle is used as a template for prediction in the next cycle. This procedure was applied to X-ray data for 215 structures released by the Protein Data Bank in a recent six-month period. In 87% of cases our procedure yielded a model with at least 50% of C α atoms matching those in the deposited models within 2 Å. Predictions from the iterative template-guided prediction procedure were more accurate than those obtained without templates. It is concluded that AlphaFold predictions obtained based on sequence information alone are usually accurate enough to solve the crystallographic phase problem with molecular replacement, and a general strategy for macromolecular structure determination that includes AI-based prediction both as a starting point and as a method of model optimization is suggested., (open access.)

Machine-learning prediction algorithms such as AlphaFold and RoseTTAFold can create remarkably accurate protein models, but these models usually have some regions that are predicted with low confidence or poor accuracy. We hypothesized that by implicitly including new experimental information such as a density map, a greater portion of a model could be predicted accurately, and that this might synergistically improve parts of the model that were not fully addressed by either machine learning or experiment alone. An iterative procedure was developed in which AlphaFold models are automatically rebuilt on the basis of experimental density maps and the rebuilt models are used as templates in new AlphaFold predictions. We show that including experimental information improves prediction beyond the improvement obtained with simple rebuilding guided by the experimental data. This procedure for AlphaFold modeling with density has been incorporated into an automated procedure for interpretation of crystallographic and electron cryo-microscopy maps., (© 2022. The Author(s).)

AlphaFold has recently become an important tool in providing models for experimental structure determination by X-ray crystallography and cryo-EM. Large parts of the predicted models typically approach the accuracy of experimentally determined structures, although there are frequently local errors and errors in the relative orientations of domains. Importantly, residues in the model of a protein predicted by AlphaFold are tagged with a predicted local distance difference test score, informing users about which regions of the structure are predicted with less confidence. AlphaFold also produces a predicted aligned error matrix indicating its confidence in the relative positions of each pair of residues in the predicted model. The phenix.process_predicted_model tool downweights or removes low-confidence residues and can break a model into confidently predicted domains in preparation for molecular replacement or cryo-EM docking. These confidence metrics are further used in ISOLDE to weight torsion and atom-atom distance restraints, allowing the complete AlphaFold model to be interactively rearranged to match the docked fragments and reducing the need for the rebuilding of connecting regions., (open access.)

Heterozygous mutations in BMPR2 (bone morphogenetic protein (BMP) receptor type II) cause pulmonary arterial hypertension. BMPRII is a receptor for over 15 BMP ligands, but why BMPR2 mutations cause lung-specific pathology is unknown. To elucidate the molecular basis of BMP:BMPRII interactions, we report crystal structures of binary and ternary BMPRII receptor complexes with BMP10, which contain an ensemble of seven different BMP10:BMPRII 1:1 complexes. BMPRII binds BMP10 at the knuckle epitope, with the A-loop and β4 strand making BMPRII-specific interactions. The BMPRII binding surface on BMP10 is dynamic, and the affinity is weaker in the ternary complex than in the binary complex. Hydrophobic core and A-loop interactions are important in BMPRII-mediated signalling. Our data reveal how BMPRII is a low affinity receptor, implying that forming a signalling complex requires high concentrations of BMPRII, hence mutations will impact on tissues with highest BMPR2 expression such as the lung vasculature., (© 2022. The Author(s).)

The process of platelet production has so far been understood to be a 2-stage process: megakaryocyte maturation from hematopoietic stem cells followed by proplatelet formation, with each phase regulating the peripheral blood platelet count. Proplatelet formation releases into the bloodstream beads-on-a-string preplatelets, which undergo fission into mature platelets. For the first time, we show that preplatelet maturation is a third, tightly regulated, critical process akin to cytokinesis that regulates platelet count. We show that deficiency in cytokine receptor-like factor 3 (CRLF3) in mice leads to an isolated and sustained 25% to 48% reduction in the platelet count without any effect on other blood cell lineages. We show that Crlf3-/- preplatelets have increased microtubule stability, possibly because of increased microtubule glutamylation via the interaction of CRLF3 with key members of the Hippo pathway. Using a mouse model of JAK2 V617F essential thrombocythemia, we show that a lack of CRLF3 leads to long-term lineage-specific normalization of the platelet count. We thereby postulate that targeting CRLF3 has therapeutic potential for treatment of thrombocythemia., (© 2022 by The American Society of Hematology. Licensed under Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International (CC BY-NC-ND 4.0), permitting only noncommercial, nonderivative use with attribution. All other rights reserved.)

The AlphaFold2 results in the 14th edition of Critical Assessment of Structure Prediction (CASP14) showed that accurate (low root-mean-square deviation) in silico models of protein structure domains are on the horizon, whether or not the protein is related to known structures through high-coverage sequence similarity. As highly accurate models become available, generated by harnessing the power of correlated mutations and deep learning, one of the aspects of structural biology to be impacted will be methods of phasing in crystallography. Here, the data from CASP14 are used to explore the prospects for changes in phasing methods, and in particular to explore the prospects for molecular-replacement phasing using in silico models., (open access.)

The assessment of CASP models for utility in molecular replacement is a measure of their use in a valuable real-world application. In CASP7, the metric for molecular replacement assessment involved full likelihood-based molecular replacement searches; however, this restricted the assessable targets to crystal structures with only one copy of the target in the asymmetric unit, and to those where the search found the correct pose. In CASP10, full molecular replacement searches were replaced by likelihood-based rigid-body refinement of models superimposed on the target using the LGA algorithm, with the metric being the refined log-likelihood-gain (LLG) score. This enabled multi-copy targets and very poor models to be evaluated, but a significant further issue remained: the requirement of diffraction data for assessment. We introduce here the relative-expected-LLG (reLLG), which is independent of diffraction data. This reLLG is also independent of any crystal form, and can be calculated regardless of the source of the target, be it X-ray, NMR or cryo-EM. We calibrate the reLLG against the LLG for targets in CASP14, showing that it is a robust measure of both model and group ranking. Like the LLG, the reLLG shows that accurate coordinate error estimates add substantial value to predicted models. We find that refinement by CASP groups can often convert an inadequate initial model into a successful MR search model. Consistent with findings from others, we show that the AlphaFold2 models are sufficiently good, and reliably so, to surpass other current model generation strategies for attempting molecular replacement phasing., (© 2021 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals LLC.)

DeepMind presented notably accurate predictions at the recent 14th Critical Assessment of Structure Prediction (CASP14) conference. We explored network architectures that incorporate related ideas and obtained the best performance with a three-track network in which information at the one-dimensional (1D) sequence level, the 2D distance map level, and the 3D coordinate level is successively transformed and integrated. The three-track network produces structure predictions with accuracies approaching those of DeepMind in CASP14, enables the rapid solution of challenging x-ray crystallography and cryo-electron microscopy structure modeling problems, and provides insights into the functions of proteins of currently unknown structure. The network also enables rapid generation of accurate protein-protein complex models from sequence information alone, short-circuiting traditional approaches that require modeling of individual subunits followed by docking. We make the method available to the scientific community to speed biological research., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works.)

SAD phasing can be challenging when the signal-to-noise ratio is low. In such cases, having an accurate estimate of the substructure content can determine whether or not the substructure of anomalous scatterer positions can successfully be determined. Here, a likelihood-based target function is proposed to accurately estimate the strength of the anomalous scattering contribution directly from the measured intensities, determining a complex correlation parameter relating the Bijvoet mates as a function of resolution. This gives a novel measure of the intrinsic anomalous signal. The SAD likelihood target function also accounts for correlated errors in the measurement of intensities from Bijvoet mates, which can arise from the effects of radiation damage. When the anomalous signal is assumed to come primarily from a substructure comprising one anomalous scatterer with a known value of f'' and when the protein composition of the crystal is estimated correctly, the refined complex correlation parameters can be interpreted in terms of the atomic content of the primary anomalous scatterer before the substructure is known. The maximum-likelihood estimation of substructure content was tested on a curated database of 357 SAD cases with useful anomalous signal. The prior estimates of substructure content are highly correlated to the content determined by phasing calculations, with a correlation coefficient (on a log-log basis) of 0.72., (open access.)

When building atomic models into weak and/or low-resolution density, a common strategy is to restrain their conformation to that of a higher resolution model of the same or similar sequence. When doing so, it is important to avoid over-restraining to the reference model in the face of disagreement with the experimental data. The most common strategy for this is the use of `top-out' potentials. These act like simple harmonic restraints within a defined range, but gradually weaken when the deviation between the model and reference grows beyond that range. In each current implementation the rate at which the potential flattens at large deviations follows a fixed form, although the form chosen varies among implementations. A restraint potential with a tuneable rate of flattening would provide greater flexibility to encode the confidence in any given restraint. Here, two new such potentials are described: a Cartesian distance restraint derived from a recent generalization of common loss functions and a periodic torsion restraint based on a renormalization of the von Mises distribution. Further, their implementation as user-adjustable/switchable restraints in ISOLDE is described and their use in some real-world examples is demonstrated., (open access.)

The angiotensin peptides that control blood pressure are released from the non-inhibitory plasma serpin, angiotensinogen, on cleavage of its extended N-terminal tail by the specific aspartyl-protease, renin. Angiotensinogen had previously been assumed to be a passive substrate, but we describe here how recent studies reveal an inherent conformational mechanism that is critical to the cleavage and release of the angiotensin peptides and consequently to the control of blood pressure. A series of crystallographic structures of angiotensinogen and its derivative forms, together with its complexes with renin show in molecular detail how the interaction with renin triggers a profound shift of the amino-terminal tail of angiotensinogen with modulation occurring at several levels. The tail of angiotensinogen is restrained by a labile disulfide bond, with changes in its redox status affecting angiotensin release, as demonstrably so in the hypertensive complication of pregnancy, pre-eclampsia. The shift of the tail also enhances the binding of renin through a tail-in-mouth allosteric mechanism. The N-terminus is now seen to insert into a pocket equivalent to the hormone-binding site on other serpins, with helix H of angiotensinogen unwinding to form key interactions with renin. The findings explain the precise species specificity of the interaction with renin and with variant carbohydrate linkages. Overall, the studies provide new insights into the physiological regulation of angiotensin release, with an ability to respond to local tissue and temperature changes, and with the opening of strategies for the development of novel agents for the treatment of hypertension., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2021 Shu, Wan, Read, Carrell and Zhou.)

Detection of translational noncrystallographic symmetry (TNCS) can be critical for success in crystallographic phasing, particularly when molecular-replacement models are poor or anomalous phasing information is weak. If the correct TNCS is detected then expected intensity factors for each reflection can be refined, so that the maximum-likelihood functions underlying molecular replacement and single-wavelength anomalous dispersion use appropriate structure-factor normalization and variance terms. Here, an analysis of a curated database of protein structures from the Protein Data Bank to investigate how TNCS manifests in the Patterson function is described. These studies informed an algorithm for the detection of TNCS, which includes a method for detecting the number of vectors involved in any commensurate modulation (the TNCS order). The algorithm generates a ranked list of possible TNCS associations in the asymmetric unit for exploration during structure solution., (open access.)

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., (open access.)

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., (open access.)

A density-modification procedure for improving maps from single-particle electron cryogenic microscopy (cryo-EM) is presented. The theoretical basis of the method is identical to that of maximum-likelihood density modification, previously used to improve maps from macromolecular X-ray crystallography. Key differences from applications in crystallography are that the errors in Fourier coefficients are largely in the phases in crystallography but in both phases and amplitudes in cryo-EM, and that half-maps with independent errors are available in cryo-EM. These differences lead to a distinct approach for combination of information from starting maps with information obtained in the density-modification process. The density-modification procedure was applied to a set of 104 datasets and improved map-model correlation and increased the visibility of details in many of the maps. The procedure requires two unmasked half-maps and a sequence file or other source of information on the volume of the macromolecule that has been imaged.

The information gained by making a measurement, termed the Kullback-Leibler divergence, assesses how much more precisely the true quantity is known after the measurement was made (the posterior probability distribution) than before (the prior probability distribution). It provides an upper bound for the contribution that an observation can make to the total likelihood score in likelihood-based crystallographic algorithms. This makes information gain a natural criterion for deciding which data can legitimately be omitted from likelihood calculations. Many existing methods use an approximation for the effects of measurement error that breaks down for very weak and poorly measured data. For such methods a different (higher) information threshold is appropriate compared with methods that account well for even large measurement errors. Concerns are raised about a current trend to deposit data that have been corrected for anisotropy, sharpened and pruned without including the original unaltered measurements. If not checked, this trend will have serious consequences for the reuse of deposited data by those who hope to repeat calculations using improved new methods., (open access.)

Good prior estimates of the effective root-mean-square deviation (r.m.s.d.) between the atomic coordinates of the model and the target optimize the signal in molecular replacement, thereby increasing the success rate in difficult cases. Previous studies using protein structures solved by X-ray crystallography as models showed that optimal error estimates (refined after structure solution) were correlated with the sequence identity between the model and target, and with the number of residues in the model. Here, this work has been extended to find additional correlations between parameters of the model and the target and hence improved prior estimates of the coordinate error. Using a graph database, a curated set of 6030 molecular-replacement calculations using models that had been solved by X-ray crystallography was analysed to consider about 120 model and target parameters. Improved estimates were achieved by replacing the sequence identity with the Gonnet score for sequence similarity, as well as by considering the resolution of the target structure and the MolProbity score of the model. This approach was extended by analysing 12 610 additional molecular-replacement calculations where the model was determined by NMR. The median r.m.s.d. between pairs of models in an ensemble was found to be correlated with the estimated r.m.s.d. to the target. For models solved by NMR, the overall coordinate error estimates were larger than for structures determined by X-ray crystallography, and were more highly correlated with the number of residues., (open access.)

To identify novel causes of hereditary thrombocytopenia, we performed a genetic association analysis of whole-genome sequencing data from 13 037 individuals enrolled in the National Institute for Health Research (NIHR) BioResource, including 233 cases with isolated thrombocytopenia. We found an association between rare variants in the transcription factor-encoding gene IKZF5 and thrombocytopenia. We report 5 causal missense variants in or near IKZF5 zinc fingers, of which 2 occurred de novo and 3 co-segregated in 3 pedigrees. A canonical DNA-zinc finger binding model predicts that 3 of the variants alter DNA recognition. Expression studies showed that chromatin binding was disrupted in mutant compared with wild-type IKZF5, and electron microscopy revealed a reduced quantity of α granules in normally sized platelets. Proplatelet formation was reduced in megakaryocytes from 7 cases relative to 6 controls. Comparison of RNA-sequencing data from platelets, monocytes, neutrophils, and CD4+ T cells from 3 cases and 14 healthy controls showed 1194 differentially expressed genes in platelets but only 4 differentially expressed genes in each of the other blood cell types. In conclusion, IKZF5 is a novel transcriptional regulator of megakaryopoiesis and the eighth transcription factor associated with dominant thrombocytopenia in humans., (© 2019 by The American Society of Hematology.)

Performance in the template-based modeling (TBM) category of CASP13 is assessed here, using a variety of metrics. Performance of the predictor groups that participated is ranked using the primary ranking score that was developed by the assessors for CASP12. This reveals that the best results are obtained by groups that include contact predictions or inter-residue distance predictions derived from deep multiple sequence alignments. In cases where there is a good homolog in the wwPDB (TBM-easy category), the best results are obtained by modifying a template. However, for cases with poorer homologs (TBM-hard), very good results can be obtained without using an explicit template, by deep learning algorithms trained on the wwPDB. Alternative metrics are introduced, to allow testing of aspects of structural models that are not addressed by traditional CASP metrics. These include comparisons to the main-chain and side-chain torsion angles of the target, and the utility of models for solving crystal structures by the molecular replacement method. The alternative metrics are poorly correlated with the traditional metrics, and it is proposed that modeling has reached a sufficient level of maturity that the best models should be expected to satisfy this wider range of criteria., (© 2019 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals, Inc.)

Performance in the model refinement category of the 13th round of Critical Assessment of Structure Prediction (CASP13) is assessed, showing that some groups consistently improve most starting models whereas the majority of participants continue to degrade the starting model on average. Using the ranking formula developed for CASP12, it is shown that only 7 of 32 groups perform better than a "naïve predictor" who just submits the starting model. Common features in their approaches include a dependence on physics-based force fields to judge alternative conformations and the use of molecular dynamics to relax models to local minima, usually with some restraints to prevent excessively large movements. In addition to the traditional CASP metrics that focus largely on the quality of the overall fold, alternative metrics are evaluated, including comparisons of the main-chain and side-chain torsion angles, and the utility of the models for solving crystal structures by the molecular replacement method. It is proposed that the introduction of these metrics, as well as consideration of the accuracy of coordinate error estimates, would improve the discrimination between good and very good models., (© 2019 The Authors. Proteins: Structure, Function, and Bioinformatics published by Wiley Periodicals, Inc.)

AMPylation is an inactivating modification that alters the activity of the major endoplasmic reticulum (ER) chaperone BiP to match 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 of residues along an inhibitory pathway linking 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. Moreover, 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., (© 2019 The Authors. Published under the terms of the CC BY 4.0 license.)

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., (open access.)

Data pathologies caused by effects such as diffraction anisotropy and translational noncrystallographic symmetry (tNCS) can dramatically complicate the solution of the crystal structures of macromolecules. Such problems were encountered in determining the structure of a mutant form of Rab27a, a member of the Rab GTPases. Mutant Rab27a constructs that crystallize in the free form were designed for use in the discovery of drugs to reduce primary tumour invasiveness and metastasis. One construct, hRab27a Mut , crystallized within 24 h and diffracted to 2.82 Å resolution, with a unit cell possessing room for a large number of protein copies. Initial efforts to solve the structure using molecular replacement by Phaser were not successful. Analysis of the data set revealed that the crystals suffered from both extreme anisotropy and strong tNCS. As a result, large numbers of reflections had estimated standard deviations that were much larger than their measured intensities and their expected intensities, revealing problems with the use of such data at the time in Phaser. By eliminating extremely weak reflections with the largest combined effects of anisotropy and tNCS, these problems could be avoided, allowing a molecular-replacement solution to be found. The lessons that were learned in solving this structure have guided improvements in the numerical analysis used in Phaser, particularly in identifying diffraction measurements that convey very little information content. The calculation of information content could also be applied as an alternative to ellipsoidal truncation. The post-mortem analysis also revealed an oversight in accounting for measurement errors in the fast rotation function. While the crystal of mutant Rab27a is not amenable to drug screening, the structure can guide new modifications to obtain more suitable crystal forms., (open access.)

The renin-angiotensin cascade is a hormone system that regulates blood pressure and fluid balance. Renin-mediated cleavage of the angiotensin I peptide from the N terminus of angiotensinogen (AGT) is the rate-limiting step of this cascade; however, the detailed molecular mechanism underlying this step is unclear. Here, we solved the crystal structures of glycosylated human AGT (2.30 Å resolution), its encounter complex with renin (2.55 Å), AGT cleaved in its reactive center loop (RCL; 2.97 Å), and spent AGT from which the N-terminal angiotensin peptide was removed (2.63 Å). These structures revealed that AGT undergoes profound conformational changes and binds renin through a tail-into-mouth allosteric mechanism that inserts the N terminus into a pocket equivalent to a hormone-binding site on other serpins. These changes fully extended the N-terminal tail, with the scissile bond for angiotensin release docked in renin's active site. Insertion of the N terminus into this pocket accompanied a complete unwinding of helix H of AGT, which, in turn, formed key interactions with renin in the complementary binding interface. Mutagenesis and kinetic analyses confirmed that renin-mediated production of angiotensin I is controlled by interactions of amino acid residues and glycan components outside renin's active-site cleft. Our findings indicate that AGT adapts unique serpin features for hormone delivery and binds renin through concerted movements in the N-terminal tail and in its main body to modulate angiotensin release. These insights provide a structural basis for the development of agents that attenuate angiotensin release by targeting AGT's hormone binding pocket., (© 2019 Yan et al.)

Angiotensinogen (AGT) is a critical protein in the renin-angiotensin-aldosterone system and may have an important role in the pathogenesis of pre-eclampsia. The disulphide linkage between cysteines 18 and 138 has a key role in the redox switch of AGT which modulates the release of angiotensin I with consequential effects on blood pressure. In this paper, we report a quantitative targeted LC-MS/MS method for the reliable measurement of the total AGT and its reduced and oxidised forms in human plasma. AGT was selectively enriched from human plasma using two-dimensional chromatography employing concanavalin A lectin affinity and reversed phase steps and then deglycosylated using PNGase F. A differential alkylation approach was coupled with targeted LC-MS/MS method to identify the two AGT forms in the plasma chymotryptic digest. An additional AGT proteolytic marker peptide was identified and used to measure total AGT levels. The developed MS workflow enabled the reproducible detection of total AGT and its two distinct forms in human plasma with analytical precision of ≤ 15%. The LC-MS/MS assay for total AGT in plasma showed a linear response (R 2  = 0.992) with a limit of quantification in the low nanomolar range. The method gave suitable validation characteristics for biomedical application to the quantification of the oxidation level and the total level of AGT in plasma samples collected from normal and pre-eclamptic patients.

The ability of bacteria to infect a host relies in part on the secretion of molecular virulence factors across the cell envelope. Pseudomonas aeruginosa, a ubiquitous environmental bacterium causing opportunistic infections in humans, employs the type II secretion system (T2SS) to transport effector proteins across its cellular envelope as part of a diverse array of virulence strategies. General secretory pathway protein L (GspL) is an essential inner-membrane component of the T2SS apparatus, and is thought to facilitate transduction of the energy from ATP hydrolysis in the cytoplasm to the periplasmic components of the system. However, our incomplete understanding of the assembly principles of the T2SS machinery prevents the mechanistic deconvolution of T2SS-mediated protein secretion. Here we show via two crystal structures that the periplasmic ferredoxin-like domain of GspL (GspL fld ) is a dimer stabilized by hydrophobic interactions, and that this interface may allow significant interdomain plasticity. The general dimerization mode of GspL fld is shared with GspL from Vibrio parahaemolyticus suggesting a conserved oligomerization mode across the GspL family. Furthermore, we identified a tetrameric form of the complete periplasmic segment of GspL (GspL peri ) which indicates that GspL may be able to adopt multiple oligomeric states as part of its dynamic role in the T2SS apparatus.

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., (open access.)

This data article makes available the informed computation of the whole Protein Data Bank (PDB) to investigate diffraction anisotropy on a large scale and to perform statistics. This data has been investigated in detail in "X-ray diffraction reveals the intrinsic difference in the physical properties of membrane and soluble proteins" [1]. Diffraction anisotropy is traditionally associated with absence of contacts in-between macromolecules within the crystals in a given direction of space. There are however many case that do not follow this empirical rule. To investigate and sort out this discrepancy, we computed diffraction anisotropy for every entry of the PDB, and put them in context of relevant metrics to compare X-ray diffraction in reciprocal space to the crystal packing in real space. These metrics were either extracted from PDB files when available (resolution, space groups, cell parameters, solvent content), or calculated using standard procedures (anisotropy, crystal contacts, presence of ligands). More specifically, we separated entries to compare soluble vs membrane proteins, and further separated the later in subcategories according to their insertion in the membrane, function, or type of crystallization (Type I vs Type II crystal packing). This informed database is being made available to investigators in the raw and curated formats that can be re-used for further downstream studies. This dataset is useful to test ideas and to ascertain hypothesis based on statistical analysis.