Nonspecific interactions in spontaneous assembly of empty versus functional single-stranded RNA viruses

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.

Mechanisms for destabilisation of RNA viruses at air-water and liquid-liquid interfaces

Understanding the interactions between viruses and surfaces or interfaces is important, as they provide the principles underpinning the cleaning and disinfection of contaminated surfaces. Yet, the physics of such interactions is currently poorly understood. For instance, there are longstanding experimental observations suggesting that the presence of air-water interfaces can generically inactivate and kill viruses, yet the mechanism underlying this phenomenon remains unknown. Here we use theory and simulations to show that electrostatics may provide one such mechanism, and that this is very general. Thus, we predict that the electrostatic free energy of an RNA virus should increase by several thousands of k_BT as the virion breaches an air-water interface. We also show that the fate of a virus approaching a generic liquid-liquid interface depends strongly on the detailed balance between interfacial and electrostatic forces, which can be tuned, for instance, by choosing different media to contact a virus-laden respiratory droplet. Tunability arises because both the electrostatic and interfacial forces scale similarly with viral size. We propose that these results can be used to design effective strategies for surface disinfection.

We know that air-water interfaces can generically inactivate viruses, but the mechanisms behind this observation are unclear. Here the authors use simulations to uncover those mechanisms and find that the electrostatic repulsive free energy of an RNA virus increases by several thousands of kBT as it approaches an air-water interface, providing a mechanism for viral destabilization which may induce inactivation.

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.

The impact of local assembly rules on RNA packaging in a T = 1 satellite plant virus

The vast majority of viruses consist of a nucleic acid surrounded by a protective icosahedral protein shell called the capsid. During viral infection of a host cell, the timing and efficiency of the assembly process is important for ensuring the production of infectious new progeny virus particles. In the class of single-stranded RNA (ssRNA) viruses, the assembly of the capsid takes place in tandem with packaging of the ssRNA genome in a highly cooperative co-assembly process. In simple ssRNA viruses such as the bacteriophage MS2 and small RNA plant viruses such as STNV, this cooperative process results from multiple interactions between the protein shell and sites in the RNA genome which have been termed packaging signals. Using a stochastic assembly algorithm which includes cooperative interactions between the protein shell and packaging signals in the RNA genome, we demonstrate that highly efficient assembly of STNV capsids arises from a set of simple local rules. Altering the local assembly rules results in different nucleation scenarios with varying assembly efficiencies, which in some cases depend strongly on interactions with RNA packaging signals. Our results provide a potential simple explanation based on local assembly rules for the ability of some ssRNA viruses to spontaneously assemble around charged polymers and other non-viral RNAs in vitro.

Virus Isoelectric Point Estimation: Theories and Methods

Much of virus fate, both in the environment and in physical/chemical treatment, is dependent on electrostatic interactions. Developing an accurate means of predicting virion isoelectric point (pI) would help to understand and anticipate virus fate and transport, especially for viruses that are not readily propagated in the lab. One simple approach to predicting pI estimates the pH at which the sum of charges from ionizable amino acids in capsid proteins approaches zero. However, predicted pIs based on capsid charges frequently deviate by several pH units from empirically measured pIs. Recently, the discrepancy between empirical and predicted pI was attributed to the electrostatic neutralization of predictable polynucleotide-binding regions (PBRs) of the capsid interior. In this paper, we review models presupposing (i) the influence of the viral polynucleotide on surface charge or (ii) the contribution of only exterior residues to surface charge. We then compare these models to the approach of excluding only PBRs and hypothesize a conceptual electrostatic model that aligns with this approach. The PBR exclusion method outperformed methods based on three-dimensional (3D) structure and accounted for major discrepancies in predicted pIs without adversely affecting pI prediction for a diverse range of viruses. In addition, the PBR exclusion method was determined to be the best available method for predicting virus pI, since (i) PBRs are predicted independently of the impact on pI, (ii) PBR prediction relies on proteome sequences rather than detailed structural models, and (iii) PBR exclusion was successfully demonstrated on a diverse set of viruses. These models apply to nonenveloped viruses only. A similar model for enveloped viruses is complicated by a lack of data on enveloped virus pI, as well as uncertainties regarding the influence of the phospholipid envelope on charge and ion gradients.

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.

