Electron cryomicroscopy and bioinformatics suggest protein fold models for rice dwarf virus

Ab initio resolution measurement for single particle structures

A computational method is described that allows the measurement of the signal-to-noise ratio and resolution of a three-dimensional structure obtained by single particle electron microscopy and reconstruction. The method does not rely on the availability of the original image data or the calculation of several structures from different parts of the data that are needed for the commonly used Fourier Shell Correlation criterion. Instead, the correlation between neighboring Fourier pixels is calculated and used to distinguish signal from noise. The new method has been conveniently implemented in a computer program called RMEASURE and is available to the microscopy community.

Ab initio random model method facilitates 3D reconstruction of icosahedral particles

Model-based, three-dimensional (3D) image reconstruction procedures require a starting model to initiate data analysis. We have designed an ab initio method, which we call the random model (RM) method, that automatically generates models to initiate structural analysis of icosahedral viruses imaged by cryo-electron microscopy. The robustness of the RM procedure was demonstrated on experimental sets of images for five representative viruses. The RM method also provides a straightforward way to generate unbiased starting models to derive independent 3D reconstructions and obtain a more reliable assessment of resolution. The fundamental scheme embodied in the RM method should be relatively easy to integrate into other icosahedral software packages.

Automated segmentation of molecular subunits in electron cryomicroscopy density maps

Electron cryomicroscopy (cryoEM) is capable of imaging large macromolecular machines composed of multiple components. However, it is currently only possible to achieve moderate resolution at which it may be possible to computationally extract the individual components in the machine. In this work, we present application details of an automated method for detecting and segmenting the components of a large machine in an experimentally determined density map. This method is applicable to object with and without symmetry and takes advantage of global and local symmetry axes if present. We have applied this segmentation algorithm to several cryoEM data sets already deposited in EMDB with various complexities, symmetries and resolutions and validated the results using manually segmented density and available structures of the components in the PDB. As such, automated segmentation could become a useful tool for the analysis of the ever-increasing number of structures of macromolecular machines derived from cryoEM.

Refinement of Protein Structures by Iterative Comparative Modeling and CryoEM Density Fitting

We developed a method for structure characterization of assembly components by iterative comparative protein structure modeling and fitting into cryo-electron microscopy (cryoEM) density maps. Specifically, we calculate a comparative model of a given component by considering many alternative alignments between the target sequence and a related template structure while optimizing the fit of a model into the corresponding density map. The method relies on the previously developed Moulder protocol that iterates over alignment, model building, and model assessment. The protocol was benchmarked using 20 varied target-template pairs of known structures with less than 30% sequence identity and corresponding simulated density maps at resolutions from 5A to 25A. Relative to the models based on the best existing sequence profile alignment methods, the percentage of C(alpha) atoms that are within 5A of the corresponding C(alpha) atoms in the superposed native structure increases on average from 52% to 66%, which is half-way between the starting models and the models from the best possible alignments (82%). The test also reveals that despite the improvements in the accuracy of the fitness function, this function is still the bottleneck in reducing the remaining errors. To demonstrate the usefulness of the protocol, we applied it to the upper domain of the P8 capsid protein of rice dwarf virus that has been studied by cryoEM at 6.8A. The C(alpha) root-mean-square deviation of the model based on the remotely related template, bluetongue virus VP7, improved from 8.7A to 6.0A, while the best possible model has a C(alpha) RMSD value of 5.3A. Moreover, the resulting model fits better into the cryoEM density map than the initial template structure. The method is being implemented in our program MODELLER for protein structure modeling by satisfaction of spatial restraints and will be applicable to the rapidly increasing number of cryoEM density maps of macromolecular assemblies.

Structural biology of cellular machines

Multi-component macromolecular machines contribute to all essential biological processes, from cell motility and signal transduction to information storage and processing. Structural analysis of assemblies at atomic resolution is emerging as the field of structural cell biology. Several recent studies, including those focused on the ribosome, the acrosomal bundle and bacterial flagella, have demonstrated the ability of a hybrid approach that combines imaging, crystallography and computational tools to generate testable atomic models of fundamental biological machines. A complete understanding of cellular and systems biology will require the detailed structural understanding of hundreds of biological machines. The realization of this goal demands a concerted effort to develop and apply new strategies for the systematic identification, isolation, structural characterization and mechanistic analysis of multi-component assemblies at all resolution ranges. The establishment of a database describing the structural and dynamic properties of protein assemblies will provide novel opportunities to define the molecular and atomic mechanisms controlling overall cell physiology.

Electron cryomicroscopy of single particles at subnanometer resolution

Electron cryomicroscopy and single-particle reconstruction have advanced substantially over the past two decades. There are now numerous examples of structures that have been solved using this technique to better than 10 A resolution. At such resolutions, direct identification of alpha helices is possible and, often, beta-sheet-containing regions can be identified. The most numerous subnanometer resolution structures are the icosahedral viruses, as higher resolution is easier to achieve with higher symmetry. Important non-icosahedral structures solved to subnanometer resolution include several ribosome structures, clathrin assemblies and, most recently, the Ca2+ release channel. There is now hope that, in the next few years, this technique will achieve resolutions approaching 4 A, permitting a complete trace of the protein backbone without reference to a crystal structure.

Combining electron microscopy and comparative protein structure modeling

