Protein structure fitting and refinement guided by cryo-EM density

Structure of adeno-associated virus-2 in complex with neutralizing monoclonal antibody A20

The use of adeno-associated virus (AAV) as a gene therapy vector is limited by the host neutralizing immune response. The cryo-electron microscopy (EM) structure at 8.5Å resolution is determined for a complex of AAV-2 with the Fab' fragment of monoclonal antibody (MAb) A20, the most extensively characterized AAV MAb. The binding footprint is determined through fitting the cryo-EM reconstruction with a homology model following sequencing of the variable domain, and provides a structural basis for integrating diverse prior epitope mappings. The footprint extends from the previously implicated plateau to the side of the spike, and into the conserved canyon, covering a larger area than anticipated. Comparison with structures of binding and non-binding serotypes indicates that recognition depends on a combination of subtle serotype-specific features. Separation of the neutralizing epitope from the heparan sulfate cell attachment site encourages attempts to develop immune-resistant vectors that can still bind to target cells.

Structural analysis of coxsackievirus A7 reveals conformational changes associated with uncoating

Coxsackievirus A7 (CAV7) is a rarely detected and poorly characterized serotype of the Enterovirus species Human enterovirus A (HEV-A) within the Picornaviridae family. The CAV7-USSR strain has caused polio-like epidemics and was originally thought to represent the fourth poliovirus type, but later evidence linked this strain to the CAV7-Parker prototype. Another isolate, CAV7-275/58, was also serologically similar to Parker but was noninfectious in a mouse model. Sequencing of the genomic region encoding the capsid proteins of the USSR and 275/58 strains and subsequent comparison with the corresponding amino acid sequences of the Parker strain revealed that the Parker and USSR strains are nearly identical, while the 275/58 strain is more distant. Using electron cryomicroscopy and three-dimensional image reconstruction, the structures of the CAV7-USSR virion and empty capsid were resolved to 8.2-Å and 6.1-Å resolutions, respectively. This is one of the first detailed structural analyses of the HEV-A species. Using homology modeling, reconstruction segmentation, and flexible fitting, we constructed a pseudoatomic T = 1 (pseudo T = 3) model incorporating the three major capsid proteins (VP1 to VP3), addressed the conformational changes of the capsid and its constituent viral proteins occurring during RNA release, and mapped the capsid proteins' variable regions to the structure. During uncoating, VP4 and RNA are released analogously to poliovirus 1, the interfaces of VP2 and VP3 are rearranged, and VP1 rotates. Variable regions in the capsid proteins were predicted to map mainly to the surface of VP1 and are thus likely to affect the tropism and pathogenicity of CAV7.

Capsid structure and its stability at the late stages of bacteriophage SPP1 assembly

The structure of the bacteriophage SPP1 capsid was determined at subnanometer resolution by cryo-electron microscopy and single-particle analysis. The icosahedral capsid is composed of the major capsid protein gp13 and the auxiliary protein gp12, which are organized in a T=7 lattice. DNA is arranged in layers with a distance of ~24.5 Å. gp12 forms spikes that are anchored at the center of gp13 hexamers. In a gp12-deficient mutant, the centers of hexamers are closed by loops of gp13 coming together to protect the SPP1 genome from the outside environment. The HK97-like fold was used to build a pseudoatomic model of gp13. Its structural organization remains unchanged upon tail binding and following DNA release. gp13 exhibits enhanced thermostability in the DNA-filled capsid. A remarkable convergence between the thermostability of the capsid and those of the other virion components was found, revealing that the overall architecture of the SPP1 infectious particle coevolved toward high robustness.

EBs recognize a nucleotide-dependent structural cap at growing microtubule ends

Growing microtubule ends serve as transient binding platforms for essential proteins that regulate microtubule dynamics and their interactions with cellular substructures. End-binding proteins (EBs) autonomously recognize an extended region at growing microtubule ends with unknown structural characteristics and then recruit other factors to the dynamic end structure. Using cryo-electron microscopy, subnanometer single-particle reconstruction, and fluorescence imaging, we present a pseudoatomic model of how the calponin homology (CH) domain of the fission yeast EB Mal3 binds to the end regions of growing microtubules. The Mal3 CH domain bridges protofilaments except at the microtubule seam. By binding close to the exchangeable GTP-binding site, the CH domain is ideally positioned to sense the microtubule's nucleotide state. The same microtubule-end region is also a stabilizing structural cap protecting the microtubule from depolymerization. This insight supports a common structural link between two important biological phenomena, microtubule dynamic instability and end tracking.

ATP-triggered conformational changes delineate substrate-binding and -folding mechanics of the GroEL chaperonin

The chaperonin GroEL assists the folding of nascent or stress-denatured polypeptides by actions of binding and encapsulation. ATP binding initiates a series of conformational changes triggering the association of the cochaperonin GroES, followed by further large movements that eject the substrate polypeptide from hydrophobic binding sites into a GroES-capped, hydrophilic folding chamber. We used cryo-electron microscopy, statistical analysis, and flexible fitting to resolve a set of distinct GroEL-ATP conformations that can be ordered into a trajectory of domain rotation and elevation. The initial conformations are likely to be the ones that capture polypeptide substrate. Then the binding domains extend radially to separate from each other but maintain their binding surfaces facing the cavity, potentially exerting mechanical force upon kinetically trapped, misfolded substrates. The extended conformation also provides a potential docking site for GroES, to trigger the final, 100° domain rotation constituting the “power stroke” that ejects substrate into the folding chamber.

Newly folded substrates inside the molecular cage of the HtrA chaperone DegQ

