Real-space Refinement Research Articles

Automated Real-Space Refinement of Protein Structures Using a Realistic Backbone Move Set

Crystals of many important biological macromolecules diffract to limited resolution, rendering accurate model building and refinement difficult and time-consuming. We present a torsional optimization protocol that is applicable to many such situations and combines Protein Data Bank-based torsional optimization with real-space refinement against the electron density derived from crystallography or cryo-electron microscopy. Our method converts moderate- to low-resolution structures at initial (e.g., backbone trace only) or late stages of refinement to structures with increased numbers of hydrogen bonds, improved crystallographic R-factors, and superior backbone geometry. This automated method is applicable to DNA-binding and membrane proteins of any size and will aid studies of structural biology by improving model quality and saving considerable effort. The method can be extended to improve NMR and other structures. Our backbone score and its sequence profile provide an additional standard tool for evaluating structural quality.

BCL::EM-Fit: Rigid body fitting of atomic structures into density maps using geometric hashing and real space refinement

Cryo-electron microscopy (cryoEM) can visualize large macromolecular assemblies at resolutions often below 10Å and recently as good as 3.8-4.5 Å. These density maps provide important insights into the biological functioning of molecular machineries such as viruses or the ribosome, in particular if atomic-resolution crystal structures or models of individual components of the assembly can be placed into the density map. The present work introduces a novel algorithm termed BCL::EM-Fit that accurately fits atomic-detail structural models into medium resolution density maps. In an initial step, a "geometric hashing" algorithm provides a short list of likely placements. In a follow up Monte Carlo/Metropolis refinement step, the initial placements are optimized by their cross correlation coefficient. The resolution of density maps for a reliable fit was determined to be 10 Å or better using tests with simulated density maps. The algorithm was applied to fitting of capsid proteins into an experimental cryoEM density map of human adenovirus at a resolution of 6.8 and 9.0 Å, and fitting of the GroEL protein at 5.4 Å. In the process, the handedness of the cryoEM density map was unambiguously identified. The BCL::EM-Fit algorithm offers an alternative to the established Fourier/Real space fitting programs. BCL::EM-Fit is free for academic use and available from a web server or as downloadable binary file at http://www.meilerlab.org.

An Improved Model for Dynamin Assembly Revealed by Cryo-EM

Dynamin is a multidomain GTPase that assembles into collar-like structures at the necks of deeply invaginated coated pits during the final stages of clathrin-mediated endocytosis (CME) and catalyzes membrane scission. Assembly of purified dynamin tetramers in vitro yields helical structures comparable to those observed in vivo. The formation of these oligomers stimulates dynamin's basal GTP hydrolysis >100-fold. Mutational analysis indicates that dynamin's stimulated GTP hydrolysis is required for CME; however, mounting evidence suggests that this activity causes disassembly of the dynamin collar rather than direct membrane severing. Despite recent structural studies showing that stimulated hydrolysis arises from the transition-dependent dimerization of dynamin's catalytic G domains, little is known about the conformational changes that precede and/or result from this interaction in the context of the polymer. Specifically, it is unclear how the G domains are properly oriented, which subunits associate, and how catalysis triggers dissociation of the pleckstrin homology (PH) domain at the membrane surface. Much of this ambiguity can be attributed to the low resolution (>20A) of previous dynamin polymer models and the absence of a complete dynamin tetramer crystal structure. To clarify these issues, we have used cryo-EM and iterative helical real space refinement to generate an 11A reconstruction of a truncated form of dynamin (ΔPRD) in the assembled, GMPPCP-bound state. This map reveals new structural characteristics including a twisted, interlacing interaction that stabilizes the middle/GED stalk and a previously uncharacterized density feature adjacent to the exterior GTPase head. Computational docking of crystallized dynamin fragments reveals the location and connectivity of different domains within the assembled polymer. Chemical crosslinking experiments also provide new insights into the architecture and organization of dynamin tetramer. These data have important implications regarding the conformational changes associated with dynamin catalyzed GTP hydrolysis and membrane fission.

