Search Models For Molecular Replacement Research Articles

High-accuracy protein structure prediction in CASP14.

The application of state-of-the-art deep-learning approaches to the protein modeling problem has expanded the "high-accuracy" category in CASP14 to encompass all targets. Building on the metrics used for high-accuracy assessment in previous CASPs, we evaluated the performance of all groups that submitted models for at least 10 targets across all difficulty classes, and judged the usefulness of those produced by AlphaFold2 (AF2) as molecular replacement search models with AMPLE. Driven by the qualitative diversity of the targets submitted to CASP, we also introduce DipDiff as a new measure for the improvement in backbone geometry provided by a model versus available templates. Although a large leap in high-accuracy is seen due to AF2, the second-best method in CASP14 out-performed the best in CASP13, illustrating the role of community-based benchmarking in the development and evolution of the protein structure prediction field.

Fragment-based determination of a proteinase K structure from MicroED data using ARCIMBOLDO_SHREDDER

Structure determination of novel biological macromolecules by X-ray crystallography can be facilitated by the use of small structural fragments, some of only a few residues in length, as effective search models for molecular replacement to overcome the phase problem. Independence from the need for a complete pre-existing model with sequence similarity to the crystallized molecule is the primary appeal of ARCIMBOLDO, a suite of programs which employs this ab initio algorithm for phase determination. Here, the use of ARCIMBOLDO is investigated to overcome the phase problem with the electron cryomicroscopy (cryoEM) method known as microcrystal electron diffraction (MicroED). The results support the use of the ARCIMBOLDO_SHREDDER pipeline to provide phasing solutions for a structure of proteinase K from 1.6 Å resolution data using model fragments derived from the structures of proteins sharing a sequence identity of as low as 20%. ARCIMBOLDO_SHREDDER identified the most accurate polyalanine fragments from a set of distantly related sequence homologues. Alternatively, such templates were extracted in spherical volumes and given internal degrees of freedom to refine towards the target structure. Both modes relied on the rotation function in Phaser to identify or refine fragment models and its translation function to place them. Model completion from the placed fragments proceeded through phase combination of partial solutions and/or density modification and main-chain autotracing using SHELXE. The combined set of fragments was sufficient to arrive at a solution that resembled that determined by conventional molecular replacement using the known target structure as a search model. This approach obviates the need for a single, complete and highly accurate search model when phasing MicroED data, and permits the evaluation of large fragment libraries for this purpose.

Utilizing guanine-coordinated Zn2+ ions to determine DNA crystal structures by single-wavelength anomalous diffraction.

The experimental phase determination of crystal structures of nucleic acids and nucleic acid-ligand complexes would benefit from a facile method. Even for double-stranded DNA, software-generated models are generally insufficiently accurate to serve as molecular replacement search models, necessitating experimental phasing. Here, it is demonstrated that Zn2+ ions coordinated to the N7 atom of guanine bases generate sufficient anomalous signal for single-wavelength anomalous diffraction (SAD) phasing of DNA crystal structures. Using zinc SAD, three crystal structures of double-stranded DNA oligomers, 5'-AGGGATCCCT-3', 5'-GGGATCCC-3' and 5'-GAGGCCTC-3', were determined. By determining the crystal structure of one of these oligomers, GAGGCCTC, in the presence of Mg2+ instead of Zn2+, it was demonstrated that Zn2+ is not structurally perturbing. These structures allowed the analysis of structural changes in the DNA on the binding of analogues of the natural product mithramycin to two of these oligomers, AGGGATCCCT and GAGGCCTC. Zinc SAD may become a routine approach for determining the crystal structures of nucleic acids and their complexes with small molecules.

Volume-ensemble search models for molecular replacement

Volume-ensemble search models for molecular replacement

Recent developments in MrBUMP: better search-model preparation, graphical interaction with search models, and solution improvement and assessment.

