Periodic Table of Virus Capsids: Implications for Natural Selection and Design

Background: For survival, most natural viruses depend upon the existence of spherical capsids: protective shells of various sizes composed of protein subunits. So far, general evolutionary pressures shaping capsid design have remained elusive, even though an understanding of such properties may help in rationally impeding the virus life cycle and designing efficient nano-assemblies. Principal Findings: This report uncovers an unprecedented and species-independent evolutionary pressure on virus capsids, based on the the notion that the simplest capsid designs (or those capsids with the lowest ‘‘hexamer complexity’’, C h) are the fittest, which was shown to be true for all available virus capsids. The theories result in a physically meaningful periodic table of virus capsids that uncovers strong and overarching evolutionary pressures, while also offering geometric explanations to other capsid properties (rigidity, pleomorphy, auxiliary requirements, etc.) that were previously considered to be unrelatable properties of the individual virus. Significance: Apart from describing a universal rule for virus capsid evolution, our work (especially the periodic table) provides a language with which highly diverse virus capsids, unified only by geometry, may be described and related to each other. Finally, the available virus structure databases and other published data reiterate the predicted geometryderive

The BackMAP Python module: how a simpler Ramachandran number can simplify the life of a protein simulator

Protein backbones occupy diverse conformations, but compact metrics to describe such conformations and transitions between them have been missing. This report re-introduces the Ramachandran number (ℛ) as a residue-level structural metric that could simply the life of anyone contending with large numbers of protein backbone conformations (e.g., ensembles from NMR and trajectories from simulations). Previously, the Ramachandran number (ℛ) was introduced using a complicated closed form, which made the Ramachandran number difficult to implement. This report discusses a much simpler closed form of ℛ that makes it much easier to calculate, thereby making it easy to implement. Additionally, this report discusses how ℛ dramatically reduces the dimensionality of the protein backbone, thereby making it ideal for simultaneously interrogating large numbers of protein structures. For example, 200 distinct conformations can easily be described in one graphic using ℛ (rather than 200 distinct Ramachandran plots). Finally, a new Python-based backbone analysis tool—BackMAP—is introduced, which reiterates how ℛ can be used as a simple and succinct descriptor of protein backbones and their dynamics

A Universal Trend among Proteomes Indicates an Oily Last Common Ancestor

