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

Hierarchical organization and assembly of the archaeal cell sheath from an amyloid-like protein

Certain archaeal cells possess external proteinaceous sheath, whose structure and organization are both unknown. By cellular cryogenic electron tomography (cryoET), here we have determined sheath organization of the prototypical archaeon, Methanospirillum hungatei. Fitting of Alphafold-predicted model of the sheath protein (SH) monomer into the 7.9 Å-resolution structure reveals that the sheath cylinder consists of axially stacked β-hoops, each of which is comprised of two to six 400 nm-diameter rings of β-strand arches (β-rings). With both similarities to and differences from amyloid cross-β fibril architecture, each β-ring contains two giant β-sheets contributed by ~ 450 SH monomers that entirely encircle the outer circumference of the cell. Tomograms of immature cells suggest models of sheath biogenesis: oligomerization of SH monomers into β-ring precursors after their membrane-proximal cytoplasmic synthesis, followed by translocation through the unplugged end of a dividing cell, and insertion of nascent β-hoops into the immature sheath cylinder at the junction of two daughter cells.

Cryo-EM reconstruction of the Cafeteria roenbergensis virus capsid suggests novel assembly pathway for giant viruses

Whereas the protein composition and overall shape of several giant virus capsids have been described, the mechanism by which these large capsids assemble remains enigmatic. Here, we present a reconstruction of the capsid of Cafeteria roenbergensis virus (CroV), one of the largest viruses analyzed by cryo-electron microscopy (cryo-EM) to date. The CroV capsid has a diameter of 3,000 Å and a Triangulation number of 499. Unlike related mimiviruses, the CroV capsid is not decorated with glycosylated surface fibers, but features 30 Å-long surface protrusions that are formed by loops of the major capsid protein. Based on the orientation of capsomers in the cryo-EM reconstruction, we propose that the capsids of CroV and related giant viruses are assembled by a newly conceived assembly pathway that initiates at a five-fold vertex and continuously proceeds outwards in a spiraling fashion.

Double-stranded RNA virus outer shell assembly by bona fide domain-swapping

Correct outer protein shell assembly is a prerequisite for virion infectivity in many multi-shelled dsRNA viruses. In the prototypic dsRNA bacteriophage φ6, the assembly reaction is promoted by calcium ions but its biomechanics remain poorly understood. Here, we describe the near-atomic resolution structure of the φ6 double-shelled particle. The outer T=13 shell protein P8 consists of two alpha-helical domains joined by a linker, which allows the trimer to adopt either a closed or an open conformation. The trimers in an open conformation swap domains with each other. Our observations allow us to propose a mechanistic model for calcium concentration regulated outer shell assembly. Furthermore, the structure provides a prime exemplar of bona fide domain-swapping. This leads us to extend the theory of domain-swapping from the level of monomeric subunits and multimers to closed spherical shells, and to hypothesize a mechanism by which closed protein shells may arise in evolution.

Double-shelled bacteriophage φ6 is a well-studied model system used to understand assembly of dsRNA viruses. Here the authors report a near-atomic resolution cryo-EM structure of φ6 and propose a model for the structural transitions occurring in the outer shell during genome packaging.

De Novo modeling in cryo-EM density maps with Pathwalking

As electron cryo-microscopy (cryo-EM) can now frequently achieve near atomic resolution, accurate interpretation of these density maps in terms of atomistic detail has become paramount in deciphering macromolecular structure and function. However, there are few software tools for modeling protein structure from cryo-EM density maps in this resolution range. Here, we present an extension of our original Pathwalking protocol, which can automatically trace a protein backbone directly from a near-atomic resolution (3-6Å) density map. The original Pathwalking approach utilized a Traveling Salesman Problem solver for backbone tracing, but manual adjustment was still required during modeling. In the new version, human intervention is minimized and we provide a more robust approach for backbone modeling. This includes iterative secondary structure identification, termini detection and the ability to model multiple subunits without prior segmentation. Overall, the new Pathwalking procedure provides a more complete and robust tool for annotating protein structure function in near-atomic resolution density maps.

The right place at the right time: chaperoning core histone variants

Histone proteins dynamically regulate chromatin structure and epigenetic signaling to maintain cell homeostasis. These processes require controlled spatial and temporal deposition and eviction of histones by their dedicated chaperones. With the evolution of histone variants, a network of functionally specific histone chaperones has emerged. Molecular details of the determinants of chaperone specificity for different histone variants are only slowly being resolved. A complete understanding of these processes is essential to shed light on the genuine biological roles of histone variants, their chaperones, and their impact on chromatin dynamics.

Sequence Conservation, Radial Distance and Packing Density in Spherical Viral Capsids

The conservation level of a residue is a useful measure about the importance of that residue in protein structure and function. Much information about sequence conservation comes from aligning homologous sequences. Profiles showing the variation of the conservation level along the sequence are usually interpreted in evolutionary terms and dictated by site similarities of a proper set of homologous sequences. Here, we report that, of the viral icosahedral capsids, the sequence conservation profile can be determined by variations in the distances between residues and the centroid of the capsid – with a direct inverse proportionality between the conservation level and the centroid distance – as well as by the spatial variations in local packing density. Examining both the centroid and the packing density models against a dataset of 51 crystal structures of nonhomologous icosahedral capsids, we found that many global patterns and minor features derived from the viral structures are consistent with those present in the sequence conservation profiles. The quantitative link between the level of conservation and structural features like centroid-distance or packing density allows us to look at residue conservation from a structural viewpoint as well as from an evolutionary viewpoint.

