Controlling viral capsid assembly with templating

Hierarchical assembly is more robust than egalitarian assembly in synthetic capsids

Self-assembly of complex and functional materials remains a grand challenge in soft material science. Efficient assembly depends on a delicate balance between thermodynamic and kinetic effects, requiring fine-tuning affinities and concentrations of subunits. By contrast, we introduce an assembly paradigm that allows large error-tolerance in the subunit affinity and helps avoid kinetic traps. Our combined experimental and computational approach uses a model system of triangular subunits programmed to assemble into T = 3 icosahedral capsids comprising 60 units. The experimental platform uses DNA origami to create monodisperse colloids whose three-dimensional geometry is controlled to nanometer precision, with two distinct bonds whose affinities are controlled to kBT precision, quantified in situ by static light scattering. The computational model uses a coarse-grained representation of subunits, short-ranged potentials, and Langevin dynamics. Experimental observations and modeling reveal that when the bond affinities are unequal, two distinct hierarchical assembly pathways occur, in which the subunits first form dimers in one case and pentamers in another. These hierarchical pathways produce complete capsids faster and are more robust against affinity variation than egalitarian pathways, in which all binding sites have equal strengths. This finding suggests that hierarchical assembly may be a general engineering principle for optimizing self-assembly of complex target structures.

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.

Crystal structure of the potato leafroll virus coat protein and implications for viral assembly

Luteoviruses, poleroviruses, and enamoviruses are insect-transmitted, agricultural pathogens that infect a wide array of plants, including staple food crops. Previous cryo-electron microscopy studies of virus-like particles show that luteovirid viral capsids are built from a structural coat protein that organizes with T = 3 icosahedral symmetry. Here, we present the crystal structure of a truncated version of the coat protein monomer from potato leafroll virus at 1.80-Å resolution. In the crystal lattice, monomers pack into flat sheets that preserve the two-fold and three-fold axes of icosahedral symmetry and show minimal structural deviations when compared to the full-length subunits of the assembled virus-like particle. These observations have important implications in viral assembly and maturation and suggest that the CP N-terminus and its interactions with RNA play an important role in generating capsid curvature.

Equilibrium mechanisms of self-limiting assembly

Self-assembly is a ubiquitous process in synthetic and biological systems, broadly defined as the spontaneous organization of multiple subunits (e.g. macromolecules, particles) into ordered multi-unit structures. The vast majority of equilibrium assembly processes give rise to two states: one consisting of dispersed disassociated subunits, and the other, a bulk-condensed state of unlimited size. This review focuses on the more specialized class of self-limiting assembly, which describes equilibrium assembly processes resulting in finite-size structures. These systems pose a generic and basic question, how do thermodynamic processes involving non-covalent interactions between identical subunits "measure" and select the size of assembled structures? In this review, we begin with an introduction to the basic statistical mechanical framework for assembly thermodynamics, and use this to highlight the key physical ingredients that ensure equilibrium assembly will terminate at finite dimensions. Then, we introduce examples of self-limiting assembly systems, and classify them within this framework based on two broad categories: self-closing assemblies and open-boundary assemblies. These include well-known cases in biology and synthetic soft matter - micellization of amphiphiles and shell/tubule formation of tapered subunits - as well as less widely known classes of assemblies, such as short-range attractive/long-range repulsive systems and geometrically-frustrated assemblies. For each of these self-limiting mechanisms, we describe the physical mechanisms that select equilibrium assembly size, as well as potential limitations of finite-size selection. Finally, we discuss alternative mechanisms for finite-size assemblies, and draw contrasts with the size-control that these can achieve relative to self-limitation in equilibrium, single-species assemblies.

Mechanisms of Scaffold-Mediated Microcompartment Assembly and Size Control

This article describes a theoretical and computational study of the dynamical assembly of a protein shell around a complex consisting of many cargo molecules and long, flexible scaffold molecules. Our study is motivated by bacterial microcompartments, which are proteinaceous organelles that assemble around a condensed droplet of enzymes and reactants. As in many examples of cytoplasmic liquid-liquid phase separation, condensation of the microcompartment interior cargo is driven by flexible scaffold proteins that have weak multivalent interactions with the cargo. Our results predict that the shell size, amount of encapsulated cargo, and assembly pathways depend sensitively on properties of the scaffold, including its length and valency of scaffold-cargo interactions. Moreover, the ability of self-assembling protein shells to change their size to accommodate scaffold molecules of different lengths depends crucially on whether the spontaneous curvature radius of the protein shell is smaller or larger than a characteristic elastic length scale of the shell. Beyond natural microcompartments, these results have important implications for synthetic biology efforts to target alternative molecules for encapsulation by microcompartments or viral shells. More broadly, the results elucidate how cells exploit coupling between self-assembly and liquid-liquid phase separation to organize their interiors.

