Energy landscapes for shells assembled from pentagonal and hexagonal pyramids

Effect of coat-protein concentration on the self-assembly of bacteriophage MS2 capsids around RNA

Self-assembly is a vital part of the life cycle of certain icosahedral RNA viruses. Furthermore, the assembly process can be harnessed to make icosahedral virus-like particles (VLPs) from coat protein and RNA in vitro. Although much previous work has explored the effects of RNA-protein interactions on the assembly products, relatively little research has explored the effects of coat-protein concentration. We mix coat protein and RNA from bacteriophage MS2, and we use a combination of gel electrophoresis, dynamic light scattering, and transmission electron microscopy to investigate the assembly products. We show that with increasing coat-protein concentration, the products transition from well-formed MS2 VLPs to "monster" particles consisting of multiple partial capsids to RNA-protein condensates consisting of large networks of RNA and partially assembled capsids. We argue that the transition from well-formed to monster particles arises because the assembly follows a nucleation-and-growth pathway in which the nucleation rate depends sensitively on the coat-protein concentration, such that at high protein concentrations, multiple nuclei can form on each RNA strand. To understand the formation of the condensates, which occurs at even higher coat-protein concentrations, we use Monte Carlo simulations with coarse-grained models of capsomers and RNA. These simulations suggest that the formation of condensates occurs by the adsorption of protein to the RNA followed by the assembly of capsids. Multiple RNA molecules can become trapped when a capsid grows from capsomers attached to two different RNA molecules or when excess protein bridges together growing capsids on different RNA molecules. Our results provide insight into an important biophysical process and could inform design rules for making VLPs for various applications.

Virtual indentation of the empty capsid of the minute virus of mice using a minimal coarse-grained model

A minimal coarse-grained model for T=1 viral capsids assembled from 20 protein rigid trimers has been designed by extending a previously proposed form of the interaction energy written as a sum of anisotropic pairwise interactions between the trimeric capsomers. The extension of the model has been performed to properly account for the coupling between two internal coordinates: the one that measures the intercapsomer distance and the other that gives the intercapsomer dihedral angle. The model has been able to fit with less than a 10% error the atomic force microscopy (AFM) indentation experimental data for the empty capsid of the minute virus of mice (MVM), providing in this way an admissible picture of the main mechanisms behind the capsid deformations. In this scenario, the bending of the intercapsomer dihedral angle is the angular internal coordinate that can support larger deformations away from its equilibrium values, determining important features of the AFM indentation experiments as the elastic constants along the three symmetry axes of the capsid and the critical indentations. From the value of one of the parameters of our model, we conclude that trimers in the MVM must be quite oblate tops, in excellent agreement with their known structure. The transition from the linear to the nonlinear regimes sampled in the indentation process appears to be an interesting topic for future research in physical virology.

Combining experiment and energy landscapes to explore anaerobic heme breakdown in multifunctional hemoproteins

To survive, many pathogens extract heme from their host organism and break down the porphyrin scaffold to sequester the Fe2+ ion via a heme oxygenase. Recent studies have revealed that certain pathogens can anaerobically degrade heme. Our own research has shown that one such pathway proceeds via NADH-dependent heme degradation, which has been identified in a family of hemoproteins from a range of bacteria. HemS, from Yersinia enterocolitica, is the main focus of this work, along with HmuS (Yersinia pestis), ChuS (Escherichia coli) and ShuS (Shigella dysenteriae). We combine experiments, Energy Landscape Theory, and a bioinformatic investigation to place these homologues within a wider phylogenetic context. A subset of these hemoproteins are known to bind certain DNA promoter regions, suggesting not only that they can catalytically degrade heme, but that they are also involved in transcriptional modulation responding to heme flux. Many of the bacterial species responsible for these hemoproteins (including those that produce HemS, ChuS and ShuS) are known to specifically target oxygen-depleted regions of the gastrointestinal tract. A deeper understanding of anaerobic heme breakdown processes exploited by these pathogens could therefore prove useful in the development of future strategies for disease prevention.

Effect of ionic strength on the assembly of simian vacuolating virus capsid protein around poly(styrene sulfonate)