Increasing sophistication in molecular-replacement (MR) software and the rapid expansion of the PDB in recent years have allowed the technique to become the dominant method for determining the phases of a target structure in macromolecular X-ray crystallography. In addition, improvements in bioinformatic techniques for finding suitable homologous structures for use as MR search models, combined with developments in refinement and model-building techniques, have pushed the applicability of MR to lower sequence identities and made weak MR solutions more amenable to refinement and improvement. MrBUMP is a CCP4 pipeline which automates all stages of the MR procedure. Its scope covers everything from the sourcing and preparation of suitable search models right through to rebuilding of the positioned search model. Recent improvements to the pipeline include the adoption of more sensitive bioinformatic tools for sourcing search models, enhanced model-preparation techniques including better ensembling of homologues, and the use of phase improvement and model building on the resulting solution. The pipeline has also been deployed as an online service through CCP4 online, which allows its users to exploit large bioinformatic databases and coarse-grained parallelism to speed up the determination of a possible solution. Finally, the molecular-graphics application CCP4mg has been combined with MrBUMP to provide an interactive visual aid to the user during the process of selecting and manipulating search models for use in MR. Here, these developments in MrBUMP are described with a case study to explore how some of the enhancements to the pipeline and to CCP4mg can help to solve a difficult case.

Re-refinement of the spliceosomal U4 snRNP core-domain structure.

The core domain of small nuclear ribonucleoprotein (snRNP), comprised of a ring of seven paralogous proteins bound around a single-stranded RNA sequence, functions as the assembly nucleus in the maturation of U1, U2, U4 and U5 spliceosomal snRNPs. The structure of the human U4 snRNP core domain was initially solved at 3.6 Å resolution by experimental phasing using data with tetartohedral twinning. Molecular replacement from this model followed by density modification using untwinned data recently led to a structure of the minimal U1 snRNP at 3.3 Å resolution. With the latter structure providing a search model for molecular replacement, the U4 core-domain structure has now been re-refined. The U4 Sm site-sequence AAUUUUU has been shown to bind to the seven Sm proteins SmF-SmE-SmG-SmD3-SmB-SmD1-SmD2 in an identical manner as the U1 Sm-site sequence AAUUUGU, except in SmD1 where the bound U replaces G. The progression from the initial to the re-refined structure exemplifies a tortuous route to accuracy: where well diffracting crystals of complex assemblies are initially unavailable, the early model errors are rectified by exploiting preliminary interpretations in further experiments involving homologous structures. New insights are obtained from the more accurate model.

Using Molecular Replacement Phasing to Study the Structure and Function of RNA.

In recent years a wide variety of RNA molecules regulating fundamental cellular processes has been discovered. Therefore, RNA structure determination is experiencing a boost and many more RNA structures are likely to be determined in the years to come. The broader availability of experimentally determined RNA structures implies that molecular replacement (MR) will be used more and more frequently as a method for phasing future crystallographic structures. In this report we describe various aspects relative to RNA structure determination by MR. First, we describe how to select and create MR search models for nucleic acids. Second, we describe how to perform MR searches on RNA using available crystallographic software. Finally, we describe how to refine and interpret the successful MR solutions. These protocols are applicable to determine novel RNA structures as well as to establish structural-functional relationships on existing RNA structures.

Exploring the speed and performance of molecular replacement with AMPLE using QUARK ab initio protein models.

AMPLE clusters and truncates ab initio protein structure predictions, producing search models for molecular replacement. Here, an interesting degree of complementarity is shown between targets solved using the different ab initio modelling programs QUARK and ROSETTA. Search models derived from either program collectively solve almost all of the all-helical targets in the test set. Initial solutions produced by Phaser after only 5 min perform surprisingly well, improving the prospects for in situ structure solution by AMPLE during synchrotron visits. Taken together, the results show the potential for AMPLE to run more quickly and successfully solve more targets than previously suspected.

Application of theAMPLEcluster-and-truncate approach to NMR structures for molecular replacement

AMPLE is a program developed for clustering and truncating ab initio protein structure predictions into search models for molecular replacement. Here, it is shown that its core cluster-and-truncate methods also work well for processing NMR ensembles into search models. Rosetta remodelling helps to extend success to NMR structures bearing low sequence identity or high structural divergence from the target protein. Potential future routes to improved performance are considered and practical, general guidelines on using AMPLE are provided.

Crystal structure of the C-terminal globular domain of the third paralog of the Archaeoglobus fulgidus oligosaccharyltransferases