Despite progresses in ancestral protein sequence reconstruction, much needs to be unraveled about the nature of the putative last common ancestral proteome that served as the prototype of all extant lifeforms. Here, we present data that indicate a steady decline (oil escape) in proteome hydrophobicity over species evolvedness (node number) evident in 272 diverse proteomes, which indicates a highly hydrophobic (oily) last common ancestor (LCA). This trend, obtained from simple considerations (free from sequence reconstruction methods), was corroborated by regression studies within homologous and orthologous protein clusters as well as phylogenetic estimates of the ancestral oil content. While indicating an inherent irreversibility in molecular evolution, oil escape also serves as a rare and universal reaction-coordinate for evolution (reinforcing Darwin's principle of Common Descent), and may prove important in matters such as (i) explaining the emergence of intrinsically disordered proteins, (ii) developing composition- and speciation-based “global” molecular clocks, and (iii) improving the statistical methods for ancestral sequence reconstruction.Chemistry and Chemical Biolog

How simple can a model of an empty viral capsid be? Charge distributions in viral capsids

We investigate and quantify salient features of the charge distributions on viral capsids. Our analysis combines the experimentally determined capsid geometry with simple models for ionization of amino acids, thus yielding the detailed description of spatial distribution for positive and negative charge across the capsid wall. The obtained data is processed in order to extract the mean radii of distributions, surface charge densities and dipole moment densities. The results are evaluated and examined in light of previously proposed models of capsid charge distributions, which are shown to have to some extent limited value when applied to real viruses.Comment: 10 pages, 10 figures; accepted for publication in Journal of Biological Physic

Figure 9: Alternative representations of the Ramachandran plot.

The Ramachandran plot is important to structural biology as it describes a peptide backbone in the context of its dominant degrees of freedom—the backbone dihedral angles φ and ψ (Ramachandran, Ramakrishnan & Sasisekharan, 1963). Since its introduction, the Ramachandran plot has been a crucial tool to characterize protein backbone features. However, the conformation or twist of a backbone as a function of φ and ψ has not been completely described for both cis and trans backbones. Additionally, little intuitive understanding is available about a peptide’s conformation simply from knowing the φ and ψ values of a peptide (e.g., is the regular peptide defined by φ = ψ =  − 100°  left-handed or right-handed?). This report provides a new metric for backbone handedness (h) based on interpreting a peptide backbone as a helix with axial displacement d and angular displacement θ, both of which are derived from a peptide backbone’s internal coordinates, especially dihedral angles φ, ψ and ω. In particular, h equals sin(θ)d∕|d|, with range [−1, 1] and negative (or positive) values indicating left(or right)-handedness. The metric h is used to characterize the handedness of every region of the Ramachandran plot for both cis (ω = 0°) and trans (ω = 180°) backbones, which provides the first exhaustive survey of twist handedness in Ramachandran (φ, ψ) space. These maps fill in the ‘dead space’ within the Ramachandran plot, which are regions that are not commonly accessed by structured proteins, but which may be accessible to intrinsically disordered proteins, short peptide fragments, and protein mimics such as peptoids. Finally, building on the work of (Zacharias & Knapp, 2013), this report presents a new plot based on d and θ that serves as a universal and intuitive alternative to the Ramachandran plot. The universality arises from the fact that the co-inhabitants of such a plot include every possible peptide backbone including cis and trans backbones. The intuitiveness arises from the fact that d and θ provide, at a glance, numerous aspects of the backbone including compactness, handedness, and planarity

Mechanical and Assembly Units of Viral Capsids Identified via Quasi-Rigid Domain Decomposition

Key steps in a viral life-cycle, such as self-assembly of a protective protein container or in some cases also subsequent maturation events, are governed by the interplay of physico-chemical mechanisms involving various spatial and temporal scales. These salient aspects of a viral life cycle are hence well described and rationalised from a mesoscopic perspective. Accordingly, various experimental and computational efforts have been directed towards identifying the fundamental building blocks that are instrumental for the mechanical response, or constitute the assembly units, of a few specific viral shells. Motivated by these earlier studies we introduce and apply a general and efficient computational scheme for identifying the stable domains of a given viral capsid. The method is based on elastic network models and quasi-rigid domain decomposition. It is first applied to a heterogeneous set of well-characterized viruses (CCMV, MS2, STNV, STMV) for which the known mechanical or assembly domains are correctly identified. The validated method is next applied to other viral particles such as L-A, Pariacoto and polyoma viruses, whose fundamental functional domains are still unknown or debated and for which we formulate verifiable predictions. The numerical code implementing the domain decomposition strategy is made freely available

Mechanisms of Size Control and Polymorphism in Viral Capsid Assembly

We simulate the assembly dynamics of icosahedral capsids from subunits that interconvert between different conformations (or quasi-equivalent states). The simulations identify mechanisms by which subunits form empty capsids with only one morphology but adaptively assemble into different icosahedral morphologies around nanoparticle cargoes with varying sizes, as seen in recent experiments with brome mosaic virus (BMV) capsid proteins. Adaptive cargo encapsidation requires moderate cargo-subunit interaction strengths; stronger interactions frustrate assembly by stabilizing intermediates with incommensurate curvature. We compare simulation results to experiments with cowpea chlorotic mottle virus empty capsids and BMV capsids assembled on functionalized nanoparticles and suggest new cargo encapsidation experiments. Finally, we find that both empty and templated capsids maintain the precise spatial ordering of subunit conformations seen in the crystal structure even if interactions that preserve this arrangement are favored by as little as the thermal energy, consistent with experimental observations that different subunit conformations are highly similar

Origination of the Protein Fold Repertoire from Oily Pluripotent Peptides

While the repertoire of protein folds that exists today underlies most of life’s capabilities, our mechanistic picture of protein fold origination is incomplete. This paper discusses a hypothetical mechanism for the emergence of the protein fold repertoire from highly dynamic and collapsed peptides, exemplified by peptides with high oil content or hydrophobicity. These peptides are called pluripotent to emphasize their capacity to evolve into numerous folds transiently available to them. As evidence, the paper will discuss previous simulation work on the superior fold evolvability of oily peptides, trace (“fossil”) evidence within proteomes seen today, and a general relationship between protein dynamism and evolvability. Aside from implications on the origination of protein folds, the hypothesis implies that the vanishing utility of a random peptide in protein origination may be relatively exaggerated, as some random peptides with a certain composition (e.g., oily) may fare better than others. In later sections, the hypothesis is discussed in the context of existing discussions regarding the spontaneous origination of biomolecules