Irreversible and reversible morphological changes in the φ6 capsid and similar viral shells: symmetry and micromechanics

Understanding the physicochemical processes occurring in viruses during their maturation is of fundamental importance since only mature viruses can infect host cells. Here we consider the irreversible and reversible morphological changes that occur with the dodecahedral φ6 procapsid during the sequential packaging of 3 RNA segments forming the viral genome. It is shown that the dodecahedral shape of all the four observed capsid states is perfectly reproduced by a sphere radially deformed by only two irreducible spherical harmonics with icosahedral symmetry and wave numbers l = 6 and l = 10. The rotation of proteins around the 3-fold axes at the Procapsid → Intermediate 1 irreversible transformation is in fact also well described with the shear field containing only two irreducible harmonics with the same two wave numbers. The high stability of the Intermediate 1 state is discussed and the shapes of the Intermediate 2 state and Capsid (reversibly transforming back to the Intermediate 1 state) are shown to be mainly due to the isotropic pressure that the encapsidated RNA segments exert on the shell walls. The hidden symmetry of the capsid and the physicochemical features of the in vitro genome extraction from the viral shell are also elucidated.

Viruses with different genome types adopt a similar strategy to pack nucleic acids based on positively charged protein domains

Capsid proteins often present a positively charged arginine-rich sequence at their terminal regions, which has a fundamental role in genome packaging and particle stability for some icosahedral viruses. These sequences show little to no conservation and are structurally dynamic such that they cannot be easily detected by common sequence or structure comparisons. As a result, the occurrence and distribution of positively charged domains across the viral universe are unknown. Based on the net charge calculation of discrete protein segments, we identified proteins containing amino acid stretches with a notably high net charge (Q > + 17), which are enriched in icosahedral viruses with a distinctive bias towards arginine over lysine. We used viral particle structural data to calculate the total electrostatic charge derived from the most positively charged protein segment of capsid proteins and correlated these values with genome charges arising from the phosphates of each nucleotide. We obtained a positive correlation (r = 0.91, p-value <0001) for a group of 17 viral families, corresponding to 40% of all families with icosahedral structures described to date. These data indicated that unrelated viruses with diverse genome types adopt a common underlying mechanism for capsid assembly based on R-arms.

How a Virus Circumvents Energy Barriers to Form Symmetric Shells

Previous self-assembly experiments on a model icosahedral plant virus have shown that, under physiological conditions, capsid proteins initially bind to the genome through an en masse mechanism and form nucleoprotein complexes in a disordered state, which raises the question as to how virions are assembled into a highly ordered structure in the host cell. Using small-angle X-ray scattering, we find out that a disorder-order transition occurs under physiological conditions upon an increase in capsid protein concentrations. Our cryo-transmission electron microscopy reveals closed spherical shells containing in vitro transcribed viral RNA even at pH 7.5, in marked contrast with the previous observations. We use Monte Carlo simulations to explain this disorder-order transition and find that, as the shell grows, the structures of disordered intermediates in which the distribution of pentamers does not belong to the icosahedral subgroups become energetically so unfavorable that the caps can easily dissociate and reassemble, overcoming the energy barriers for the formation of perfect icosahedral shells. In addition, we monitor the growth of capsids under the condition that the nucleation and growth is the dominant pathway and show that the key for the disorder-order transition in both en masse and nucleation and growth pathways lies in the strength of elastic energy compared to the other forces in the system including protein-protein interactions and the chemical potential of free subunits. Our findings explain, at least in part, why perfect virions with icosahedral order form under different conditions including physiological ones.

A dimorphism shift of hepatitis B virus capsids in response to ionic conditions