Performance of phased rotation, conformation and translation function: accurate protein model building with tripeptidic and tetrapeptidic fragments

The automatic building of protein structures with tripeptidic and tetrapeptidic fragments was investigated. The oligopeptidic conformers were positioned in the electron-density map by a phased rotation, conformation and translation function and refined by a real-space refinement. The number of successfully located fragments lay within the interval 75-95% depending on the resolution and phase quality. The overlaps of partially located fragments were analyzed. The correctly positioned fragments were connected into chains. Chains formed in this way were extended directly into the electron density and a sequence was assigned. In the initial stage of the model building the number of located fragments was between 60% and 95%, but this number could be increased by several cycles of reciprocal-space refinement and automatic model rebuilding. A nearly complete structure can be obtained on the condition that the resolution is reasonable. Computer graphics will only be needed for a final check and small corrections.

Der p 5 Crystal Structure Provides Insight into the Group 5 Dust Mite Allergens

Group 5 allergens from house dust mites elicit strong IgE antibody binding in mite-allergic patients. The structure of Der p 5 was determined by x-ray crystallography to better understand the IgE epitopes, to investigate the biologic function in mites, and to compare with the conflicting published Blo t 5 structures, designated 2JMH and 2JRK in the Protein Data Bank. Der p 5 is a three-helical bundle similar to Blo t 5, but the interactions of the helices are more similar to 2JMH than 2JRK. The crystallographic asymmetric unit contains three dimers of Der p 5 that are not exactly alike. Solution scattering techniques were used to assess the multimeric state of Der p 5 in vitro and showed that the predominant state was monomeric, similar to Blo t 5, but larger multimeric species are also present. In the crystal, the formation of the Der p 5 dimer creates a large hydrophobic cavity of approximately 3000 A(3) that could be a ligand-binding site. Many allergens are known to bind hydrophobic ligands, which are thought to stimulate the innate immune system and have adjuvant-like effects on IgE-mediated inflammatory responses.

Features and development of Coot

Coot is a molecular-graphics application for model building and validation of biological macromolecules. The program displays electron-density maps and atomic models and allows model manipulations such as idealization, real-space refinement, manual rotation/translation, rigid-body fitting, ligand search, solvation, mutations, rotamers and Ramachandran idealization. Furthermore, tools are provided for model validation as well as interfaces to external programs for refinement, validation and graphics. The software is designed to be easy to learn for novice users, which is achieved by ensuring that tools for common tasks are 'discoverable' through familiar user-interface elements (menus and toolbars) or by intuitive behaviour (mouse controls). Recent developments have focused on providing tools for expert users, with customisable key bindings, extensions and an extensive scripting interface. The software is under rapid development, but has already achieved very widespread use within the crystallographic community. The current state of the software is presented, with a description of the facilities available and of some of the underlying methods employed.

Combining density functional theory (DFT) and pair distribution function (PDF) analysis to solve the structure of metastable materials: the case of metakaolin

Understanding the atomic structure of complex metastable (including glassy) materials is of great importance in research and industry, however, such materials resist solution by most standard techniques. Here, a novel technique combining thermodynamics and local structure is presented to solve the structure of the metastable aluminosilicate material metakaolin (calcined kaolinite) without the use of chemical constraints. The structure is elucidated by iterating between least-squares real-space refinement using neutron pair distribution function data, and geometry optimisation using density functional modelling. The resulting structural representation is both energetically feasible and in excellent agreement with experimental data. This accurate structural representation of metakaolin provides new insight into the local environment of the aluminium atoms, with evidence of the existence of tri-coordinated aluminium. By the availability of this detailed chemically feasible atomic description, without the need to artificially impose constraints during the refinement process, there exists the opportunity to tailor chemical and mechanical processes involving metakaolin and other complex metastable materials at the atomic level to obtain optimal performance at the macro-scale.