Frustration and Direct-Coupling Analyses to Predict Formation and Function of Adeno-Associated Virus

Adeno-associated virus (AAV) is a promising gene therapy vector because of its efficient gene delivery and relatively mild immunogenicity. To improve delivery target specificity, researchers use combinatorial and rational library design strategies to generate novel AAV capsid variants. These approaches frequently propose high proportions of nonforming or noninfective capsid protein sequences that reduce the effective depth of synthesized vector DNA libraries, thereby raising the discovery cost of novel vectors. We evaluated two computational techniques for their ability to estimate the impact of residue mutations on AAV capsid protein-protein interactions and thus predict changes in vector fitness, reasoning that these approaches might inform the design of functionally enriched AAV libraries and accelerate therapeutic candidate identification. The Frustratometer computes an energy function derived from the energy landscape theory of protein folding. Direct-coupling analysis (DCA) is a statistical framework that captures residue coevolution within proteins. We applied the Frustratometer to select candidate protein residues predicted to favor assembled or disassembled capsid states, then predicted mutation effects at these sites using the Frustratometer and DCA. Capsid mutants were experimentally assessed for changes in virus formation, stability, and transduction ability. The Frustratometer-based metric showed a counterintuitive correlation with viral stability, whereas a DCA-derived metric was highly correlated with virus transduction ability in the small population of residues studied. Our results suggest that coevolutionary models may be able to elucidate complex capsid residue-residue interaction networks essential for viral function, but further study is needed to understand the relationship between protein energy simulations and viral capsid metastability.

Spike Proteins of SARS-CoV and SARS-CoV-2 Utilize Different Mechanisms to Bind With Human ACE2

The ongoing outbreak of COVID-19 has been a serious threat to human health worldwide. The virus SARS-CoV-2 initiates its infection to the human body via the interaction of its spike (S) protein with the human Angiotensin-Converting Enzyme 2 (ACE2) of the host cells. Therefore, understanding the fundamental mechanisms of how SARS-CoV-2 S protein receptor binding domain (RBD) binds to ACE2 is highly demanded for developing treatments for COVID-19. Here we implemented multi-scale computational approaches to study the binding mechanisms of human ACE2 and S proteins of both SARS-CoV and SARS-CoV-2. Electrostatic features, including electrostatic potential, electric field lines, and electrostatic forces of SARS-CoV and SARS-CoV-2 were calculated and compared in detail. The results demonstrate that SARS-CoV and SARS-CoV-2 S proteins are both attractive to ACE2 by electrostatic forces even at different distances. However, the residues contributing to the electrostatic features are quite different due to the mutations between SARS-CoV S protein and SARS-CoV-2 S protein. Such differences are analyzed comprehensively. Compared to SARS-CoV, the SARS-CoV-2 binds with ACE2 using a more robust strategy: The electric field line related residues are distributed quite differently, which results in a more robust binding strategy of SARS-CoV-2. Also, SARS-CoV-2 has a higher electric field line density than that of SARS-CoV, which indicates stronger interaction between SARS-CoV-2 and ACE2, compared to that of SARS-CoV. Key residues involved in salt bridges and hydrogen bonds are identified in this study, which may help the future drug design against COVID-19.

Using computational approaches to study dengue virus capsid assembly

Dengue viral capsid plays a significant role in viral life cycle of dengue, especially in vial genome protection and virus-cell fusion. Revealing mechanisms of the viral capsid protein assembly may lead to the discovery of anti-viral drugs that inhibit the assembly of the viral capsid. The E and M-proteins are arranged into heterotetramers, which consists of two copies of E and M-protein. The heterotetramers are assembled into a highly ordered capsid. While many investigations of the interactions between E and M-proteins have been performed, there are very few studies on the interactions between the heterotetramers and their roles in capsid assembly. Utilizing a series of computational approaches, this study focuses on the assembly mechanism of the heterotetramers. Our electrostatic analyses lead to the identification of four binding modes between each two dengue heterotetramers that repeat periodically throughout the virus capsid. Among these four binding modes, heterotetramers in binding modes I, II and IV are attractive. But in the binding mode III the heterotetramers repel each other, making mode III a suitable target for drug design. Furthermore, MD simulations were performed following by salt bridges analysis. This study demonstrates that using computational approaches is a promising direction to study the dengue virus.