The dimorphism of HBV capsids (coexistence of T = 3 and T = 4 capsids) was found to be regulatable by controlling the rate of capsid nucleation using cations such as K+ or Ca2+: a quick addition of highly concentrated monovalent and/or multivalent counter-cations resulted in a morphism transition from a thermodynamically more stable, T = 4 capsid-dominant state (>80% of total capsids) to a new state containing ∼1 : 1 amounts of T = 3 and T = 4 capsids. These results suggested that the salts with strong charge screening ability could narrow the difference in nucleation energy barriers between the two states, which were not inter-convertible once formed. The effect of salts was more significant than other factors such as pH or protein concentration in achieving such a dimorphism shift. The general mechanism of HBV capsid dimorphism described here provides a new perspective in understanding the virus assembly during infection and directing the design of non-infectious capsids for nanotechnology applications.

Electrostatics-Driven Inflation of Elastic Icosahedral Shells as a Model for Swelling of Viruses

We develop a clear theoretical description of radial swelling in virus-like particles that delineates the importance of electrostatic contributions to swelling in the absence of any conformational changes. The model couples the elastic parameters of the capsid-represented as a continuous elastic shell-to the electrostatic pressure acting on it. We show that different modifications of the electrostatic interactions brought about by, for instance, changes in pH or solution ionic strength are often sufficient to achieve the experimentally observed swelling (∼10% of the capsid radius). Additionally, we derive analytical expressions for the electrostatics-driven radial swelling of virus-like particles that enable one to quickly estimate the magnitudes of physical quantities involved.

Systematic analysis of biological roles of charged amino acid residues located throughout the structured inner wall of a virus capsid

Structure-based mutational analysis of viruses is providing many insights into the relationship between structure and biological function of macromolecular complexes. We have systematically investigated the individual biological roles of charged residues located throughout the structured capsid inner wall (outside disordered peptide segments) of a model spherical virus, the minute virus of mice (MVM). The functional effects of point mutations that altered the electrical charge at 16 different positions at the capsid inner wall were analyzed. The results revealed that MVM capsid self-assembly is rather tolerant to point mutations that alter the number and distribution of charged residues at the capsid inner wall. However, mutations that either increased or decreased the number of positive charges around capsid-bound DNA segments reduced the thermal resistance of the virion. Moreover, mutations that either removed or changed the positions of negatively charged carboxylates in rings of acidic residues around capsid pores were deleterious by precluding a capsid conformational transition associated to through-pore translocation events. The results suggest that number, distribution and specific position of electrically charged residues across the inner wall of a spherical virus may have been selected through evolution as a compromise between several different biological requirements.

Adsorption and encapsulation of flexible polyelectrolytes in charged spherical vesicles

We present a theory of adsorption of flexible polyelectrolytes on the interior and exterior surfaces of a charged vesicle in an electrolyte solution. The criteria for adsorption and the density profiles of the adsorbed polymer chain are derived in terms of various characteristics of the polymer, vesicle, and medium, such as the charge density and length of the polymer, charge density and size of the vesicle, electrolyte concentration and dielectric constant of the medium. For adsorption inside the vesicle, the competition between the loss of conformational entropy and gain in adsorption energy results in two kinds of encapsulated states, depending on the strength of the polymer-vesicle interaction. By considering also the adsorption from outside the vesicle, we derive the entropic and energy contributions to the free energy change to transfer an adsorbed chain in the interior to an adsorbed chain on the exterior. In this paper, we have used the Wentzel-Kramers-Brillouin (WKB) method to solve the equation for the probability distribution function of the chain. The present WKB results are compared with the previous results based on variational methods. The WKB and variational results are in good agreement for both the interior and exterior states of adsorption, except in the zero-salt limit for adsorption in the exterior region. The adsorption criteria and density profiles for both the interior and exterior states are presented in terms of various experimentally controllable variables. Calculation of the dependencies of free energy change to transfer an adsorbed chain from the interior to the exterior surface on salt concentration and vesicle radius shows that the free energy penalty to expel a chain from a vesicle is only of the order of thermal energy.

Self consistent field theory of virus assembly