The HtrA protein family combines chaperone and protease activities and is essential for protein quality control in many organisms. Whereas the mechanisms underlying the proteolytic function of HtrA proteins are well characterized, their chaperone activity remains poorly understood. Here we describe cryo-EM structures of Escherichia coli DegQ in its 12- and 24-mer states in complex with model substrates, providing a structural model of HtrA chaperone action. Up to six lysozyme substrates bind inside the DegQ 12-mer cage and are visualized in a close-to-native state. An asymmetric reconstruction reveals the binding of a well-ordered lysozyme to four DegQ protomers. DegQ PDZ domains are located adjacent to substrate density and their presence is required for chaperone activity. The substrate-interacting regions appear conserved in 12- and 24-mer cages, suggesting a common mechanism of chaperone function.

Reconstructing virus structures from nanometer to near-atomic resolutions with cryo-electron microscopy and tomography

The past few decades have seen tremendous advances in single-particle electron -cryo-microscopy (cryo-EM). The field has matured to the point that near-atomic resolution density maps can be generated for icosahedral viruses without the need for crystallization. In parallel, substantial progress has been made in determining the structures of nonicosahedrally arranged proteins in viruses by employing either single-particle cryo-EM or cryo-electron tomography (cryo-ET). Implicit in this course have been the availability of a new generation of electron cryo-microscopes and the development of the computational tools that are essential for generating these maps and models. This methodology has enabled structural biologists to analyze structures in increasing detail for virus particles that are in different morphogenetic states. Furthermore, electron imaging of frozen, hydrated cells, in the process of being infected by viruses, has also opened up a new avenue for studying virus structures "in situ". Here we present the common techniques used to acquire and process cryo-EM and cryo-ET data and discuss their implications for structural virology both now and in the future.

Applications and limitations of in silico models in drug discovery

Drug discovery in the late twentieth and early twenty-first century has witnessed a myriad of changes that were adopted to predict whether a compound is likely to be successful, or conversely enable identification of molecules with liabilities as early as possible. These changes include integration of in silico strategies for lead design and optimization that perform complementary roles to that of the traditional in vitro and in vivo approaches. The in silico models are facilitated by the availability of large datasets associated with high-throughput screening, bioinformatics algorithms to mine and annotate the data from a target perspective, and chemoinformatics methods to integrate chemistry methods into lead design process. This chapter highlights the applications of some of these methods and their limitations. We hope this serves as an introduction to in silico drug discovery.

Macromolecular assembly structures by comparative modeling and electron microscopy

Advances in electron microscopy allow for structure determination of large biological machines at increasingly higher resolutions. A key step in this process is fitting component structures into the electron microscopy-derived density map of their assembly. Comparative modeling can contribute by providing atomic models of the components, via fold assignment, sequence-structure alignment, model building, and model assessment. All four stages of comparative modeling can also benefit from consideration of the density map. In this chapter, we describe numerous types of modeling problems restrained by a density map and available protocols for finding solutions. In particular, we provide detailed instructions for density map-guided modeling using the Integrative Modeling Platform (IMP), MODELLER, and UCSF Chimera.

Putting structure into context: fitting of atomic models into electron microscopic and electron tomographic reconstructions

A complete understanding of complex dynamic cellular processes such as cell migration or cell adhesion requires the integration of atomic level structural information into the larger cellular context. While direct atomic-level information at the cellular level remains inaccessible, electron microscopy, electron tomography and their associated computational image processing approaches have now matured to a point where sub-cellular structures can be imaged in three dimensions at the nanometer scale. Atomic-resolution information obtained by other means can be combined with this data to obtain three-dimensional models of large macromolecular assemblies in their cellular context. This article summarizes some recent advances in this field.

RANKING VALID TOPOLOGIES OF THE SECONDARY STRUCTURE ELEMENTS USING A CONSTRAINT GRAPH

Electron cryo-microscopy is a fast advancing biophysical technique to derive three-dimensional structures of large protein complexes. Using this technique, many density maps have been generated at intermediate resolution such as 6–10 Å resolution. Although it is challenging to derive the backbone of the protein directly from such density maps, secondary structure elements such as helices and β-sheets can be computationally detected. Our work in this paper provides an approach to enumerate the top-ranked possible topologies instead of enumerating the entire population of the topologies. This approach is particularly practical for large proteins. We developed a directed weighted graph, the topology graph, to represent the secondary structure assignment problem. We prove that the problem of finding the valid topology with the minimum cost is NP hard. We developed an O(N2 2N) dynamic programming algorithm to identify the topology with the minimum cost. The test of 15 proteins suggests that our dynamic programming approach is feasible to work with proteins of much larger size than we could before. The largest protein in the test contains 18 helical sticks detected from the density map out of 33 helices in the protein.

Consensus among flexible fitting approaches improves the interpretation of cryo-EM data

Cryo-elecron microscopy (cryo-EM) can provide important structural information of large macromolecular assemblies in different conformational states. Recent years have seen an increase in structures deposited in the Protein Data Bank (PDB) by fitting a high-resolution structure into its low-resolution cryo-EM map. A commonly used protocol for accommodating the conformational changes between the X-ray structure and the cryo-EM map is rigid body fitting of individual domains. With the emergence of different flexible fitting approaches, there is a need to compare and revise these different protocols for the fitting. We have applied three diverse automated flexible fitting approaches on a protein dataset for which rigid domain fitting (RDF) models have been deposited in the PDB. In general, a consensus is observed in the conformations, which indicates a convergence from these theoretically different approaches to the most probable solution corresponding to the cryo-EM map. However, the result shows that the convergence might not be observed for proteins with complex conformational changes or with missing densities in cryo-EM map. In contrast, RDF structures deposited in the PDB can represent conformations that not only differ from the consensus obtained by flexible fitting but also from X-ray crystallography. Thus, this study emphasizes that a "consensus" achieved by the use of several automated flexible fitting approaches can provide a higher level of confidence in the modeled configurations. Following this protocol not only increases the confidence level of fitting, but also highlights protein regions with uncertain fitting. Hence, this protocol can lead to better interpretation of cryo-EM data.