Analysis Of Ribosomal Dynamics As Revealed By Cryo-EM And Flexible Fitting

Ribosomes are the molecular machines that translate the genetic message into nascent peptides, through a complex dynamics interplay with mRNAs, tRNAs, and various protein factors. A prominent example for ribosomal dynamics is the rotation of the small ribosomal subunit with respect to the large subunit, characterized as the “ratchet motion”, which is triggered by the binding of several translation factors. Based on density maps of ribosomal complexes obtained by cryo-EM, we analyzed two kinds of ribosomal ratchet motions, induced by the binding of EF-G and RF3, respectively. By using the flexible fitting technique (real-space refinement) (1) and an RNA secondary structure display tool (coloRNA) (2), quasi-atomic models of the ribosome were obtained in these ratchet-motion-related functional states. The observed differences in rRNA were further mapped onto the highly conserved RNA secondary structure diagram. Comparisons between the two sets of ratchet motions revealed that, while the overall patterns of the RNA displacement are very similar, several local regions stand out in their differential behavior, including the highly conserved GAC (GTPase-associated-center) region. We postulate that these regions are important in modulating the general ratchet motion and bestowing it with the dynamic characteristics required for the specific function.(1) Gao H, Sengupta J, Valle M, Korostelev A, Eswar N, Stagg SM, Van Roey P, Agrawal RK, Harvey SC, Sali A, Chapman MS, Frank J, 2003. Study of the structural dynamics of the E. coli 70S ribosome using real space refinement. Cell 113: 789-801.(2) LeBarron J, Mitra K, Frank J, 2007. Displaying 3D data on RNA secondary structures: coloRNA. J Struct Biol 157: 262-70.

Autofix for backward-fit sidechains: using MolProbity and real-space refinement to put misfits in their place

Misfit sidechains in protein crystal structures are a stumbling block in using those structures to direct further scientific inference. Problems due to surface disorder and poor electron density are very difficult to address, but a large class of systematic errors are quite common even in well-ordered regions, resulting in sidechains fit backwards into local density in predictable ways. The MolProbity web site is effective at diagnosing such errors, and can perform reliable automated correction of a few special cases such as 180° flips of Asn or Gln sidechain amides, using all-atom contacts and H-bond networks. However, most at-risk residues involve tetrahedral geometry, and their valid correction requires rigorous evaluation of sidechain movement and sometimes backbone shift. The current work extends the benefits of robust automated correction to more sidechain types. The Autofix method identifies candidate systematic, flipped-over errors in Leu, Thr, Val, and Arg using MolProbity quality statistics, proposes a corrected position using real-space refinement with rotamer selection in Coot, and accepts or rejects the correction based on improvement in MolProbity criteria and on χ angle change. Criteria are chosen conservatively, after examining many individual results, to ensure valid correction. To test this method, Autofix was run and analyzed for 945 representative PDB files and on the 50S ribosomal subunit of file 1YHQ. Over 40% of Leu, Val, and Thr outliers and 15% of Arg outliers were successfully corrected, resulting in a total of 3,679 corrected sidechains, or 4 per structure on average. Summary Sentences: A common class of misfit sidechains in protein crystal structures is due to systematic errors that place the sidechain backwards into the local electron density. A fully automated method called “Autofix” identifies such errors for Leu, Val, Thr, and Arg and corrects over one third of them, using MolProbity validation criteria and Coot real-space refinement of rotamers.Electronic supplementary materialThe online version of this article (doi:10.1007/s10969-008-9045-8) contains supplementary material, which is available to authorized users.

Cryo-EM Structure of the DNA-Dependent Protein Kinase Catalytic Subunit at Subnanometer Resolution Reveals α Helices and Insight into DNA Binding