Virus-like particles (VLPs) are noninfectious nanocapsules that can be used for drug delivery or vaccine applications. VLPs can be assembled from virus capsid proteins around a condensing agent, such as RNA, DNA, or a charged polymer. Electrostatic interactions play an important role in the assembly reaction. VLPs assemble from many copies of capsid protein, with a combinatorial number of intermediates. Hence, the mechanism of the reaction is poorly understood. In this paper, we combined solution small-angle X-ray scattering (SAXS), cryo-transmission electron microscopy (TEM), and computational modeling to determine the effect of ionic strength on the assembly of Simian Vacuolating Virus 40 (SV40)-like particles. We mixed poly(styrene sulfonate) with SV40 capsid protein pentamers at different ionic strengths. We then characterized the assembly product by SAXS and cryo-TEM. To analyze the data, we performed Langevin dynamics simulations using a coarse-grained model that revealed incomplete, asymmetric VLP structures consistent with the experimental data. We found that close to physiological ionic strength, [Formula: see text] VLPs coexisted with VP1 pentamers. At lower or higher ionic strengths, incomplete particles coexisted with pentamers and [Formula: see text] particles. Including the simulated structures was essential to explain the SAXS data in a manner that is consistent with the cryo-TEM images.

Microcompartment assembly around multicomponent fluid cargoes

This article describes dynamical simulations of the assembly of an icosahedral protein shell around a bicomponent fluid cargo. Our simulations are motivated by bacterial microcompartments, which are protein shells found in bacteria that assemble around a complex of enzymes and other components involved in certain metabolic processes. The simulations demonstrate that the relative interaction strengths among the different cargo species play a key role in determining the amount of each species that is encapsulated, their spatial organization, and the nature of the shell assembly pathways. However, the shell protein–shell protein and shell protein–cargo component interactions that help drive assembly and encapsulation also influence cargo composition within certain parameter regimes. These behaviors are governed by a combination of thermodynamic and kinetic effects. In addition to elucidating how natural microcompartments encapsulate multiple components involved within reaction cascades, these results have implications for efforts in synthetic biology to colocalize alternative sets of molecules within microcompartments to accelerate specific reactions. More broadly, the results suggest that coupling between self-assembly and multicomponent liquid–liquid phase separation may play a role in the organization of the cellular cytoplasm.

Minimal Design Principles for Icosahedral Virus Capsids

The geometrical structures of single- and multiple-shell icosahedral virus capsids are reproduced as the targets that minimize the cost corresponding to relatively simple design functions. Capsid subunits are first identified as building blocks at a given coarse-grained scale and then represented in these functions as point particles located on an appropriate number of concentric spherical surfaces. Minimal design cost is assigned to optimal spherical packings of the particles. The cost functions are inspired by the packings favored for the Thomson problem, which minimize the electrostatic potential energy between identical charged particles. In some cases, icosahedral symmetry constraints are incorporated as external fields acting on the particles. The simplest cost functions can be obtained by separating particles in disjoint nonequivalent sets with distinct interactions, or by introducing interacting holes (the absence of particles). These functions can be adapted to reproduce any capsid structure found in real viruses. Structures absent in Nature require significantly more complex designs. Measures of information content and complexity are assigned to both the cost functions and the capsid geometries. In terms of these measures, icosahedral structures and the corresponding cost functions are the simplest solutions.

Assembly and Stability of Simian Virus 40 Polymorphs

Understanding viral assembly pathways is of critical importance to biology, medicine, and nanotechology. Here, we study the assembly path of a system with various structures, the simian vacuolating virus 40 (SV40) polymorphs. We simulate the templated assembly process of VP1 pentamers, which are the constituents of SV40, into icosahedal shells made of N = 12 pentamers (T = 1). The simulations include connections formed between pentamers by C-terminal flexible lateral units, termed here "C-terminal ligands", which are shown to control assembly behavior and shell dynamics. The model also incorporates electrostatic attractions between the N-terminal peptide strands (ligands) and the negatively charged cargo, allowing for agreement with experiments of RNA templated assembly at various pH and ionic conditions. During viral assembly, pentamers bound to any template increase its effective size due to the length and flexibility of the C-terminal ligands, which can connect to other VP1 pentamers and recruit them to a partially completed capsid. All closed shells formed other than the T = 1 feature the ability to dynamically rearrange and are thus termed "pseudo-closed". The N = 13 shell can even spontaneously "self-correct" by losing a pentamer and become a T = 1 capsid when the template size fluctuates. Bound pentamers recruiting additional pentamers to dynamically rearranging capsids allow closed shells to continue growing via the pseudo-closed growth mechanism, for which experimental evidence already exists. Overall, we show that the C-terminal ligands control the dynamic assembly paths of SV40 polymorphs.

