A general method to quantify quasi-equivalence in icosahedral viruses

Integrated cryoEM structure of a spumaretrovirus reveals cross-kingdom evolutionary relationships and the molecular basis for assembly and virus entry

Foamy viruses (FVs) are an ancient lineage of retroviruses, with an evolutionary history spanning over 450 million years. Vector systems based on Prototype Foamy Virus (PFV) are promising candidates for gene and oncolytic therapies. Structural studies of PFV contribute to the understanding of the mechanisms of FV replication, cell entry and infection, and retroviral evolution. Here we combine cryoEM and cryoET to determine high-resolution in situ structures of the PFV icosahedral capsid (CA) and envelope glycoprotein (Env), including its type III transmembrane anchor and membrane-proximal external region (MPER), and show how they are organized in an integrated structure of assembled PFV particles. The atomic models reveal an ancient retroviral capsid architecture and an unexpected relationship between Env and other class 1 fusion proteins of the Mononegavirales. Our results represent the de novo structure determination of an assembled retrovirus particle.

Inside and outside of virus-like particles HBc and HBc/4M2e: A comprehensive study of the structure

Hepatitis B virus core antigen (HBc) with the insertion of four external domains of the influenza A M2 protein (HBc/4M2e) form virus-like particles whose structure was studied using a combination of molecular modeling and cryo-electron microscopy (cryo-EM). It was also shown that self-assembling of the particles occurs inside bacterial cells, but despite the big inner volume of the core shell particle, purified HBc/4M2e contain an insignificant amount of bacterial proteins. It was shown that a fragment of the M2e corresponding to 4M2e insertion is prone to formation of amyloid-like fibrils. However, as the part of the immunodominant loop, M2e insertion does not show a tendency to intermolecular interaction. A full-atomic HBc-4M2e model with the resolution of about 3 Å (3.13 Å for particles of Т = 4 symmetry, 3.7 Å for particles of Т = 3 symmetry) was obtained by molecular modeling methods based on cryo-EM data.

Spindle-shaped archaeal viruses evolved from rod-shaped ancestors to package a larger genome

Spindle- or lemon-shaped viruses infect archaea in diverse environments. Due to the highly pleomorphic nature of these virions, which can be found with cylindrical tails emanating from the spindle-shaped body, structural studies of these capsids have been challenging. We have determined the atomic structure of the capsid of Sulfolobus monocaudavirus 1, a virus that infects hosts living in nearly boiling acid. A highly hydrophobic protein, likely integrated into the host membrane before the virions assemble, forms 7 strands that slide past each other in both the tails and the spindle body. We observe the discrete steps that occur as the tail tubes expand, and these are due to highly conserved quasiequivalent interactions with neighboring subunits maintained despite significant diameter changes. Our results show how helical assemblies can vary their diameters, becoming nearly spherical to package a larger genome and suggest how all spindle-shaped viruses have evolved from archaeal rod-like viruses.

VIPERdb: A Tool for Virus Research

