Resolution and Probabilistic Models of Components in CryoEM Maps of Mature P22 Bacteriophage

CryoEM continues to produce density maps of larger and more complex assemblies with multiple protein components of mixed symmetries. Resolution is not always uniform throughout a cryoEM map, and it can be useful to estimate the resolution in specific molecular components of a large assembly. In this study, we present procedures to 1) estimate the resolution in subcomponents by gold-standard Fourier shell correlation (FSC); 2) validate modeling procedures, particularly at medium resolutions, which can include loop modeling and flexible fitting; and 3) build probabilistic models that combine high-accuracy priors (such as crystallographic structures) with medium-resolution cryoEM densities. As an example, we apply these methods to new cryoEM maps of the mature bacteriophage P22, reconstructed without imposing icosahedral symmetry. Resolution estimates based on gold-standard FSC show the highest resolution in the coat region (7.6 Å), whereas other components are at slightly lower resolutions: portal (9.2 Å), hub (8.5 Å), tailspike (10.9 Å), and needle (10.5 Å). These differences are indicative of inherent structural heterogeneity and/or reconstruction accuracy in different subcomponents of the map. Probabilistic models for these subcomponents provide new insights, to our knowledge, and structural information when taking into account uncertainty given the limitations of the observed density.     

Visualization of the externalized VP2 N termini of infectious human parvovirus B19

The structures of infectious human parvovirus B19 and empty wild-type particles were determined by cryoelectron microscopy (cryoEM) to 7.5-Å and 11.3-Å resolution, respectively, assuming icosahedral symmetry. Both of these, DNA filled and empty, wild-type particles contain a few copies of the minor capsid protein VP1. Comparison of wild-type B19 with the crystal structure and cryoEM reconstruction of recombinant B19 particles consisting of only the major capsid protein VP2 showed structural differences in the vicinity of the icosahedral fivefold axes. Although the unique N-terminal region of VP1 could not be visualized in the icosahedrally averaged maps, the N terminus of VP2 was shown to be exposed on the viral surface adjacent to the fivefold β-cylinder. The conserved glycine-rich region is positioned between two neighboring, fivefold-symmetrically related VP subunits and not in the fivefold channel as observed for other parvoviruses.     

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.     

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.     

Advances in structure analysis using small-angle scattering in solution

The resolution and reliability of solution scattering models have been significantly improved by ab initio shape and domain structure determination, and by detailed modelling of macromolecular complexes using rigid-body refinement. Substantial progress has also been made in the quantitative analysis and modelling of assembly and folding processes, and intermolecular interactions.     

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.     

Atomic Model of Rabbit Hemorrhagic Disease Virus by Cryo-Electron Microscopy and Crystallography

Rabbit hemorrhagic disease, first described in China in 1984, causes hemorrhagic necrosis of the liver. Its etiological agent, rabbit hemorrhagic disease virus (RHDV), belongs to the Lagovirus genus in the family Caliciviridae. The detailed molecular structure of any lagovirus capsid has yet to be determined. Here, we report a cryo-electron microscopic (cryoEM) reconstruction of wild-type RHDV at 6.5 Å resolution and the crystal structures of the shell (S) and protruding (P) domains of its major capsid protein, VP60, each at 2.0 Å resolution. From these data we built a complete atomic model of the RHDV capsid. VP60 has a conserved S domain and a specific P2 sub-domain that differs from those found in other caliciviruses. As seen in the shell portion of the RHDV cryoEM map, which was resolved to ∼5.5 Å, the N-terminal arm domain of VP60 folds back onto its cognate S domain. Sequence alignments of VP60 from six groups of RHDV isolates revealed seven regions of high variation that could be mapped onto the surface of the P2 sub-domain and suggested three putative pockets might be responsible for binding to histo-blood group antigens. A flexible loop in one of these regions was shown to interact with rabbit tissue cells and contains an important epitope for anti-RHDV antibody production. Our study provides a reliable, pseudo-atomic model of a Lagovirus and suggests a new candidate for an efficient vaccine that can be used to protect rabbits from RHDV infection.     

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

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

Gorgon and pathwalking: macromolecular modeling tools for subnanometer resolution density maps

The complex interplay of proteins and other molecules, often in the form of large transitory assemblies, are critical to cellular function. Today, X-ray crystallography and electron cryo-microscopy (cryo-EM) are routinely used to image these macromolecular complexes, though often at limited resolutions. Despite the rapidly growing number of macromolecular structures, few tools exist for modeling and annotating structures in the range of 3-10 ? resolution. To address this need, we have developed a number of utilities specifically targeting subnanometer resolution density maps. As part of the 2010 Cryo-EM Modeling Challenge, we demonstrated two of our latest de novo modeling tools, Pathwalking and Gorgon, as well as a tool for secondary structure identification (SSEHunter) and a new rigid-body/flexible fitting tool in Gorgon. In total, we submitted 30 structural models from ten different subnanometer resolution data sets in four of the six challenge categories. Each of our utlities produced accurate structural models and annotations across the various density maps. In the end, the utilities that we present here offer users a robust toolkit for analyzing and modeling protein structure in macromolecular assemblies at non-atomic resolutions.     