A minimal coarse-grained model for the low-frequency normal mode analysis of icosahedral viral capsids

The main goal of this work is the design of a coarse-grained theoretical model of minimal resolution for the study of the physical properties of icosahedral virus capsids within the linear-response regime. In this model the capsid is represented as an interacting many-body system whose composing elements are capsid subunits (capsomers), which are treated as three-dimensional rigid bodies. The total interaction potential energy is written as a sum of pairwise capsomer-capsomer interactions. Based on previous work [Gomez Llorente et al., Soft Matter, 2014, 10, 3560], a minimal and complete anisotropic binary interaction that includes a full Hessian matrix of independent force constants is proposed. In this interaction model, capsomers have rotational symmetry around an axis of order n > 2. The full coarse-grained model is applied to analyse the low-frequency normal-mode spectrum of icosahedral T = 1 capsids. The model performance is evaluated by fitting its predicted spectrum to the full-atom results for the Satellite Tobacco Necrosis Virus (STNV) capsid [Dykeman and Sankey, Phys. Rev. Lett., 2008, 100, 028101]. Two capsomer choices that are compatible with the capsid icosahedral symmetry are checked, namely pentamers (n = 5) and trimers (n = 3). Both subunit types provide fair fits, from which the magnitude of the coarse-grained force constants for a real virus is obtained. The model is able to uncover latent instabilities whose analysis is fully consistent with the current knowledge about the STNV capsid, which does not self-assemble in the absence of RNA and is thermally unstable. The straightforward generalisability of the model beyond the linear regime and its completeness make it a promising tool to theoretically interpret many experimental data such as those provided by the atomic force microscopy or even to better understand processes far from equilibrium such as the capsid self-assembly.

Kinetics of empty viral capsid assembly in a minimal model

The efficient construction of a protective protein shell or capsid is one of the most crucial steps in the replication cycle of a virus. The formation of the simplest capsid typically proceeds by the spontaneous assembly of identical building blocks. This process can also be achieved in vitro even in the absence of genetic material, thus opening the door to the production of artificial viral cages for a myriad of applications. In this work, we analyze the efficiency and the kinetic peculiarities of this self-assembly process using Brownian Dynamics simulations. We use a minimal model that considers identical assembly units and is able to reproduce successfully the correct final architecture of spherical capsids. The selection of a specific size and structure is achieved by changing a single parameter that imposes an angular anisotropy on the interaction. We analyze how the geometrical constraints of the interaction affect the efficiency of the assembly. We find that the optimal conditions for an efficient assembly from a kinetic point of view strongly depart from the lowest capsid energy corresponding to the minimum of the potential energy landscape. Our work illustrates the important differences between the equilibrium and dynamic characteristics of viral self-assembly, and provides important insights on how to design specific interactions for a successful assembly of artificial viral cages.

The role of the encapsulated cargo in microcompartment assembly

Bacterial microcompartments are large, roughly icosahedral shells that assemble around enzymes and reactants involved in certain metabolic pathways in bacteria. Motivated by microcompartment assembly, we use coarse-grained computational and theoretical modeling to study the factors that control the size and morphology of a protein shell assembling around hundreds to thousands of molecules. We perform dynamical simulations of shell assembly in the presence and absence of cargo over a range of interaction strengths, subunit and cargo stoichiometries, and the shell spontaneous curvature. Depending on these parameters, we find that the presence of a cargo can either increase or decrease the size of a shell relative to its intrinsic spontaneous curvature, as seen in recent experiments. These features are controlled by a balance of kinetic and thermodynamic effects, and the shell size is assembly pathway dependent. We discuss implications of these results for synthetic biology efforts to target new enzymes to microcompartment interiors.

Multiple Pathways in Capsid Assembly