Allosteric pathway selection in templated assembly

Assembling large numbers of molecular building blocks into functional nanostructures is no trivial task. It relies on guiding building blocks through complex energy landscapes shaped by synergistic and antagonistic supramolecular interactions. In nature, the use of molecular templates is a potent strategy to navigate the process to the desired structure with high fidelity. Yet, nature's templating strategy remains to be fully exploited in man-made nanomaterials. Designing effective template-guided self-assembling systems can only be realized through precise insight into how the chemical design of building blocks and the resulting balance of repulsive and attractive forces give rise to pathway selection and suppression of trapped states. We develop a minimal model to unravel the kinetic pathways and pathway selection of the templated assembly of molecular building blocks on a template. We show how allosteric activation of the associative interactions can suppress undesired solution-aggregation pathways and gives rise to a true template-assembly path.

Packaging of Genomic RNA in Positive-Sense Single-Stranded RNA Viruses: A Complex Story

The packaging of genomic RNA in positive-sense single-stranded RNA viruses is a key part of the viral infectious cycle, yet this step is not fully understood. Unlike double-stranded DNA and RNA viruses, this process is coupled with nucleocapsid assembly. The specificity of RNA packaging depends on multiple factors: (i) one or more packaging signals, (ii) RNA replication, (iii) translation, (iv) viral factories, and (v) the physical properties of the RNA. The relative contribution of each of these factors to packaging specificity is different for every virus. In vitro and in vivo data show that there are different packaging mechanisms that control selective packaging of the genomic RNA during nucleocapsid assembly. The goals of this article are to explain some of the key experiments that support the contribution of these factors to packaging selectivity and to draw a general scenario that could help us move towards a better understanding of this step of the viral infectious cycle.

Assembly and Characterization of HBc Derived Virus-like Particles with Magnetic Core

Core-virus like particles (VLPs) assembly is a kinetically complex cascade of interactions between viral proteins, nanoparticle's surface and an ionic environment. Despite many in silico simulations regarding this process, there is still a lack of experimental data. The main goal of this study was to investigate the capsid protein of hepatitis B virus (HBc) assembly into virus-like particles with superparamagnetic iron oxide nanoparticles (SPIONs) as a magnetic core in relation to their characteristics. The native form of HBc was obtained via agroinfection of Nicotiana benthamiana with pEAQ-HBc plasmid. SPIONs of diameter of 15 nm were synthesized and functionalized with two ligands, providing variety in ζ-potential and hydrodynamic diameter. The antigenic potential of the assembled core-VLPs was assessed with enzyme-linked immunosorbent assay (ELISA). Morphology of SPIONs and core-VLPs was evaluated via transmission electron microscopy (TEM). The most successful core-VLPs assembly was obtained for SPIONs functionalized with dihexadecyl phosphate (DHP) at SPIONs/HBc ratio of 0.2/0.05 mg/mL. ELISA results indicate significant decrease of antigenicity concomitant with core-VLPs assembly. In summary, this study provides an experimental assessment of the crucial parameters guiding SPION-HBc VLPs assembly and evaluates the antigenicity of the obtained structures.

Self-assembly of convex particles on spherocylindrical surfaces

The precise control of assembly and packing of proteins and colloids on curved surfaces has fundamental implications in nanotechnology. In this paper, we describe dynamical simulations of the self-assembly of conical subunits around a spherocylindrical template, and a continuum theory for the bending energy of a triangular lattice with spontaneous curvature on a surface with arbitrary curvature. We find that assembly depends sensitively on mismatches between subunit spontaneous curvature and the mean curvature of the template, as well as anisotropic curvature of the template (mismatch between the two principal curvatures). Our simulations predict assembly morphologies that closely resemble those observed in experiments in which virus capsid proteins self-assemble around metal nanorods. Below a threshold curvature mismatch, our simulations identify a regime of optimal assembly leading to complete, symmetrical particles. Outside of this regime we observe defective particles, whose morphologies depend on the degree of curvature mismatch. To learn how assembly is affected by the nonuniform curvature of a spherocylinder, we also study the simpler cases of assembly around spherical and cylindrical cores. Our results show that both the intrinsic (Gaussian) and extrinsic (mean) curvatures of a template play significant roles in guiding the assembly of anisotropic subunits, providing a rich design space for the formation of nanoscale materials.