UCSF Chimera, MODELLER, and IMP: an integrated modeling system

Structural modeling of macromolecular complexes greatly benefits from interactive visualization capabilities. Here we present the integration of several modeling tools into UCSF Chimera. These include comparative modeling by MODELLER, simultaneous fitting of multiple components into electron microscopy density maps by IMP MultiFit, computing of small-angle X-ray scattering profiles and fitting of the corresponding experimental profile by IMP FoXS, and assessment of amino acid sidechain conformations based on rotamer probabilities and local interactions by Chimera.

The Cryo-EM structure of a complete 30S translation initiation complex from Escherichia coli

Formation of the 30S initiation complex (30S IC) is an important checkpoint in regulation of gene expression. The selection of mRNA, correct start codon, and the initiator fMet-tRNA^fMet requires the presence of three initiation factors (IF1, IF2, IF3) of which IF3 and IF1 control the fidelity of the process, while IF2 recruits fMet-tRNA^fMet. Here we present a cryo-EM reconstruction of the complete 30S IC, containing mRNA, fMet-tRNA^fMet, IF1, IF2, and IF3. In the 30S IC, IF2 contacts IF1, the 30S subunit shoulder, and the CCA end of fMet-tRNA^fMet, which occupies a novel P/I position (P/I1). The N-terminal domain of IF3 contacts the tRNA, whereas the C-terminal domain is bound to the platform of the 30S subunit. Binding of initiation factors and fMet-tRNA^fMet induces a rotation of the head relative to the body of the 30S subunit, which is likely to prevail through 50S subunit joining until GTP hydrolysis and dissociation of IF2 take place. The structure provides insights into the mechanism of mRNA selection during translation initiation.

Translation is the process by which a ribosome converts the sequence of a messenger RNA (mRNA)—produced from a gene—into the sequence of amino acids that comprise a protein. Bacterial ribosomes each have one large and one small subunit: the 50S and 30S subunits. Initiation of translation entails selection of an mRNA, identification of the correct starting point from which to read its code, and engagement of the initial amino acid carrier (tRNA). These events take place in the 30S subunit and require the presence of three initiation factors (IF1, IF2, IF3). Formation of this 30S initiation complex precedes joining with the 50S subunit to assemble the functional ribosome. By using a cryo-electron microscopy approach to visualize the structures without fixation or staining, we have determined the structure of a complete 30S initiation complex and identified the positions and orientations of the tRNA and all three initiation factors. We found that the presence of the initiation factors and tRNA induces rotation of the head relative to the body of the 30S subunit, which may be essential for rapid binding to the 50S subunit and for regulating selection of the mRNA. IF3 had not been seen previously in the context of the 30S structure and its visualization gives insight into a potential role in preventing association of the two ribosomal subunits. These findings are important for understanding how the interplay of elements during the early stages of translation selects the mRNA and regulates formation of functional ribosomes.

Automated minimization of steric clashes in protein structures

Molecular modeling of proteins including homology modeling, structure determination, and knowledge-based protein design requires tools to evaluate and refine three-dimensional protein structures. Steric clash is one of the artifacts prevalent in low-resolution structures and homology models. Steric clashes arise due to the unnatural overlap of any two nonbonding atoms in a protein structure. Usually, removal of severe steric clashes in some structures is challenging since many existing refinement programs do not accept structures with severe steric clashes. Here, we present a quantitative approach of identifying steric clashes in proteins by defining clashes based on the Van der Waals repulsion energy of the clashing atoms. We also define a metric for quantitative estimation of the severity of clashes in proteins by performing statistical analysis of clashes in high-resolution protein structures. We describe a rapid, automated, and robust protocol, Chiron, which efficiently resolves severe clashes in low-resolution structures and homology models with minimal perturbation in the protein backbone. Benchmark studies highlight the efficiency and robustness of Chiron compared with other widely used methods. We provide Chiron as an automated web server to evaluate and resolve clashes in protein structures that can be further used for more accurate protein design.

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.

Accurate flexible fitting of high-resolution protein structures into cryo-electron microscopy maps using coarse-grained pseudo-energy minimization

Cryo-electron microscopy (cryo-EM) has been widely used to explore conformational states of large biomolecular assemblies. The detailed interpretation of cryo-EM data requires the flexible fitting of a known high-resolution protein structure into a low-resolution cryo-EM map. To this end, we have developed what we believe is a new method based on a two-bead-per-residue protein representation, and a modified form of the elastic network model that allows large-scale conformational changes while maintaining pseudobonds and secondary structures. Our method minimizes a pseudo-energy which linearly combines various terms of the modified elastic network model energy with a cryo-EM-fitting score and a collision energy that penalizes steric collisions. Unlike previous flexible fitting efforts using the lowest few normal modes, our method effectively utilizes all normal modes so that both global and local structural changes can be fully modeled. We have validated our method for a diverse set of 10 pairs of protein structures using simulated cryo-EM maps with a range of resolutions and in the absence/presence of random noise. We have shown that our method is both accurate and efficient compared with alternative techniques, and its performance is robust to the addition of random noise. Our method is also shown to be useful for the flexible fitting of three experimental cryo-EM maps.

Structure of the Drosophila apoptosome at 6.9 å resolution