For a three-dimensional structure to spontaneously self-assemble from many identical components, the steps on the pathway must be kinetically accessible. Many virus capsids are icosahedral and assembled from hundreds of identical proteins, but how they navigate the assembly process is poorly understood. Capsid assembly is thought to involve stepwise addition of subunits to a growing capsid fragment. Coarse-grained models suggest that the reaction occurs on a downhill energy landscape, so intermediates are expected to be fleeting. In this work, charge detection mass spectrometry (CDMS) has been used to track assembly of the hepatitis B virus (HBV) capsid in real time. The icosahedral T = 4 capsid of HBV is assembled from 120 capsid protein dimers. Our results indicate that there are multiple pathways for assembly. Under conditions that favor a modest association energy there is no accumulation of large intermediates, which indicates that available pathways include ones on a downhill energy surface. Under higher salt conditions, where subunit interactions are strengthened, around half of the products of the initial assembly reaction have masses close to the T = 4 capsid and the other half are stalled intermediates which emerge abruptly at around 90 dimers, indicating a bifurcation in the ensemble of assembly paths. When incubated at room temperature, the 90-dimer intermediates accumulate dimers and gradually shift to higher mass and merge with the capsid peak. Though free subunits are present in solution, the stalled intermediates indicate the presence of a local minima on the energy landscape. Some intermediates may result from hole closure, where the growing capsid distorts to close the hole due to the missing capsid proteins or from a species where subsequent additions are particularly labile.

Theoretical approaches for dynamical ordering of biomolecular systems

BACKGROUND

Living systems are characterized by the dynamic assembly and disassembly of biomolecules. The dynamical ordering mechanism of these biomolecules has been investigated both experimentally and theoretically. The main theoretical approaches include quantum mechanical (QM) calculation, all-atom (AA) modeling, and coarse-grained (CG) modeling. The selected approach depends on the size of the target system (which differs among electrons, atoms, molecules, and molecular assemblies). These hierarchal approaches can be combined with molecular dynamics (MD) simulation and/or integral equation theories for liquids, which cover all size hierarchies.

SCOPE OF REVIEW

We review the framework of quantum mechanical/molecular mechanical (QM/MM) calculations, AA MD simulations, CG modeling, and integral equation theories. Applications of these methods to the dynamical ordering of biomolecular systems are also exemplified.

MAJOR CONCLUSIONS

The QM/MM calculation enables the study of chemical reactions. The AA MD simulation, which omits the QM calculation, can follow longer time-scale phenomena. By reducing the number of degrees of freedom and the computational cost, CG modeling can follow much longer time-scale phenomena than AA modeling. Integral equation theories for liquids elucidate the liquid structure, for example, whether the liquid follows a radial distribution function.

GENERAL SIGNIFICANCE

These theoretical approaches can analyze the dynamic behaviors of biomolecular systems. They also provide useful tools for exploring the dynamic ordering systems of biomolecules, such as self-assembly. This article is part of a Special Issue entitled "Biophysical Exploration of Dynamical Ordering of Biomolecular Systems" edited by Dr. Koichi Kato.

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.

Physical Ingredients Controlling Stability and Structural Selection of Empty Viral Capsids

One of the crucial steps in the viral replication cycle is the self-assembly of its protein shell. Typically, each native virus adopts a unique architecture, but the coat proteins of many viruses have the capability to self-assemble in vitro into different structures by changing the assembly conditions. However, the mechanisms determining which of the possible capsid shapes and structures is selected by a virus are still not well-known. We present a coarse-grained model to analyze and understand the physical mechanisms controlling the size and structure selection in the assembly of empty viral capsids. Using this model and Monte Carlo simulations, we have characterized the phase diagram and stability of T = 1,3,4,7 and snub cube shells. In addition, we have studied the tolerance of different shells to changes in physical parameters related to ambient conditions, identifying possible strategies to induce misassembly or failure. Finally, we discuss the factors that select the shape of a capsid as spherical, faceted, elongated, or decapsidated. Our model sheds important light on the ingredients that control the assembly and stability of viral shells. This knowledge is essential to get capsids with well-defined size and structure that could be used for promising applications in medicine or bionanotechnology.

Many-molecule encapsulation by an icosahedral shell