Defects and Chirality in the Nanoparticle-Directed Assembly of Spherocylindrical Shells of Virus Coat Proteins

Virus coat proteins of small isometric plant viruses readily assemble into symmetric, icosahedral cages encapsulating noncognate cargo, provided the cargo meets a minimal set of chemical and physical requirements. While this capability has been intensely explored for certain virus-enabled nanotechnologies, additional applications require lower symmetry than that of an icosahedron. Here, we show that the coat proteins of an icosahedral virus can efficiently assemble around metal nanorods into spherocylindrical closed shells with hexagonally close-packed bodies and icosahedral caps. Comparison of chiral angles and packing defects observed by in situ atomic force microscopy with those obtained from molecular dynamics models offers insight into the mechanism of growth, and the influence of stresses associated with intrinsic curvature and assembly pathways.

Oligonucleotide Length-Dependent Formation of Virus-Like Particles

Understanding the assembly pathway of viruses can contribute to creating monodisperse virus-based materials. In this study, the cowpea chlorotic mottle virus (CCMV) is used to determine the interactions between the capsid proteins of viruses and their cargo. The assembly of the capsid proteins in the presence of different lengths of short, single-stranded (ss) DNA is studied at neutral pH, at which the protein-protein interactions are weak. Chromatography, electrophoresis, microscopy, and light scattering data show that the assembly efficiency and speed of the particles increase with increasing length of oligonucleotides. The minimal length required for assembly under the conditions used herein is 14 nucleotides. Assembly of particles containing such short strands of ssDNA can take almost a month. This slow assembly process enabled the study of intermediate states, which confirmed a low cooperative assembly for CCMV and allowed for further expansion of current assembly theories.

Computational support for a scaffolding mechanism of centriole assembly

Centrioles are essential for forming cilia, flagella and centrosomes. Successful centriole assembly requires proteins of the SAS-6 family, which can form oligomeric ring structures with ninefold symmetry in vitro. While important progress has been made in understanding SAS-6 protein biophysics, the mechanisms enabling ring formation in vivo remain elusive. Likewise, the mechanisms by which a nascent centriole forms near-orthogonal to an existing one are not known. Here, we investigate possible mechanisms of centriole assembly using coarse-grained Brownian dynamics computer simulations in combination with a rate equation approach. Our results suggest that without any external factors, strong stabilization associated with ring closure would be needed to enable efficient ring formation. Strikingly, our simulations reveal that a scaffold-assisted assembly mechanism can trigger robust ring formation owing to local cooperativity, and that this mechanism can also impart orthogonalilty to centriole assembly. Overall, our findings provide novel insights into the organizing principles governing the assembly of this important organelle.

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

Allosteric Control of Icosahedral Capsid Assembly

During the life cycle of a virus, viral proteins and other components self-assemble to form an ordered protein shell called a capsid. This assembly process is subject to multiple competing constraints, including the need to form a thermostable shell while avoiding kinetic traps. It has been proposed that viral assembly satisfies these constraints through allosteric regulation, including the interconversion of capsid proteins among conformations with different propensities for assembly. In this article, we use computational and theoretical modeling to explore how such allostery affects the assembly of icosahedral shells. We simulate assembly under a wide range of protein concentrations, protein binding affinities, and two different mechanisms of allosteric control. We find that above a threshold strength of allosteric control, assembly becomes robust over a broad range of subunit binding affinities and concentrations, allowing the formation of highly thermostable capsids. Our results suggest that allostery can significantly shift the range of protein binding affinities that lead to successful assembly and thus should be taken into account in models that are used to estimate interaction parameters from experimental data.

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.

Topological effects on capsomer-polyion co-assembly

On the basis of a T = 1 icosahedral capsid model, the capsomer-polyion co-assembly process has been investigated by molecular dynamics simulations using capsomers with different net charge and charge distribution as well as linear, branched, and hyper-branched polyions. The assembly process was characterized in terms of the time-dependent cluster size probabilities, averaged cluster size, encapsulation efficiency, and polyion extension. The kinetics of the capsid formation displayed a two-step process. The first one comprised adsorption of capsomers on the polyion, driven by their electrostatic attraction, whereas the second one involved a relocation and/or reorientation of adsorbed capsomers, which rate is reduced upon increasing electrostatic interaction. We found that increased polyion branching facilitated a more rapid encapsulation process towards a higher yield. Moreover, the hyper-branched polyions were entirely encapsulated at all polyion-capsid charge ratios considered.