BackgroundProtein N-glycosylation occurs in the three domains of life. Oligosaccharyltransferase (OST) transfers an oligosaccharide chain to the asparagine residue in the N-glycosylation sequons. The catalytic subunits of the OST enzyme are STT3 in eukaryotes, AglB in archaea and PglB in eubacteria. The genome of a hyperthermophilic archaeon, Archaeoglobus fulgidus, encodes three paralogous AglB proteins. We previously solved the crystal structures of the C-terminal globular domains of two paralogs, AglB-Short 1 and AglB-Short 2.ResultsWe determined the crystal structure of the C-terminal globular domain of the third AglB paralog, AglB-Long, at 1.9 Å resolutions. The crystallization of the fusion protein with maltose binding protein (MBP) afforded high quality protein crystals. Two MBP-AglB-L molecules formed a swapped dimer in the crystal. Since the fusion protein behaved as a monomer upon gel filtration, we reconstituted the monomer structure from the swapped dimer by exchanging the swapped segments. The C-terminal domain of A. fulgidus AglB-L includes a structural unit common to AglB-S1 and AglB-S2. This structural unit contains the evolutionally conserved WWDYG and DK motifs. The present structure revealed that A. fulgidus AglB-L contained a variant type of the DK motif with a short insertion, and confirmed that the second signature residue, Lys, of the DK motif participates in the formation of a pocket that binds to the serine and threonine residues at the +2 position of the N-glycosylation sequon.ConclusionsThe structure of A. fulgidus AglB-L, together with the two previously solved structures of AglB-S1 and AglB-S2, provides a complete overview of the three AglB paralogs encoded in the A. fulgidus genome. All three AglBs contain a variant type of the DK motif. This finding supports a previously proposed rule: The STT3/AglB/PglB paralogs in one organism always contain the same type of Ser/Thr-binding pocket. The present structure will be useful as a search model for molecular replacement in the structural determination of the full-length A. fulgidus AglB-L.

The interplay between RPGR, PDEδ and Arl2/3 regulate the ciliary targeting of farnesylated cargo

Defects in primary cilia result in human diseases known as ciliopathies. The retinitis pigmentosa GTPase regulator (RPGR), mutated in the most severe form of the eye disease, is located at the transition zone of the ciliary organelle. The RPGR-interacting partner PDEδ is involved in trafficking of farnesylated ciliary cargo, but the significance of this interaction is unknown. The crystal structure of the propeller domain of RPGR shows the location of patient mutations and how they perturb the structure. The RPGR·PDEδ complex structure shows PDEδ on a highly conserved surface patch of RPGR. Biochemical experiments and structural considerations show that RPGR can bind with high affinity to cargo-loaded PDEδ and exposes the Arl2/Arl3-binding site on PDEδ. On the basis of these results, we propose a model where RPGR is acting as a scaffold protein recruiting cargo-loaded PDEδ and Arl3 to release lipidated cargo into cilia.

Comparison of the Crystal and Solution Structures of a Cyanobacterial Hemoglobin

GlbN is a monomeric 2/2 hemoglobin found in several cyanobacteria. Physiological data suggest that, in Synechococcus sp. PCC 7002, GlbN protects the cell from damage associated with nitrate metabolism [1]. Under reducing conditions, GlbN undergoes a facile posttranslational cross-linking between His117 and the heme 2-vinyl. In a study of the enzymatic and chemical properties of GlbN, the solution structure of the protein in the ferric, bis-histidine state with heme posttranslational modification was determined by NMR methods. As this protein is paramagnetic (S = 1/2) and can be prepared in the reduced state (S = 0), heme methyl hyperfine chemical shifts were used to determine the orientation of the ligating histidines, and proton pseudocontact shifts were applied in the refinement. The NMR model contains regions that appear to sample multiple conformations, specifically the A helix and EF loop. Here we present the X-ray structure of GlbN. Crystals were obtained by the vapor diffusion hanging drop method, and data were collected under cryogenic conditions (National Synchrotron Light Source, BNL). The structure of the related Synechocystis sp. PCC 6803 GlbN (1RTX) was used as the search model for molecular replacement. Iterative model building yielded a final structure at 1.65 A resolution. The X-ray data revealed a ruffled heme, provided a reliable representation of its environment in solution, and allowed for an assessment of the paramagnetic NMR approach. The X-ray data also gave an indication of the effect of crystal packing on the flexible regions of GlbN. Comparison with the Synechocystis GlbN structure highlighted the resilience of the 2/2 globin fold as well as subtle differences in the solvent accessibility of the heme active site.[1] Scott et al. (2010) Biochemistry 49:7000Supported by NSF MCB 0843439