The Drosophila Apaf-1 related killer forms an apoptosome in the intrinsic cell death pathway. In this study we show that Dark forms a single ring when initiator procaspases are bound. This Dark-Dronc complex cleaves DrICE efficiently; hence, a single ring represents the Drosophila apoptosome. We then determined the 3D structure of a double ring at ∼6.9 Å resolution and created a model of the apoptosome. Subunit interactions in the Dark complex are similar to those in Apaf-1 and CED-4 apoptosomes, but there are significant differences. In particular, Dark has "lost" a loop in the nucleotide-binding pocket, which opens a path for possible dATP exchange in the apoptosome. In addition, caspase recruitment domains (CARDs) form a crown on the central hub of the Dark apoptosome. This CARD geometry suggests that conformational changes will be required to form active Dark-Dronc complexes. When taken together, these data provide insights into apoptosome structure, function, and evolution.

Cryo-EM analysis reveals new insights into the mechanism of action of pyruvate carboxylase

Pyruvate carboxylase (PC) is a conserved multifunctional enzyme linked to important metabolic diseases. PC homotetramer is arranged in two layers with two opposing monomers per layer. Cryo-EM explores the conformational variability of PC in the presence of different substrates. The results demonstrate that the biotin-carboxyl carrier protein (BCCP) domain localizes near the biotin carboxylase (BC) domain of its own monomer and travels to the carboxyltransferase (CT) domain of the opposite monomer. All density maps show noticeable conformational differences between layers, mainly for the BCCP and BC domains. This asymmetry may be indicative of a coordination mechanism where monomers from different layers catalyze the BC and CT reactions consecutively. A conformational change of the PC tetramerization (PT) domain suggests a new functional role in communication. A long-range communication pathway between subunits in different layers, via interacting PT-PT and BC-BC domains, may be responsible for the cooperativity of PC from Staphylococcus aureus.

Exploring the spatial and temporal organization of a cell's proteome

To increase our current understanding of cellular processes, such as cell signaling and division, knowledge is needed about the spatial and temporal organization of the proteome at different organizational levels. These levels cover a wide range of length and time scales: from the atomic structures of macromolecules for inferring their molecular function, to the quantitative description of their abundance, and spatial distribution in the cell. Emerging new experimental technologies are greatly increasing the availability of such spatial information on the molecular organization in living cells. This review addresses three fields that have significantly contributed to our understanding of the proteome's spatial and temporal organization: first, methods for the structure determination of individual macromolecular assemblies, specifically the fitting of atomic structures into density maps generated from electron microscopy techniques; second, research that visualizes the spatial distributions of these complexes within the cellular context using cryo electron tomography techniques combined with computational image processing; and third, methods for the spatial modeling of the dynamic organization of the proteome, specifically those methods for simulating reaction and diffusion of proteins and complexes in crowded intracellular fluids. The long-term goal is to integrate the varied data about a proteome's organization into a spatially explicit, predictive model of cellular processes.

Determining macromolecular assembly structures by molecular docking and fitting into an electron density map

Structural models of macromolecular assemblies are instrumental for gaining a mechanistic understanding of cellular processes. Determining these structures is a major challenge for experimental techniques, such as X-ray crystallography, NMR spectroscopy and electron microscopy (EM). Thus, computational modeling techniques, including molecular docking, are required. The development of most molecular docking methods has so far been focused on modeling of binary complexes. We have recently introduced the MultiFit method for modeling the structure of a multisubunit complex by simultaneously optimizing the fit of the model into an EM density map of the entire complex and the shape complementarity between interacting subunits. Here, we report algorithmic advances of the MultiFit method that result in an efficient and accurate assembly of the input subunits into their density map. The successful predictions and the increasing number of complexes being characterized by EM suggests that the CAPRI challenge could be extended to include docking-based modeling of macromolecular assemblies guided by EM.

Modeling protein structure at near atomic resolutions with Gorgon

Electron cryo-microscopy (cryo-EM) has played an increasingly important role in elucidating the structure and function of macromolecular assemblies in near native solution conditions. Typically, however, only non-atomic resolution reconstructions have been obtained for these large complexes, necessitating computational tools for integrating and extracting structural details. With recent advances in cryo-EM, maps at near-atomic resolutions have been achieved for several macromolecular assemblies from which models have been manually constructed. In this work, we describe a new interactive modeling toolkit called Gorgon targeted at intermediate to near-atomic resolution density maps (10-3.5 Å), particularly from cryo-EM. Gorgon's de novo modeling procedure couples sequence-based secondary structure prediction with feature detection and geometric modeling techniques to generate initial protein backbone models. Beyond model building, Gorgon is an extensible interactive visualization platform with a variety of computational tools for annotating a wide variety of 3D volumes. Examples from cryo-EM maps of Rotavirus and Rice Dwarf Virus are used to demonstrate its applicability to modeling protein structure.

Three-dimensional asymmetric reconstruction of tailed bacteriophage

A universal goal in studying the structures of macromolecules and macromolecular complexes by means of electron cryo-microscopy (cryo-TEM) and three-dimensional (3D) image reconstruction is the derivation of a reliable atomic or pseudoatomic model. Such a model provides the foundation for exploring in detail the mechanisms by which biomolecules function. Though a variety of highly ordered, symmetric specimens such as 2D crystals, helices, and icosahedral virus capsids have been studied by these methods at near-atomic resolution, until recently, numerous challenges have made it difficult to achieve sub-nanometer resolution with large (≥~500Å), asymmetric molecules such as the tailed bacteriophages. After briefly reviewing some of the history behind the development of asymmetric virus reconstructions, we use recent structural studies of the prolate phage ϕ29 as an example to illustrate the step-by-step procedures used to compute an asymmetric reconstruction at sub-nanometer resolution. In contrast to methods that have been employed to study other asymmetric complexes, we demonstrate how symmetries in the head and tail components of the phage can be exploited to obtain the structure of the entire phage in an expedited, stepwise process. Prospects for future enhancements to the procedures currently employed are noted in the concluding section.