A molecular thermodynamic model for the stability of hepatitis B capsids

Self-assembly of capsid proteins and genome encapsidation are two critical steps in the life cycle of most plant and animal viruses. A theoretical description of such processes from a physiochemical perspective may help better understand viral replication and morphogenesis thus provide fresh insights into the experimental studies of antiviral strategies. In this work, we propose a molecular thermodynamic model for predicting the stability of Hepatitis B virus (HBV) capsids either with or without loading nucleic materials. With the key components represented by coarse-grained thermodynamic models, the theoretical predictions are in excellent agreement with experimental data for the formation free energies of empty T4 capsids over a broad range of temperature and ion concentrations. The theoretical model predicts T3/T4 dimorphism also in good agreement with the capsid formation at in vivo and in vitro conditions. In addition, we have studied the stability of the viral particles in response to physiological cellular conditions with the explicit consideration of the hydrophobic association of capsid subunits, electrostatic interactions, molecular excluded volume effects, entropy of mixing, and conformational changes of the biomolecular species. The course-grained model captures the essential features of the HBV nucleocapsid stability revealed by recent experiments.

Mechanisms of virus assembly

Viruses are nanoscale entities containing a nucleic acid genome encased in a protein shell called a capsid and in some cases are surrounded by a lipid bilayer membrane. This review summarizes the physics that govern the processes by which capsids assemble within their host cells and in vitro. We describe the thermodynamics and kinetics for the assembly of protein subunits into icosahedral capsid shells and how these are modified in cases in which the capsid assembles around a nucleic acid or on a lipid bilayer. We present experimental and theoretical techniques used to characterize capsid assembly, and we highlight aspects of virus assembly that are likely to receive significant attention in the near future.

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.

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.

Using Markov state models to study self-assembly

Markov state models (MSMs) have been demonstrated to be a powerful method for computationally studying intramolecular processes such as protein folding and macromolecular conformational changes. In this article, we present a new approach to construct MSMs that is applicable to modeling a broad class of multi-molecular assembly reactions. Distinct structures formed during assembly are distinguished by their undirected graphs, which are defined by strong subunit interactions. Spatial inhomogeneities of free subunits are accounted for using a recently developed Gaussian-based signature. Simplifications to this state identification are also investigated. The feasibility of this approach is demonstrated on two different coarse-grained models for virus self-assembly. We find good agreement between the dynamics predicted by the MSMs and long, unbiased simulations, and that the MSMs can reduce overall simulation time by orders of magnitude.

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.

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.

To build a virus on a nucleic acid substrate

Many viruses package their genomes concomitant with assembly. Here, we show that this reaction can be described by three coefficients: association of capsid protein (CP) to nucleic acid (NA), KNA; CP-CP interaction, ω; and α, proportional to the work required to package NA. The value of α can vary as NA is packaged. A phase diagram of average lnα versus lnω identifies conditions where assembly is likely to fail or succeed. NA morphology can favor (lnα > 0) or impede (lnα < 0) assembly. As lnω becomes larger, capsids become more stable and assembly becomes more cooperative. Where (lnα + lnω) < 0, the CP is unable to contain the NA, so that assembly results in aberrant particles. This phase diagram is consistent with quantitative studies of cowpea chlorotic mottle virus, hepatitis B virus, and simian virus 40 assembling on ssRNA and dsDNA substrates. Thus, the formalism we develop is suitable for describing and predicting behavior of experimental studies of CP assembly on NA.

Budding pathway in the templated assembly of viruslike particles

A new pathway for the assembly of viral capsid protein around inorganic nanoparticle cores was observed by time-course light scattering and cryo-electron tomography. Gold nanoparticles with an average diameter of 11.3 nm have been used as a template for the assembly of Brome mosaic virus (BMV) capsid protein at different concentrations. At least at low protein concentrations the kinetic features of the scattering and extinction measurements are consistent with the initial rapid formation of large nanoparticle-protein clusters, which subsequently separate into individual viruslike particles (VLPs). The occurrence of multiparticle clusters at short times after mixing nanoparticles and proteins was confirmed by cryo-EM. Cryo-electron tomography of the multiparticle clusters yielded an average surface-to-surface interparticle distance of ∼7.5 nm, equivalent to ∼1.5 times the thickness of a protein shell. We propose a scenario in which VLP generation may take place through monomer exchange between aggregated particles with defect-ridden or incomplete shells, leading to the formation of stable icosahedral shells, which eventually bud off the aggregate. Together with results from previous works, the findings highlight the astonishing versatility of plant virus capsid protein assembly. This previously unknown mechanism for VLP formation has features that may have relevance for the crowded environment characterizing virus factories in the cell.

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