The ground state dominance approximation (GSDA) has been extensively used to study the assembly of viral shells. In this work we employ the self-consistent field theory (SCFT) to investigate the adsorption of RNA onto positively charged spherical viral shells and examine the conditions when GSDA does not apply and SCFT has to be used to obtain a reliable solution. We find that there are two regimes in which GSDA does work. First, when the genomic RNA length is long enough compared to the capsid radius, and second, when the interaction between the genome and capsid is so strong that the genome is basically localized next to the wall. We find that for the case in which RNA is more or less distributed uniformly in the shell, regardless of the length of RNA, GSDA is not a good approximation. We observe that as the polymer-shell interaction becomes stronger, the energy gap between the ground state and first excited state increases and thus GSDA becomes a better approximation. We also present our results corresponding to the genome persistence length obtained through the tangent-tangent correlation length and show that it is zero in case of GSDA but is equal to the inverse of the energy gap when using SCFT.

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.

Varieties of charge distributions in coat proteins of ssRNA+ viruses

A major part of the interactions involved in the assembly and stability of icosahedral, positive-sense single-stranded RNA (ssRNA+) viruses is electrostatic in nature, as can be inferred from the strong pH- and salt-dependence of their assembly phase diagrams. Electrostatic interactions do not act only between the capsid coat proteins (CPs), but just as often provide a significant contribution to the interactions of the CPs with the genomic RNA, mediated to a large extent by positively charged, flexible N-terminal tails of the CPs. In this work, we provide two clear and complementary definitions of an N-terminal tail of a protein, and use them to extract the tail sequences of a large number of CPs of ssRNA+  viruses. We examine the pH-dependent interplay of charge on both tails and CPs alike, and show that-in contrast to the charge on the CPs-the net positive charge on the N-tails persists even to very basic pH values. In addition, we note a limit to the length of the wild-type genomes of those viruses which utilize positively charged tails, when compared to viruses without charged tails and similar capsid size. At the same time, we observe no clear connection between the charge on the N-tails and the genome lengths of the viruses included in our study.

50th Anniversary Perspective: A Perspective on Polyelectrolyte Solutions

From the beginning of life with the information-containing polymers until the present era of a plethora of water-based materials in health care industry and biotechnology, polyelectrolytes are ubiquitous with a broad range of structural and functional properties. The main attribute of polyelectrolyte solutions is that all molecules are strongly correlated both topologically and electrostatically in their neutralizing background of charged ions in highly polarizable solvent. These strong correlations and the necessary use of numerous variables in experiments on polyelectrolytes have presented immense challenges toward fundamental understanding of the various behaviors of charged polymeric systems. This Perspective presents the author's subjective summary of several conceptual advances and the remaining persistent challenges in the contexts of charge and size of polymers, structures in homogeneous solutions, thermodynamic instability and phase transitions, structural evolution with oppositely charged polymers, dynamics in polyelectrolyte solutions, kinetics of phase separation, mobility of charged macromolecules between compartments, and implications to biological systems.

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.

Single Particle Observation of SV40 VP1 Polyanion-Induced Assembly Shows That Substrate Size and Structure Modulate Capsid Geometry

Simian virus 40 capsid protein (VP1) is a unique system for studying substrate-dependent assembly of a nanoparticle. Here, we investigate a simplest case of this system where 12 VP1 pentamers and a single polyanion, e.g., RNA, form a T = 1 particle. To test the roles of polyanion substrate length and structure during assembly, we characterized the assembly products with size exclusion chromatography, transmission electron microscopy, and single-particle resistive-pulse sensing. We found that 500 and 600 nt RNAs had the optimal length and structure for assembly of uniform T = 1 particles. Longer 800 nt RNA, shorter 300 nt RNA, and a linear 600 unit poly(styrene sulfonate) (PSS) polyelectrolyte produced heterogeneous populations of products. This result was surprising as the 600mer PSS and 500-600 nt RNA have similar mass and charge. Like ssRNA, PSS also has a short 4 nm persistence length, but unlike RNA, PSS lacks a compact tertiary structure. These data indicate that even for flexible substrates, shape as well as size affect assembly and are consistent with the hypothesis that work, derived from protein-protein and protein-substrate interactions, is used to compact the substrate.

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.

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.