A newly isolated reovirus has the simplest genomic and structural organization of any reovirus

A total of 2,691 mosquitoes representing 17 species was collected from eight locations in southwest Cameroon and screened for pathogenic viruses. Ten isolates of a novel reovirus (genus Dinovernavirus) were detected by culturing mosquito pools on Aedes albopictus (C6/36) cell cultures. A virus that caused overt cytopathic effects was isolated, but it did not infect vertebrate cells or produce detectable disease in infant mice after intracerebral inoculation. The virus, tentatively designated Fako virus (FAKV), represents the first 9-segment, double-stranded RNA (dsRNA) virus to be isolated in nature. FAKV appears to have a broad mosquito host range, and its detection in male specimens suggests mosquito-to-mosquito transmission in nature. The structure of the T=1 FAKV virion, determined to subnanometer resolution by cryoelectron microscopy (cryo-EM), showed only four proteins per icosahedral asymmetric unit: a dimer of the major capsid protein, one turret protein, and one clamp protein. While all other turreted reoviruses of known structures have at least two copies of the clamp protein per asymmetric unit, FAKV's clamp protein bound at only one conformer of the major capsid protein. The FAKV capsid architecture and genome organization represent the most simplified reovirus described to date, and phylogenetic analysis suggests that it arose from a more complex ancestor by serial loss-of-function events.

IMPORTANCE

We describe the detection, genetic, phenotypic, and structural characteristics of a novel Dinovernavirus species isolated from mosquitoes collected in Cameroon. The virus, tentatively designated Fako virus (FAKV), is related to both single-shelled and partially double-shelled viruses. The only other described virus in this genus was isolated from cultured mosquito cells. It was previously unclear whether the phenotypic characteristics of that virus were reflective of this genus in nature or were altered during serial passaging in the chronically infected cell line. FAKV is a naturally occurring single-shelled reovirus with a unique virion architecture that lacks several key structural elements thought to stabilize a single-shelled reovirus virion, suggesting what may be the minimal number of proteins needed to form a viable reovirus particle. FAKV evolved from more complex ancestors by losing a genome segment and several virion proteins.

Four levels of hierarchical organization, including noncovalent chainmail, brace the mature tumor herpesvirus capsid against pressurization

Like many double-stranded DNA viruses, tumor gammaherpesviruses Epstein-Barr virus and Kaposi's sarcoma-associated herpesvirus withstand high internal pressure. Bacteriophage HK97 uses covalent chainmail for this purpose, but how this is achieved noncovalently in the much larger gammaherpesvirus capsid is unknown. Our cryoelectron microscopy structure of a gammaherpesvirus capsid reveals a hierarchy of four levels of organization: (1) Within a hexon capsomer, each monomer of the major capsid protein (MCP), 1,378 amino acids and six domains, interacts with its neighboring MCPs at four sites. (2) Neighboring capsomers are linked in pairs by MCP dimerization domains and in groups of three by heterotrimeric triplex proteins. (3) Small (∼280 amino acids) HK97-like domains in MCP monomers alternate with triplex heterotrimers to form a belt that encircles each capsomer. (4) One hundred sixty-two belts concatenate to form noncovalent chainmail. The triplex heterotrimer orchestrates all four levels and likely drives maturation to an angular capsid that can withstand pressurization.

Transcriptome profiling confirmed correlations between symptoms and transcriptional changes in RDV infected rice and revealed nucleolus as a possible target of RDV manipulation

BACKGROUND

Rice dwarf virus (RDV) is the causal agent of rice dwarf disease, which limits rice production in many areas of south East Asia. Transcriptional changes of rice in response to RDV infection have been characterized by Shimizu et al. and Satoh et al.. Both studies found induction of defense related genes and correlations between transcriptional changes and symptom development in RDV-infected rice. However, the same rice cultivar, namely Nipponbare belonging to the Japonic subspecies of rice was used in both studies.

METHODS

Gene expression changes of the indica subspecies of rice, namely Oryza sativa L. ssp. indica cv Yixiang2292 that show moderate resistance to RDV, in response to RDV infection were characterized using an Affymetrix Rice Genome Array. Differentially expressed genes (DEGs) were classified according to their Gene Ontology (GO) annotation. The effects of transient expression of Pns11 in Nicotiana benthaminana on the expression of nucleolar genes were studied using real-time PCR (RT-PCR).

RESULTS

856 genes involved in defense or other physiological processes were identified to be DEGs, most of which showed up-regulation. Ribosome- and nucleolus related genes were significantly enriched in the DEGs. Representative genes related to nucleolar function exhibited altered expression in N. benthaminana plants transiently expressing Pns11 of RDV.

CONCLUSIONS

Induction of defense related genes is common for rice infected with RDV. There is a co-relation between symptom severity and transcriptional alteration in RDV infected rice. Besides ribosome, RDV may also target nucleolus to manipulate the translation machinery of rice. Given the tight links between nucleolus and ribosome, it is intriguing to speculate that RDV may enhance expression of ribosomal genes by targeting nucleolus through Pns11.

3D structures of fungal partitiviruses