Encapsulation of gold nanoparticles into self-assembling protein nanoparticles

BACKGROUND

Gold nanoparticles are useful tools for biological applications due to their attractive physical and chemical properties. Their applications can be further expanded when they are functionalized with biological molecules. The biological molecules not only provide the interfaces for interactions between nanoparticles and biological environment, but also contribute their biological functions to the nanoparticles. Therefore, we used self-assembling protein nanoparticles (SAPNs) to encapsulate gold nanoparticles. The protein nanoparticles are formed upon self-assembly of a protein chain that is composed of a pentameric coiled-coil domain at the N-terminus and trimeric coiled-coil domain at the C-terminus. The self-assembling protein nanoparticles form a central cavity of about 10 nm in size, which is ideal for the encapsulation of gold nanoparticles with similar sizes.

RESULTS

We have used SAPNs to encapsulate several commercially available gold nanoparticles. The hydrodynamic size and the surface coating of gold nanoparticles are two important factors influencing successful encapsulation by the SAPNs. Gold nanoparticles with a hydrodynamic size of less than 15 nm can successfully be encapsulated. Gold nanoparticles with citrate coating appear to have stronger interactions with the proteins, which can interfere with the formation of regular protein nanoparticles. Upon encapsulation gold nanoparticles with polymer coating interfere less strongly with the ability of the SAPNs to assemble into nanoparticles. Although the central cavity of the SAPNs carries an overall charge, the electrostatic interaction appears to be less critical for the efficient encapsulation of gold nanoparticles into the protein nanoparticles.

CONCLUSIONS

The SAPNs can be used to encapsulate gold nanoparticles. The SAPNs can be further functionalized by engineering functional peptides or proteins to either their N- or C-termini. Therefore encapsulation of gold nanoparticles into SAPNs can provide a useful platform to generate a multifunctional biodevices.

Unlocking internal prestress from protein nanoshells

The capsids of icosahedral viruses are closed shells assembled from a hexagonal lattice of proteins with fivefold angular defects located at the icosahedral vertices. Elasticity theory predicts that these disclinations are subject to an internal compressive prestress, which provides an explanation for the link between size and shape of capsids. Using a combination of experiment and elasticity theory we investigate the question of whether macromolecular assemblies are subject to residual prestress, due to basic geometric incompatibility of the subunits. Here we report the first direct experimental test of the theory: by controlled removal of protein pentamers from the icosahedral vertices, we measure the mechanical response of so-called "whiffle ball" capsids of herpes simplex virus, and demonstrate the signature of internal prestress locked into wild-type capsids during assembly.

Early stages of the HIV-1 capsid protein lattice formation

The early stages in the formation of the HIV-1 capsid (CA) protein lattice are investigated. The underlying coarse-grained (CG) model is parameterized directly from experimental data and examined under various native contact interaction strengths, CA dimer interfacial configurations, and local surface curvatures. The mechanism of early contiguous mature-style CA p6 lattice formation is explored, and a trimer-of-dimers structure is found to be crucial for CA lattice production. Quasi-equivalent generation of both the pentamer and hexamer components of the HIV-1 viral CA is also demonstrated, and the formation of pentamers is shown to be highly sensitive to local curvature, supporting the view that such inclusions in high-curvature regions allow closure of the viral CA surface. The complicated behavior of CA lattice self-assembly is shown to be reducible to a relatively simple function of the trimer-of-dimers behavior.

Langevin dynamics simulation of polymer-assisted virus-like assembly

Starting from a coarse grained representation of the building units of the minute virus of mice and a flexible polyelectrolyte molecule, we have explored the mechanism of assembly into icosahedral structures with the help of Langevin dynamics simulations and the parallel tempering technique. Regular icosahedra with appropriate symmetry form only in a narrow range of temperature and polymer length. Within this region of parameters where successful assembly would proceed, we have systematically investigated the growth kinetics. The assembly of icosahedra is found to follow the classical nucleation and growth mechanism in the absence of the polymer, with the three regimes of nucleation, linear growth, and slowing down in the later stage. The calculated average nucleation time obeys the laws expected from the classical nucleation theory. The linear growth rate is found to obey the laws of secondary nucleation as in the case of lamellar growth in polymer crystallization. The same mechanism is seen in the simulations of the assembly of icosahedra in the presence of the polymer as well. The polymer reduces the nucleation barrier significantly by enhancing the local concentration of subunits via adsorbing them on their backbone. The details of growth in the presence of the polymer are also found to be consistent with the classical nucleation theory, despite the smallness of the assembled structures.