Recently, advances have been made in methods and applications that integrate electron microscopy density maps and comparative modeling to produce atomic structures of macromolecular assemblies. Electron microscopy can benefit from comparative modeling through the fitting of comparative models into electron microscopy density maps. Also, comparative modeling can benefit from electron microscopy through the use of intermediate-resolution density maps in fold recognition, template selection and sequence-structure alignment.

Features of reovirus outer capsid protein mu1 revealed by electron cryomicroscopy and image reconstruction of the virion at 7.0 Angstrom resolution

Reovirus is a useful model for addressing the molecular basis of membrane penetration by one of the larger nonenveloped animal viruses. We now report the structure of the reovirus virion at approximately 7.0 A resolution as obtained by electron cryomicroscopy and three-dimensional image reconstruction. Several features of the myristoylated outer capsid protein mu1, not seen in a previous X-ray crystal structure of the mu1-sigma3 heterohexamer, are evident in the virion. These features appear to be important for stabilizing the outer capsid, regulating the conformational changes in mu1 that accompany perforation of target membranes, and contributing directly to membrane penetration during cell entry.

Fast rotational matching of single-particle images

The presence of noise and absence of contrast in electron micrographs lead to a reduced resolution of the final 3D reconstruction, due to the inherent limitations of single-particle image alignment. The fast rotational matching (FRM) algorithm was introduced recently for an accurate alignment of 2D images under such challenging conditions. Here, we implemented this algorithm for the first time in a standard 3D reconstruction package used in electron microscopy. This allowed us to carry out exhaustive tests of the robustness and reliability in iterative orientation determination, classification, and 3D reconstruction on simulated and experimental image data. A classification test on GroEL chaperonin images demonstrates that FRM assigns up to 13% more images to their correct reference orientation, compared to the classical self-correlation function method. Moreover, at sub-nanometer resolution, GroEL and rice dwarf virus reconstructions exhibit a remarkable resolution gain of 10-20% that is attributed to the novel image alignment kernel.

Automatic ultrastructure segmentation of reconstructed cryoEM maps of icosahedral viruses

We present an automatic algorithm to segment all the local and global asymmetric units of a three-dimensional density map of icosahedral viruses. This approach is readily applicable to the structural analysis of a broad range of virus structures that are reconstructed using cryo-electron microscopy (cryo-EM) technique. Our algorithm includes three major steps operating on the three dimensional density map: the detection of critical points of the volumetric density function, the detection of global and local symmetry axes, and, finally, the boundary segmentation of all the asymmetric units. We demonstrate the efficacy of our algorithm and report our results on several experimental volumetric datasets, consisting of both reconstructed cryo-EM molecular density maps taken from the European Bioinformatics Institute archive, as well our own synthetically generated (blurred) maps calculated from X-ray resolution molecular structural data taken from the Protein Data Bank.

Compressed representations of macromolecular structures and properties

We introduce a new and unified, compressed volumetric representation for macromolecular structures at varying feature resolutions, as well as for many computed associated properties. Important caveats of this compressed representation are fast random data access and decompression operations. Many computational tasks for manipulating large structures, including those requiring interactivity such as real-time visualization, are greatly enhanced by utilizing this compact representation. The compression scheme is obtained by using a custom designed hierarchical wavelet basis construction. Due to the continuity offered by these wavelets, we retain very good accuracy of molecular surfaces, at very high compression ratios, for macromolecular structures at multiple resolutions.

Determining Protein Topology from Skeletons of Secondary Structures

We report a novel computational procedure for determining protein native topology, or fold, by defining loop connectivity based on skeletons of secondary structures that can usually be obtained from low to intermediate-resolution density maps. The procedure primarily involves a knowledge-based geometry filter followed by an energetics-based evaluation. It was tested on a large set of skeletons covering a wide range of protein architecture, including one modeled from an experimentally determined 7.6A cryo-electron microscopy (cryo-EM) density map. The results showed that the new procedure could effectively deduce protein folds without high-resolution structural data, a feature that could also be used to recognize native fold in structure prediction and to interpret data in fields like structure genomics. Most importantly, in the energetics-based evaluation, it was revealed that, despite the inevitable errors in the artificially constructed structures and limited accuracy of knowledge-based potential functions, the average energy of an ensemble of structures with slightly different configurations around the native skeleton is a much more robust parameter for marking native topology than the energy of individual structures in the ensemble. This result implies that, among all the possible topology candidates for a given skeleton, evolution has selected the native topology as the one that can accommodate the largest structural variations, not the one rigidly trapped in a deep, but narrow, conformational energy well.

The structural basis of recognition and removal of cellular mRNA 7-methyl G 'caps' by a viral capsid protein: a unique viral response to host defense

The single segment, double-stranded RNA genome of the L-A virus (L-A) of yeast encodes two proteins: the major coat protein Gag (76 kDa) and the Gag-Pol fusion protein (180 kDa). The icosahedral L-A capsid is formed by 120 copies of Gag and has architecture similar to that seen in the reovirus, blue tongue virus and rice dwarf virus inner protein shells. Gag chemically removes the m7GMP caps from host cellular mRNAs. Previously we identified a trench on the outer surface of Gag that included His154, to which caps are covalently attached. Here we report the refined L-A coordinates at 3.4 angstroms resolution with additional structural features and the structure of L-A with bound m7GDP at 6.5 angstroms resolution, which shows the conformational change of the virus upon ligand binding. Based on site-directed mutations, residues in or adjacent to the trench that are essential (or dispensable) for the decapping reaction are described here. Along with His154, the reaction requires a cluster of positive charge adjoining the trench and residues Tyr 452, Tyr150 and either Tyr or Phe at position 538. A tentative mechanism for decapping is proposed.