Partitiviruses constitute one of the nine currently recognized families of viruses with encapsidated, double-stranded (ds)RNA genomes. The partitivirus genome is bisegmented, and each genome segment is packaged inside a separate viral capsid. Different partitiviruses infect plants, fungi, or protozoa. Recent studies have shed light on the three-dimensional structures of the virions of three representative fungal partitiviruses. These structures include a number of distinctive features, allowing informative comparisons with the structures of dsRNA viruses from other families. The results and comparisons suggest several new conclusions about the functions, assembly, and evolution of these viruses.

FOLD-EM: automated fold recognition in medium- and low-resolution (4-15 Å) electron density maps

MOTIVATION

Owing to the size and complexity of large multi-component biological assemblies, the most tractable approach to determining their atomic structure is often to fit high-resolution radiographic or nuclear magnetic resonance structures of isolated components into lower resolution electron density maps of the larger assembly obtained using cryo-electron microscopy (cryo-EM). This hybrid approach to structure determination requires that an atomic resolution structure of each component, or a suitable homolog, is available. If neither is available, then the amount of structural information regarding that component is limited by the resolution of the cryo-EM map. However, even if a suitable homolog cannot be identified using sequence analysis, a search for structural homologs should still be performed because structural homology often persists throughout evolution even when sequence homology is undetectable, As macromolecules can often be described as a collection of independently folded domains, one way of searching for structural homologs would be to systematically fit representative domain structures from a protein domain database into the medium/low resolution cryo-EM map and return the best fits. Taken together, the best fitting non-overlapping structures would constitute a 'mosaic' backbone model of the assembly that could aid map interpretation and illuminate biological function.

RESULT

Using the computational principles of the Scale-Invariant Feature Transform (SIFT), we have developed FOLD-EM-a computational tool that can identify folded macromolecular domains in medium to low resolution (4-15 Å) electron density maps and return a model of the constituent polypeptides in a fully automated fashion. As a by-product, FOLD-EM can also do flexible multi-domain fitting that may provide insight into conformational changes that occur in macromolecular assemblies.

Macromolecular structure modeling from 3D EM using VolRover 2.0

We review tools for structure identification and model-based refinement from three-dimensional electron microscopy implemented in our in-house software package, VOLROVER 2.0. For viral density maps with icosahedral symmetry, we segment the capsid, polymeric, and monomeric subunits using techniques based on automatic symmetry detection and multidomain fast marching. For large biomolecules without symmetry information, we again use our multidomain fast-marching method with manual or fit-based multiseeding to segment meaningful substructures. In either case, we subject the resulting segmented subunit to secondary structure detection when the EM resolution is sufficiently high, and rigid-body structure fitting when the corresponding X-ray structure is available. Secondary structure elements are identified by three techniques: our earlier volume-based and boundary-based skeletonization methods as well as a new method, currently in development, based on solving the grassfire flow equation. For rigid-body fitting, we adapt our earlier fast Fourier-based correlation scheme F2Dock. Our reported segmentation, secondary structure elements identification, and rigid-body fitting techniques, implemented in VOLROVER 2.0 are applied to the PSB 2011 cryo-EM modeling challenge data, and our results are briefly compared to similar results submitted from other research groups. The comparisons show that our techniques are equally capable of segmenting relatively accurate subunits from a viral or protein assembly, and that high segmentation quality leads in turn to higher-quality results of secondary structure elements identification and correlation-based rigid-body fitting. © 2012 Wiley Periodicals, Inc. Biopolymers 97: 709-731, 2012.

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.

Detection of secondary and supersecondary structures of proteins from cryo-electron microscopy

Recent advances in three-dimensional electron microscopy (3D EM) have enabled the quantitative visualization of the structural building blocks of proteins at improved resolutions. We provide algorithms to detect the secondary structures (α-helices and β-sheets) from proteins for which the volumetric maps are reconstructed at 6-10Å resolution. Additionally, we show that when the resolution is coarser than 10Å, some of the supersecondary structures can be detected from 3D EM maps. For both these algorithms, we employ tools from computational geometry and differential topology, specifically the computation of stable/unstable manifolds of certain critical points of the distance function induced by the molecular surface. Our results connect mathematically well-defined constructions with bio-chemically induced structures observed in proteins.

DERIVING TOPOLOGY AND SEQUENCE ALIGNMENT FOR THE HELIX SKELETON IN LOW-RESOLUTION PROTEIN DENSITY MAPS

Cryoelectron microscopy (cryoEM) is an experimental technique to determine the three-dimensional (3D) structure of large protein complexes. Currently, this technique is able to generate protein density maps at 6–9 Å resolution, at which the skeleton of the structure (which is composed of α-helices and β-sheets) can be visualized. As a step towards predicting the entire backbone of the protein from the protein density map, we developed a method to predict the topology and sequence alignment for the skeleton helices. Our method combines the geometrical information of the skeleton helices with the Rosetta ab initio structure prediction method to derive a consensus topology and sequence alignment for the skeleton helices. We tested the method with 60 proteins. For 45 proteins, the majority of the skeleton helices were assigned a correct topology from one of our top ten predictions. The offsets of the alignment for most of the assigned helices were within ±2 amino acids in the sequence. We also analyzed the use of the skeleton helices as a clustering tool for the decoy structures generated by Rosetta. Our comparison suggests that the topology clustering is a better method than a general overlap clustering method to enrich the ranking of decoys, particularly when the decoy pool is small.

