Geometric considerations in virus capsid size specificity, auxiliary requirements, and buckling

Modeling reveals the strength of weak interactions in stacked-ring assembly

Cells employ many large macromolecular machines for the execution and regulation of processes that are vital for cell and organismal viability. Interestingly, cells cannot synthesize these machines as functioning units. Instead, cells synthesize the molecular parts that must then assemble into the functional complex. Many important machines, including chaperones such as GroEL and proteases such as the proteasome, comprise protein rings that are stacked on top of one another. While there is some experimental data regarding how stacked-ring complexes such as the proteasome self-assemble, a comprehensive understanding of the dynamics of stacked-ring assembly is currently lacking. Here, we developed a mathematical model of stacked-trimer assembly and performed an analysis of the assembly of the stacked homomeric trimer, which is the simplest stacked-ring architecture. We found that stacked rings are particularly susceptible to a form of kinetic trapping that we term "deadlock," in which the system gets stuck in a state where there are many large intermediates that are not the fully assembled structure but that cannot productively react. When interaction affinities are uniformly strong, deadlock severely limits assembly yield. We thus predicted that stacked rings would avoid situations where all interfaces in the structure have high affinity. Analysis of available crystal structures indicated that indeed the majority-if not all-of stacked trimers do not contain uniformly strong interactions. Finally, to better understand the origins of deadlock, we developed a formal pathway analysis and showed that, when all the binding affinities are strong, many of the possible pathways are utilized. In contrast, optimal assembly strategies utilize only a small number of pathways. Our work suggests that deadlock is a critical factor influencing the evolution of macromolecular machines and provides general principles for understanding the self-assembly efficiency of existing machines.

α-Carboxysome Size Is Controlled by the Disordered Scaffold Protein CsoS2

Carboxysomes are protein microcompartments that function in the bacterial CO2 concentrating mechanism (CCM) to facilitate CO2 assimilation. To do so, carboxysomes assemble from thousands of constituent proteins into an icosahedral shell, which encapsulates the enzymes Rubisco and carbonic anhydrase to form structures typically > 100 nm and > 300 megadaltons. Although many of the protein interactions driving the assembly process have been determined, it remains unknown how size and composition are precisely controlled. Here, we show that the size of α-carboxysomes is controlled by the disordered scaffolding protein CsoS2. CsoS2 contains two classes of related peptide repeats that bind to the shell in a distinct fashion, and our data indicate that size is controlled by the relative number of these interactions. We propose an energetic and structural model wherein the two repeat classes bind at the junction of shell hexamers but differ in their preferences for the shell contact angles, and thus the local curvature. In total, this model suggests that a set of specific and repeated interactions between CsoS2 and shell proteins collectively achieve the large size and monodispersity of α-carboxysomes.

Principles for enhancing virus capsid capacity and stability from a thermophilic virus capsid structure

The capsids of double-stranded DNA viruses protect the viral genome from the harsh extracellular environment, while maintaining stability against the high internal pressure of packaged DNA. To elucidate how capsids maintain stability in an extreme environment, we use cryoelectron microscopy to determine the capsid structure of thermostable phage P74-26 to 2.8-Å resolution. We find P74-26 capsids exhibit an overall architecture very similar to those of other tailed bacteriophages, allowing us to directly compare structures to derive the structural basis for enhanced stability. Our structure reveals lasso-like interactions that appear to function like catch bonds. This architecture allows the capsid to expand during genome packaging, yet maintain structural stability. The P74-26 capsid has T = 7 geometry despite being twice as large as mesophilic homologs. Capsid capacity is increased with a larger, flatter major capsid protein. Given these results, we predict decreased icosahedral complexity (i.e. T ≤ 7) leads to a more stable capsid assembly.

Viral capsids need to protect the genome against harsh environmental conditions and cope with high internal pressure from the packaged genome. Here, the authors determine the structure of the thermostable phage P74-26 capsid at 2.8-Å resolution and identify features underlying enhanced capsid capacity and structural stability.