Physics of shell assembly: line tension, hole implosion, and closure catastrophe

The self-assembly of perfectly ordered closed shells is a challenging process involved in many biological and nanoscale systems. However, most of the aspects that determine their formation are still unknown. Here we investigate the growth of shells by simulating the assembly of spherical structures made of N identical subunits. Remarkably, we show that the formation and energetics of partially assembled shells are dominated by an effective line-tension that can be described in simple thermodynamic terms. In addition, we unveil two mechanisms that can prevent the correct formation of defect-free structures: "hole implosion," which leads to a premature closure of the shell; and "closure catastrophe," which causes a dramatic production of structural disorder during the later stages of the growth of big shells.

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.

Exploring the paths of (virus) assembly

Assembly of viruses that have hundreds of subunits or folding of proteins that have hundreds of amino acids-complex biological reactions-are often spontaneous and rapid. Here, we examine the complete set of intermediates available for the assembly of a hypothetical viruslike particle and the connectivity between these intermediates in a graph-theory-inspired study. Using a build-up procedure, assuming ideal geometry, we enumerated the complete set of 2,423,313 species for formation of an icosahedron from 30 dimeric subunits. Stability of each n-subunit intermediate was defined by the number of contacts between subunits. The probability of forming an intermediate was based on the number of paths to it from its precedecessors. When defining population subsets predicted to have the greatest impact on assembly, both stability- and probability-based criteria select a small group of compact and degenerate species; ergo, only a few hundred intermediates make a measurable contribution to assembly. Though the number of possible intermediates grows combinatorially with the number of subunits in the capsid, the number of intermediates that make a significant contribution to the reaction grows by a much smaller function, a result that may contribute to our understanding of assembly and folding reactions.

Encapsulation of a polymer by an icosahedral virus

The coat proteins of many viruses spontaneously form icosahedral capsids around nucleic acids or other polymers. Elucidating the role of the packaged polymer in capsid formation could promote biomedical efforts to block viral replication and enable use of capsids in nanomaterials applications. To this end, we perform Brownian dynamics on a coarse-grained model that describes the dynamics of icosahedral capsid assembly around a flexible polymer. We identify several mechanisms by which the polymer plays an active role in its encapsulation, including cooperative polymer-protein motions. These mechanisms are related to experimentally controllable parameters such as polymer length, protein concentration and solution conditions. Furthermore, the simulations demonstrate that assembly mechanisms are correlated with encapsulation efficiency, and we present a phase diagram that predicts assembly outcomes as a function of experimental parameters. We anticipate that our simulation results will provide a framework for designing in vitro assembly experiments on single-stranded RNA virus capsids.

Mechanisms of capsid assembly around a polymer

Capsids of many viruses assemble around nucleic acids or other polymers. Understanding how the properties of the packaged polymer affect the assembly process could promote biomedical efforts to prevent viral assembly or nanomaterials applications that exploit assembly. To this end, we simulate on a lattice the dynamical assembly of closed, hollow shells composed of several hundred to 1000 subunits, around a flexible polymer. We find that assembly is most efficient at an optimum polymer length that scales with the surface area of the capsid; polymers that are significantly longer than optimal often lead to partial-capsids with unpackaged polymer "tails" or a competition between multiple partial-capsids attached to a single polymer. These predictions can be tested with bulk experiments in which capsid proteins assemble around homopolymeric RNA or synthetic polyelectrolytes. We also find that the polymer can increase the net rate of subunit accretion to a growing capsid both by stabilizing the addition of new subunits and by enhancing the incoming flux of subunits; the effects of these processes may be distinguishable with experiments that monitor the assembly of individual capsids.

Viral coat proteins as flexible nano-building-blocks for nanoparticle encapsulation

Viral capsid-nanoparticle hybrid structures offer new opportunities for nanobiotechnology. We previously generated virus-based nanoparticles (VNPs) of simian virus 40 (SV40) containing quantum dots (QDs) for cellular imaging. However, as an interesting issue of nano-bio interfaces, the mechanism of nanoparticle (NP) encapsulation by viral coat proteins remains unclear. Here, four kinds of QDs with the same core/shell but different surface coatings are tested for encapsulation. All the QDs can be encapsulated efficiently and there is no correlation between the encapsulation efficiency and the surface charge of the QDs. All the SV40 VNPs encapsulating differently modified QDs show similar structures, fluorescence properties, and activity in entering living cells. These results demonstrate the flexibility of SV40 major capsid protein VP1 in NP encapsulation and provide new clues to the mechanism of NP packaging by viral shells.