The DNA-dependent protein kinase catalytic subunit (DNA-PKcs) regulates the nonhomologous end joining pathway for repair of double-stranded DNA (dsDNA) breaks. Here, we present a 7A resolution structure of DNA-PKcs determined by cryo-electron microscopy single-particle reconstruction. This structure is composed of density rods throughout the molecule that are indicative of alpha helices and reveals structural features not observed in lower resolution EM structures. Docking of homology models into the DNA-PKcs structure demonstrates that up to eight helical HEAT repeat motifs fit well within the density. Surprisingly, models for the kinase domain can be docked into either the crown or base of the molecule at this resolution, although real space refinement suggests that the base location is the best fit. We propose a model for the interaction of DNA with DNA-PKcs in which one turn of dsDNA enters the central channel and interacts with a resolved alpha-helical protrusion.

Effect of Calcium on Calmodulin Bound to the IQ Motifs of Myosin V

The long neck of unconventional myosin V is composed of six tandem "IQ motifs," which are fully occupied by calmodulin (CaM) in the absence of calcium. Calcium regulates the activity, the folded-to-extended conformational transition, and the processive run length of myosin V, and thus, it is important to understand how calcium affects CaM binding to the IQ motifs. Here we used electron cryomicroscopy together with computer-based docking of crystal structures into three-dimensional reconstructions of actin decorated with a motor domain-two IQ complex to provide an atomic model of myosin V in the presence of calcium. Calcium causes a major rearrangement of the bound CaMs, dissociation of CaM bound to IQ motif 2, and propagated changes in the motor domain. Tryptophan fluorescence spectroscopy showed that calcium-CaM binds to IQ motifs 1, 3, and 5 in a different conformation than apoCaM. Proteolytic cleavage was consistent with CaM preferentially dissociating from the second IQ motif. The enzymatic and mechanical functions of myosin V can, therefore, be modulated both by calcium-dependent conformational changes of bound CaM as well as by CaM dissociation.

Real-space refinement of single-crystal electron diffuse scattering and its application toBi2Ru2O7−δ

A real-space atomistic refinement approach to the analysis of experimental electrondiffraction patterns is described. The method employs the reverse Monte Carlo algorithmto produce atomistic configurations capable of qualitatively reproducing diffuse electronscattering patterns. Its implementation in the program EDRMC is described indetail, together with a number of additional constraints/restraints that can beused to guide the refinement process. In particular, appropriate restraints ensurethe individual atomic displacements introduced to model the diffuse scatteringpatterns are simultaneously consistent with the known average structure. Theapproach is then used to interpret electron diffraction patterns measured forBi2Ru2O7−δ. The diffuse scattering patterns observed are shown to arise primarily from concertedtranslations of Bi atoms. These translations can be interpreted in terms of rotations of[O′Bi4 ] tetrahedra correlated along the crystal axes and uncorrelated along orthogonal directions.

PDFfit2 and PDFgui: computer programs for studying nanostructure in crystals

PDFfit2 is a program as well as a library for real-space refinement of crystal structures. Itis capable of fitting a theoretical three-dimensional (3D) structure to atomic pairdistribution function data and is ideal for nanoscale investigations. The fit system accountsfor lattice constants, atomic positions and anisotropic atomic displacement parameters,correlated atomic motion, and experimental factors that may affect the data.The atomic positions and thermal coefficients can be constrained to follow thesymmetry requirements of an arbitrary space group. The PDFfit2 engine is written inC++ and is accessible via Python, allowing it to inter-operate with other Pythonprograms. PDFgui is a graphical interface built on the PDFfit2 engine. PDFguiorganizes fits and simplifies many data analysis tasks, such as configuring andplotting multiple fits. PDFfit2 and PDFgui are freely available via the Internet.

High-resolution Electron Microscopy of Helical Specimens: A Fresh Look at Tobacco Mosaic Virus