Controlled reversible buckling of polydopamine spherical microcapsules: revealing the hidden rich phenomena of post-buckling of spherical polymeric shells

Under external pressure compression, various kinds of artificial microcapsules can undergo buckling induced deformation and catastrophic rupturing failure, which needs to be understood for their diverse practical applications. For this, many theories and numerical simulations have recently emerged, leading to some intriguing but often debatable predictions and scaling laws. However, experimental testing of these predictions is very limited, due to challenges in realizing prescribed buckling pathways and in situ monitoring of the buckling procedure. Herein, we report the buckling behaviors of well-defined spherical polydopamine (PDA) capsules with tunable sizes and homogeneous nanoscale shells. Simple but controlled solvent evaporation was implemented inside a home-made optical chamber to induce buckling of PDA capsules by following a prescribed pathway toward targeted shapes that are only dictated by the inherent material properties of the capsules. In addition, the buckling speed was slowed down to the timescale of minutes, which can prevent buckling from being trapped at some metastable intermediate states as well as facilitating in situ optical monitoring of the whole buckling procedure in slow motion. In this way, several classic buckling behaviors were clearly observed, including the sudden appearance of spinodal-like dimples above critical pressures, transition of the indentation rim from the axisymmetric to polygonal shape, and evolution of multi-indented buckling into single indented buckling following Ostwald ripening. These observations are qualitatively comparable with recent predictions from numerical results. Furthermore, some novel buckling phenomena have been reported for the first time, which might stimulate further theories and numerical simulations.

Stepwise reversible nanomechanical buckling in a viral capsid

Viruses are nanoscale infectious agents constructed of a proteinaceous capsid that protects the packaged genomic material. Nanoindentation experiments using atomic force microscopy have, in recent years, provided unprecedented insight into the elastic properties, structural stability and maturation-dependent mechanical changes in viruses. However, the dynamics of capsid behavior are still unresolved. Here we used high-resolution nanoindentation experiments on mature, DNA-filled T7 bacteriophage particles. The elastic regime of the nanoindentation force trace contained discrete, stepwise transitions that cause buckling of the T7 capsid with magnitudes that are integer multiples of ∼0.6 nm. Remarkably, the transitions are reversible and contribute to the rapid consolidation of the capsid structure against a force during cantilever retraction. The stepwise transitions were present even following the removal of the genomic DNA by heat treatment, indicating that they are related to the structure and dynamics of the capsomeric proteins. Dynamic force spectroscopy experiments revealed that the thermally activated consolidation step is ∼10⁴ times faster than spontaneous buckling, suggesting that the capsid stability is under strong dynamic control. Capsid structural dynamics may play an important role in protecting the genomic material from harsh environmental impacts. The nanomechanics approach employed here may be used to investigate the structural dynamics of other viruses and nanoscale containers as well.

Hidden symmetry of small spherical viruses and organization principles in "anomalous" and double-shelled capsid nanoassemblies

We propose the principles of structural organization in spherical nanoassemblies with icosahedral symmetry constituted by asymmetric protein molecules. The approach modifies the paradigmatic geometrical Caspar and Klug (CK) model of icosahedral viral capsids and demonstrates the common origin of both the "anomalous" and conventional capsid structures. In contrast to all previous models of "anomalous" viral capsids the proposed modified model conserves the basic structural principles of the CK approach and reveals the common hidden symmetry underlying all small viral shells. We demonstrate the common genesis of the "anomalous" and conventional capsids and explain their structures in the same frame. The organization principles are derived from the group theory analysis of the positional order on the spherical surface. The relationship between the modified CK geometrical model and the theory of two-dimensional spherical crystallization is discussed. We also apply the proposed approach to complex double-shelled capsids and capsids with protruding knob-like proteins. The introduced notion of commensurability for the concentric nanoshells explains the peculiarities of their organization and helps to predict analogous, but yet undiscovered, double-shelled viral capsid nanostructures.