Collapse and hybridization of RNA: view from replica technique approach

The replica technique method is applied to investigate the kinetic behavior of the coarse-grained model for the RNA molecule. A non-equilibrium phase transition of second order between the glassy phase and the ensemble of freely fluctuating structures has been observed. The non-equilibrium steady state is investigated as well and the thermodynamic characteristics of the system have been evaluated. The non-equilibrium behavior of the specific heat is discussed. Based on our analysis, we point out the state in the kinetic pathway in which the RNA molecule is most prone to hybridization.

The Globular State of the Single-Stranded RNA: Effect of the Secondary Structure Rearrangements

The mutual influence of the slow rearrangements of secondary structure and fast collapse of the long single-stranded RNA (ssRNA) in approximation of coarse-grained model is studied with analytic calculations. It is assumed that the characteristic time of the secondary structure rearrangement is much longer than that for the formation of the tertiary structure. A nonequilibrium phase transition of the 2nd order has been observed.

Inverted critical adsorption of polyelectrolytes in confinement

What are the fundamental laws for the adsorption of charged polymers onto oppositely charged surfaces, for convex, planar, and concave geometries? This question is at the heart of surface coating applications, various complex formation phenomena, as well as in the context of cellular and viral biophysics. It has been a long-standing challenge in theoretical polymer physics; for realistic systems the quantitative understanding is however often achievable only by computer simulations. In this study, we present the findings of such extensive Monte-Carlo in silico experiments for polymer-surface adsorption in confined domains. We study the inverted critical adsorption of finite-length polyelectrolytes in three fundamental geometries: planar slit, cylindrical pore, and spherical cavity. The scaling relations extracted from simulations for the critical surface charge density σc-defining the adsorption-desorption transition-are in excellent agreement with our analytical calculations based on the ground-state analysis of the Edwards equation. In particular, we confirm the magnitude and scaling of σc for the concave interfaces versus the Debye screening length 1/κ and the extent of confinement a for these three interfaces for small κa values. For large κa the critical adsorption condition approaches the known planar limit. The transition between the two regimes takes place when the radius of surface curvature or half of the slit thickness a is of the order of 1/κ. We also rationalize how σc(κ) dependence gets modified for semi-flexible versus flexible chains under external confinement. We examine the implications of the chain length for critical adsorption-the effect often hard to tackle theoretically-putting an emphasis on polymers inside attractive spherical cavities. The applications of our findings to some biological systems are discussed, for instance the adsorption of nucleic acids onto the inner surfaces of cylindrical and spherical viral capsids.

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.

Convexity and Stiffness in Energy Functions for Electrostatic Simulations

We study the properties of convex functionals which have been proposed for the simulation of charged molecular systems within the Poisson-Boltzmann approximation. We consider the extent to which the functionals reproduce the true fluctuations of electrolytes and thus the one-loop correction to mean field theory-including the Debye-Hückel correction to the free energy of ionic solutions. We also compare the functionals for use in numerical optimization of a mean field model of a charged polymer and show that different functionals have very different stiffnesses leading to substantial differences in accuracy and speed.

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.

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.

Multivalent ion effects on electrostatic stability of virus-like nano-shells

Electrostatic properties and stability of charged virus-like nano-shells are examined in ionic solutions with monovalent and multivalent ions. A theoretical model based on a thin charged spherical shell and multivalent ions within the "dressed multivalent ion" approximation, yielding their distribution across the shell and the corresponding electrostatic (osmotic) pressure acting on the shell, is compared with extensive implicit Monte-Carlo simulations. It is found to be accurate for positive or low negative surface charge densities of the shell and for sufficiently high (low) monovalent (multivalent) salt concentrations. Phase diagrams involving electrostatic pressure exhibit positive and negative values, corresponding to an outward and an inward facing force on the shell, respectively. This provides an explanation for the high sensitivity of viral shell stability and self-assembly of viral capsid shells on the ionic environment.