Crystal Structure of a Two-Subunit TrkA Octameric Gating Ring Assembly

The TM1088 locus of T. maritima codes for two proteins designated TM1088A and TM1088B, which combine to form the cytosolic portion of a putative Trk K⁺ transporter. We report the crystal structure of this assembly to a resolution of 3.45 Å. The high resolution crystal structures of the components of the assembly, TM1088A and TM1088B, were also determined independently to 1.50 Å and 1.55 Å, respectively. The TM1088 proteins are structurally homologous to each other and to other K⁺ transporter proteins, such as TrkA. These proteins form a cytosolic gating ring assembly that controls the flow of K⁺ ions across the membrane. TM1088 represents the first structure of a two-subunit Trk assembly. Despite the atypical genetics and chain organization of the TM1088 assembly, it shares significant structural homology and an overall quaternary organization with other single-subunit K⁺ gating ring assemblies. This structure provides the first structural insights into what may be an evolutionary ancestor of more modern single-subunit K⁺ gating ring assemblies.     

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.     

Structural polymorphism of the major capsid protein of rotavirus

Rotaviruses are important human pathogens with a triple-layered icosahedral capsid. The major capsid protein VP6 is shown here to self-assemble into spherical or helical particles mainly depending upon pH. Assembly is inhibited either by low pH (<3.0) or by a high concentration (>100 mM) of divalent cations (Ca(2+) and Zn(2+)). The structures of two types of helical tubes were determined by electron cryomicroscopy and image analysis to a resolution of 2.0 and 2.5 nm. In both reconstructions, the molecular envelope of VP6 fits the atomic model determined by X-ray crystallography remarkably well. The 3-fold symmetry of the VP6 trimer, being incompatible with the helical symmetry, is broken at the level of the trimer contacts. One type of contact is maintained within all VP6 particles (tubes and virus), strongly suggesting that VP6 assemblies arise from different packings of a unique dimer of trimers. Our data show that the protonation state and thus the charge distribution are important switches governing the assembly of macromolecular assemblies.     

Scoring functions for cryoEM density fitting

In fitting atomic structures into cryoEM density maps of macromolecular assemblies, the cross-correlation function (CCF) is the most prevalent method of scoring the goodness-of-fit. However, there are still many possible, less studied ways of scoring fits. In this paper, we introduce four scores new to cryoEM fitting and compare their performance to three known scores. Our benchmark consists of (a) 4 protein assemblies with simulated maps at 5-20 ? resolution, including the heptameric ring of GroEL; and (b) 4 experimental maps of GroEL at ～6-23 ? resolution with corresponding fitted atomic models. We perturb each fit 1000 times and assess each new fit with each score. The correlation between a score and the Cα RMSD of each fit from the "correctly" fitted structure shows that the CCF is one of the best scores, but in certain situations could be augmented or even replaced by other scores. For instance, our implementation of a score based on mutual information outperforms or is comparable to the CCF in almost all test cases, and our new "envelope score" works as well as the CCF at sub-nanometer resolution but is an order of magnitude faster to calculate. The results also suggest that the width of the Gaussian function used to blur the atomic structure into a density map can significantly affect the fitting process. Finally, we show that our score-testing method, when combined with the Laplacian CCF or the mutual information scores, can be used as a statistical tool for improving cryoEM density fitting.     

Joint use of small-angle X-ray and neutron scattering to study biological macromolecules in solution

Novel techniques for simultaneous analysis of X-ray and neutron scattering patterns from macromolecular complexes in solution are presented. They include ab initio shape and internal structure determination of multicomponent particles and more detailed rigid body modeling of complexes using high resolution structures of subunits. The methods fit simultaneously X-ray and neutron scattering curves including contrast variation data sets from selectively deuterated complexes. Biochemically sound interconnected models without steric clashes between the components displaying a pre-defined symmetry are generated. For rigid body modeling, distance restraints between specified residues/nucleotides or their ranges are taken into account. The efficiency of the methods is demonstrated in model examples, and potential sources of ambiguity are discussed.     

Combined approaches to flexible fitting and assessment in virus capsids undergoing conformational change

Fitting of atomic components into electron cryo-microscopy (cryoEM) density maps is routinely used to understand the structure and function of macromolecular machines. Many fitting methods have been developed, but a standard protocol for successful fitting and assessment of fitted models has yet to be agreed upon among the experts in the field. Here, we created and tested a protocol that highlights important issues related to homology modelling, density map segmentation, rigid and flexible fitting, as well as the assessment of fits. As part of it, we use two different flexible fitting methods (Flex-EM and iMODfit) and demonstrate how combining the analysis of multiple fits and model assessment could result in an improved model. The protocol is applied to the case of the mature and empty capsids of Coxsackievirus A7 (CAV7) by flexibly fitting homology models into the corresponding cryoEM density maps at 8.2 and 6.1 Å resolution. As a result, and due to the improved homology models (derived from recently solved crystal structures of a close homolog – EV71 capsid – in mature and empty forms), the final models present an improvement over previously published models. In close agreement with the capsid expansion observed in the EV71 structures, the new CAV7 models reveal that the expansion is accompanied by ∼5° counterclockwise rotation of the asymmetric unit, predominantly contributed by the capsid protein VP1. The protocol could be applied not only to viral capsids but also to many other complexes characterised by a combination of atomic structure modelling and cryoEM density fitting.     