Role of surface charge density in nanoparticle-templated assembly of bromovirus protein cages

Self-assembling icosahedral protein cages have potentially useful physical and chemical characteristics for a variety of nanotechnology applications, ranging from therapeutic or diagnostic vectors to building blocks for hierarchical materials. For application-specific functional control of protein cage assemblies, a deeper understanding of the interaction between the protein cage and its payload is necessary. Protein-cage encapsulated nanoparticles, with their well-defined surface chemistry, allow for systematic control over key parameters of encapsulation such as the surface charge, hydrophobicity, and size. Independent control over these variables allows experimental testing of different assembly mechanism models. Previous studies done with Brome mosaic virus capsids and negatively charged gold nanoparticles indicated that the result of the self-assembly process depends on the diameter of the particle. However, in these experiments, the surface-ligand density was maintained at saturation levels, while the total charge and the radius of curvature remained coupled variables, making the interpretation of the observed dependence on the core size difficult. The current work furnishes evidence of a critical surface charge density for assembly through an analysis aimed at decoupling the surface charge and the core size.

Understanding the concentration dependence of viral capsid assembly kinetics--the origin of the lag time and identifying the critical nucleus size

The kinetics for the assembly of viral proteins into a population of capsids can be measured in vitro with size exclusion chromatography or dynamic light scattering, but extracting mechanistic information from these studies is challenging. For example, it is not straightforward to determine the critical nucleus size or the elongation time (the time required for a nucleus to grow to completion). In this work, we study theoretical and computational models for capsid assembly to show that the critical nucleus size can be determined from the concentration dependence of the assembly half-life and that the elongation time is revealed by the length of the lag phase. Furthermore, we find that the system becomes kinetically trapped when nucleation becomes fast compared to elongation. Implications of this constraint for determining elongation mechanisms from experimental assembly data are discussed.

Thermodynamics of nanospheres encapsulated in virus capsids

We investigate the thermodynamics of complexation of functionalized charged nanospheres with viral proteins. The physics of this problem is governed not only by electrostatic interaction between the proteins and the nanosphere cores (screened by salt ions), but also by configurational degrees of freedom of the charged protein N tails. We approach the problem by constructing an appropriate complexation free-energy functional. On the basis of both numerical and analytical studies of this functional we construct the phase diagram for the assembly which contains the information on the assembled structures that appear in the thermodynamical equilibrium, depending on the size and surface charge density of the nanosphere cores. We show that both the nanosphere core charge and its radius determine the size of the capsid that forms around the core.

Optimal architectures of elongated viruses

Many viruses protect their genetic material by a closed elongated protein shell. Unlike spherical viruses, the structure of these prolates is not yet well understood, and only a few of them have been fully characterized. We present the results of a simple phenomenological model, which describes the remarkable structures of prolate or bacilliform viral shells. Surprisingly, we find that the special well-defined geometry of these elongated viruses arises just as a consequence of free-energy minimization of a generic interaction between the structural units of the capsid. Hemispherical T-number caps centered along the 5-, 3-, and 2-fold axes with hexagonally ordered cylindrical bodies are found to be local energy minima, thus justifying their occurrence as optimal viral structures. Moreover, closed elongated viruses show a sequence of magic numbers for the end-caps, leading to strict selection rules for the length and structure of the body as well as for the number of capsomers and proteins of the capsid. The model reproduces the architecture of spherical and bacilliform viruses, both in vivo and in vitro, and constitutes an important step towards understanding viral assembly and its potential control for biological and nanotechnological applications.

Self-assembly of monodisperse clusters: Dependence on target geometry

We apply a simple model system of patchy particles to study monodisperse self-assembly using the Platonic solids as target structures. We find marked differences between the assembly behaviors of the different systems. Tetrahedra, octahedral, and icosahedra assemble easily, while cubes are more challenging and dodecahedra do not assemble. We relate these differences to the kinetics and thermodynamics of assembly, with the formation of large disordered aggregates a particular important competitor to correct assembly. In particular, the free energy landscapes of those targets that are easy to assemble are funnel-like, whereas for the dodecahedral system the landscape is relatively flat with little driving force to facilitate escape from disordered aggregates.