The treatment of helical objects as a string of single particles has become an established technique to resolve their three-dimensional (3D) structure using electron cryo-microscopy. It can be applied to a wide range of helical particles such as viruses, microtubules and helical filaments. We have made improvements to this approach using Tobacco Mosaic Virus (TMV) as a test specimen and obtained a map from 210,000 asymmetric units at a resolution better than 5 Å. This was made possible by performing a full correction of the contrast transfer function of the microscope. Alignment of helical segments was helped by constraints derived from the helical symmetry of the virus. Furthermore, symmetrization was implemented by multiple inclusions of symmetry-related views in the 3D reconstruction. We used the density map to build an atomic model of TMV. The model was refined using a real-space refinement strategy that accommodates multiple conformers. The atomic model shows significant deviations from the deposited model for the helical form of TMV at the lower-radius region (residues 88 to 109). This region appears more ordered with well-defined secondary structure, compared with the earlier helical structure. The RNA phosphate backbone is sandwiched between two arginine side-chains, stabilizing the interaction between RNA and coat protein. A cluster of two or three carboxylates is buried in a hydrophobic environment isolating it from neighboring subunits. These carboxylates may represent the so-called Caspar carboxylates that form a metastable switch for viral disassembly. Overall, the observed differences suggest that the new model represents a different, more stable state of the virus, compared with the earlier published model.

Structures of modified eEF2·80S ribosome complexes reveal the role of GTP hydrolysis in translocation

On the basis of kinetic data on ribosome protein synthesis, the mechanical energy for translocation of the mRNA-tRNA complex is thought to be provided by GTP hydrolysis of an elongation factor (eEF2 in eukaryotes, EF-G in bacteria). We have obtained cryo-EM reconstructions of eukaryotic ribosomes complexed with ADP-ribosylated eEF2 (ADPR-eEF2), before and after GTP hydrolysis, providing a structural basis for analyzing the GTPase-coupled mechanism of translocation. Using the ADP-ribosyl group as a distinct marker, we observe conformational changes of ADPR-eEF2 that are due strictly to GTP hydrolysis. These movements are likely representative of native eEF2 motions in a physiological context and are sufficient to uncouple the mRNA-tRNA complex from two universally conserved bases in the ribosomal decoding center (A1492 and A1493 in Escherichia coli) during translocation. Interpretation of these data provides a detailed two-step model of translocation that begins with the eEF2/EF-G binding-induced ratcheting motion of the small ribosomal subunit. GTP hydrolysis then uncouples the mRNA-tRNA complex from the decoding center so translocation of the mRNA-tRNA moiety may be completed by a head rotation of the small subunit.

Improving amino-acid identification, fit and Cαprediction using the Simplex method in automated model building

Automated methods for protein model building in X-ray crystallography typically use a two-phased approach that involves first modeling the protein backbone followed by building in the side chains. The latter phase requires the identification of the amino-acid side-chain type as well as fitting of the side-chain model into the observed electron density. While mistakes in identification of individual side chains are common for a number of reasons, sequence alignment can sometimes be used to correct errors by mapping fragments into the true (expected) amino-acid sequence and exploiting contiguity constraints among neighbors. However, side chains cannot always be confidently aligned; this depends on having sufficient accuracy in the initial calls. The recognition of amino-acid side-chains based on the surrounding pattern of electron density, whether by features, density correlation or free atoms, can be sensitive to inaccuracies in the coordinates of the predicted backbone C(alpha) atoms to which they are anchored. By incorporating a Nelder-Mead Simplex search into the side-chain identification and model-building routines of TEXTAL, it is demonstrated that this form of residue-by-residue rigid-body real-space refinement (in which the C(alpha) itself is allowed to shift) can improve the initial accuracy of side-chain selection by over 25% on average (from 25% average identity to 32% on a test set of five representative proteins, without corrections by sequence alignment). This improvement in amino-acid selection accuracy in TEXTAL is often sufficient to bring the pairwise amino-acid identity of chains in the model out of the so-called ;twilight zone' for sequence-alignment methods. When coupled with sequence alignment, use of the Simplex search yielded improvements in side-chain accuracy on average by over 13 percentage points (from 64 to 77%) and up to 38 percentage points (from 40 to 78%) in one case compared with using sequence alignment alone.