Atomic model of CPV reveals the mechanism used by this single-shelled virus to economically carry out functions conserved in multishelled reoviruses

Unlike the multishelled viruses in the Reoviridae, cytoplasmic polyhedrosis virus (CPV) is single shelled, yet stable and fully capable of carrying out functions conserved within Reoviridae. Here, we report a 3.1 Å resolution cryo electron microscopy structure of CPV and derive its atomic model, consisting of 60 turret proteins (TPs), 120 each of capsid shell proteins (CSPs) and large protrusion proteins (LPPs). Two unique segments of CSP contribute to CPV's stability: an inserted protrusion domain interacting with neighboring proteins, and an N-anchor tying up CSPs together through strong interactions such as β sheet augmentation. Without the need to interact with outer shell proteins, LPP retains only the N-terminal two-third region containing a conserved helix-barrel core and interacts exclusively with CSP. TP is also simplified, containing only domains involved in RNA capping. Our results illustrate how CPV proteins have evolved in a coordinative manner to economically carry out their conserved functions.

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.

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.

Fingerprint-based structure retrieval using electron density

We present a computational approach that can quickly search a large protein structural database to identify structures that fit a given electron density, such as determined by cryo-electron microscopy. We use geometric invariants (fingerprints) constructed using 3D Zernike moments to describe the electron density, and reduce the problem of fitting of the structure to the electron density to simple fingerprint comparison. Using this approach, we are able to screen the entire Protein Data Bank and identify structures that fit two experimental electron densities determined by cryo-electron microscopy.

Atomic resolution cryo electron microscopy of macromolecular complexes

Single-particle cryo electron microscopy (cryoEM) is a technique for determining three-dimensional (3D) structures from projection images of molecular complexes preserved in their "native," noncrystalline state. Recently, atomic or near-atomic resolution structures of several viruses and protein assemblies have been determined by single-particle cryoEM, allowing ab initio atomic model building by following the amino acid side chains or nucleic acid bases identifiable in their cryoEM density maps. In particular, these cryoEM structures have revealed extended arms contributing to molecular interactions that are otherwise not resolved by the conventional structural method of X-ray crystallography at similar resolutions. High-resolution cryoEM requires careful consideration of a number of factors, including proper sample preparation to ensure structural homogeneity, optimal configuration of electron imaging conditions to record high-resolution cryoEM images, accurate determination of image parameters to correct image distortions, efficient refinement and computation to reconstruct a 3D density map, and finally appropriate choice of modeling tools to construct atomic models for functional interpretation. This progress illustrates the power of cryoEM and ushers it into the arsenal of structural biology, alongside conventional techniques of X-ray crystallography and NMR, as a major tool (and sometimes the preferred one) for the studies of molecular interactions in supramolecular assemblies or machines.

Bluetongue virus coat protein VP2 contains sialic acid-binding domains, and VP5 resembles enveloped virus fusion proteins

Bluetongue virus (BTV) is transmitted by blood-feeding insects (Culicoides sp.) and causes hemorrhagic diseases in livestock. BTV is a nonenveloped, double-stranded RNA (dsRNA) virus with two capsids: a well-studied, stable core enclosing the dsRNA genome and a highly unstable, poorly studied coat responsible for host cell attachment and entry. Here, based on cryo-electron microscopy (cryoEM), we report a 7-A resolution structure of the infectious BTV virion, including the coat proteins. We show that unlike other dsRNA viruses, the VP2 attachment trimer has a triskelion shape composed of three tip domains branching from a central hub domain. We identify three putative sialic acid-binding pockets in the hub and present supporting biochemical data indicating sugar moiety binding is important for BTV infection. Despite being a nonenveloped virus, the putative VP5 membrane penetration trimer, located slightly inward of the VP2 attachment trimer, has a central coiled-coil alpha-helical bundle, similar to the fusion proteins of many enveloped viruses (e.g., HIV, herpesviruses, vesicular stomatitis virus, and influenza virus). Moreover, mapping of the amino acid sequence of VP5 to the secondary structural elements identified by cryoEM locates 15 amphipathic alpha-helical regions on the external surface of each VP5 trimer. The cryoEM density map also reveals few, weak interactions between the VP5 trimer and both the outer-coat VP2 trimer and the underlying core VP7 trimer, suggesting that the surface of VP5 could unfurl like an umbrella during penetration and shedding of the coat to release the transcriptionally active core particle.

Quantitative analysis of cryo-EM density map segmentation by watershed and scale-space filtering, and fitting of structures by alignment to regions

Cryo-electron microscopy produces 3D density maps of molecular machines, which consist of various molecular components such as proteins and RNA. Segmentation of individual components in such maps is a challenging task, and is mostly accomplished interactively. We present an approach based on the immersive watershed method and grouping of the resulting regions using progressively smoothed maps. The method requires only three parameters: the segmentation threshold, a smoothing step size, and the number of smoothing steps. We first apply the method to maps generated from molecular structures and use a quantitative metric to measure the segmentation accuracy. The method does not attain perfect accuracy, however it produces single or small groups of regions that roughly match individual proteins or subunits. We also present two methods for fitting of structures into density maps, based on aligning the structures with single regions or small groups of regions. The first method aligns centers and principal axes, whereas the second aligns centers and then rotates the structure to find the best fit. We describe both interactive and automated ways of using these two methods. Finally, we show segmentation and fitting results for several experimentally-obtained 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.