The VIrus Particle ExploreR database (VIPERdb) ( http://viperdb.scripps.edu ) is a database and web portal for primarily icosahedral virus capsid structures that integrates structure-derived information with visualization and analysis tools accessed through a set of web interfaces. Our aim in developing VIPERdb is to provide comprehensive structure-derived information on viruses comprising simple to detailed attributes such as size (diameter), architecture ( T number), genome type, taxonomy, intersubunit association energies, and surface-accessible residues. In addition, a number of web-based tools are provided to enable users to interact with the structures and compare and contrast structure-derived properties between different viruses. Recently, we have constructed a series of data visualizations using modern JavaScript charting libraries such as Google Charts that allow users to explore trends and gain insights based on the various data available in the database. Furthermore, we now include helical viruses and nonicosahedral capsids by implementing modified procedures for data curation and analysis. This article provides an up-to-date overview of VIPERdb, describing various data and tools that are currently available and how to use them to facilitate structure-based bioinformatics analysis of virus capsids.

Viral Capsid Assembly: A Quantified Uncertainty Approach

Most of the existing research in assembly pathway prediction/analysis of viral capsids makes the simplifying assumption that the configuration of the intermediate states can be extracted directly from the final configuration of the entire capsid. This assumption does not take into account the conformational changes of the constituent proteins as well as minor changes to the binding interfaces that continue throughout the assembly process until stabilization. This article presents a statistical-ensemble-based approach that samples the configurational space for each monomer with the relative local orientation between monomers, to capture the uncertainties in binding and conformations. Further, instead of using larger capsomers (trimers, pentamers) as building blocks, we allow all possible subassemblies to bind in all possible combinations. We represent the resulting assembly graph in two different ways: First, we use the Wilcoxon signed-rank measure to compare the distributions of binding free energy computed on the sampled conformations to predict likely pathways. Second, we represent chemical equilibrium aspects of the transitions as a Bayesian Factor graph where both associations and dissociations are modeled based on concentrations and the binding free energies. We applied these protocols on the feline panleukopenia virus and the Nudaurelia capensis virus. Results from these experiments showed a significant departure from those that one would obtain if only the static configurations of the proteins were considered. Hence, we establish the importance of an uncertainty-aware protocol for pathway analysis, and we provide a statistical framework as an important first step toward assembly pathway prediction with high statistical confidence.

Uncertainty Quantified Computational Analysis of the Energetics of Virus Capsid Assembly

Most of the existing research in assembly pathway prediction/analysis of viral capsids makes the simplifying assumption that the configuration of the intermediate states can be extracted directly from the final configuration of the entire capsid. This assumption does not take into account the conformational changes of the constituent proteins as well as minor changes to the binding interfaces that continue throughout the assembly process until stabilization. This paper presents a statistical-ensemble based approach which samples the configurational space for each monomer with the relative local orientation between monomers, to capture the uncertainties in binding and conformations. Furthermore, instead of using larger capsomers (trimers, pentamers) as building blocks, we allow all possible subassemblies to bind in all possible combinations. We represent the resulting assembly graph in two different ways: First, we use the Wilcoxon signed rank measure to compare the distributions of binding free energy computed on the sampled conformations to predict likely pathways. Second, we represent chemical equilibrium aspects of the transitions as a Bayesian Factor graph where both associations and dissociations are modeled based on concentrations and the binding free energies. We applied these protocols on the feline panleukopenia virus and the Nudaurelia capensis virus. Results from these experiments showed significant departure from those one would obtain if only the static configurations of the proteins were considered. Hence, we establish the importance of an uncertainty-aware protocol for pathway analysis, and provide a statistical framework as an important first step towards assembly pathway prediction with high statistical confidence.

Structure based sequence analysis of viral and cellular protein assemblies

It is well accepted that, in general, protein structural similarity is strongly related to the amino acid sequence identity. To analyze in great detail the correlation, distribution and variation levels of conserved residues in the protein structure, we analyzed all available high-resolution structural data of 5245 cellular complex-forming proteins and 293 spherical virus capsid proteins (VCPs). We categorized and compare them in terms of protein structural regions. In all cases, the buried core residues are the most conserved, followed by the residues at the protein-protein interfaces. The solvent-exposed surface shows greater sequence variations. Our results provide evidence that cellular monomers and VCPs could be two extremes in the quaternary structural space, with cellular dimers and oligomers in between. Moreover, based on statistical analysis, we detected a distinct group of icosahedral virus families whose capsid proteins seem to evolve much slower than the rest of the protein complexes analyzed in this work.

Nodavirus-based biological container for targeted delivery system

Biological containers such as virus-like particles (VLPs) have gained increasing interest in the fields of gene therapy and vaccine development. Several virus-based materials have been studied, but the toxicity, biodistribution, and immunology of these systems still require extensive investigation. The specific goal of this review is to provide information about nodaviruses, which are causative infectious agents of insects and aquatic animals, but not humans. By understanding the structure and biophysical properties of such viruses, further chemical or genetic modification for novel nanocarriers could be developed. Therefore, their application for therapeutic purposes, particularly in humans, is of great interest.

CapsidMaps: protein-protein interaction pattern discovery platform for the structural analysis of virus capsids using Google Maps

Structural analysis and visualization of protein-protein interactions is a challenging task since it is difficult to appreciate easily the extent of all contacts made by the residues forming the interfaces. In the case of viruses, structural analysis becomes even more demanding because several interfaces coexist and, in most cases, these are formed by hundreds of contacting residues that belong to multiple interacting coat proteins. CapsidMaps is an interactive analysis and visualization tool that is designed to benefit the structural virology community. Developed as an improved extension of the φ-ψ Explorer, here we describe the details of its design and implementation. We present results of analysis of a spherical virus to showcase the features and utility of the new tool. CapsidMaps also facilitates the comparison of quaternary interactions between two spherical virus particles by computing a similarity (S)-score. The tool can also be used to identify residues that are solvent exposed and in the process of locating antigenic epitope regions as well as residues forming the inside surface of the capsid that interact with the nucleic acid genome. CapsidMaps is part of the VIPERdb Science Gateway, and is freely available as a web-based and cross-browser compliant application at http://viperdb.scripps.edu.

PCalign: a method to quantify physicochemical similarity of protein-protein interfaces

Background

Structural comparison of protein-protein interfaces provides valuable insights into the functional relationship between proteins, which may not solely arise from shared evolutionary origin. A few methods that exist for such comparative studies have focused on structural models determined at atomic resolution, and may miss out interesting patterns present in large macromolecular complexes that are typically solved by low-resolution techniques.

Results

We developed a coarse-grained method, PCalign, to quantitatively evaluate physicochemical similarities between a given pair of protein-protein interfaces. This method uses an order-independent algorithm, geometric hashing, to superimpose the backbone atoms of a given pair of interfaces, and provides a normalized scoring function, PC-score, to account for the extent of overlap in terms of both geometric and chemical characteristics. We demonstrate that PCalign outperforms existing methods, and additionally facilitates comparative studies across models of different resolutions, which are not accommodated by existing methods. Furthermore, we illustrate potential application of our method to recognize interesting biological relationships masked by apparent lack of structural similarity.

Conclusions

PCalign is a useful method in recognizing shared chemical and spatial patterns among protein-protein interfaces. It outperforms existing methods for high-quality data, and additionally facilitates comparison across structural models with different levels of details with proven robustness against noise.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-015-0471-x) contains supplementary material, which is available to authorized users.