Unified data resource for cryo-EM

Three-dimensional (3D) cryoelectron microscopy reconstruction methods are uniquely able to reveal structures of many important macromolecules and macromolecular complexes. EMDataBank.org, a joint effort of the Protein Databank in Europe (PDBe), the Research Collaboratory for Structural Bioinformatics (RCSB), and the National Center for Macromolecular Imaging (NCMI), is a "one-stop shop" resource for global deposition and retrieval of cryo-EM map, model, and associated metadata. The resource unifies public access to the two major EM Structural Data archives: EM Data Bank (EMDB) and Protein Data Bank (PDB), and facilitates use of EM structural data of macromolecules and macromolecular complexes by the wider scientific community.

Modeling of proteins and their assemblies with the integrative modeling platform

To understand the workings of the living cell, we need to characterize protein assemblies that constitute the cell (for example, the ribosome, 26S proteasome, and the nuclear pore complex). A reliable high-resolution structural characterization of these assemblies is frequently beyond the reach of current experimental methods, such as X-ray crystallography, NMR spectroscopy, electron microscopy, footprinting, chemical cross-linking, FRET spectroscopy, small-angle X-ray scattering, and proteomics. However, the information garnered from different methods can be combined and used to build computational models of the assembly structures that are consistent with all of the available datasets. Here, we describe a protocol for this integration, whereby the information is converted to a set of spatial restraints and a variety of optimization procedures can be used to generate models that satisfy the restraints as much as possible. These generated models can then potentially inform about the precision and accuracy of structure determination, the accuracy of the input datasets, and further data generation. We also demonstrate the Integrative Modeling Platform (IMP) software, which provides the necessary computational framework to implement this protocol, and several applications for specific-use cases.

Improved protein surface comparison and application to low-resolution protein structure data

Background

Recent advancements of experimental techniques for determining protein tertiary structures raise significant challenges for protein bioinformatics. With the number of known structures of unknown function expanding at a rapid pace, an urgent task is to provide reliable clues to their biological function on a large scale. Conventional approaches for structure comparison are not suitable for a real-time database search due to their slow speed. Moreover, a new challenge has arisen from recent techniques such as electron microscopy (EM), which provide low-resolution structure data. Previously, we have introduced a method for protein surface shape representation using the 3D Zernike descriptors (3DZDs). The 3DZD enables fast structure database searches, taking advantage of its rotation invariance and compact representation. The search results of protein surface represented with the 3DZD has showngood agreement with the existing structure classifications, but some discrepancies were also observed.

Results

The three new surface representations of backbone atoms, originally devised all-atom-surface representation, and the combination of all-atom surface with the backbone representation are examined. All representations are encoded with the 3DZD. Also, we have investigated the applicability of the 3DZD for searching protein EM density maps of varying resolutions. The surface representations are evaluated on structure retrieval using two existing classifications, SCOP and the CE-based classification.

Conclusions

Overall, the 3DZDs representing backbone atoms show better retrieval performance than the original all-atom surface representation. The performance further improved when the two representations are combined. Moreover, we observed that the 3DZD is also powerful in comparing low-resolution structures obtained by electron microscopy.

RNA channelling by the eukaryotic exosome

The eukaryotic exosome is a key nuclease for the degradation, processing and quality control of a wide variety of RNAs. Here, we report electron microscopic reconstructions and pseudo-atomic models of the ten-subunit Saccharomyces cerevisiae exosome in the unbound and RNA-bound states. In the RNA-bound structures, extra density that is visible at the entry and exit sites of the exosome channel indicates that a substrate-threading mechanism is used by the eukaryotic exosome. This channelling mechanism seems to be conserved in exosome-like complexes from all domains of life, and might have been present in the most recent common ancestor.

Template-free 13-protofilament microtubule-MAP assembly visualized at 8 A resolution

Microtubule-associated proteins (MAPs) are essential for regulating and organizing cellular microtubules (MTs). However, our mechanistic understanding of MAP function is limited by a lack of detailed structural information. Using cryo-electron microscopy and single particle algorithms, we solved the 8 Å structure of doublecortin (DCX)-stabilized MTs. Because of DCX's unusual ability to specifically nucleate and stabilize 13-protofilament MTs, our reconstruction provides unprecedented insight into the structure of MTs with an in vivo architecture, and in the absence of a stabilizing drug. DCX specifically recognizes the corner of four tubulin dimers, a binding mode ideally suited to stabilizing both lateral and longitudinal lattice contacts. A striking consequence of this is that DCX does not bind the MT seam. DCX binding on the MT surface indirectly stabilizes conserved tubulin-tubulin lateral contacts in the MT lumen, operating independently of the nucleotide bound to tubulin. DCX's exquisite binding selectivity uncovers important insights into regulation of cellular MTs.

Molecular rearrangements involved in the capsid shell maturation of bacteriophage T7

Maturation of dsDNA bacteriophages involves assembling the virus prohead from a limited set of structural components followed by rearrangements required for the stability that is necessary for infecting a host under challenging environmental conditions. Here, we determine the mature capsid structure of T7 at 1 nm resolution by cryo-electron microscopy and compare it with the prohead to reveal the molecular basis of T7 shell maturation. The mature capsid presents an expanded and thinner shell, with a drastic rearrangement of the major protein monomers that increases in their interacting surfaces, in turn resulting in a new bonding lattice. The rearrangements include tilting, in-plane rotation, and radial expansion of the subunits, as well as a relative bending of the A- and P-domains of each subunit. The unique features of this shell transformation, which does not employ the accessory proteins, inserted domains, or molecular interactions observed in other phages, suggest a simple capsid assembling strategy that may have appeared early in the evolution of these viruses.