An automated procedure for detecting protein folds from sub-nanometer resolution electron density

The use of sub-nanometer resolution electron density as spatial constraints for de novo and ab initio structure prediction requires knowledge of protein boundaries to accurately segment the electron density for the prediction algorithms. Here we present a procedure where even poorly segmented density can be used to determine the fold of the protein. The method is automated, fast, capable of searching for multiple copies of a protein fold, and accessible to densities encompassing more than a thousand residues. The automation is particularly powerful as it allows the procedure to take full advantage of the expanding repository in the Protein Data Bank. We have tested the method on nine segmented sub-nanometer image reconstruction electron densities. The method successfully identifies the correct fold for the six densities for which an atomic structure is known, identifies a fold that agrees with prior structural data, a fold that agrees with predictions from the Fold & Function Assignment server, and a fold that correlates with secondary structure prediction. The identified folds in the last three examples can be used as templates for comparative modeling of the bacteriophage P22 tail-machine (a 3MDa complex composed of 39 protein subunits).

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.

Multiresolution phase extension of a trypsin inhibitor structure from 5 A to 2 A based on diffraction amplitudes alone

For a specific structure of trypsin inhibitor, starting from a 5 A phase set, about 80% of the 2 A phases are correctly determined within an error of 18 degrees by applying a multiresolution refinement procedure. The refinement proceeds both in real and reciprocal spaces. In extending the structure from 5 to 2 A, the amplitudes of the reflections are the only requirements for this procedure. In contrast to the conventional phase extension scheme, the amplitudes are used not only in the Fourier synthesis but also in the real-space density modifications.

Refinement of protein structures into low-resolution density maps using rosetta

We describe a method based on Rosetta structure refinement for generating high-resolution, all-atom protein models from electron cryomicroscopy density maps. A local measure of the fit of a model to the density is used to directly guide structure refinement and to identify regions incompatible with the density that are then targeted for extensive rebuilding. Over a range of test cases using both simulated and experimentally generated data, the method consistently increases the accuracy of starting models generated either by comparative modeling or by hand-tracing the density. The method can achieve near-atomic resolution starting from density maps at 4-6 A resolution.

A structural model of the pore-forming region of the skeletal muscle ryanodine receptor (RyR1)

Ryanodine receptors (RyRs) are ion channels that regulate muscle contraction by releasing calcium ions from intracellular stores into the cytoplasm. Mutations in skeletal muscle RyR (RyR1) give rise to congenital diseases such as central core disease. The absence of high-resolution structures of RyR1 has limited our understanding of channel function and disease mechanisms at the molecular level. Here, we report a structural model of the pore-forming region of RyR1. Molecular dynamics simulations show high ion binding to putative pore residues D4899, E4900, D4938, and D4945, which are experimentally known to be critical for channel conductance and selectivity. We also observe preferential localization of Ca(2+) over K(+) in the selectivity filter of RyR1. Simulations of RyR1-D4899Q mutant show a loss of preference to Ca(2+) in the selectivity filter as seen experimentally. Electrophysiological experiments on a central core disease mutant, RyR1-G4898R, show constitutively open channels that conduct K(+) but not Ca(2+). Our simulations with G4898R likewise show a decrease in the preference of Ca(2+) over K(+) in the selectivity filter. Together, the computational and experimental results shed light on ion conductance and selectivity of RyR1 at an atomistic level.

Hybrid approaches: applying computational methods in cryo-electron microscopy

Recent advances in cryo-electron microscopy have led to an increasing number of high (3-5A) to medium (5-10A) resolution cryoEM density maps. These density maps contain valuable information about the protein structure but frequently require computational algorithms to aid their structural interpretation. It is these hybrid approaches between cryoEM and computational protein structure prediction algorithms that will shape protein structure elucidation from density maps.

Estimating contrast transfer function and associated parameters by constrained non-linear optimization

The three-dimensional reconstruction of macromolecules from two-dimensional single-particle electron images requires determination and correction of the contrast transfer function (CTF) and envelope function. A computational algorithm based on constrained non-linear optimization is developed to estimate the essential parameters in the CTF and envelope function model simultaneously and automatically. The application of this estimation method is demonstrated with focal series images of amorphous carbon film as well as images of ice-embedded icosahedral virus particles suspended across holes.

Atomic structure reveals the unique capsid organization of a dsRNA virus

For most dsRNA viruses, the genome-enclosing capsid comprises 120 copies of a single capsid protein (CP) organized into 60 icosahedrally equivalent dimers, generally identified as 2 nonsymmetricallyinteracting CP molecules with extensive lateral contacts. The crystal structure of a partitivirus, Penicillium stoloniferum virus F (PsV-F), reveals a different organization, in which the CP dimer is related by almost-perfect local 2-fold symmetry, forms prominent surface arches, and includes extensive structure swapping between the 2 subunits. An electron cryomicroscopy map of PsV-F shows that the disordered N terminus of each CP molecule interacts with the dsRNA genome and probably participates in its packaging or transcription. Intact PsV-F particles mediate semiconservative transcription, and transcripts are likely to exit through negatively charged channels at the icosahedral 5-fold axes. Other findings suggest that the PsV-F capsid is assembled from dimers of CP dimers, with an arrangement similar to flavivirus E glycoproteins.