Electron cryomicroscopy of biological machines at subnanometer resolution

Advances in electron cryomicroscopy (cryo-EM) have made possible the structural determination of large biological machines in the resolution range of 6-9 angstroms. Rice dwarf virus and the acrosomal bundle represent two distinct types of machines amenable to cryo-EM investigations at subnanometer resolutions. However, calculating the density map is only the first step, and much analysis remains to extract structural insights and the mechanism of action in these machines. This paper will review the computational and visualization methodologies necessary for analysis (structure mining) of the computed cryo-EM maps of these machines. These steps include component segmentation, averaging based on local symmetry among components, density connectivity trace, incorporation of bioinformatics analysis, and fitting of high-resolution component data, if available. The consequences of these analyses can not only identify accurately some of the secondary structure elements of the molecular components in machines but also suggest structural mechanisms related to their biological functions.

Seeing GroEL at 6 A resolution by single particle electron cryomicroscopy

We present a reconstruction of native GroEL by electron cryomicroscopy (cryo-EM) and single particle analysis at 6 A resolution. alpha helices are clearly visible and beta sheet density is also visible at this resolution. While the overall conformation of this structure is quite consistent with the published X-ray data, a measurable shift in the positions of three alpha helices in the intermediate domain is observed, not consistent with any of the 7 monomeric structures in the Protein Data Bank model (1OEL). In addition, there is evidence for slight rearrangement or flexibility in parts of the apical domain. The 6 A resolution cryo-EM GroEL structure clearly demonstrates the veracity and expanding scope of cryo-EM and the single particle reconstruction technique for macromolecular machines.

Structural characterization of components of protein assemblies by comparative modeling and electron cryo-microscopy

We explore structural characterization of protein assemblies by a combination of electron cryo-microscopy (cryoEM) and comparative protein structure modeling. Specifically, our method finds an optimal atomic model of a given assembly subunit and its position within an assembly by fitting alternative comparative models into a cryoEM map. The alternative models are calculated by MODELLER [J. Mol. Biol. 234 (1993) 313] from different sequence alignments between the modeled protein and its template structures. The fitting of these models into a cryoEM density map is performed either by FOLDHUNTER [J. Mol. Biol. 308 (2001) 1033] or by a new density fitting module of MODELLER (Mod-EM). Identification of the most accurate model is based on the correlation between the model accuracy and the quality of fit into the cryoEM density map. To quantify this correlation, we created a benchmark consisting of eight proteins of different structural folds with corresponding density maps simulated at five resolutions from 5 to 15 angstroms, with three noise levels each. Each of the proteins in the set was modeled based on 300 different alignments to their remotely related templates (12-32% sequence identity), spanning the range from entirely inaccurate to essentially accurate alignments. The benchmark revealed that one of the most accurate models can usually be identified by the quality of its fit into the cryoEM density map, even for noisy maps at 15 angstroms resolution. Therefore, a cryoEM density map can be helpful in improving the accuracy of a comparative model. Moreover, a pseudo-atomic model of a component in an assembly may be built better with comparative models of the native subunit sequences than with experimentally determined structures of their homologs.

Noise bias in the refinement of structures derived from single particles

One of the main goals in the determination of three-dimensional macromolecular structures from electron microscope images of individual molecules and complexes (single particles) is a sufficiently high spatial resolution, about 4 A, at which the interpretation with an atomic model becomes possible. To reach high resolution, an iterative refinement procedure using an expectation maximization algorithm is often used that leads to a more accurate alignment of the positional and orientational parameters for each particle. We show here the results of refinement algorithms that use a phase residual, a linear correlation coefficient, or a weighted correlation coefficient to align individual particles. The algorithms were applied to computer-generated data sets that contained projections from model structures, as well as noise. The algorithms show different degrees of over-fitting, especially at high resolution where the signal is weak. We demonstrate that the degree of over-fitting is reduced with a weighting scheme that depends on the signal-to-noise ratio in the data. The weighting also improves the accuracy of resolution measurement by the commonly used Fourier shell correlation. The performance of the refinement algorithms is compared to that using a maximum likelihood approach. The weighted correlation coefficient was implemented in the computer program FREALIGN.

A structural perspective on protein-protein interactions

Structures of macromolecular complexes are necessary for a mechanistic description of biochemical and cellular processes. They can be solved by experimental methods, such as X-ray crystallography, NMR spectroscopy and electron microscopy, as well as by computational protein structure prediction, docking and bioinformatics. Recent advances and applications of these methods emphasize the need for hybrid approaches that combine a variety of data to achieve better efficiency, accuracy, resolution and completeness.

UCSF Chimera--a visualization system for exploratory research and analysis