EMDataBank.org: unified data resource for CryoEM

Cryo-electron microscopy reconstruction methods are uniquely able to reveal structures of many important macromolecules and macromolecular complexes. EMDataBank.org, a joint effort of the Protein Data Bank in Europe (PDBe), the Research Collaboratory for Structural Bioinformatics (RCSB) and the National Center for Macromolecular Imaging (NCMI), is a global ‘one-stop shop’ resource for deposition and retrieval of cryoEM maps, models and associated metadata. The resource unifies public access to the two major archives containing EM-based structural data: EM Data Bank (EMDB) and Protein Data Bank (PDB), and facilitates use of EM structural data of macromolecules and macromolecular complexes by the wider scientific community.

Direct visualization of secondary structures of F-actin by electron cryomicroscopy

F-actin is a helical assembly of actin, which is a component of muscle fibres essential for contraction and has a crucial role in numerous cellular processes, such as the formation of lamellipodia and filopodia, as the most abundant component and regulator of cytoskeletons by dynamic assembly and disassembly (from G-actin to F-actin and vice versa). Actin is a ubiquitous protein and is involved in important biological functions, but the definitive high-resolution structure of F-actin remains unknown. Although a recent atomic model well reproduced X-ray fibre diffraction intensity data from a highly oriented liquid-crystalline sol specimen, its refinement without experimental phase information has certain limitations. Direct visualization of the structure by electron cryomicroscopy, however, has been difficult because it is relatively thin and flexible. Here we report the F-actin structure at 6.6 Å resolution, made obtainable by recent advances in electron cryomicroscopy. The density map clearly resolves all the secondary structures of G-actin, such as α-helices, β-structures and loops, and makes unambiguous modelling and refinement possible. Complex domain motions that open the nucleotide-binding pocket on F-actin formation, specific D-loop and terminal conformations, and relatively tight axial but markedly loose interprotofilament interactions hydrophilic in nature are revealed in the F-actin model, and all seem to be important for dynamic functions of actin.

A fast mathematical programming procedure for simultaneous fitting of assembly components into cryoEM density maps

MOTIVATION

Single-particle cryo electron microscopy (cryoEM) typically produces density maps of macromolecular assemblies at intermediate to low resolution (approximately 5-30 A). By fitting high-resolution structures of assembly components into these maps, pseudo-atomic models can be obtained. Optimizing the quality-of-fit of all components simultaneously is challenging due to the large search space that makes the exhaustive search over all possible component configurations computationally unfeasible.

RESULTS

We developed an efficient mathematical programming algorithm that simultaneously fits all component structures into an assembly density map. The fitting is formulated as a point set matching problem involving several point sets that represent component and assembly densities at a reduced complexity level. In contrast to other point matching algorithms, our algorithm is able to match multiple point sets simultaneously and not only based on their geometrical equivalence, but also based on the similarity of the density in the immediate point neighborhood. In addition, we present an efficient refinement method based on the Iterative Closest Point registration algorithm. The integer quadratic programming method generates an assembly configuration in a few seconds. This efficiency allows the generation of an ensemble of candidate solutions that can be assessed by an independent scoring function. We benchmarked the method using simulated density maps of 11 protein assemblies at 20 A, and an experimental cryoEM map at 23.5 A resolution. Our method was able to generate assembly structures with root-mean-square errors <6.5 A, which have been further reduced to <1.8 A by the local refinement procedure.

AVAILABILITY

The program is available upon request as a Matlab code package.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics Online.

Integrative structure modeling of macromolecular assemblies from proteomics data

Proteomics techniques have been used to generate comprehensive lists of protein interactions in a number of species. However, relatively little is known about how these interactions result in functional multiprotein complexes. This gap can be bridged by combining data from proteomics experiments with data from established structure determination techniques. Correspondingly, integrative computational methods are being developed to provide descriptions of protein complexes at varying levels of accuracy and resolution, ranging from complex compositions to detailed atomic structures.

Mechanism of eIF6-mediated inhibition of ribosomal subunit joining

During the process of ribosomal assembly, the essential eukaryotic translation initiation factor 6 (eIF6) is known to act as a ribosomal anti-association factor. However, a molecular understanding of the anti-association activity of eIF6 is still missing. Here we present the cryo-electron microscopy reconstruction of a complex of the large ribosomal subunit with eukaryotic eIF6 from Saccharomyces cerevisiae. The structure reveals that the eIF6 binding site involves mainly rpL23 (L14p in Escherichia coli). Based on our structural data, we propose that the mechanism of the anti-association activity of eIF6 is based on steric hindrance of intersubunit bridge formation around the dynamic bridge B6.

Structure of an apoptosome-procaspase-9 CARD complex

Apaf-1 coassembles with cytochrome c to form the apoptosome, which then binds and activates procaspase-9 (pc-9). We removed pc-9 catalytic domains from the holoapoptosome by site-directed thrombinolysis. A structure of the resulting apoptosome-pc-9 CARD complex was then determined at approximately 9.5 A resolution. In our model, the central hub is constructed like other AAA+ protein rings but also contains novel features. At higher radius, the regulatory region of each Apaf-1 is comprised of tandem seven and eight blade beta-propellers with cytochrome c docked between them. Remarkably, Apaf-1 CARDs are disordered in the ground state. During activation, each Apaf-1 CARD interacts with a pc-9 CARD and these heterodimers form a flexibly tethered "disk" that sits above the central hub. When taken together, the data reveal conformational changes during Apaf-1 assembly that allow pc-9 activation. The model also provides a plausible explanation for the effects of NOD mutations that have been mapped onto the central hub.