Reduction of the secondary structure topological space through direct estimation of the contact energy formed by the secondary structures

Background

Electron cryomicroscopy is a fast developing technique aiming at the determination of the 3-dimensional structures of large protein complexes. Using this technique, protein density maps can be generated with 6 to 10 Å resolution. At such resolutions, the secondary structure elements such as helices and β-strands appear to be skeletons and can be computationally detected. However, it is not known which segment of the protein sequence corresponds to which of the skeletons. The topology in this paper refers to the linear order and the directionality of the secondary structures. For a protein with N helices and M strands, there are (N!2^N)(M!2^M) different topologies, each of which maps N helix segments and M strand segments on the protein sequence to N helix and M strand skeletons. Since the backbone position is not available in the skeleton, each topology of the skeletons corresponds to additional freedom to position the atoms in the skeletons.

Results

We have developed a method to construct the possible atomic structures for the helix skeletons by sampling the solution space of all the possible topologies of the skeletons. Our method also ranks the possible structures based on the contact energy formed by the secondary structures, rather than the entire chain. If we assume that the backbone atomic positions are known for the skeletons, then the native topology of the secondary structures can be found in the top 30% of the ranked list of all possible topologies for all the 30 proteins tested, and within the top 5% for most of the 30 proteins. Without assuming the backbone location of the skeletons, the possible atomic structures of the skeletons can be constructed using the axis of the skeleton and the sequence segments. The best constructed structure for the skeletons has RMSD to native between 4 and 5 Å for the four tested α-proteins. These best constructed structures were ranked the 17^th, 31^st, 16^th and 5^th respectively for the four proteins out of 32066, 391833, 98755 and 192935 possible assignments in the pool.

Conclusion

Our work suggested that the direct estimation of the contact energy formed by the secondary structures is quite effective in reducing the topological space to a small subset that includes a near native structure for the skeletons.

Rotavirus architecture at subnanometer resolution

Rotavirus, a nonturreted member of the Reoviridae, is the causative agent of severe infantile diarrhea. The double-stranded RNA genome encodes six structural proteins that make up the triple-layer particle. X-ray crystallography has elucidated the structure of one of these capsid proteins, VP6, and two domains from VP4, the spike protein. Complementing this work, electron cryomicroscopy (cryoEM) has provided relatively low-resolution structures for the triple-layer capsid in several biochemical states. However, a complete, high-resolution structural model of rotavirus remains unresolved. Combining new structural analysis techniques with the subnanometer-resolution cryoEM structure of rotavirus, we now provide a more detailed structural model for the major capsid proteins and their interactions within the triple-layer particle. Through a series of intersubunit interactions, the spike protein (VP4) adopts a dimeric appearance above the capsid surface, while forming a trimeric base anchored inside one of the three types of aqueous channels between VP7 and VP6 capsid layers. While the trimeric base suggests the presence of three VP4 molecules in one spike, only hints of the third molecule are observed above the capsid surface. Beyond their interactions with VP4, the interactions between VP6 and VP7 subunits could also be readily identified. In the innermost T=1 layer composed of VP2, visualization of the secondary structure elements allowed us to identify the polypeptide fold for VP2 and examine the complex network of interactions between this layer and the T=13 VP6 layer. This integrated structural approach has resulted in a relatively high-resolution structural model for the complete, infectious structure of rotavirus, as well as revealing the subtle nuances required for maintaining interactions in such a large macromolecular assembly.

Computational approaches for automatic structural analysis of large biomolecular complexes

We present computational solutions to two problems of macromolecular structure interpretation from reconstructed three-dimensional electron microscopy (3D-EM) maps of large bio-molecular complexes at intermediate resolution (5A-15 A). The two problems addressed are: 1) 3D structural alignment (matching) between identified and segmented 3D maps of structure units (e.g. trimeric configuration of proteins), and 2) the secondary structure identification of a segmented protein 3D map (i.e.locations of alpha-helices, beta-sheets). For problem 1, we present an efficient algorithm to correlate spatially (and structurally) two 3D maps of structure units. Besides providing a similarity score between structure units, the algorithm yields an effective technique for resolution refinement of repeated structure units, by 3D alignment and averaging. For problem 2, we present an efficient algorithm to compute eigenvalues and link eigenvectors of a Gaussian convoluted structure tensor derived from the protein 3D Map, thereby identifying and locating secondary structural motifs of proteins. The efficiency and performance of our approach is demonstrated on several experimentally reconstructed 3D maps of virus capsid shells from single-particle cryo-electron microscopy (cryo-EM), as well as computationally simulated protein structure density 3D maps generated from protein model entries in the Protein Data Bank.

Macromolecular modeling with rosetta

Advances over the past few years have begun to enable prediction and design of macromolecular structures at near-atomic accuracy. Progress has stemmed from the development of reasonably accurate and efficiently computed all-atom potential functions as well as effective conformational sampling strategies appropriate for searching a highly rugged energy landscape, both driven by feedback from structure prediction and design tests. A unified energetic and kinematic framework in the Rosetta program allows a wide range of molecular modeling problems, from fibril structure prediction to RNA folding to the design of new protein interfaces, to be readily investigated and highlights areas for improvement. The methodology enables the creation of novel molecules with useful functions and holds promise for accelerating experimental structural inference. Emerging connections to crystallographic phasing, NMR modeling, and lower-resolution approaches are described and critically assessed.