We computationally study how an icosahedral shell assembles around hundreds of molecules. Such a process occurs during the formation of the carboxysome, a bacterial microcompartment that assembles around many copies of the enzymes ribulose 1,5-bisphosphate carboxylase/ oxygenase and carbonic anhydrase to facilitate carbon fixation in cyanobacteria. Our simulations identify two classes of assembly pathways leading to encapsulation of many-molecule cargoes. In one, shell assembly proceeds concomitantly with cargo condensation. In the other, the cargo first forms a dense globule; then, shell proteins assemble around and bud from the condensed cargo complex. Although the model is simplified, the simulations predict intermediates and closure mechanisms not accessible in experiments, and show how assembly can be tuned between these two pathways by modulating protein interactions. In addition to elucidating assembly pathways and critical control parameters for microcompartment assembly, our results may guide the reengineering of viruses as nanoreactors that self-assemble around their reactants.

DOI:http://dx.doi.org/10.7554/eLife.14078.001

Bacterial microcompartments are protein shells that are found inside bacteria and enclose enzymes and other chemicals required for certain biological reactions. For example, the carboxysome is a type of microcompartment that enables the bacteria to convert the products of photosynthesis into sugars. During the formation of a microcompartment, the outer protein shell assembles around hundreds of enzymes and chemicals. This formation process is tightly controlled and involves multiple interactions between the shell proteins and the cargo – the enzymes and other reaction ingredients – they will enclose. Understanding how to control which enzymes are encapsulated within microcompartments could help researchers to re-engineer the microcompartments so that they contain drugs or other useful products.

Recent studies have used microscopy to visualize how microcompartments are assembled. However, most of the intermediate structures that form during assembly are too small and short-lived to be seen. It has therefore not been possible to explore in detail how shell proteins collect the necessary cargo and then assemble into an ordered shell with the cargo on the inside. Experiments alone are probably not enough to understand the process, especially since microcompartment assembly can currently only be studied within live cells or cellular extract. Within these complex environments it is difficult to determine the effect of any individual factor on the overall assembly process.

Perlmutter, Mohajerani and Hagan have now taken a different approach by developing computational and theoretical models to explore how microcompartments assemble. Computer simulations showed that microcompartments could assemble by two pathways. In one pathway, the protein shell and cargo coalesce at the same time. In the other pathway, the cargo molecules first assemble into a large disordered complex, with the shell proteins attached on the outside. The shell proteins then assemble, carving out a piece of the cargo complex. The simulations showed that many factors affect how the shell assembles, such as the strengths of the interactions between the shell proteins and the cargo. They also identified a factor that controls how much cargo ends up inside the assembled shell.

Perlmutter, Mohajerani and Hagan found that, in addition to revealing how microcompartments may assemble within their natural setting, the simulations provided guidance on how to re-engineer microcompartments to assemble around other components. This would enable researchers to create customizable compartments that self-assemble within bacteria or other host organisms, for example to carry out carbon fixation or make biofuels.

A future challenge will be to investigate other aspects of microcompartment assembly, such as the factors that control the size of these compartments.

DOI:http://dx.doi.org/10.7554/eLife.14078.002

Design strategies for self-assembly of discrete targets

Both biological and artificial self-assembly processes can take place by a range of different schemes, from the successive addition of identical building blocks to hierarchical sequences of intermediates, all the way to the fully addressable limit in which each component is unique. In this paper, we introduce an idealized model of cubic particles with patterned faces that allows self-assembly strategies to be compared and tested. We consider a simple octameric target, starting with the minimal requirements for successful self-assembly and comparing the benefits and limitations of more sophisticated hierarchical and addressable schemes. Simulations are performed using a hybrid dynamical Monte Carlo protocol that allows self-assembling clusters to rearrange internally while still providing Stokes-Einstein-like diffusion of aggregates of different sizes. Our simulations explicitly capture the thermodynamic, dynamic, and steric challenges typically faced by self-assembly processes, including competition between multiple partially completed structures. Self-assembly pathways are extracted from the simulation trajectories by a fully extendable scheme for identifying structural fragments, which are then assembled into history diagrams for successfully completed target structures. For the simple target, a one-component assembly scheme is most efficient and robust overall, but hierarchical and addressable strategies can have an advantage under some conditions if high yield is a priority.