The design, implementation, and capabilities of an extensible visualization system, UCSF Chimera, are discussed. Chimera is segmented into a core that provides basic services and visualization, and extensions that provide most higher level functionality. This architecture ensures that the extension mechanism satisfies the demands of outside developers who wish to incorporate new features. Two unusual extensions are presented: Multiscale, which adds the ability to visualize large-scale molecular assemblies such as viral coats, and Collaboratory, which allows researchers to share a Chimera session interactively despite being at separate locales. Other extensions include Multalign Viewer, for showing multiple sequence alignments and associated structures; ViewDock, for screening docked ligand orientations; Movie, for replaying molecular dynamics trajectories; and Volume Viewer, for display and analysis of volumetric data. A discussion of the usage of Chimera in real-world situations is given, along with anticipated future directions. Chimera includes full user documentation, is free to academic and nonprofit users, and is available for Microsoft Windows, Linux, Apple Mac OS X, SGI IRIX, and HP Tru64 Unix from http://www.cgl.ucsf.edu/chimera/.

Interactions between the inner and outer capsids of bluetongue virus

Bluetongue virus is a large and structurally complex virus composed of three concentric capsid layers that surround 10 segments of a double-stranded RNA genome. X-ray crystallographic analysis of the particles without the outer capsid layer has provided atomic structural details of VP3 and VP7, which form the inner two layers. However, limited structural information is available on the other five proteins in the virion-two of which are important for receptor recognition, hemagglutination, and membrane interaction-are in the outer layer, and the others, important for endogenous transcriptase activity are internal. Here we report the electron cryomicroscopy (cryo-EM) reconstruction of the mature particle, which shows that the outer layer has a unique non-T = 13 icosahedral organization consisting of two distinct triskelion and globular motifs interacting extensively with the underlying T = 13 layer. Comparative cryo-EM analysis of the recombinant corelike particles has shown that VP1 (viral polymerase) and VP4 (capping enzyme) together form a flower-shaped structure attached to the underside of VP3, directly beneath the fivefold axis. The structural data have been substantiated by biochemical studies demonstrating the interactions between the individual outer and inner capsid proteins.

Conformational changes in Sindbis virions resulting from exposure to low pH and interactions with cells suggest that cell penetration may occur at the cell surface in the absence of membrane fusion

Alphaviruses have the ability to induce cell-cell fusion after exposure to acid pH. This observation has served as an article of proof that these membrane-containing viruses infect cells by fusion of the virus membrane with a host cell membrane upon exposure to acid pH after incorporation into a cell endosome. We have investigated the requirements for the induction of virus-mediated, low pH-induced cell-cell fusion and cell-virus fusion. We have correlated the pH requirements for this process to structural changes they produce in the virus by electron cryo-microscopy. We found that exposure to acid pH was required to establish conditions for membrane fusion but that membrane fusion did not occur until return to neutral pH. Electron cryo-microscopy revealed dramatic changes in the structure of the virion as it was moved to acid pH and then returned to neutral pH. None of these treatments resulted in the disassembly of the virus protein icosahedral shell that is a requisite for the process of virus membrane-cell membrane fusion. The appearance of a prominent protruding structure upon exposure to acid pH and its disappearance upon return to neutral pH suggested that the production of a "pore"-like structure at the fivefold axis may facilitate cell penetration as has been proposed for polio (J. Virol. 74 (2000) 1342) and human rhino virus (Mol. Cell 10 (2002) 317). This transient structural change also provided an explanation for how membrane fusion occurs after return to neutral pH. Examination of virus-cell complexes at neutral pH supported the contention that infection occurs at the cell surface at neutral pH by the production of a virus structure that breaches the plasma membrane bilayer. These data suggest an alternative route of infection for Sindbis virus that occurs by a process that does not involve membrane fusion and does not require disassembly of the virus protein shell.

Using surface envelopes for discrimination of molecular models

Shape information about macromolecules is increasingly available but is difficult to use in modeling efforts. We demonstrate that shape information alone can often distinguish structural models of biological macromolecules. By using a data structure called a surface envelope (SE) to represent the shape of the molecule, we propose a method that generates a fitness score for the shape of a particular molecular model. This score correlates well with root mean squared deviation (RMSD) of the model to the known test structures and can be used to filter models in decoy sets. The scoring method requires both alignment of the model to the SE in three-dimensional space and assessment of the degree to which atoms in the model fill the SE. Alignment combines a hybrid algorithm using principal components and a previously published iterated closest point algorithm. We test our method against models generated from random atom perturbation from crystal structures, published decoy sets used in structure prediction, and models created from the trajectories of atoms in molecular modeling runs. We also test our alignment algorithm against experimental electron microscopic data from rice dwarf virus. The alignment performance is reliable, and we show a high correlation between model RMSD and score function. This correlation is stronger for molecular models with greater oblong character (as measured by the ratio of largest to smallest principal component).

A Structural-informatics approach for tracing beta-sheets: building pseudo-C(alpha) traces for beta-strands in intermediate-resolution density maps

We report the development of two computational methods to assist density map interpretation at intermediate resolutions: sheettracer for building pseudo-C(alpha) models of beta-sheets, and a deconvolution method for enhancing features attributed to major secondary structural elements. Sheettracer is tightly coupled with sheetminer, which was developed to locate sheet densities in intermediate-resolution density maps. The results from sheetminer are used as inputs to sheettracer, which employs a multi-step ad hoc morphological analysis of sheet densities to trace individual strands of beta-sheets. The methods were tested on simulated density maps from 12 protein crystal structures that represent a reasonably complete sampling of sheet morphology. The sheet-tracing results were quantitatively assessed in terms of sensitivity, specificity and rms deviations. Furthermore, sheettracer and the deconvolution method were rigorously tested on experimental maps of the lambda2 protein of reovirus at resolutions of 7.6A and 11.8A. Our results clearly demonstrate the capability of sheettracer in building pseudo-C(alpha) models of beta-sheets in intermediate-resolution density maps and the power of the deconvolution method in enhancing the performance of sheettracer. These computational methods, along with other related ones, should facilitate recognition and analysis of folding motifs from experimental data at intermediate resolutions.