Elongation Arrest by SecM via a Cascade of Ribosomal RNA Rearrangements

In E. coli, the SecM nascent polypeptide causes elongation arrest, while interacting with 23S RNA bases A2058 and A749-753 in the exit tunnel of the large ribosomal subunit. We compared atomic models fitted by real-space refinement into cryo-electron microscopy reconstructions of a pretranslocational and SecM-stalled E. coli ribosome complex. A cascade of RNA rearrangements propagates from the exit tunnel throughout the large subunit, affecting intersubunit bridges and tRNA positions, which in turn reorient small subunit RNA elements. Elongation arrest could result from the inhibition of mRNA.(tRNAs) translocation, E site tRNA egress, and perhaps translation factor activation at the GTPase-associated center. Our study suggests that the specific secondary and tertiary arrangement of ribosomal RNA provides the basis for internal signal transduction within the ribosome. Thus, the ribosome may itself have the ability to regulate its progression through translation by modulating its structure and consequently its receptivity to activation by cofactors.

Molding Atomic Structures into Intermediate- Resolution Cryo-EM Density Maps of Ribosomal Complexes Using Real-Space Refinement

Real-space refinement has been previously introduced as a flexible fitting method to interpret medium-resolution cryo-EM density maps in terms of atomic structures. In this way, conformational changes related to functional processes can be analyzed on the molecular level. In the application of the technique to the ribosome, quasiatomic models have been derived that have advanced our understanding of translocation. In this article, the choice of parameters for the fitting procedure is discussed. The quality of the fitting depends critically on the number of rigid pieces into which the model is divided. Suitable quality indicators are crosscorrelation, R factor, and density residual, all of which can also be locally applied. The example of the ribosome may provide some guidelines for general applications of real-space refinement to flexible fitting problems.

Fitting of High-Resolution Structures into Electron Microscopy Reconstruction Images

Dynamic macromolecular assemblies, such as ribosomes, viruses, and muscle protein complexes, are often more amenable to visualization by electron microscopy than by high-resolution X-ray crystallography or NMR. When high-resolution structures of component structures are available, it is possible to build an atomic model that gives information about the molecular interactions at greater detail than the experimental resolution, due to constraints of modeling placed upon the interpretation. There are now several competing computational methods to search systematically for orientations and positions of components that match the experimental image density, and continuing developments will be reviewed. Attention is now also moving toward the related task of optimization, with flexible and/or multifragment models and sometimes with stereochemically restrained refinement methods. This paper will review the various approaches and describe advances in the authors' methods and applications of real-space refinement.

Real-space protein-model completion: an inverse-kinematics approach

Rapid protein-structure determination relies greatly on software that can automatically build a protein model into an experimental electron-density map. In favorable circumstances, various software systems are capable of building over 90% of the final model. However, completeness falls off rapidly with the resolution of the diffraction data. Manual completion of these partial models is usually feasible, but is time-consuming and prone to subjective interpretation. Except for the N- and C-termini of the chain, the end points of each missing fragment are known from the initial model. Hence, fitting fragments reduces to an inverse-kinematics problem. A method has been developed that combines fast inverse-kinematics algorithms with a real-space torsion-angle refinement procedure in a two-stage approach to fit missing main-chain fragments into the electron density between two anchor points. The first stage samples a large number of closing conformations, guided by the electron density. These candidates are ranked according to density fit. In a subsequent refinement stage, optimization steps are projected onto a carefully chosen subspace of conformation space to preserve rigid geometry and closure. Experimental results show that fitted fragments are in excellent agreement with the final refined structure for lengths of up to 12-15 residues in areas of weak or ambiguous electron density, even at medium to low resolution.