Kinetics versus Thermodynamics in Virus Capsid Polymorphism

Virus coat proteins spontaneously self-assemble into empty shells in aqueous solution under the appropriate physicochemical conditions, driven by an interaction free energy per bond on the order of 2-5 times the thermal energy kBT. For this seemingly modest interaction strength, each protein building block nonetheless gains a very large binding free energy, between 10 and 20 kBT. Because of this, there is debate about whether the assembly process is reversible or irreversible. Here we discuss capsid polymorphism observed in in vitro experiments from the perspective of nucleation theory and of the thermodynamics of mass action. We specifically consider the potential contribution of a curvature free energy term to the effective interaction potential between the proteins. From these models, we propose experiments that may conclusively reveal whether virus capsid assembly into a mixture of polymorphs is a reversible or an irreversible process.

Hierarchical assembly may be a way to make large information-rich structures

Self-assembly in the laboratory can now yield 'information-rich' nanostructures in which each component is of a distinct type and has a defined spatial position. Ensuring the thermodynamic stability of such structures requires inter-component interaction energies to increase logarithmically with structure size, in order to counter the entropy gained upon mixing component types in solution. However, self-assembly in the presence of strong interactions results in general in kinetic trapping, so suggesting a limit to the size of an (equilibrium) structure that can be self-assembled from distinguishable components. Here we study numerically a two-dimensional hierarchical assembly scheme already considered in experiment. We show that this scheme is immune to the kinetic traps associated with strong 'native' interactions (interactions designed to stabilize the intended structure), and so, in principle, offers a way to make large information-rich structures. In this scheme the size of an assembled structure scales exponentially with the stage of assembly, and assembly can continue as long as random motion is able to bring structures into contact. The resulting superstructure could provide a template for building in the third dimension. The chief drawback of this scheme is that it is particularly susceptible to kinetic traps that result from 'non-native' interactions (interactions not required to stabilize the intended structure); the scale on which such a scheme can be realized therefore depends upon how effectively this latter kind of interaction can be suppressed.

Protein-Protein Interfaces in Viral Capsids Are Structurally Unique

Viral capsids exhibit elaborate and symmetrical architectures of defined sizes and remarkable mechanical properties not seen with cellular macromolecular complexes. Given the uniqueness of the higher-order organization of viral capsid proteins in the virosphere, we explored the question of whether the patterns of protein-protein interactions within viral capsids are distinct from those in generic protein complexes. Our comparative analysis involving a non-redundant set of 551 inter-subunit interfaces in viral capsids from VIPERdb and 20,014 protein-protein interfaces in non-capsid protein complexes from the Protein Data Bank found 418 generic protein-protein interfaces that share similar physicochemical patterns with some protein-protein interfaces in the capsid set, using the program PCalign we developed for comparing protein-protein interfaces. This overlap in the structural space of protein-protein interfaces is significantly small, with a p-value <0.0001, based on a permutation test on the total set of protein-protein interfaces. Furthermore, the generic protein-protein interfaces that bear similarity in their spatial and chemical arrangement with capsid ones are mostly small in size with fewer than 20 interfacial residues, which results from the relatively limited choices of natural design for small interfaces rather than having significant biological implications in terms of functional relationships. We conclude based on this study that protein-protein interfaces in viral capsids are non-representative of patterns in the smaller, more compact cellular protein complexes. Our finding highlights the design principle of building large biological containers from repeated, self-assembling units and provides insights into specific targets for antiviral drug design for improved efficacy.

Modeling Viral Capsid Assembly

I present a review of the theoretical and computational methodologies that have been used to model the assembly of viral capsids. I discuss the capabilities and limitations of approaches ranging from equilibrium continuum theories to molecular dynamics simulations, and I give an overview of some of the important conclusions about virus assembly that have resulted from these modeling efforts. Topics include the assembly of empty viral shells, assembly around single-stranded nucleic acids to form viral particles, and assembly around synthetic polymers or charged nanoparticles for nanotechnology or biomedical applications. I present some examples in which modeling efforts have promoted experimental breakthroughs, as well as directions in which the connection between modeling and experiment can be strengthened.