Building and refining protein models within cryo-electron microscopy density maps based on homology modeling and multiscale structure refinement

Automatic modeling methods using cryoelectron microscopy (cryoEM) density maps as constraints are promising approaches to building atomic models of individual proteins or protein domains. However, their application to large macromolecular assemblies has not been possible largely due to computational limitations inherent to such unsupervised methods. Here we describe a new method, EM-IMO (electron microscopy-iterative modular optimization), for building, modifying and refining local structures of protein models using cryoEM maps as a constraint. As a supervised refinement method, EM-IMO allows users to specify parameters derived from inspections so as to guide, and as a consequence, significantly speed up the refinement. An EM-IMO-based refinement protocol is first benchmarked on a data set of 50 homology models using simulated density maps. A multiscale refinement strategy that combines EM-IMO-based and molecular dynamics-based refinement is then applied to build backbone models for the seven conformers of the five capsid proteins in our near-atomic-resolution cryoEM map of the grass carp reovirus virion, a member of the Aquareovirus genus of the Reoviridae family. The refined models allow us to reconstruct a backbone model of the entire grass carp reovirus capsid and provide valuable functional insights that are described in the accompanying publication [Cheng, L., Zhu, J., Hui, W. H., Zhang, X., Honig, B., Fang, Q. & Zhou, Z. H. (2010). Backbone model of an aquareovirus virion by cryo-electron microscopy and bioinformatics. J. Mol. Biol. (this issue). doi:10.1016/j.jmb.2009.12.027.]. Our study demonstrates that the integrated use of homology modeling and a multiscale refinement protocol that combines supervised and automated structure refinement offers a practical strategy for building atomic models based on medium- to high-resolution cryoEM density maps.

Analyses of subnanometer resolution cryo-EM density maps

Today, electron cryomicroscopy (cryo-EM) can routinely achieve subnanometer resolutions of complex macromolecular assemblies. From a density map, one can extract key structural and functional information using a variety of computational analysis tools. At subnanometer resolution, these tools make it possible to isolate individual subunits, identify secondary structures, and accurately fit atomic models. With several cryo-EM studies achieving resolutions beyond 5Å, computational modeling and feature recognition tools have been employed to construct backbone and atomic models of the protein components directly from a density map. In this chapter, we describe several common classes of computational tools that can be used to analyze and model subnanometer resolution reconstructions from cryo-EM. A general protocol for analyzing subnanometer resolution density maps is presented along with a full description of steps used in analyzing the 4.3Å resolution structure of Mm-cpn.

From electron microscopy maps to atomic structures using normal mode-based fitting

Electron microscopy (EM) has made possible to solve the structure of many proteins. However, the resolution of some of the EM maps is too low for interpretation at the atomic level, which is particularly important to describe function. We describe methods that combine low-resolution EM data with atomic structures for different conformations of the same protein in order to produce atomic models compatible with the EM map.We illustrate these methods with EM data from decavanadate-induced tubular crystals of a pseudo-phosphorylated intermediate of Ca-ATPase and the various atomic structures of other intermediates available in the Protein Data Bank (PDB). Determination of atomic structure permits not only to analyse protein-protein interactions in the crystals, but also to localize residues in the proximity of the crystallizing agent both within Ca-ATPase and between Ca-ATPase molecules.

Semi-isometric registration of line features for flexible fitting of protein structures

In this paper, we study a registration problem that is motivated by a practical biology problem - fitting protein structures to low-resolution density maps. We consider registration between two sets of lines features (e.g., helices in the proteins) that have undergone not a single, but multiple isometric transformations (e.g., hinge-motions). The problem is further complicated by the presence of symmetry in each set. We formulate the problem as a clique-finding problem in a product graph, and propose a heuristic solution that includes a fast clique-finding algorithm unique to the structure of this graph. When tested on a suite of real protein structures, the algorithm achieved high accuracy even for very large inputs containing hundreds of helices.

Comprehensive molecular structure of the eukaryotic ribosome

Despite the emergence of a large number of X-ray crystallographic models of the bacterial 70S ribosome over the past decade, an accurate atomic model of the eukaryotic 80S ribosome is still not available. Eukaryotic ribosomes possess more ribosomal proteins and ribosomal RNA than do bacterial ribosomes, which are implicated in extraribosomal functions in the eukaryotic cells. By combining cryo-EM with RNA and protein homology modeling, we obtained an atomic model of the yeast 80S ribosome complete with all ribosomal RNA expansion segments and all ribosomal proteins for which a structural homolog can be identified. Mutation or deletion of 80S ribosomal proteins can abrogate maturation of the ribosome, leading to several human diseases. We have localized one such protein unique to eukaryotes, rpS19e, whose mutations are associated with Diamond-Blackfan anemia in humans. Additionally, we characterize crucial interactions between the dynamic stalk base of the ribosome with eukaryotic elongation factor 2.

EM-fold: De novo folding of alpha-helical proteins guided by intermediate-resolution electron microscopy density maps