RNA topology remolds electrostatic stabilization of viruses

Simple RNA viruses efficiently encapsulate their genome into a nano-sized protein shell: the capsid. Spontaneous coassembly of the genome and the capsid proteins is driven predominantly by electrostatic interactions between the negatively charged RNA and the positively charged inner capsid wall. Using field theoretic formulation we show that the inherently branched RNA secondary structure allows viruses to maximize the amount of encapsulated genome and make assembly more efficient, allowing viral RNAs to out-compete cellular RNAs during replication in infected host cells.

Mechanical and assembly units of viral capsids identified via quasi-rigid domain decomposition

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

The genetic material of viruses is packaged inside capsids constituted from a few tens to thousands of proteins. The latter can organize in multimers that serve as fundamental blocks for the viral shell assembly or that control the capsid conformational transitions and response to mechanical stress. In this work, we introduce and apply a computational scheme that identifies the fundamental protein blocks from the structural fluctuations of the capsids in thermal equilibrium. These can be derived from phenomenological elastic network models with minimal computational expenditure. Accordingly, the basic functional protein units of a capsid can be obtained from the sole input of the capsid crystal structure. The method is applied to a heterogeneous set of viruses of various size and geometries. These include well-characterised instances for validation purposes, as well as debated ones for which predictions are formulated.

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.

Hepatitis virus capsid polymorph stability depends on encapsulated cargo size

Protein cages providing a controlled environment to encapsulated cargo are a ubiquitous presence in any biological system. Well-known examples are capsids, the regular protein shells of viruses, which protect and deliver the viral genome. Since some virus capsids can be loaded with nongenomic cargoes, they are interesting for a variety of applications ranging from biomedical delivery to energy harvesting. A question of vital importance for such applications is how does capsid stability depend on the size of the cargo? A nanoparticle-templated assembly approach was employed here to determine how different polymorphs of the Hepatitis B virus icosahedral capsid respond to a gradual change in the encapsulated cargo size. It was found that assembly into complete virus-like particles occurs cooperatively around a variety of core diameters, albeit the degree of cooperativity varies. Among these virus-like particles, it was found that those of an outer diameter corresponding to an icosahedral array of 240 proteins (T = 4) are able to accommodate the widest range of cargo sizes.

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

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

Single-stranded RNAs of simple viruses seem to be topologically more compact than other types of single-stranded RNA. It has been suggested that this has an evolutionary purpose: more compact structures are more easily encapsulated in the limited space that the cavity of the virus capsid offers. We employ a simple Flory theory to calculate the optimal amount of polymers confined in a viral shell. We find that the free energy gain or more specifically the efficiency of RNA encapsidation increases substantially with topological compactness. We also find that the optimal length of RNA encapsidated in a capsid increases with the degree of branching of the genome even though this effect is very weak. Further, we show that if the structure of the branching of the polymer is allowed to anneal, the optimal loading increases substantially.

Effect of capsid confinement on the chromatin organization of the SV40 minichromosome

Using small-angle X-ray scattering, we determined the three-dimensional packing architecture of the minichromosome confined within the SV40 virus. In solution, the minichromosome, composed of closed circular dsDNA complexed in nucleosomes, was shown to be structurally similar to cellular chromatin. In contrast, we find a unique organization of the nanometrically encapsidated chromatin, whereby minichromosomal density is somewhat higher at the center of the capsid and decreases towards the walls. This organization is in excellent agreement with a coarse-grained computer model, accounting for tethered nucleosomal interactions under viral capsid confinement. With analogy to confined liquid crystals, but contrary to the solenoid structure of cellular chromatin, our simulations indicate that the nucleosomes within the capsid lack orientational order. Nucleosomes in the layer adjacent to the capsid wall, however, align with the boundary, thereby inducing a 'molten droplet' state of the chromatin. These findings indicate that nucleosomal interactions suffice to predict the genome organization in polyomavirus capsids and underscore the adaptable nature of the eukaryotic chromatin architecture to nanoscale confinement.

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.

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

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

