Impact of the topology of viral RNAs on their encapsulation by virus coat proteins

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.

Long ssRNA undergoes continuous compaction in the presence of polyvalent cations

In the presence of polyvalent cations, long double-stranded DNA (dsDNA) in dilute solution undergoes a single-molecule, first-order, phase transition ("condensation"), a phenomenon that has been documented and analyzed by many years of experimental and theoretical studies. There has been no systematic effort, however, to determine whether long single-stranded RNA (ssRNA) shows an analogous behavior. In this study, using dynamic light scattering, analytical ultracentrifugation, and gel electrophoresis, we examine the effects of increasing polyvalent cation concentrations on the effective size of long ssRNAs ranging from 3000 to 12,000 nucleotides. Our results indicate that ssRNA does not undergo a discontinuous condensation as does dsDNA but rather a "continuous" decrease in size with increasing polyvalent cation concentration. And, instead of the 10-fold decrease in size shown by long dsDNA, we document a 50% decrease, as demonstrated for a range of lengths and sequences of ssRNA.

Relaxational dynamics of the T-number conversion of virus capsids

We extend a recently proposed kinetic theory of virus capsid assembly based on Model A kinetics and study the dynamics of the interconversion of virus capsids of different sizes triggered by a quench, that is, by sudden changes in the solution conditions. The work is inspired by in vitro experiments on functionalized coat proteins of the plant virus cowpea chlorotic mottle virus, which undergo a reversible transition between two different shell sizes (T = 1 and T = 3) upon changing the acidity and salinity of the solution. We find that the relaxation dynamics are governed by two time scales that, in almost all cases, can be identified as two distinct processes. Initially, the monomers and one of the two types of capsids respond to the quench. Subsequently, the monomer concentration remains essentially constant, and the conversion between the two capsid species completes. In the intermediate stages, a long-lived metastable steady state may present itself, where the thermodynamically less stable species predominate. We conclude that a Model A based relaxational model can reasonably describe the early and intermediate stages of the conversion experiments. However, it fails to provide a good representation of the time evolution of the state of assembly of the coat proteins in the very late stages of equilibration when one of the two species disappears from the solution. It appears that explicitly incorporating the nucleation barriers to assembly and disassembly is crucial for an accurate description of the experimental findings, at least under conditions where these barriers are sufficiently large.

Electrostatic Theory of the Acidity of the Solution in the Lumina of Viruses and Virus-Like Particles

Recently, Maassen et al. measured an appreciable pH difference between the bulk solution and the solution in the lumen of virus-like particles, self-assembled in an aqueous buffer solution containing the coat proteins of a simple plant virus and polyanions (Maassen, S. J.; et al. Small 2018, 14, 1802081). They attribute this to the Donnan effect, caused by an imbalance between the number of negative charges on the encapsulated polyelectrolyte molecules and the number of positive charges on the RNA binding domains of the coat proteins that make up the virus shell or capsid. By applying Poisson-Boltzmann theory, we confirm this conclusion and show that simple Donnan theory is accurate even for the smallest of viruses and virus-like particles. This, in part, is due to the additional screening caused by the presence of a large number of immobile charges in the cavity of the shell. The presence of a net charge on the outer surface of the capsid we find in practice to not have a large effect on the pH shift. Hence, Donnan theory can indeed be applied to connect the local pH and the amount of encapsulated material. The large shifts up to a full pH unit that we predict must have consequences for applications of virus capsids as nanocontainers in bionanotechnology and artificial cell organelles.

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.

Relationships between RNA topology and nucleocapsid structure in a model icosahedral virus

The process of genome packaging in most of viruses is poorly understood, notably the role of the genome itself in the nucleocapsid structure. For simple icosahedral single-stranded RNA viruses, the branched topology due to the RNA secondary structure is thought to lower the free energy required to complete a virion. We investigate the structure of nucleocapsids packaging RNA segments with various degrees of compactness by small-angle x-ray scattering and cryotransmission electron microscopy. The structural differences are mild even though compact RNA segments lead on average to better-ordered and more uniform particles across the sample. Numerical calculations confirm that the free energy is lowered for the RNA segments displaying the larger number of branch points. The effect is, however, opposite with synthetic polyelectrolytes, in which a star topology gives rise to more disorder in the capsids than a linear topology. If RNA compactness and size account in part for the proper assembly of the nucleocapsid and the genome selectivity, other factors most likely related to the host cell environment during viral assembly must come into play as well.

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.

Effect of the charge distribution of virus coat proteins on the length of packaged RNAs