The Role of Packaging Sites in Efficient and Specific Virus Assembly

During the life cycle of many single-stranded RNA viruses, including many human pathogens, a protein shell called the capsid spontaneously assembles around the viral genome. Understanding the mechanisms by which capsid proteins selectively assemble around the viral RNA amidst diverse host RNAs is a key question in virology. In one proposed mechanism, short sequences (packaging sites) within the genomic RNA promote rapid and efficient assembly through specific interactions with the capsid proteins. In this work, we develop a coarse-grained particle-based computational model for capsid proteins and RNA that represents protein-RNA interactions arising both from nonspecific electrostatics and from specific packaging site interactions. Using Brownian dynamics simulations, we explore how the efficiency and specificity of assembly depend on solution conditions (which control protein-protein and nonspecific protein-RNA interactions) and the strength and number of packaging sites. We identify distinct regions in parameter space in which packaging sites lead to highly specific assembly via different mechanisms and others in which packaging sites lead to kinetic traps. We relate these computational predictions to in vitro assays for specificity in which cognate viral RNAs compete against non-cognate RNAs for assembly by capsid proteins.

Simulations show that virus assembly and budding are facilitated by membrane microdomains

For many viruses, assembly and budding occur simultaneously during virion formation. Understanding the mechanisms underlying this process could promote biomedical efforts to block viral propagation and enable use of capsids in nanomaterials applications. To this end, we have performed molecular dynamics simulations on a coarse-grained model that describes virus assembly on a fluctuating lipid membrane. Our simulations show that the membrane can promote association of adsorbed subunits through dimensional reduction, but it also introduces thermodynamic and kinetic effects that can inhibit complete assembly. We find several mechanisms by which membrane microdomains, such as lipid rafts, reduce these effects, and thus, enhance assembly. We show how these predicted mechanisms can be experimentally tested. Furthermore, the simulations demonstrate that assembly and budding depend crucially on the system dynamics via multiple timescales related to membrane deformation, protein diffusion, association, and adsorption onto the membrane.

Biomolecular engineering of virus-like particles aided by computational chemistry methods

Multi-scale investigation of VLP self-assembly aided by computational methods is facilitating the design, redesign, and modification of functionalized VLPs.

Pathways for virus assembly around nucleic acids

Understanding the pathways by which viral capsid proteins assemble around their genomes could identify key intermediates as potential drug targets. In this work, we use computer simulations to characterize assembly over a wide range of capsid protein-protein interaction strengths and solution ionic strengths. We find that assembly pathways can be categorized into two classes, in which intermediates are either predominantly ordered or disordered. Our results suggest that estimating the protein-protein and the protein-genome binding affinities may be sufficient to predict which pathway occurs. Furthermore, the calculated phase diagrams suggest that knowledge of the dominant assembly pathway and its relationship to control parameters could identify optimal strategies to thwart or redirect assembly to block infection. Finally, analysis of simulation trajectories suggests that the two classes of assembly pathways can be distinguished in single-molecule fluorescence correlation spectroscopy or bulk time-resolved small-angle X-ray scattering experiments.

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.

Icosahedral capsid formation by capsomers and short polyions

Kinetical and structural aspects of the capsomer-polyion co-assembly into icosahedral viruses have been simulated by molecular dynamics using a coarse-grained model comprising cationic capsomers and short anionic polyions. Conditions were found at which the presence of polyions of a minimum length was necessary for capsomer formation. The largest yield of correctly formed capsids was obtained at which the driving force for capsid formation was relatively weak. Relatively stronger driving forces, i.e., stronger capsomer-capsomer short-range attraction and∕or stronger electrostatic interaction, lead to larger fraction of kinetically trapped structures and aberrant capsids. The intermediate formation was investigated and different evolving scenarios were found by just varying the polyion length.

Viral genome structures are optimal for capsid assembly