DNA poised for release in bacteriophage phi29

We present here the first asymmetric, three-dimensional reconstruction of a tailed dsDNA virus, the mature bacteriophage phi29, at subnanometer resolution. This structure reveals the rich detail of the asymmetric interactions and conformational dynamics of the phi29 protein and DNA components, and provides novel insight into the mechanics of virus assembly. For example, the dodecameric head-tail connector protein undergoes significant rearrangement upon assembly into the virion. Specific interactions occur between the tightly packed dsDNA and the proteins of the head and tail. Of particular interest and novelty, an approximately 60A diameter toroid of dsDNA was observed in the connector-lower collar cavity. The extreme deformation that occurs over a small stretch of DNA is likely a consequence of the high pressure of the packaged genome. This toroid structure may help retain the DNA inside the capsid prior to its injection into the bacterial host.

Subnanometer-resolution structures of the grass carp reovirus core and virion

Grass carp reovirus (GCRV) is a member of the Aquareovirus genus of the family Reoviridae, a large family of double-stranded RNA (dsRNA) viruses infecting plants, insects, fishes and mammals. We report the first subnanometer-resolution three-dimensional structures of both GCRV core and virion by cryoelectron microscopy. These structures have allowed the delineation of interactions among the over 1000 molecules in this enormous macromolecular machine and a detailed comparison with other dsRNA viruses at the secondary-structure level. The GCRV core structure shows that the inner proteins have strong structural similarities with those of orthoreoviruses even at the level of secondary-structure elements, indicating that the structures involved in viral dsRNA interaction and transcription are highly conserved. In contrast, the level of similarity in structures decreases in the proteins situated in the outer layers of the virion. The proteins involved in host recognition and attachment exhibit the least similarities to other members of Reoviridae. Furthermore, in GCRV, the RNA-translocating turrets are in an open state and lack a counterpart for the sigma1 protein situated on top of the close turrets observed in mammalian orthoreovirus. Interestingly, the distribution and the organization of GCRV core proteins resemble those of the cytoplasmic polyhedrosis virus, a cypovirus and the structurally simplest member of the Reoviridae family. Our results suggest that GCRV occupies a unique structure niche between the simpler cypoviruses and the considerably more complex mammalian orthoreovirus, thus providing an important model for understanding the structural and functional conservation and diversity of this enormous family of dsRNA viruses.

Towards atomic resolution structural determination by single-particle cryo-electron microscopy

Recent advances in cryo-electron microscopy and single-particle reconstruction (collectively referred to as 'cryoEM') have made it possible to determine the three-dimensional (3D) structures of several macromolecular complexes at near-atomic resolution ( approximately 3.8-4.5A). These achievements were accomplished by overcoming the challenges in sample handling, instrumentation, image processing, and model building. At near-atomic resolution, many detailed structural features can be resolved, such as the turns and deep grooves of helices, strand separation in beta sheets, and densities for loops and bulky amino acid side chains. Such structural data of the cytoplasmic polyhedrosis virus (CPV), the Epsilon 15 bacteriophage and the GroEL complex have provided valuable constraints for atomic model building using integrative tools, thus significantly enhancing the value of the cryoEM structures. The CPV structure revealed a drastic conformational change from a helix to a beta hairpin associated with RNA packaging and replication, coupling of RNA processing and release, and the long sought-after polyhedrin-binding domain. These latest advances in single-particle cryoEM provide exciting opportunities for the 3D structural determination of viruses and macromolecular complexes that are either too large or too heterogeneous to be investigated by conventional X-ray crystallography or nuclear magnetic resonance (NMR) methods.

Integrating diverse data for structure determination of macromolecular assemblies

To understand the cell, we need to determine the macromolecular assembly structures, which may consist of tens to hundreds of components. First, we review the varied experimental data that characterize the assemblies at several levels of resolution. We then describe computational methods for generating the structures using these data. To maximize completeness, resolution, accuracy, precision, and efficiency of the structure determination, a computational approach is required that uses spatial information from a variety of experimental methods. We propose such an approach, defined by its three main components: a hierarchical representation of the assembly, a scoring function consisting of spatial restraints derived from experimental data, and an optimization method that generates structures consistent with the data. This approach is illustrated by determining the configuration of the 456 proteins in the nuclear pore complex (NPC) from baker's yeast. With these tools, we are poised to integrate structural information gathered at multiple levels of the biological hierarchy--from atoms to cells--into a common framework.

Symmetry: A guide to its application in 2D electron crystallography

A defining property of a crystal is its symmetry. This mini-review sets out to summarize all aspects that define 2D crystallographic symmetry as applied to the study of macromolecular structure. It begins by defining molecular point symmetries, before covering crystallographic symmetry operations in 2D, common notation, a summary of crystallographic plane groups and theoretical methods and important considerations for the identification and application of symmetry in 2D crystal images for 3D structure determination. While many of the concepts covered here may be equally applicable to point symmetry and space group symmetry in 3D, this review has been written from the perspective of 2D electron crystallography and deals specifically with symmetry operations and crystallographic space groups in 2D crystal projection images.