Three-Dimensional Electron Microscopy at Molecular Resolution

Emerging methods in cryo-electron microscopy allow determination of the three-dimensional architectures of objects ranging in size from small proteins to large eukaryotic cells, spanning a size range of more than 12 orders of magnitude. Advances in determining structures by "single particle" microscopy and by "electron tomography" provide exciting opportunities to describe the structures of subcellular assemblies that are either too large or too heterogeneous to be investigated by conventional crystallographic methods. Here, we review selected aspects of progress in structure determination by cryo-electron microscopy at molecular resolution, with a particular emphasis on topics at the interface of single particle and tomographic approaches. The rapid pace of development in this field suggests that comprehensive descriptions of the structures of whole cells and organelles in terms of the spatial arrangements of their molecular components may soon become routine.

The atomic structure of rice dwarf virus reveals the self-assembly mechanism of component proteins

Rice dwarf virus (RDV), the causal agent of rice dwarf disease, is a member of the genus Phytoreovirus in the family Reoviridae. RDV is a double-shelled virus with a molecular mass of approximately 70 million Dalton. This virus is widely prevalent and is one of the viruses that cause the most economic damage in many Asian countries. The atomic structure of RDV was determined at 3.5 A resolution by X-ray crystallography. The double-shelled structure consists of two different proteins, the core protein P3 and the outer shell protein P8. The atomic structure shows structural and electrostatic complementarities between both homologous (P3-P3 and P8-P8) and heterologous (P3-P8) interactions, as well as overall conformational changes found in P3-P3 dimer caused by the insertion of amino-terminal loop regions of one of the P3 protein into the other. These interactions suggest how the 900 protein components are built into a higher-ordered virus core structure.

Emerging themes in rotavirus cell entry, genome organization, transcription and replication

Rotaviruses, causative agents of gastroenteritis in young animals and humans, are large icosahedral viruses with a complex architecture. The double-stranded RNA (dsRNA) genome composed of 11 segments, which codes for 6 structural and 6 non-structural proteins, is enclosed within three concentric capsid layers. In addition to facilitating host-specific interactions, the design of the capsid architecture in rotaviruses as in other dsRNA viruses should also be conducive to the requirement of transcribing the enclosed genome segments repeatedly and simultaneously within the capsid interior. Several non-structural proteins facilitate the subsequent processes of genome replication and packaging. Electron cryomicroscopy studies of intact virions, recombinant virus-like particles, functional complexes, together with recent X-ray crystallographic studies on rotavirus proteins have provided structural insights into the capsid architecture, genome organization, antibody interaction, cell entry, trypsin-enhanced infectivity, endogenous transcription and replication. These studies underscore contrasting features and unifying themes between rotavirus and other dsRNA viruses.

Applications of a bilateral denoising filter in biological electron microscopy

Due to the sensitivity of biological sample to the radiation damage, the low dose imaging conditions used for electron microscopy result in extremely noisy images. The processes of digitization, image alignment, and 3D reconstruction also introduce additional sources of noise in the final 3D structure. In this paper, we investigate the effectiveness of a bilateral denoising filter in various biological electron microscopy applications. In contrast to the conventional low pass filters, which inevitably smooth out both noise and structural features simultaneously, we found that bilateral filter holds a distinct advantage in being capable of effectively suppressing noise without blurring the high resolution details. In as much, we have applied this technique to individual micrographs, entire 3D reconstructions, segmented proteins, and tomographic reconstructions.

Web portal to an image database for high-resolution three-dimensional reconstruction

The exponential increase of image data in high-resolution reconstructions by electron cryomicroscopy (cryoEM) has posed a need for efficient data management solutions in addition to powerful data processing procedures. Although relational databases and web portals are commonly used to manage sequences and structures in biological research, their application in cryoEM has been limited due to the complexity in accomplishing the dual tasks of interacting with proprietary software and simultaneously providing data access to users without database knowledge. Here, we report our results in developing web portal to SQL image databases used by the Image Management and Icosahedral Reconstruction System (IMIRS) to manage cryoEM images for subnanometer-resolution reconstructions. Fundamental issues related to the design and deployment of web portals to image databases are described. A web browser-based user interface was designed to accomplish data reporting and other database-related services, including user authentication, data entry, graph-based data mining, and various query and reporting tasks with interactive image manipulation capabilities. With an integrated web portal, IMIRS represents the first cryoEM application that incorporates both web-based data reporting tools and a complete set of data processing modules. Our examples should thus provide general guidelines applicable to other cryoEM technology development efforts.

FEMME database: topologic and geometric information of macromolecules