Single-stranded RNA viruses efficiently encapsulate their genome into a protein shell called the capsid. Electrostatic interactions between the positive charges in the capsid protein's N-terminal tail and the negatively charged genome have been postulated as the main driving force for virus assembly. Recent experimental results indicate that the N-terminal tail with the same number of charges and same lengths packages different amounts of RNA, which reveals that electrostatics alone cannot explain all the observed outcomes of the RNA self-assembly experiments. Using a mean-field theory, we show that the combined effect of genome configurational entropy and electrostatics can explain to some extent the amount of packaged RNA with mutant proteins where the location and number of charges on the tails are altered. Understanding the factors contributing to the virus assembly could promote the attempt to block viral infections or to build capsids for gene therapy applications.

Physical virology: From virus self-assembly to particle mechanics

Viruses are highly ordered supramolecular complexes that have evolved to propagate by hijacking the host cell's machinery. Although viruses are very diverse, spreading through cells of all kingdoms of life, they share common functions and properties. Next to the general interest in virology, fundamental viral mechanisms are of growing importance in other disciplines such as biomedicine and (bio)nanotechnology. However, in order to optimally make use of viruses and virus-like particles, for instance as vehicle for targeted drug delivery or as building blocks in electronics, it is essential to understand their basic chemical and physical properties and characteristics. In this context, the number of studies addressing the mechanisms governing viral properties and processes has recently grown drastically. This review summarizes a specific part of these scientific achievements, particularly addressing physical virology approaches aimed to understand the self-assembly of viruses and the mechanical properties of viral particles. Using a physicochemical perspective, we have focused on fundamental studies providing an overview of the molecular basis governing these key aspects of viral systems. This article is categorized under: Biology-Inspired Nanomaterials > Protein and Virus-Based Structures Nanotechnology Approaches to Biology > Nanoscale Systems in Biology.

A New Model System for Exploring Assembly Mechanisms of the HIV-1 Immature Capsid In Vivo

The assembly of the HIV-1 immature capsid (HIC) is an essential step in the virus life cycle. In vivo, the HIC is composed of [Formula: see text] hexameric building blocks, and it takes 5-6 min to complete the assembly process. The involvement of numerous building blocks and the rapid timecourse makes it difficult to understand the HIC assembly process. In this work, we study HIC assembly in vivo by using differential equations. We first obtain a full model with 420 differential equations. Then, we reduce six addition reactions for separate building blocks to a single complex reaction. This strategy reduces the full model to 70 equations. Subsequently, the theoretical analysis of the reduced model shows that it might not be an effective way to decrease the HIC concentration at the equilibrium state by decreasing the microscopic on-rate constants. Based on experimental data, we estimate that the nucleating structure is much smaller than the HIC. We also estimate that the microscopic on-rate constant for nucleation reactions is far less than that for elongation reactions. The parametric collinearity investigation testifies the reliability of these two characteristics, which might explain why free building blocks do not readily polymerize into higher-order polymers until their concentration reaches a threshold value. These results can provide further insight into the assembly mechanisms of the HIC in vivo.

The different faces of mass action in virus assembly

The spontaneous encapsulation of genomic and non-genomic polyanions by coat proteins of simple icosahedral viruses is driven, in the first instance, by electrostatic interactions with polycationic RNA binding domains on these proteins. The efficiency with which the polyanions can be encapsulated in vitro, and presumably also in vivo, must in addition be governed by the loss of translational and mixing entropy associated with co-assembly, at least if this co-assembly constitutes a reversible process. These forms of entropy counteract the impact of attractive interactions between the constituents and hence they counteract complexation. By invoking mass action-type arguments and a simple model describing electrostatic interactions, we show how these forms of entropy might settle the competition between negatively charged polymers of different molecular weights for co-assembly with the coat proteins. In direct competition, mass action turns out to strongly work against the encapsulation of RNAs that are significantly shorter, which is typically the case for non-viral (host) RNAs. We also find that coat proteins favor forming virus particles over nonspecific binding to other proteins in the cytosol even if these are present in vast excess. Our results rationalize a number of recent in vitro co-assembly experiments showing that short polyanions are less effective at attracting virus coat proteins to form virus-like particles than long ones do, even if both are present at equal weight concentrations in the assembly mixture.

The effect of RNA stiffness on the self-assembly of virus particles