Quasi-atomic model of bacteriophage t7 procapsid shell: insights into the structure and evolution of a basic fold

The existence of similar folds among major structural subunits of viral capsids has shown unexpected evolutionary relationships suggesting common origins irrespective of the capsids' host life domain. Tailed bacteriophages are emerging as one such family, and we have studied the possible existence of the HK97-like fold in bacteriophage T7. The procapsid structure at approximately 10 A resolution was used to obtain a quasi-atomic model by fitting a homology model of the T7 capsid protein gp10 that was based on the atomic structure of the HK97 capsid protein. A number of fold similarities, such as the fitting of domains A and P into the L-shaped procapsid subunit, are evident between both viral systems. A different feature is related to the presence of the amino-terminal domain of gp10 found at the inner surface of the capsid that might play an important role in the interaction of capsid and scaffolding proteins.

Averaging tens to hundreds of icosahedral particle images to resolve protein secondary structure elements using a Multi-Path Simulated Annealing optimization algorithm

Accurately determining a cryoEM particle's alignment parameters is crucial to high resolution single particle 3-D reconstruction. We developed Multi-Path Simulated Annealing, a Monte-Carlo type of optimization algorithm, for globally aligning the center and orientation of a particle simultaneously. A consistency criterion was developed to ensure the alignment parameters are correct and to remove some bad particles from a large pool of images of icosahedral particles. Without using any a priori model, this procedure is able to reconstruct a structure from a random initial model. Combining the procedure above with a new empirical double threshold particle selection method, we are able to pick tens of best quality particles to reconstruct a subnanometer resolution map from scratch. Using the best 62 particles of rice dwarf virus, the reconstruction reached 9.6A resolution at which four helices of the P3A subunit of RDV are resolved. Furthermore, with the 284 best particles, the reconstruction is improved to 7.9A resolution, and 21 of 22 helices and six of seven beta sheets are resolved.

The pore structure of the closed RyR1 channel

Using single particle electron cryomicroscopy, several helices in the membrane-spanning region of RyR1, including an inner transmembrane helix, a short pore helix, and a helix parallel to the membrane on the cytoplasmic side, have been clearly resolved. Our model places a highly conserved glycine (G4934) at the hinge position of the bent inner helix and two rings of negative charges at the luminal and cytoplasmic mouths of the pore. The kinked inner helix closely resembles the inner helix of the open MthK channel, suggesting that kinking alone does not open RyR1, as proposed for K+ channels.

3D electron microscopy of biological nanomachines: principles and applications

Transmission electron microscopy is a powerful technique for studying the three-dimensional (3D) structure of a wide range of biological specimens. Knowledge of this structure is crucial for fully understanding complex relationships among macromolecular complexes and organelles in living cells. In this paper, we present the principles and main application domains of 3D transmission electron microscopy in structural biology. Moreover, we survey current developments needed in this field, and discuss the close relationship of 3D transmission electron microscopy with other experimental techniques aimed at obtaining structural and dynamical information from the scale of whole living cells to atomic structure of macromolecular complexes.

Identification of secondary structure elements in intermediate-resolution density maps

An increasing number of structural studies of large macromolecular complexes, both in X-ray crystallography and cryo-electron microscopy, have resulted in intermediate-resolution (5-10 A) density maps. Despite being limited in resolution, significant structural and functional information may be extractable from these maps. To aid in the analysis and annotation of these complexes, we have developed SSEhunter, a tool for the quantitative detection of alpha helices and beta sheets. Based on density skeletonization, local geometry calculations, and a template-based search, SSEhunter has been tested and validated on a variety of simulated and authentic subnanometer-resolution density maps. The result is a robust, user-friendly approach that allows users to quickly visualize, assess, and annotate intermediate-resolution density maps. Beyond secondary structure element identification, the skeletonization algorithm in SSEhunter provides secondary structure topology, which is potentially useful in leading to structural models of individual molecular components directly from the density.

Cryoelectron microscopy of icosahedral virus particles

With the rapid progresses in both instrumentation and computing, it is increasingly straightforward and routine to determine the structures of icosahedral viruses to subnanometer resolutions (6-10 A) by cryoelectron microscopy and image reconstruction. In this resolution range, secondary structure elements of protein subunits can be clearly discerned. Combining the three-dimensional density map and bioinformatics of the protein components, the folds of the virus capsid shell proteins can be derived. This chapter will describe the experimental and computational procedures that lead to subnanometer resolution structural determinations of icosahedral virus particles. In addition, we will describe how to extract useful structural information from the three-dimensional maps.

Rotavirus proteins: structure and assembly

Rotavirus is a major pathogen of infantile gastroenteritis. It is a large and complex virus with a multilayered capsid organization that integrates the determinants of host specificity, cell entry, and the enzymatic functions necessary for endogenous transcription of the genome that consists of 11 dsRNA segments. These segments encode six structural and six nonstructural proteins. In the last few years, there has been substantial progress in our understanding of both the structural and functional aspects of a variety of molecular processes involved in the replication of this virus. Studies leading to this progress using of a variety of structural and biochemical techniques including the recent application of RNA interference technology have uncovered several unique and intriguing features related to viral morphogenesis. This review focuses on our current understanding of the structural basis of the molecular processes that govern the replication of rotavirus.