Atomic model of an infectious rotavirus particle

Non‐enveloped viruses of different types have evolved distinct mechanisms for penetrating a cellular membrane during infection. Rotavirus penetration appears to occur by a process resembling enveloped‐virus fusion: membrane distortion linked to conformational changes in a viral protein. Evidence for such a mechanism comes from crystallographic analyses of fragments of VP4, the rotavirus‐penetration protein, and infectivity analyses of structure‐based VP4 mutants. We describe here the structure of an infectious rotavirus particle determined by electron cryomicroscopy (cryoEM) and single‐particle analysis at about 4.3 Å resolution. The cryoEM image reconstruction permits a nearly complete trace of the VP4 polypeptide chain, including the positions of most side chains. It shows how the two subfragments of VP4 (VP8^* and VP5^*) retain their association after proteolytic cleavage, reveals multiple structural roles for the β‐barrel domain of VP5^*, and specifies interactions of VP4 with other capsid proteins. The virion model allows us to integrate structural and functional information into a coherent mechanism for rotavirus entry.     

Symmetry-adapted spherical harmonics method for high-resolution 3D single-particle reconstructions

Three-dimensional (3D) reconstruction is the last and an essential step toward high-resolution structural determination in single-particle cryo-electron microscopy (cryoEM). We have implemented a new algorithm for reconstructing 3D structures of macromolecular complexes with icosahedral symmetry from cryoEM images. Icosahedral symmetry-adapted functions (ISAFs) are used to interpolate structural factors in the reciprocal space to generate a 3D reconstruction in spherical coordinates. In our implementation, we introduced a recursive method for deriving higher order ISAFs from three lower order seed functions. We demonstrate improvements of our new method in both the noise suppression and the effective resolution in 3D reconstruction over the commonly used Fourier-Bessel synthesis method introduced by Crowther et al. three decades ago. Our 3D reconstruction method can be extended to macromolecular complexes with other symmetry types and is thus likely to impact future high-resolution cryoEM single-particle reconstruction efforts in general.     

Crystal Structures of a Piscine Betanodavirus: Mechanisms of Capsid Assembly and Viral Infection

Betanodaviruses cause massive mortality in marine fish species with viral nervous necrosis. The structure of a T = 3 Grouper nervous necrosis virus-like particle (GNNV-LP) is determined by the ab initio method with non-crystallographic symmetry averaging at 3.6 Å resolution. Each capsid protein (CP) shows three major domains: (i) the N-terminal arm, an inter-subunit extension at the inner surface; (ii) the shell domain (S-domain), a jelly-roll structure; and (iii) the protrusion domain (P-domain) formed by three-fold trimeric protrusions. In addition, we have determined structures of the T = 1 subviral particles (SVPs) of (i) the delta-P-domain mutant (residues 35−217) at 3.1 Å resolution; and (ii) the N-ARM deletion mutant (residues 35−338) at 7 Å resolution; and (iii) the structure of the individual P-domain (residues 214−338) at 1.2 Å resolution. The P-domain reveals a novel DxD motif asymmetrically coordinating two Ca²⁺ ions, and seems to play a prominent role in the calcium-mediated trimerization of the GNNV CPs during the initial capsid assembly process. The flexible N-ARM (N-terminal arginine-rich motif) appears to serve as a molecular switch for T = 1 or T = 3 assembly. Finally, we find that polyethylene glycol, which is incorporated into the P-domain during the crystallization process, enhances GNNV infection. The present structural studies together with the biological assays enhance our understanding of the role of the P-domain of GNNV in the capsid assembly and viral infection by this betanodavirus.     

Structural analysis of flexible proteins in solution by small angle X-ray scattering combined with crystallography

In the last few years, SAXS of biological materials has been rapidly evolving and promises to move structural analysis to a new level. Recent innovations in SAXS data analysis allow ab initio shape predictions of proteins in solution. Furthermore, experimental scattering data can be compared to calculated scattering curves from the growing data base of solved structures and also identify aggregation and unfolded proteins. Combining SAXS results with atomic resolution structures enables detailed characterizations in solution of mass, radius, conformations, assembly, and shape changes associated with protein folding and functions. SAXS can efficiently reveal the spatial organization of protein domains, including domains missing from or disordered in known crystal structures, and establish cofactor or substrate-induced conformational changes. For flexible domains or unstructured regions that are not amenable for study by many other structural techniques, SAXS provides a unique technology. Here, we present SAXS shape predictions for PCNA that accurately predict a trimeric ring assembly and for a full-length DNA repair glycosylase with a large unstructured region. These new results in combination with illustrative published data show how SAXS combined with high resolution crystal structures efficiently establishes architectures, assemblies, conformations, and unstructured regions for proteins and protein complexes in solution.