Crystallographic and NMR Evidence for Flexibility in Oligosaccharyltransferases and Its Catalytic Significance

Oligosaccharyltransferase (OST) is a membrane-bound enzyme that catalyzes the transfer of an oligosaccharide to an asparagine residue in glycoproteins. It possesses a binding pocket that recognizes Ser and Thr residues at the+2 position in theN-glycosylation consensus, Asn-X-Ser/Thr. We determined the crystal structures of the C-terminal globular domains of the catalytic subunits of two archaeal OSTs. A comparison with previously determined structures identified a segment with remarkable conformational plasticity, induced by crystal contact effects. We characterized its dynamic properties in solution by (15)N NMR relaxation analyses. Intriguingly, the mobile region contains the+2 Ser/Thr-binding pocket. In agreement, the flexibility restriction forced by an engineered disulfide crosslink abolished the enzymatic activity, and its cleavage fully restored activity. These results suggest the necessity of multiple conformational states in the reaction. The dynamic nature of the Ser/Thr pocket could facilitate the efficient scanning of N-glycosylation sequons along nascent polypeptide chains.

Nearest-cell: a fast and easy tool for locating crystal matches in the PDB.

When embarking upon X-ray diffraction data collection from a potentially novel macromolecular crystal form, it can be useful to ascertain whether the measured data reflect a crystal form that is already recorded in the Protein Data Bank and, if so, whether it is part of a large family of related structures. Providing such information to crystallographers conveniently and quickly, as soon as the first images have been recorded and the unit cell characterized at an X-ray beamline, has the potential to save time and effort as well as pointing to possible search models for molecular replacement. Given an input unit cell, and optionally a space group, Nearest-cell rapidly scans the Protein Data Bank and retrieves near-matches.

The Atypical Response Regulator Protein ChxR Has Structural Characteristics and Dimer Interface Interactions That Are Unique within the OmpR/PhoB Subfamily

Typically as a result of phosphorylation, OmpR/PhoB response regulators form homodimers through a receiver domain as an integral step in transcriptional activation. Phosphorylation stabilizes the ionic and hydrophobic interactions between monomers. Recent studies have shown that some response regulators retain functional activity in the absence of phosphorylation and are termed atypical response regulators. The two currently available receiver domain structures of atypical response regulators are very similar to their phospho-accepting homologs, and their propensity to form homodimers is generally retained. An atypical response regulator, ChxR, from Chlamydia trachomatis, was previously reported to form homodimers; however, the residues critical to this interaction have not been elucidated. We hypothesize that the intra- and intermolecular interactions involved in forming a transcriptionally competent ChxR are distinct from the canonical phosphorylation (activation) paradigm in the OmpR/PhoB response regulator subfamily. To test this hypothesis, structural and functional studies were performed on the receiver domain of ChxR. Two crystal structures of the receiver domain were solved with the recently developed method using triiodo compound I3C. These structures revealed many characteristics unique to OmpR/PhoB subfamily members: typical or atypical. Included was the absence of two α-helices present in all other OmpR/PhoB response regulators. Functional studies on various dimer interface residues demonstrated that ChxR forms relatively stable homodimers through hydrophobic interactions, and disruption of these can be accomplished with the introduction of a charged residue within the dimer interface. A gel shift study with monomeric ChxR supports that dimerization through the receiver domain is critical for interaction with DNA.

Improved Technologies Now Routinely Provide Protein NMR Structures Useful for Molecular Replacement

Molecular replacement (MR) is widely used for addressing the phase problem in X-ray crystallography. Historically, crystallographers have had limited success using NMR structures as MR search models. Here, we report a comprehensive investigation of the utility of protein NMR ensembles as MR search models, using data for 25 pairs of X-ray and NMR structures solved and refined using modern NMR methods. Starting from NMR ensembles prepared by an improved protocol, FindCore, correct MR solutions were obtained for 22 targets. Based on these solutions, automatic model rebuilding could be done successfully. Rosetta refinement of NMR structures provided MR solutions for another two proteins. We also demonstrate that such properly prepared NMR ensembles and X-ray crystal structures have similar performance when used as MR search models for homologous structures, particularly for targets with sequence identity >40%.

Evaluating the solution from MrBUMP and BALBES.