Under many in vitro conditions, some small viruses spontaneously encapsidate a single stranded (ss) RNA into a protein shell called the capsid. While viral RNAs are found to be compact and highly branched because of long distance base-pairing between nucleotides, recent experiments reveal that in a head-to-head competition between an ssRNA with no secondary or higher order structure and a viral RNA, the capsid proteins preferentially encapsulate the linear polymer! In this paper, we study the impact of genome stiffness on the encapsidation free energy of the complex of RNA and capsid proteins. We show that an increase in effective chain stiffness because of base-pairing could be the reason why under certain conditions linear chains have an advantage over branched chains when it comes to encapsidation efficiency. While branching makes the genome more compact, RNA base-pairing increases the effective Kuhn length of the RNA molecule, which could result in an increase of the free energy of RNA confinement, that is, the work required to encapsidate RNA, and thus less efficient packaging.

RNA Base Pairing Determines the Conformations of RNA Inside Spherical Viruses

Many simple RNA viruses enclose their genetic material by a protein shell called the capsid. While the capsid structures are well characterized for most viruses, the structure of RNA inside the shells and the factors contributing to it remain poorly understood. We study the impact of base pairing on the conformations of RNA and find that it undergoes a swollen coil to globule continuous transition as a function of the strength of the pairing interaction. We also observe a first order transition and kink profile as a function of RNA length. All these transitions could explain the different RNA profiles observed inside viral shells.

Impact of a nonuniform charge distribution on virus assembly

Many spherical viruses encapsulate their genomes in protein shells with icosahedral symmetry. This process is spontaneous and driven by electrostatic interactions between positive domains on the virus coat proteins and the negative genomes. We model the effect of the nonuniform icosahedral charge distribution from the protein shell instead using a mean-field theory. We find that this nonuniform charge distribution strongly affects the optimal genome length and that it can explain the experimentally observed phenomenon of overcharging of virus and viruslike particles.

Combining Protein Cages and Polymers: from Understanding Self-Assembly to Functional Materials

Protein cages, such as viruses, are well-defined biological nanostructures which are highly symmetrical and monodisperse. They are found in various shapes and sizes and can encapsulate or template non-native materials. Furthermore, the proteins can be chemically or genetically modified giving them new properties. For these reasons, these protein structures have received increasing attention in the field of polymer-protein hybrid materials over the past years, however, advances are still to be made. This Viewpoint highlights the different ways polymers and protein cages or their subunits have been combined to understand self-assembly and create functional materials.

Effects of RNA branching on the electrostatic stabilization of viruses

Many single-stranded (ss) ribonucleic acid (RNA) viruses self-assemble from capsid protein subunits and the nucleic acid to form an infectious virion. It is believed that the electrostatic interactions between the negatively charged RNA and the positively charged viral capsid proteins drive the encapsidation, although there is growing evidence that the sequence of the viral RNA also plays a role in packaging. In particular, the sequence will determine the possible secondary structures that the ssRNA will take in solution. In this work, we use a mean-field theory to investigate how the secondary structure of the RNA combined with electrostatic interactions affects the efficiency of assembly and stability of the assembled virions. We show that the secondary structure of RNA may result in negative osmotic pressures while a linear polymer causes positive osmotic pressures for the same conditions. This may suggest that the branched structure makes the RNA more effectively packaged and the virion more stable.

A Thermodynamic Model for Genome Packaging in Hepatitis B Virus

Understanding the fundamentals of genome packaging in viral capsids is important for finding effective antiviral strategies and for utilizing benign viral particles for gene therapy. While the structure of encapsidated genomic materials has been routinely characterized with experimental techniques such as cryo-electron microscopy and x-ray diffraction, much less is known about the molecular driving forces underlying genome assembly in an intracellular environment and its in vivo interactions with the capsid proteins. Here we study the thermodynamic basis of the pregenomic RNA encapsidation in human Hepatitis B virus in vivo using a coarse-grained molecular model that captures the essential components of nonspecific intermolecular interactions. The thermodynamic model is used to examine how the electrostatic interaction between the packaged RNA and the highly charged C-terminal domains (CTD) of capsid proteins regulate the nucleocapsid formation. The theoretical model predicts optimal RNA content in Hepatitis B virus nucleocapsids with different CTD lengths in good agreement with mutagenesis measurements, confirming the predominant role of electrostatic interactions and molecular excluded-volume effects in genome packaging. We find that the amount of encapsidated RNA is not linearly correlated with the net charge of CTD tails as suggested by earlier theoretical studies. Our thermodynamic analysis of the nucleocapsid structure and stability indicates that ∼10% of the CTD residues are free from complexation with RNA, resulting in partially exposed CTD tails. The thermodynamic model also predicts the free energy of complex formation between macromolecules, which corroborates experimental results for the impact of CTD truncation on the nucleocapsid stability.

Sequence Dependence of Viral RNA Encapsidation