In medium-resolution (7-10 A) cryo-electron microscopy (cryo-EM) density maps, alpha helices can be identified as density rods whereas beta-strand or loop regions are not as easily discerned. We are proposing a computational protein structure prediction algorithm "EM-Fold" that resolves the density rod connectivity ambiguity by placing predicted alpha helices into the density rods and adding missing backbone coordinates in loop regions. In a benchmark of 11 mainly alpha-helical proteins of known structure a native-like model is identified in eight cases (rmsd 3.9-7.9 A). The three failures can be attributed to inaccuracies in the secondary structure prediction step that precedes EM-Fold. EM-Fold has been applied to the approximately 6 A resolution cryo-EM density map of protein IIIa from human adenovirus. We report the first topological model for the alpha-helical 400 residue N-terminal region of protein IIIa. EM-Fold also has the potential to interpret medium-resolution density maps in X-ray crystallography.

Native secondary structure topology has near minimum contact energy among all possible geometrically constrained topologies

Secondary structure topology in this article refers to the order and the direction of the secondary structures, such as helices and strands, with respect to the protein sequence. Even when the locations of the secondary structure Calpha atoms are known, there are still (N!2(N))(M!2(M)) different possible topologies for a protein with N helices and M strands. This work explored the question if the native topology is likely to be identified among a large set of all possible geometrically constrained topologies through an evaluation of the residue contact energy formed by the secondary structures, instead of the entire chain. We developed a contact pair specific and distance specific multiwell function based on the statistical characterization of the side chain distances of 413 proteins in the Protein Data Bank. The multiwell function has specific parameters to each of the 210 pairs of residue contacts. We illustrated a general mathematical method to extend a single well function to a multiwell function to represent the statistical data. We have performed a mutation analysis using 50 proteins to generate all the possible geometrically constrained topologies of the secondary structures. The result shows that the native topology is within the top 25% of the list ranked by the effective contact energies of the secondary structures for all the 50 proteins, and is within the top 5% for 34 proteins. As an application, the method was used to derive the structure of the skeletons from a low resolution density map that can be obtained through electron cryomicroscopy.

Molecular dynamics flexible fitting: a practical guide to combine cryo-electron microscopy and X-ray crystallography

Hybrid computational methods for combining structural data from different sources and resolutions are becoming an essential part of structural biology, especially as the field moves toward the study of large macromolecular assemblies. We have developed the molecular dynamics flexible fitting (MDFF) method for combining high-resolution atomic structures with cryo-electron microscopy (cryo-EM) maps, that results in atomic models representing the conformational state captured by cryo-EM. The method has been applied successfully to the ribosome, a ribonucleoprotein complex responsible for protein synthesis. MDFF involves a molecular dynamics simulation in which a guiding potential, based on the cryo-EM map, is added to the standard force field. Forces proportional to the gradient of the density map guide an atomic structure, available from X-ray crystallography, into high-density regions of a cryo-EM map. In this paper we describe the necessary steps to set up, run, and analyze MDFF simulations and the software packages that implement the corresponding functionalities.

The P22 tail machine at subnanometer resolution reveals the architecture of an infection conduit

The portal channel is a key component in the life cycle of bacteriophages and herpesviruses. The bacteriophage P22 portal is a 1 megadalton dodecameric oligomer of gp1 that plays key roles in capsid assembly, DNA packaging, assembly of the infection machinery, and DNA ejection. The portal is the nucleation site for the assembly of 39 additional subunits generated from multiple copies of four gene products (gp4, gp10, gp9, and gp26), which together form the multifunctional tail machine. These components are organized with a combination of 12-fold (gp1, gp4), 6-fold (gp10, trimers of gp9), and 3-fold (gp26, gp9) symmetry. Here we present the 3-dimensional structures of the P22 assembly-naive portal formed from expressed subunits (gp1) and the intact tail machine purified from infectious virions. The assembly-naive portal structure exhibits a striking structural similarity to the structures of the portal proteins of SPP1 and phi29 derived from X-ray crystallography.

An atomic model AAA-ATPase/20S core particle sub-complex of the 26S proteasome

The 26S proteasome is the most downstream element of the ubiquitin-proteasome pathway of protein degradation. It is composed of the 20S core particle (CP) and the 19S regulatory particle (RP). The RP consists of 6 AAA-ATPases and at least 13 non-ATPase subunits. Based on a cryo-EM map of the 26S proteasome, structures of homologs, and physical protein-protein interactions we derive an atomic model of the AAA-ATPase-CP sub-complex. The ATPase order in our model (Rpt1/Rpt2/Rpt6/Rpt3/Rpt4/Rpt5) is in excellent agreement with the recently identified base-precursor complexes formed during the assembly of the RP. Furthermore, the atomic CP-AAA-ATPase model suggests that the assembly chaperone Nas6 facilitates CP-RP association by enhancing the shape complementarity between Rpt3 and its binding CP alpha subunits partners.

Single-particle reconstruction of biological macromolecules in electron microscopy--30 years

This essay gives the autho's personal account on the development of concepts underlying single-particle reconstruction, a technique in electron microscopy of macromolecular assemblies with a remarkable record of achievements as of late. The ribosome proved to be an ideal testing ground for the development of specimen preparation methods, cryo-EM techniques, and algorithms, with discoveries along the way as a rich reward. Increasingly, cryo-EM and single-particle reconstruction, in combination with classification techniques, is revealing dynamic information on functional molecular machines uninhibited by molecular contacts.

The advent of near-atomic resolution in single-particle electron microscopy

Single-particle electron microscopy (EM) can provide structural information for a large variety of biological molecules, ranging from small proteins to large macromolecular assemblies, without the need to produce crystals. The year 2008 has become a landmark year for single-particle EM as for the first time density maps have been produced at a resolution that made it possible to trace protein backbones or even to build atomic models. In this review, we highlight some of the recent successes achieved by single-particle EM and describe the individual steps involved in producing a density map by this technique. We also discuss some of the remaining challenges and areas, in which further advances would have a great impact on the results that can be achieved by single-particle EM.