Assembly, stability and dynamics of virus capsids

Most viruses use a hollow protein shell, the capsid, to enclose the viral genome. Virus capsids are large, symmetric oligomers made of many copies of one or a few types of protein subunits. Self-assembly of a viral capsid is a complex oligomerization process that proceeds along a pathway regulated by ordered interactions between the participating protein subunits, and that involves a series of (usually transient) assembly intermediates. Assembly of many virus capsids requires the assistance of scaffolding proteins or the viral nucleic acid, which interact with the capsid subunits to promote and direct the process. Once assembled, many capsids undergo a maturation reaction that involves covalent modification and/or conformational rearrangements, which may increase the stability of the particle. The final, mature capsid is a relatively robust protein complex able to protect the viral genome from physicochemical aggressions; however, it is also a metastable, dynamic structure poised to undergo controlled conformational transitions required to perform biologically critical functions during virus entry into cells, intracellular trafficking, and viral genome uncoating. This article provides an updated general overview on structural, biophysical and biochemical aspects of the assembly, stability and dynamics of virus capsids.

A novel method to map and compare protein-protein interactions in spherical viral capsids

Viral capsids are composed of multiple copies of one or a few chemically distinct capsid proteins and are mostly stabilized by inter subunit protein-protein interactions. There have been efforts to identify and analyze these protein-protein interactions, in terms of their extent and similarity, between the subunit interfaces related by quasi- and icosahedral symmetry. Here, we describe a new method to map quaternary interactions in spherical virus capsids onto polar angle space with respect to the icosahedral symmetry axes using azimuthal orthographic diagrams. This approach enables one to map the nonredundant interactions in a spherical virus capsid, irrespective of its size or triangulation number (T), onto the reference icosahedral asymmetric unit space. The resultant diagrams represent characteristic fingerprints of quaternary interactions of the respective capsids. Hence, they can be used as road maps of the protein-protein interactions to visualize the distribution and the density of the interactions. In addition, unlike the previous studies, the fingerprints of different capsids, when represented in a matrix form, can be compared with one another to quantitatively evaluate the similarity (S-score) in the subunit environments and the associated protein-protein interactions. The S-score selectively distinguishes the similarity, or lack of it, in the locations of the quaternary interactions as opposed to other well-known structural similarity metrics (e.g., RMSD, TM-score). Application of this method on a subset of T = 1 and T = 3 capsids suggests that S-score values range between 1 and 0.6 for capsids that belong to the same virus family/genus; 0.6-0.3 for capsids from different families with the same T-number and similar subunit fold; and <0.3 for comparisons of the dissimilar capsids that display different quaternary architectures (T-numbers). Finally, the sequence conserved interface residues within a virus family, whose spatial locations were also conserved have been hypothesized as the essential residues for self-assembly of the member virus capsids.