FEMME (Feature Extraction in a Multi-resolution Macromolecular Environment: http://www.biocomp.cnb.uam.es/FEMME/) database version 1.0 is a new bioinformatics data resource that collects topologic and geometric information obtained from macromolecular structures solved by three-dimensional electron microscopy (3D-EM). Although the FEMME database is focused on medium resolution data, the methodology employed (based on the so-called alpha-shape theory) is applicable to atomic resolution data as well. The alpha-shape representation allows the automatic extraction of structural features from 3D-EM volumes and their subsequent characterisation. FEMME is being populated with 3D-EM data stored in the electron microscopy database EMD-DB (http://www.ebi.ac.uk/msd/). However, and since the number of entries in EMD-DB is still relatively small, FEMME is also being populated in this initial phase with structural data from PDB and PQS databases (http://www.rcsb.org/pdb/ and pqs.ebi.ac.uk/, respectively) whose resolution has been lowered accordingly. Each FEMME entry contains macromolecular geometry and topology information with a detailed description of its structural features. Moreover, FEMME data have facilitated the study and development of a method to retrieve macromolecular structures by their structural content based on the combined use of spin images and neural networks with encouraging results. Therefore, the FEMME database constitutes a powerful tool that provides a uniform and automatic way of analysing volumes coming from 3D-EM that will hopefully help the scientific community to perform wide structural comparisons.

A structural-informatics approach for mining beta-sheets: locating sheets in intermediate-resolution density maps

Here, we report a new computational method, called sheetminer, for mining beta-sheets in the density maps at intermediate resolutions of 6 to 10A. The method employs a multi-step ad hoc morphological analysis of density maps to identify the unique characteristics of beta-sheets. It was tested on density maps from 12 protein crystal structures that were artificially blurred to intermediate resolutions. There are a total of 35 independent beta-sheets with a wide distribution of morphology. The method successfully located 34 of them and missed only one. The method was also applied to an experimental 9A electron cryomicroscopic structure and an 8A X-ray density map. In both cases, the sheet-searching results were found to agree very well with known high-resolution crystal structures. Collectively, these results demonstrate clearly the robustness of sheetminer in locating the regions belonging to beta-sheets in the intermediate-resolution density maps. Furthermore, sheetminer is completely complementary to all other existing computational methods, including helixhunter and threading algorithms. Their combined usage has the potential to significantly enhance the computational modeling capacity for a much more complete interpretation of structural data at intermediate resolutions, from which extraction of functional information would be more effective. This is particularly important in the field of structural genomics, in which the fast screening approach may not always yield crystals that diffract to atomic resolution. An exciting future application of sheetminer is as a valuable tool for revealing the structures of amyloid fibrils that are rich in beta-motifs.

Dehydron: a structurally encoded signal for protein interaction

We introduce a quantifiable structural motif, called dehydron, that is shown to be central to protein-protein interactions. A dehydron is a defectively packed backbone hydrogen bond suggesting preformed monomeric structure whose Coulomb energy is highly sensitive to binding-induced water exclusion. Such preformed hydrogen bonds are effectively adhesive, since water removal from their vicinity contributes to their stability. At the structural level, a significant correlation is established between dehydrons and sites for protein complexation, with the HIV-1 capsid protein P24 complexed with antibody light-chain FAB25.3 providing the most dramatic correlation. Furthermore, the number of dehydrons in homologous similar-fold proteins from different species is shown to be a signature of proteomic complexity. The techniques are then applied to higher levels of organization: The formation of the capsid and its organization in picornaviruses correlates strongly with the distribution of dehydrons on the rim of the virus unit. Furthermore, antibody contacts and crystal contacts may be assigned to dehydrons still prevalent after the capsid has been assembled. The implications of the dehydron as an encoded signal in proteomics, bioinformatics, and inhibitor drug design are emphasized.

Architecture of the herpes simplex virus major capsid protein derived from structural bioinformatics

The dispositions of 39 alpha helices of greater than 2.5 turns and four beta sheets in the major capsid protein (VP5, 149 kDa) of herpes simplex virus type 1 were identified by computational and visualization analysis from the 8.5A electron cryomicroscopy structure of the whole capsid. The assignment of helices in the VP5 upper domain was validated by comparison with the recently determined crystal structure of this region. Analysis of the spatial arrangement of helices in the middle domain of VP5 revealed that the organization of a tightly associated bundle of ten helices closely resembled that of a domain fold found in the annexin family of proteins. Structure-based sequence searches suggested that sequences in both the N and C-terminal portions of the VP5 sequence contribute to this domain. The long helices seen in the floor domain of VP5 form an interconnected network within and across capsomeres. The combined structural and sequence-based informatics has led to an architectural model of VP5. This model placed in the context of the capsid provides insights into the strategies used to achieve viral capsid stability.

The 3D structure of the fusion primed Sendai F-protein determined by electron cryomicroscopy

The three dimensional (3D) structure of the ectodomain of the entire fusion mediating F protein from Sendai virus [MW (trimer) approximately 177 kDa] has been determined by cryoelectron microscopy of single molecules and subsequent 3D reconstruction at a resolution of approximately 16 A. The reconstruction, which has been obtained from the native, proteolytic processed fusion primed F1+F2 form, shows the protein protruding approximately 170 A out of the membrane in a homotrimeric association. It consists of a defined approximately 65 A wide distal head and an adjacent neck, which is connected to an 70 A elongated stalk. Although the overall shape appears to be similar to the recently reported X-ray structure of the Newcastle disease virus F protein, a closer comparison reveals structural differences suggesting that the investigated Sendai F structure represents an advanced state towards the fusion active conformation.

Structural studies on adenoviruses

The adenovirus genome encodes more than 40 proteins, of which 11 combine with the viral DNA to form an icosahedral capsid of approximately 150 MDa molecular weight and approximately 900 A in diameter. This chapter reviews the information that structural biology techniques have provided about the adenovirus proteins and capsid. The structures of two capsid proteins (hexon and fiber) and two non-structural polypeptides (DNA-binding protein and protease) have been solved by X-ray crystallography. Fiber and its knob have been the focus of the latest structural studies, due to their role in host recognition and consequently in virus targeting for human gene therapy. The current model for the large capsid comes from a combination of electron microscopy and crystallography. The resultant images have revealed a surprising similarity between adenovirus and a bacterial virus, which suggests their common evolutionary origin.

Accurate determination of local defocus and specimen tilt in electron microscopy

Accurate knowledge of defocus and tilt parameters is essential for the determination of three-dimensional protein structures at high resolution using electron microscopy. We present two computer programs, CTFFIND3 and CTFTILT, which determine defocus parameters from images of untilted specimens, as well as defocus and tilt parameters from images of tilted specimens, respectively. Both programs use a simple algorithm that fits the amplitude modulations visible in a power spectrum with a calculated contrast transfer function (CTF). The background present in the power spectrum is calculated using a low-pass filter. The background is then subtracted from the original power spectrum, allowing the fitting of only the oscillatory component of the CTF. CTFTILT determines specimen tilt parameters by measuring the defocus at a series of locations on the image while constraining them to a single plane. We tested the algorithm on images of two-dimensional crystals by comparing the results with those obtained using crystallographic methods. The images also contained contrast from carbon support film that added to the visibility of the CTF oscillations. The tests suggest that the fitting procedure is able to determine the image defocus with an error of about 10nm, whereas tilt axis and tilt angle are determined with an error of about 2 degrees and 1 degrees, respectively. Further tests were performed on images of single protein particles embedded in ice that were recorded from untilted or slightly tilted specimens. The visibility of the CTF oscillations from these images was reduced due to the lack of a carbon support film. Nevertheless, the test results suggest that the fitting procedure is able to determine image defocus and tilt angle with errors of about 100 nm and 6 degrees, respectively.

Cytoplasmic polyhedrosis virus structure at 8 A by electron cryomicroscopy: structural basis of capsid stability and mRNA processing regulation

The single-shelled cytoplasmic polyhedrosis virus (CPV) is a unique member of the Reoviridae. Despite lacking protective outer shells, it exhibits striking capsid stability and is capable of endogenous RNA transcription and processing. The 8 A three-dimensional structure of CPV by electron cryomicroscopy reveals secondary structure elements present in the capsid proteins CSP, LPP, and TP, which have alpha+beta folds. The extensive nonequivalent interactions between CSP and LPP, the unique CSP protrusion domain, and the perfect inter-CSP surface complementarities may account for the enhanced capsid stability. The slanted disposition of TP functional domains and the stacking of channel constrictions suggest an iris diaphragm-like mechanism for opening/closing capsid pores and turret channels in regulating the highly coordinated steps of mRNA transcription, processing, and release.

Studying Large Viruses

In this article we have attempted to describe some structural aspects of large viruses. Although this may seem a straightforward task, it is complicated by the fact that large viruses do not represent a distinctive class of organisms and any grouping under this heading will include a range of unrelated viruses with different structures, replication strategies, and host types. To simplify matters we limited our definition to dsDNA viruses with genomes of 100 kbp or larger. However, even this restricted grouping includes viruses with diverse and seemingly unrelated structures. Furthermore, few if any structural features are exclusive to large viruses and most of what appears distinctive about their structure or assembly can also be found in smaller, and usually better characterized, viruses. Therefore we have not attempted to provide a comprehensive catalog of the properties of large viruses but have tried to illustrate particular structural points with examples from a few of the better known forms, notably herpes simplex virus (HSV) and phage T4. The two techniques used to provide rigorous analyses of virus structures are X-ray crystallography and electron cryomicroscopy with computer-assisted reconstruction. To date, X-ray crystallography has been successful only with smaller viruses, and what is known about the structures of these large viruses has come primarily from electron cryomicroscopy. However, with the notable exception of the HSV capsid, such studies have been limited in extent and of relatively low resolution, and the information obtained has been confined largely to describing the spatial distributions and relationships between the subunits. Nevertheless, these studies have given us our clearest insights into the biology of these complex particles and increases in resolution promise to extend these insights by bridging the gap between gross and atomic structures, as exemplified by the identification and mapping of secondary structural elements in the HSV capsid.

Three-dimensional structure of the type 1 inositol 1,4,5-trisphosphate receptor at 24 A resolution

We report here the first three-dimensional structure of the type 1 inositol 1,4,5-trisphosphate receptor (IP(3)R). From cryo-electron microscopic images of purified receptors embedded in vitreous ice, a three-dimensional structure was determined by use of standard single particle reconstruction techniques. The structure is strikingly different from that of the ryanodine receptor at similar resolution despite molecular similarities between these two calcium release channels. The 24 A resolution structure of the IP(3)R takes the shape of an uneven dumbbell, and is approximately 170 A tall. Its larger end is bulky, with four arms protruding laterally by approximately 50 A and, in comparison with the receptor topology, probably corresponds to the cytoplasmic domain of the receptor. The lateral dimension at the height of the protruding arms is approximately 155 A. The smaller end, whose lateral dimension is approximately 100 A, has structural features indicative of the membrane-spanning domain. A central opening in this domain, which is occluded on the cytoplasmic half, outlines a pathway for calcium flow in the open state of the channel.

How to describe protein motion without amino acid sequence and atomic coordinates

This paper reports a computational method, the quantized elastic deformational model, that can reliably describe the conformational flexibility of a protein in the absence of the amino acid sequence and atomic coordinates. The essence of this method lies in the fact that, in modeling the functionally important conformational changes such as domain movements, it is possible to abandon the traditional concepts of protein structure (bonds, angles, dihedrals, etc.) and treat the protein as an elastic object. The shape and mass distribution of the object are described by the electron density maps, at various resolutions, from methods such as x-ray diffraction or cryo-electron microscopy. The amplitudes and directionality of the elastic deformational modes of a protein, whose patterns match the biologically relevant conformational changes, can then be derived solely based on the electron density map. The method yields an accurate description of protein dynamics over a wide range of resolutions even as low as 15-20 A at which there is nearly no visually distinguishable internal structures. Therefore, this method dramatically enhances the capability of studying protein motions in structural biology. It is also expected to have ample applications in related fields such as bioinformatics, structural genomics, and proteomics, in which one's ability to extract functional information from the not-so-well-defined structural models is vitally important.

Molecular interactions and viral stability revealed by structural analyses of chemically treated cypovirus capsids

Cytoplasmic polyhedrosis virus (CPV, genus Cypovirus) is a unique member of the family Reoviridae which lacks the outer protective shells that exist in all other members, yet exhibits unusual stability. We have analyzed the effects of different acidic, basic, detergent, and urea treatments on CPV capsids. The integrity of the CPV capsids was unaffected under high-pH conditions that disrupted the orthoreovirus inner core, consistent with its ability to maintain structural integrity in extremely alkaline environments during infection. However, it was sensitive to low pH, detergents, and urea, similarly to other viruses in this family. The three-dimensional structure comparisons by electron cryomicroscopy of the intact empty CPV capsid with the "spikeless" capsid whose turrets were removed by chemical treatments revealed the interaction footprint of the turret on the capsid shell. The observed structural changes associated with the removal of the turret suggest critical structural roles of the turret in maintaining capsid integrity in addition to its enzymatic activities.

Deriving folds of macromolecular complexes through electron cryomicroscopy and bioinformatics approaches

Intermediate-resolution (7-9A) structures of large macromolecular complexes can be obtained by electron cryomicroscopy. This structural information, combined with bioinformatics data for the individual protein components or domains, can lead to a fold model for the entire complex. Such approaches have been demonstrated with the 6.8 A structure of the rice dwarf virus to derive models for the major capsid shell proteins.

High resolution structural studies of complex icosahedral viruses: a brief overview

Structural descriptions of viral particles are key to our understanding of their assembly mechanisms and properties. We will describe the application of X-ray crystallography and electron cryomicroscopy to the structural determination of the bluetongue virus core and the herpesvirus capsid. These represent the highest resolution structural studies carried out by these techniques on such complex and large icosahedral virus particles. The bluetongue virus core consists of two layers of distinct proteins with different protein packing symmetries, while the herpes virus capsid is made up of four types of proteins with 3.3 MDa per asymmetric unit. The structural results reveal that each of these proteins has distinct folds and they are packed uniquely to form stable particles.

Semi-automated icosahedral particle reconstruction at sub-nanometer resolution

Electron cryomicroscopy of large macromolecular complexes is becoming an increasingly powerful tool for revealing three-dimensional structures without the need for crystallization. The execution of image processing, however, requires experience and is error-prone due to the need for a human operator to carry out interactive and repetitive processes. We have designed an approach which is intended to make image processing simple and rapid, both for experts and for novice users. We demonstrate this approach using the well-established reconstruction scheme for icosahedral particles. Finally, we implement semi-automated virus reconstruction (SAVR), an expert system that integrates the most CPU intensive and iterative steps using the scripting language Python. SAVR is portable across platforms and has been parallelized to run on both shared and distributed memory platforms. SAVR also allows the incorporation of new algorithms and facilitates the management of the increasingly large data sets needed to achieve higher resolution reconstructions. The package has been successfully applied to several data sets and shown capable of generating icosahedral reconstructions to sub-nanometer resolutions (7-10 A ).

Combining electron microscopic with x-ray crystallographic structures

Analgorithm has been developed for placing three-dimensional atomic structures into appropriately scaled cryoelectron microscopy maps. The first stage in this process is to conduct a three-dimensional angular search in which the center of gravity of an X-ray crystallographically determined structure is placed on a selected position in the cryoelectron microscopy map. The quality of the fit is measured by the sum of the density at each atomic position. The second stage is to refine the three angles and three translational parameters for the best (usually 25 to 100) fits. Useful criteria for this refinement include the sum of densities at atomic sites, the lack of atoms in negative or low density, the absence of atomic clashes between symmetry-related positions of the atomic structure, and the distances between identifiable features in the map and their positions on the fitted atomic structure. These refinements generally lead to a convergence of the originally chosen, top scoring fits to just a few (about 3 to 8) acceptable possibilities. Usually, the best remaining fit is clearly superior to any of the others.