DNA buckling in bacteriophage cavities as a mechanism to aid virus assembly

While relatively simple biologically, bacteriophages are sophisticated biochemical machines that execute a precise sequence of events during virus assembly, DNA packaging, and ejection. These stages of the viral life cycle require intricate coordination of viral components whose structures are being revealed by single molecule experiments and high resolution (cryo-electron microscopy) reconstructions. For example, during packaging, bacteriophages employ some of the strongest known molecular motors to package DNA against increasing pressure within the viral capsid shell. Located upstream of the motor is an elaborate portal system through which DNA is threaded. A high resolution reconstruction of the portal system for bacteriophage ϕ29 reveals that DNA buckles inside a small cavity under large compressive forces. In this study, we demonstrate that DNA can also buckle in other bacteriophages including T7 and P22. Using a computational rod model for DNA, we demonstrate that a DNA buckle can initiate and grow within the small confines of a cavity under biologically-attainable force levels. The forces of DNA-cavity contact and DNA-DNA electrostatic repulsion ultimately limit cavity filling. Despite conforming to very different cavity geometries, the buckled DNA within T7 and P22 exhibits near equal volumetric energy density (∼1kT/nm(3)) and energetic cost of packaging (∼22kT). We hypothesize that a DNA buckle creates large forces on the cavity interior to signal the conformational changes to end packaging. In addition, a DNA buckle may help retain the genome prior to tail assembly through significantly increased contact area with the portal.

Recognition-driven assembly of self-limiting supramolecular protein nanoparticles displaying enzymatic activity

We report the recognition-driven assembly of self-limiting protein nanoparticles displaying enzymatic activity.

Building polyhedra by self-assembly: theory and experiment

We investigate the utility of a mathematical framework based on discrete geometry to model biological and synthetic self-assembly. Our primary biological example is the self-assembly of icosahedral viruses; our synthetic example is surface-tension-driven self-folding polyhedra. In both instances, the process of self-assembly is modeled by decomposing the polyhedron into a set of partially formed intermediate states. The set of all intermediates is called the configuration space, pathways of assembly are modeled as paths in the configuration space, and the kinetics and yield of assembly are modeled by rate equations, Markov chains, or cost functions on the configuration space. We review an interesting interplay between biological function and mathematical structure in viruses in light of this framework. We discuss in particular: (i) tiling theory as a coarse-grained description of all-atom models; (ii) the building game-a growth model for the formation of polyhedra; and (iii) the application of these models to the self-assembly of the bacteriophage MS2. We then use a similar framework to model self-folding polyhedra. We use a discrete folding algorithm to compute a configuration space that idealizes surface-tension-driven self-folding and analyze pathways of assembly and dominant intermediates. These computations are then compared with experimental observations of a self-folding dodecahedron with side 300 μm. In both models, despite a combinatorial explosion in the size of the configuration space, a few pathways and intermediates dominate self-assembly. For self-folding polyhedra, the dominant intermediates have fewer degrees of freedom than comparable intermediates, and are thus more rigid. The concentration of assembly pathways on a few intermediates with distinguished geometric properties is biologically and physically important, and suggests deeper mathematical structure.

A minimal representation of the self-assembly of virus capsids

Viruses are biological nanosystems with a capsid of protein-made capsomer units that encloses and protects the genetic material responsible for their replication. Here we show how the geometrical constraints of the capsomer-capsomer interaction in icosahedral capsids and the requirement of low frustration fix the form of the shortest and universal truncated multipolar expansion of the two-body interaction between capsomers. The structures of many of the icosahedral and related virus capsids are located as single lowest energy states of a potential energy surface built from this interaction. Our minimalist representation is consistent with other models known to produce a controllable and efficient self-assembly, and unveils relevant features of the natural design of the capsids. It promises to be very useful in physical virology and may also be of interest in fields of nanoscience and nanotechnology where similar hollow convex structures are relevant.