We develop a Flory mean-field theory for viral RNA (vRNA) molecules that extends the current RNA folding algorithms to include interactions between different sections of the secondary structure. The theory is applied to sequence-selective vRNA encapsidation. The dependence on sequence enters through a single parameter: the largest eigenvalue of the Kramers matrix of the branched polymer obtained by coarse graining the secondary structure. Differences between the work of encapsidation of vRNA molecules and of randomized isomers are found to be in the range of 20 kBT, more than sufficient to provide a strong bias in favor of vRNA encapsidation. The method is applied to a packaging competition experiment where large vRNA molecules compete for encapsidation with two smaller RNA species that together have the same nucleotide sequence as the large molecule. We encounter a substantial, generic free energy bias, that also is of the order of 20 kBT, in favor of encapsidating the single large RNA molecule. The bias is mainly the consequence of the fact that dividing up a large vRNA molecule involves the release of stored elastic energy. This provides an important, nonspecific mechanism for preferential encapsidation of single larger vRNA molecules over multiple smaller mRNA molecules with the same total number of nucleotides. The result is also consistent with recent RNA packaging competition experiments by Comas-Garcia et al.1 Finally, the Flory method leads to the result that when two RNA molecules are copackaged, they are expected to remain segregated inside the capsid.

Modeling Effects of RNA on Capsid Assembly Pathways via Coarse-Grained Stochastic Simulation

The environment of a living cell is vastly different from that of an in vitro reaction system, an issue that presents great challenges to the use of in vitro models, or computer simulations based on them, for understanding biochemistry in vivo. Virus capsids make an excellent model system for such questions because they typically have few distinct components, making them amenable to in vitro and modeling studies, yet their assembly can involve complex networks of possible reactions that cannot be resolved in detail by any current experimental technology. We previously fit kinetic simulation parameters to bulk in vitro assembly data to yield a close match between simulated and real data, and then used the simulations to study features of assembly that cannot be monitored experimentally. The present work seeks to project how assembly in these simulations fit to in vitro data would be altered by computationally adding features of the cellular environment to the system, specifically the presence of nucleic acid about which many capsids assemble. The major challenge of such work is computational: simulating fine-scale assembly pathways on the scale and in the parameter domains of real viruses is far too computationally costly to allow for explicit models of nucleic acid interaction. We bypass that limitation by applying analytical models of nucleic acid effects to adjust kinetic rate parameters learned from in vitro data to see how these adjustments, singly or in combination, might affect fine-scale assembly progress. The resulting simulations exhibit surprising behavioral complexity, with distinct effects often acting synergistically to drive efficient assembly and alter pathways relative to the in vitro model. The work demonstrates how computer simulations can help us understand how assembly might differ between the in vitro and in vivo environments and what features of the cellular environment account for these differences.

Role of Genome in the Formation of Conical Retroviral Shells

Human immunodeficiency virus (HIV) capsid proteins spontaneously assemble around the genome into a protective protein shell called the capsid, which can take on a variety of shapes broadly classified as conical, cylindrical, and irregular. The majority of capsids seen in in vivo studies are conical in shape, while in vitro experiments have shown a preference for cylindrical capsids. The factors involved in the selection of the unique shape of HIV capsids are not well understood, and in particular the impact of RNA on the formation of the capsid is not known. In this work, we study the role of the genome and its interaction with the capsid protein by modeling the genomic RNA through a mean-field theory. Our results show that the confinement free energy for a homopolymeric model genome confined in a conical capsid is lower than that in a cylindrical capsid, at least when the genome does not interact with the capsid, which seems to be the case in in vivo experiments. Conversely, the confinement free energy for the cylinder is lower than that for a conical capsid if the genome is attracted to the capsid proteins as the in vitro experiments. Understanding the factors that contribute to the formation of conical capsids may shed light on the infectivity of HIV particles.

Adsorption of annealed branched polymers on curved surfaces

The behavior of annealed branched polymers near adsorbing surfaces plays a fundamental role in many biological and industrial processes. Most importantly single stranded RNA in solution tends to fold up and self-bind to form a highly branched structure. Using a mean field theory, we both perturbatively and numerically examine the adsorption of branched polymers on surfaces of several different geometries in a good solvent. Independent of the geometry of the wall, we observe that as branching density increases, surface tension decreases. However, we find a coupling between the branching density and curvature in that a further lowering of surface tension occurs when the wall curves towards the polymer, but the amount of lowering of surface tension decreases when the wall curves away from the polymer. We find that for branched polymers confined into spherical cavities, most of branch-points are located in the vicinity of the interior wall and the surface tension is minimized for a critical cavity radius. For branch polymers next to sinusoidal surfaces, we find that branch-points accumulate at the valleys while end-points on the peaks.

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.

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.