Molecular replacement is one of the key methods used to solve the problem of determining the phases of structure factors in protein structure solution from X-ray image diffraction data. Its success rate has been steadily improving with the development of improved software methods and the increasing number of structures available in the PDB for use as search models. Despite this, in cases where there is low sequence identity between the target-structure sequence and that of its set of possible homologues it can be a difficult and time-consuming chore to isolate and prepare the best search model for molecular replacement. MrBUMP and BALBES are two recent developments from CCP4 that have been designed to automate and speed up the process of determining and preparing the best search models and putting them through molecular replacement. Their intention is to provide the user with a broad set of results using many search models and to highlight the best of these for further processing. An overview of both programs is presented along with a description of how best to use them, citing case studies and the results of large-scale testing of the software.

PVMR: assembling small helix fragments as structural solutions for molecular replacement

A new real-space implementation of the molecular-replacement method is described. The method locates the search model in the target crystal by maximizing the matching between the search-model vectors and the Patterson self and cross vectors. In previous work, a new rotation function was introduced for the molecular-replacement method [Jiang (2008). Acta Cryst. D64, 561-566]. This rotation function is calculated by matching the search model directly with both the Patterson self and cross vectors in real space. All the matches are summed and averaged to enhance the overall signal-to-noise ratio for a given orientation of the search model. Recently, to avoid the dependence of the weights derived from the linear regression on the properties of the search model and the target crystal structure, such as secondary structures, space groups and cell parameters, a dynamic correlation coefficient has been designed and used as the total rotation function score [Jiang & Ding (2010). Chin. Phys. B, 19, 106101]. This work further extends this idea to the implementation of translation search. A new real- or direct-space translation function has been implemented by matching the cross vectors between the symmetry mates of the search model to the Patterson cross vectors. This method enables effective searching for small helix fragments in the target crystal. Although the solution model assembled by using multiple fragments of helix is insufficient to start ab initio phasing of the target crystal, it can be used to identify the known protein folds in the Protein Data Bank that are homologous to the target structure. It can also be combined with other experimental and theoretical models to screen and select for better search models for molecular replacement.

Crystallization and preliminary crystallographic analysis of recombinant VSP1 from Arabidopsis thaliana.

VSP1 is a defence protein in Arabidopsis thaliana that may also be involved in control of plant development. The recombinant protein has been overexpressed in Escherichia coli, purified and crystallized using the sitting-drop vapour-diffusion method. The crystal diffracted to 1.9 A resolution and a complete X-ray data set was collected at 100 K using Cu Kalpha radiation from a rotating-anode X-ray source. The crystals belonged to space group C2. As there are no related structures that could be used as a search model for molecular replacement, work is in progress on experimental phasing using heavy-atom derivatives and selenomethionine derivatives.

Techniques and tactics used in determining the structure of the trimeric ebolavirus glycoprotein.

The trimeric membrane-anchored ebolavirus envelope glycoprotein (GP) is responsible for viral attachment, fusion and entry. Knowledge of its structure is important both for understanding ebolavirus entry and for the development of medical interventions. Crystal structures of viral glycoproteins, especially those in their metastable prefusion oligomeric states, can be difficult to achieve given the challenges in production, purification, crystallization and diffraction that are inherent in the heavily glycosylated flexible nature of these types of proteins. The crystal structure of ebolavirus GP in its trimeric prefusion conformation in complex with a human antibody derived from a survivor of the 1995 Kikwit outbreak has now been determined [Lee et al. (2008), Nature (London), 454, 177-182]. Here, the techniques, tactics and strategies used to overcome a series of technical roadblocks in crystallization and phasing are described. Glycoproteins were produced in human embryonic kidney 293T cells, which allowed rapid screening of constructs and expression of protein in milligram quantities. Complexes of GP with an antibody fragment (Fab) promoted crystallization and a series of deglycosylation strategies, including sugar mutants, enzymatic deglycosylation, insect-cell expression and glycan anabolic pathway inhibitors, were attempted to improve the weakly diffracting glycoprotein crystals. The signal-to-noise ratio of the search model for molecular replacement was improved by determining the structure of the uncomplexed Fab. Phase combination with Fab model phases and a selenium anomalous signal, followed by NCS-averaged density modification, resulted in a clear interpretable electron-density map. Model building was assisted by the use of B-value-sharpened electron-density maps and the proper sequence register was confirmed by building alternate sequences using N-linked glycan sites as anchors and secondary-structural predictions.