Recent Developments in Molecular Simulation Approaches to Study Spherical Virus Capsids

Viruses are a particularly challenging systems to study via molecular simulation methods. Virus capsids typically consist of over 100 subunit proteins and reach dimensions of over 100 nm; solvated viruses capsid systems can be over 1 million atoms in size. In this review, I will present recent developments which have attempted to overcome the significant computational expense to perform simulations which can inform experimental studies, make useful predictions about biological phenomena and calculate material properties relevant to nanotechnology design efforts.

On the morphology of viral capsids: elastic properties and buckling transitions

The morphology of icosahedral viruses ranges from highly spherical to highly faceted, and for some viruses a shape transition occurs during the viral life cycle. This phenomena is predicted from continuum elasticity, via the buckling transition theory by Nelson (Phys. Rev. E 2003, 68, 051910), in which the shape is dependent on the Foppl-von Kármán number (γ), which is a ratio of the two-dimensional Young's modulus (Y) and the bending modulus (κ). However, until now, no direct calculations have been performed on atomic-level capsid structures to test the predictions of the theory. In this study, we employ a previously described multiscale method by May and Brooks (Phys. Rev. Lett. 2011, 106, 188101) to calculate Y and κ for the bacteriophage HK97, which undergoes a spherical to faceted transition during its viral life cycle. We observe a change in γ consistent with the buckling transition theory and also a significant reduction in κ, which facilitates formation of the faceted state. We go on to examine many capsids from the T = 3 and 7 classes using only elastic network models, which allows us to calculate the ratio Y/κ, without the expense of all-atom molecular dynamics. We observe for the T = 7 capsids, there is strong correlation between the shape of the capsid and γ; however, there is no such correlation for the smaller T = 3 viruses.

Exploring the symmetry and mechanism of virus capsid maturation via an ensemble of pathways

Many icosahedral viruses undergo large-scale conformational transitions between icosahedrally symmetric conformations during their life cycles. However, whether icosahedral symmetry is maintained along the transition pathways for this process is unknown. By employing a simplified and directed structure-based potential we compute an ensemble of transition pathways for the maturation transition of bacteriophage HK97. We observe localized symmetry-breaking events, but find that the large-scale displacements are dominated by icosahedrally symmetric deformation modes. We find that all pathways obey a common mechanism characterized by formation of pentameric contacts early in the transition.

Simulated self-assembly of the HIV-1 capsid: protein shape and native contacts are sufficient for two-dimensional lattice formation

We report Monte Carlo simulations of the initial stages of self-assembly of the HIV-1 capsid protein (CA), using a coarse-grained representation that mimics the CA backbone structure and intermolecular contacts observed experimentally. A simple representation of N-terminal domain/N-terminal domain and N-terminal domain/C-terminal domain interactions, coupled with the correct protein shape, is sufficient to drive formation of an ordered lattice with the correct hexagonal symmetry in two dimensions. We derive an approximate concentration/temperature phase diagram for lattice formation, and we investigate the pathway by which the lattice develops from initially separated CA dimers. Within this model, lattice formation occurs in two stages: 1), condensation of CA dimers into disordered clusters; and 2), nucleation of the lattice by the appearance of one hexamer unit within a cluster. Trimers of CA dimers are important early intermediates, and pentamers are metastable within clusters. Introduction of a preformed hexamer at the beginning of a Monte Carlo run does not directly seed lattice formation, but does facilitate the formation of large clusters. We discuss possible connections between these simulations and experimental observations concerning CA assembly within HIV-1 and in vitro.

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 geometry-derived rules, reinforcing the role of geometry in the natural selection and design of virus capsids.