Macromolecular ab initio phasing enforcing secondary and tertiary structure.

Ab initio phasing of macromolecular structures, from the native intensities alone with no experimental phase information or previous particular structural knowledge, has been the object of a long quest, limited by two main barriers: structure size and resolution of the data. Current approaches to extend the scope of ab initio phasing include use of the Patterson function, density modification and data extrapolation. The authors' approach relies on the combination of locating model fragments such as polyalanine α-helices with the program PHASER and density modification with the program SHELXE. Given the difficulties in discriminating correct small substructures, many putative groups of fragments have to be tested in parallel; thus calculations are performed in a grid or supercomputer. The method has been named after the Italian painter Arcimboldo, who used to compose portraits out of fruit and vegetables. With ARCIMBOLDO, most collections of fragments remain a 'still-life', but some are correct enough for density modification and main-chain tracing to reveal the protein's true portrait. Beyond α-helices, other fragments can be exploited in an analogous way: libraries of helices with modelled side chains, β-strands, predictable fragments such as DNA-binding folds or fragments selected from distant homologues up to libraries of small local folds that are used to enforce nonspecific tertiary structure; thus restoring the ab initio nature of the method. Using these methods, a number of unknown macromolecules with a few thousand atoms and resolutions around 2 Å have been solved. In the 2014 release, use of the program has been simplified. The software mediates the use of massive computing to automate the grid access required in difficult cases but may also run on a single multicore workstation (http://chango.ibmb.csic.es/ARCIMBOLDO_LITE) to solve straightforward cases

Implications of AlphaFold2 for crystallographic phasing by molecular replacement.

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

Evaluation of model refinement in CASP13.

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.Wellcome Trust Marie Sklowdowska-Curie grain for EU Horizon 202

Evaluation of template-based modeling in CASP13.

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

ARCIMBOLDO on coiled coils.

ARCIMBOLDO solves the phase problem by combining the location of small model fragments using Phaser with density modification and autotracing using SHELXE. Mainly helical structures constitute favourable cases, which can be solved using polyalanine helical fragments as search models. Nevertheless, the solution of coiled-coil structures is often complicated by their anisotropic diffraction and apparent translational noncrystallographic symmetry. Long, straight helices have internal translational symmetry and their alignment in preferential directions gives rise to systematic overlap of Patterson vectors. This situation has to be differentiated from the translational symmetry relating different monomers. ARCIMBOLDO_LITE has been run on single workstations on a test pool of 150 coiled-coil structures with 15-635 amino acids per asymmetric unit and with diffraction data resolutions of between 0.9 and 3.0 Å. The results have been used to identify and address specific issues when solving this class of structures using ARCIMBOLDO. Features from Phaser v.2.7 onwards are essential to correct anisotropy and produce translation solutions that will pass the packing filters. As the resolution becomes worse than 2.3 Å, the helix direction may be reversed in the placed fragments. Differentiation between true solutions and pseudo-solutions, in which helix fragments were correctly positioned but in a reverse orientation, was found to be problematic at resolutions worse than 2.3 Å. Therefore, after every new fragment-placement round, complete or sparse combinations of helices in alternative directions are generated and evaluated. The final solution is once again probed by helix reversal, refinement and extension. To conclude, density modification and SHELXE autotracing incorporating helical constraints is also exploited to extend the resolution limit in the case of coiled coils and to enhance the identification of correct solutions. This study resulted in a specialized mode within ARCIMBOLDO for the solution of coiled-coil structures, which overrides the resolution limit and can be invoked from the command line (keyword coiled_coil) or ARCIMBOLDO_LITE task interface in CCP4i

ARCIMBOLDO on coiled coils.

ARCIMBOLDO solves the phase problem by combining the location of small model fragments using Phaser with density modification and autotracing using SHELXE. Mainly helical structures constitute favourable cases, which can be solved using polyalanine helical fragments as search models. Nevertheless, the solution of coiled-coil structures is often complicated by their anisotropic diffraction and apparent translational noncrystallographic symmetry. Long, straight helices have internal translational symmetry and their alignment in preferential directions gives rise to systematic overlap of Patterson vectors. This situation has to be differentiated from the translational symmetry relating different monomers. ARCIMBOLDO_LITE has been run on single workstations on a test pool of 150 coiled-coil structures with 15-635 amino acids per asymmetric unit and with diffraction data resolutions of between 0.9 and 3.0 Å. The results have been used to identify and address specific issues when solving this class of structures using ARCIMBOLDO. Features from Phaser v.2.7 onwards are essential to correct anisotropy and produce translation solutions that will pass the packing filters. As the resolution becomes worse than 2.3 Å, the helix direction may be reversed in the placed fragments. Differentiation between true solutions and pseudo-solutions, in which helix fragments were correctly positioned but in a reverse orientation, was found to be problematic at resolutions worse than 2.3 Å. Therefore, after every new fragment-placement round, complete or sparse combinations of helices in alternative directions are generated and evaluated. The final solution is once again probed by helix reversal, refinement and extension. To conclude, density modification and SHELXE autotracing incorporating helical constraints is also exploited to extend the resolution limit in the case of coiled coils and to enhance the identification of correct solutions. This study resulted in a specialized mode within ARCIMBOLDO for the solution of coiled-coil structures, which overrides the resolution limit and can be invoked from the command line (keyword coiled_coil) or ARCIMBOLDO_LITE task interface in CCP4i

Structure of a 13-fold superhelix (almost) determined from first principles.

Nuclear hormone receptors are cytoplasm-based transcription factors that bind a ligand, translate to the nucleus and initiate gene transcription in complex with a co-activator such as TIF2 (transcriptional intermediary factor 2). For structural studies the co-activator is usually mimicked by a peptide of circa 13 residues, which for the largest part forms an α-helix when bound to the receptor. The aim was to co-crystallize the glucocorticoid receptor in complex with a ligand and the TIF2 co-activator peptide. The 1.82 Å resolution diffraction data obtained from the crystal could not be phased by molecular replacement using the known receptor structures. HPLC analysis of the crystals revealed the absence of the receptor and indicated that only the co-activator peptide was present. The self-rotation function displayed 13-fold rotational symmetry, which initiated an exhaustive but unsuccessful molecular-replacement approach using motifs of 13-fold symmetry such as α- and β-barrels in various geometries. The structure was ultimately determined by using a single α-helix and the software ARCIMBOLDO, which assembles fragments placed by PHASER before using them as seeds for density modification model building in SHELXE. Systematic variation of the helix length revealed upper and lower size limits for successful structure determination. A beautiful but unanticipated structure was obtained that forms superhelices with left-handed twist throughout the crystal, stabilized by ligand interactions. Together with the increasing diversity of structural elements in the Protein Data Bank the results from TIF2 confirm the potential of fragment-based molecular replacement to significantly accelerate the phasing step for native diffraction data at around 2 Å resolution

ALEPH: a network-oriented approach for the generation of fragment-based libraries and for structure interpretation.

The analysis of large structural databases reveals general features and relationships among proteins, providing useful insight. A different approach is required to characterize ubiquitous secondary-structure elements, where flexibility is essential in order to capture small local differences. The ALEPH software is optimized for the analysis and the extraction of small protein folds by relying on their geometry rather than on their sequence. The annotation of the structural variability of a given fold provides valuable information for fragment-based molecular-replacement methods, in which testing alternative model hypotheses can succeed in solving difficult structures when no homology models are available or are successful. ARCIMBOLDO_BORGES combines the use of composite secondary-structure elements as a search model with density modification and tracing to reveal the rest of the structure when both steps are successful. This phasing method relies on general fold libraries describing variations around a given pattern of β-sheets and helices extracted using ALEPH. The program introduces characteristic vectors defined from the main-chain atoms as a way to describe the geometrical properties of the structure. ALEPH encodes structural properties in a graph network, the exploration of which allows secondary-structure annotation, decomposition of a structure into small compact folds, generation of libraries of models representing a variation of a given fold and finally superposition of these folds onto a target structure. These functions are available through a graphical interface designed to interactively show the results of structure manipulation, annotation, fold decomposition, clustering and library generation. ALEPH can produce pictures of the graphs, structures and folds for publication purposes

Factors influencing estimates of coordinate error for molecular replacement.

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