Understanding how virus capsids assemble around their nucleic acid (NA) genomes could promote efforts to block viral propagation or to reengineer capsids for gene therapy applications. We develop a coarse-grained model of capsid proteins and NAs with which we investigate assembly dynamics and thermodynamics. In contrast to recent theoretical models, we find that capsids spontaneously ‘overcharge’; that is, the negative charge of the NA exceeds the positive charge on capsid. When applied to specific viruses, the optimal NA lengths closely correspond to the natural genome lengths. Calculations based on linear polyelectrolytes rather than base-paired NAs underpredict the optimal length, demonstrating the importance of NA structure to capsid assembly. These results suggest that electrostatics, excluded volume, and NA tertiary structure are sufficient to predict assembly thermodynamics and that the ability of viruses to selectively encapsidate their genomic NAs can be explained, at least in part, on a thermodynamic basis.

DOI:http://dx.doi.org/10.7554/eLife.00632.001

Viruses are infectious agents made up of proteins and a genome made of DNA or RNA. Upon infecting a host cell, viruses hijack the cell’s gene expression machinery and force it to produce copies of the viral genome and proteins, which then assemble into new viruses that can eventually infect other host cells. Because assembly is an essential step in the viral life cycle, understanding how this process occurs could significantly advance the fight against viral diseases.

In many viral families, a protein shell called a capsid forms around the viral genome during the assembly process. However, capsids can also assemble around nucleic acids in solution, indicating that a host cell is not required for their formation. Since capsid proteins are positively charged, and nucleic acids are negatively charged, electrostatic interactions between the two are thought to have an important role in capsid assembly. However, it is unclear how structural features of the viral genome affect assembly, and why the negative charge on viral genomes is actually far greater than the positive charge on capsids. These questions are difficult to address experimentally because most of the intermediates that form during virus assembly are too short-lived to be imaged.

Here, Perlmutter et al. have used state of the art computational methods and advances in graphical processing units (GPUs) to produce the most realistic model of capsid assembly to date. They showed that the stability of the complex formed between the nucleic acid and the capsid depends on the length of the viral genome. Yield was highest for genomes within a certain range of lengths, and capsids that assembled around longer or shorter genomes tended to be malformed.

Perlmutter et al. also explored how structural features of the virus—including base-pairing between viral nucleic acids, and the size and charge of the capsid—determine the optimal length of the viral genome. When they included structural data from real viruses in their simulations and predicted the optimal lengths for the viral genome, the results were very similar to those seen in existing viruses. This indicates that the structure of the viral genome has been optimized to promote packaging into capsids. Understanding this relationship between structure and packaging will make it easier to develop antiviral agents that thwart or misdirect virus assembly, and could aid the redesign of viruses for use in gene therapy and drug delivery.

DOI:http://dx.doi.org/10.7554/eLife.00632.002

Theoretical studies on assembly, physical stability and dynamics of viruses

All matter has to obey the general laws of physics and living matter is not an exception. Viruses have not only learnt how to cope with them, but have managed to use them for their own survival. In this chapter we will review some of the exciting physics behind viruses and discuss simple physical models that can shed some light on different aspects of the viral life cycle and viral properties. In particular, we will focus on how the structure and shape of the capsid, its assembly and stability, and the entry and exit of viral particles and their genomes can be understood using fundamental physics theories.

Stochastic dynamics of virus capsid formation: direct versus hierarchical self-assembly

BACKGROUND

In order to replicate within their cellular host, many viruses have developed self-assembly strategies for their capsids which are sufficiently robust as to be reconstituted in vitro. Mathematical models for virus self-assembly usually assume that the bonds leading to cluster formation have constant reactivity over the time course of assembly (direct assembly). In some cases, however, binding sites between the capsomers have been reported to be activated during the self-assembly process (hierarchical assembly).

RESULTS

In order to study possible advantages of such hierarchical schemes for icosahedral virus capsid assembly, we use Brownian dynamics simulations of a patchy particle model that allows us to switch binding sites on and off during assembly. For T1 viruses, we implement a hierarchical assembly scheme where inter-capsomer bonds become active only if a complete pentamer has been assembled. We find direct assembly to be favorable for reversible bonds allowing for repeated structural reorganizations, while hierarchical assembly is favorable for strong bonds with small dissociation rate, as this situation is less prone to kinetic trapping. However, at the same time it is more vulnerable to monomer starvation during the final phase. Increasing the number of initial monomers does have only a weak effect on these general features. The differences between the two assembly schemes become more pronounced for more complex virus geometries, as shown here for T3 viruses, which assemble through homogeneous pentamers and heterogeneous hexamers in the hierarchical scheme. In order to complement the simulations for this more complicated case, we introduce a master equation approach that agrees well with the simulation results.