Density functional theory for encapsidated polyelectrolytes: a comparison with Monte Carlo simulation

Genome packaging inside viral capsids is strongly influenced by the molecular size and the backbone structure of RNA∕DNA chains and their electrostatic affinity with the capsid proteins. Coarse-grained models are able to capture the generic features of non-specific interactions and provide a useful testing ground for theoretical developments. In this work, we use the classical density functional theory (DFT) within the framework of an extended primitive model for electrolyte solutions to investigate the self-organization of flexible and semi-flexible linear polyelectrolytes in spherical capsids that are permeable to small ions but not polymer segments. We compare the DFT predictions with Monte Carlo (MC) simulation for the density distributions of polymer segments and small ions at different backbone flexibilities and several solution conditions. In general, the agreement between DFT and MC is near quantitative except when the simulation results are noticeably influenced by the boundary effects. The numerical efficiency of the DFT calculations makes it promising as a useful tool for quantification of the structural and thermodynamic properties of viral nucleocapsids in vivo and at conditions pertinent to experiments.

Interaction of virions with membrane glycolipids

Cellular membranes contain various lipids including glycolipids (GLs). The hydrophilic head groups of GLs extend from the membrane into the aqueous environment outside the cell where they act as recognition sites for specific interactions. The first steps of interaction of virions with cells often include contacts with GLs. To clarify the details of such contacts, we have used the total internal reflection fluorescence microscopy to explore the interaction of individual unlabelled virus-like particles (or, more specifically, norovirus protein capsids), which are firmly bound to a lipid bilayer, and fluorescent vesicles containing glycosphingolipids (these lipids form a subclass of GLs). The corresponding binding kinetics were earlier found to be kinetically limited, while the detachment kinetics were logarithmic over a wide range of time. Here, the detachment rate is observed to dramatically decrease with increasing concentration of glycosphingolipids from 1% to 8%. This effect has been analytically explained by using a generic model describing the statistics of bonds in the contact area between a virion and a lipid membrane. Among other factors, the model takes the formation of GL domains into account. Our analysis indicates that in the system under consideration, such domains, if present, have a characteristic size smaller than the contact area between the vesicle and the virus-like particle.

An examination of the electrostatic interactions between the N-terminal tail of the Brome Mosaic Virus coat protein and encapsidated RNAs

The coat protein of positive-stranded RNA viruses often contains a positively charged tail that extends toward the center of the capsid and interacts with the viral genome. Electrostatic interaction between the tail and the RNA has been postulated as a major force in virus assembly and stabilization. The goal of this work is to examine the correlation between electrostatic interaction and amount of RNA packaged in the tripartite Brome Mosaic Virus (BMV). Nanoindentation experiment using atomic force microscopy showed that the stiffness of BMV virions with different RNAs varied by a range that is 10-fold higher than that would be predicted by electrostatics. BMV mutants with decreased positive charges encapsidated lower amounts of RNA while mutants with increased positive charges packaged additional RNAs up to ∼900 nt. However, the extra RNAs included truncated BMV RNAs, an additional copy of RNA4, potential cellular RNAs, or a combination of the three, indicating that change in the charge of the capsid could result in several different outcomes in RNA encapsidation. In addition, mutant with specific arginines changed to lysines in the capsid also exhibited defects in the specific encapsidation of BMV RNA4. The experimental results indicate that electrostatics is a major component in RNA encapsidation but was unable to account for all of the observed effects on RNA encapsidation. Thermodynamic modeling incorporating the electrostatics was able to predict the approximate length of the RNA to be encapsidated for the majority of mutant virions, but not for a mutant with extreme clustered positive charges. Cryo-electron microscopy of virions that encapsidated an additional copy of RNA4 revealed that, despite the increase in RNA encapsidated, the capsid structure was minimally changed. These results experimentally demonstrated the impact of electrostatics and additional restraints in the encapsidation of BMV RNAs, which could be applicable to other viruses.