VIPERdb2: an enhanced and web API enabled relational database for structural virology

VIPERdb (http://viperdb.scripps.edu) is a relational database and a web portal for icosahedral virus capsid structures. Our aim is to provide a comprehensive resource specific to the needs of the virology community, with an emphasis on the description and comparison of derived data from structural and computational analyses of the virus capsids. In the current release, VIPERdb(2), we implemented a useful and novel method to represent capsid protein residues in the icosahedral asymmetric unit (IAU) using azimuthal polar orthographic projections, otherwise known as Phi-Psi (Phi-Psi) diagrams. In conjunction with a new Application Programming Interface (API), these diagrams can be used as a dynamic interface to the database to map residues (categorized as surface, interface and core residues) and identify family wide conserved residues including hotspots at the interfaces. Additionally, we enhanced the interactivity with the database by interfacing with web-based tools. In particular, the applications Jmol and STRAP were implemented to visualize and interact with the virus molecular structures and provide sequence-structure alignment capabilities. Together with extended curation practices that maintain data uniformity, a relational database implementation based on a schema for macromolecular structures and the APIs provided will greatly enhance the ability to do structural bioinformatics analysis of virus capsids.

Features of the mammalian orthoreovirus 3 Dearing l1 single-stranded RNA that direct packaging and serotype restriction

A series of recombinant mammalian orthoreoviruses (mammalian orthoreovirus 3 Dearing, MRV-3De) were generated that express an MRV-3De lambda3-CAT fusion protein. Individual viruses contain L1CAT double-stranded (ds) RNAs that range in length from a minimum of 1020 bp to 4616 bp. The engineered dsRNAs were generated from in vitro-transcribed single-stranded (ss) RNAs and incorporated into infectious virus particles by using reverse genetics. In addition to defining the sequences required for these ssRNAs to be 'identified' as l1 ssRNAs, the individual nucleotides in these regions that 'mark' each ssRNA as originating from mammalian orthoreovirus 1 Lang (MRV-1La), mammalian orthoreovirus 2 D5/Jones (MRV-2Jo) or MRV-3De have been identified. A C at position 81 in the MRV-1La 5' 129 nt sequence was able to be replaced with a U, as normally present in MRV-3De; this toggled the activity of the MRV-1La ssRNA to that of an MRV-3De 5' l1. RNA secondary-structure predictions for the 5' 129 nt of both the biologically active MRV-3De l1 ssRNA and the U81-MRV-3De-restored MRV-1La 5' ssRNA predicted a common structure.

Extent of protein-protein interactions and quasi-equivalence in viral capsids

Viral capsids are composed of multiple copies of one or a few gene products that self-assemble on their own or in the presence of the viral genome and/or auxiliary proteins into closed shells (capsids). We have analyzed 75 high-resolution virus capsid structures by calculating the average fraction of the solvent-accessible surface area of the coat protein subunits buried in the viral capsids. This fraction ranges from 0 to 1 and represents a normalized protein-protein interaction (PPI) index and is a measure of the extent of protein-protein interactions. The PPI indices were used to compare the extent of association of subunits among different capsids. We further examined the variation of the PPI indices as a function of the molecular weight of the coat protein subunit and the capsid diameter. Our results suggest that the PPI indices in T=1 and pseudo-T=3 capsids vary linearly with the molecular weight of the subunit and capsid size. This is in contrast to quasi-equivalent capsids with T>or=3, where the extent of protein-protein interactions is relatively independent of the subunit and capsid sizes. The striking outcome of this analysis is the distinctive clustering of the "T=2" capsids, which are distinguished by higher subunit molecular weights and a much lower degree of protein-protein interactions. Furthermore, the calculated residual (R(sym)) of the fraction buried surface areas of the structurally unique subunits in capsids with T>1 was used to calculate the quasi-equivalence of different subunit environments.

VIPERdb: a relational database for structural virology

VIPERdb () is a database for icosahedral virus capsid structures. Our aim is to provide a comprehensive resource specific to the needs of the structural virology community, with an emphasis on the description and comparison of derived data from structural and energetic analyses of capsids. A relational database implementation based on a schema for macromolecular structure makes the data highly accessible to the user, allowing detailed queries at the atomic level. Together with curation practices that maintain data uniformity, this will facilitate structural bioinformatics studies of virus capsids. User friendly search, visualization and educational tools on the website allow both structural and derived data to be examined easily and extensively. Links to relevant literature, sequence and taxonomy databases are provided for each entry.

Theoretical aspects of virus capsid assembly

A virus capsid is constructed from many copies of the same protein(s). Molecular recognition is central to capsid assembly. The capsid protein must polymerize in order to create a three-dimensional protein polymer. More than structure is required to understand this self-assembly reaction: one must understand how the pieces come together in solution.

Exploring icosahedral virus structures with VIPER

Virus structures are megadalton nucleoprotein complexes with an exceptional variety of protein-protein and protein-nucleic-acid interactions. Three-dimensional crystal structures of over 70 virus capsids, from more than 20 families and 30 different genera of viruses, have been solved to near-atomic resolution. The enormous amount of information contained in these structures is difficult to access, even for scientists trained in structural biology. Virus Particle Explorer (VIPER) is a web-based catalogue of structural information that describes the icosahedral virus particles. In addition to high-resolution crystal structures, VIPER has expanded to include virus structures obtained by cryo-electron microscopy (EM) techniques. The VIPER database is a powerful resource for virologists, microbiologists, virus crystallographers and EM researchers. This review describes how to use VIPER, using several examples to show the power of this resource for research and educational purposes.

Pseudo-atomic models of swollen CCMV from cryo-electron microscopy data

The capsid of cowpea chlorotic mottle virus (CCMV) can reversibly switch between two forms that are contingent on the charge of acidic residues that are clustered at the quasi-threefold axes of the T=3 icosahedral particle. The quaternary structure conformations are dependent on divalent metal ions and pH and were previously analyzed by crystallography in the native, compact form, and by cryo-electron microscopy in the compact and swollen forms (Speir et al., 1995). In this report we use the atomic models of the three structurally unique viral subunits determined by crystallography for a detailed interpretation of the 28-A-resolution electron density of the swollen form and the production of a pseudo-atomic model of this particle. The model of the quaternary structure conforms with high fidelity to conventional geometric constraints, quasi-equivalence, intersubunit association energies, and the electron density. It was derived by conserving the pentamers and hexamers of subunits whose associated electron densities are strikingly similar in the two forms of the particles. Treating these as rigid units in the modeling implies that the particle flexibility is accommodated primarily by changes in dimer interactions, an observation that is consistent with the flexible C-terminal polypeptide extensions that stabilize this contact in the crystal structure. Because the hexamers and pentamers were incrementally translated and rotated in a screw motion, with energy minimization at each of 28 steps, a path for the expansion is also implied.