CONCLUSIONS

Our analysis shows for which molecular parameters hierarchical assembly schemes can outperform direct ones and suggests that viruses with high bond stability might prefer hierarchical assembly schemes. These insights increase our physical understanding of an essential biological process, with many interesting potential applications in medicine and materials science.

Coarse-graining of multiprotein assemblies

Multiscale models are important tools to elucidate how small changes in local subunit conformations may propagate to affect the properties of macromolecular complexes. We review recent advances in coarse-graining methods for poly-protein assemblies, systems that are composed of many copies of relatively few components, with a particular focus on viral capsids and cytoskeletal filaments. These methods are grouped into two broad categories-mapping methods, which use information from one scale of representation to parameterize a lower resolution model, and bridging methods, which repeatedly connect different scales during simulation-and we provide examples of both classes at different levels of complexity. Collectively, these models illustrate the numerous approaches to information transfer between scales and demonstrate that the complexity required of the model depends in general on the nature of the information sought.

Formation and assembly-disassembly processes of ZnO hexagonal pyramids driven by dipolar and excluded volume interactions

ZnO hexagonal pyramids were obtained in hydrophilic media without any traditional stabilizers (capping agents). The absence of a thick organic shell reducing the anisotropy of nanoparticle (NP) interactions, oxide nature of the materials, and new geometry of the nanocrystals makes possible the observation of new self-organization phenomena. Several new features not present in the previous cases of NP self-organization were identified and discussed. The formation of ZnO pyramids involved recrystallization of larger amorphous NPs followed by the multistage disassembly of intermediate aggregates into individual virtually perfectly shaped nanocrystals. The evolution of NPs begins with crystallization of clustered plates within the original amorphous spherical colloids, and then agglomerated truncated pyramids are formed. These agglomerates further transform into chained pyramids, which eventually separate from each other. The crystallization and disassembly processes can be associated with the decrease of potential and anisotropy of the attractive force field around the crystallites represented in part by dipole moments. The reassembly of the pyramids can still be attained via engaging excluded volume interaction after adding similarly charged polymer. Overall, in this system, we see the first examples of (1) coupled crystallization and disassembly process; (2) induced assembly of nanoscale particles using excluded volume interactions, which were previously used only for aggregation of microscale colloids; and (3) nanoparticle assemblies with variable and experimentally verifiable relative orientation of dipoles including head-to-tail, tail-to-tail pairs, and antiparallel chains. Described assemblies of ZnO pyramids with collective behavior of individual building blocks as well as distinct and experimentally controlled stages of assembly and disassembly present a fundamentally interesting nanoparticle system with rich dynamic behavior.

Modelling the self-assembly of virus capsids

We use computer simulations to study a model, first proposed by Wales (2005 Phil. Trans. R. Soc. A 363 357), for the reversible and monodisperse self-assembly of simple icosahedral virus capsid structures. The success and efficiency of assembly as a function of thermodynamic and geometric factors can be qualitatively related to the potential energy landscape structure of the assembling system. Even though the model is strongly coarse-grained, it exhibits a number of features also observed in experiments, such as sigmoidal assembly dynamics, hysteresis in capsid formation and numerous kinetic traps. We also investigate the effect of macromolecular crowding on the assembly dynamics. Crowding agents generally reduce capsid yields at optimal conditions for non-crowded assembly, but may increase yields for parameter regimes away from the optimum. Finally, we generalize the model to a larger triangulation number T = 3, and observe assembly dynamics more complex than that seen for the original T = 1 model.

Emergent complexity from simple anisotropic building blocks: shells, tubes, and spirals

We describe a remarkably simple, generic, coarse-grained model involving anisotropic interactions, and characterize the global minima for clusters as a function of various parameters. Appropriate choices for the anisotropic interactions can reproduce a wide variety of complex morphologies as global minima, including spheroidal shells, tubular, helical and even head-tail morphologies, elucidating the physical principles that drive the assembly of these mesoscopic structures. Our model captures several experimental observations, such as the existence of competing morphologies, capsid polymorphism, and the effect of scaffolding proteins on capsid assembly.