Radial distribution of RNA genome packaged inside spherical viruses

The problem of RNA genomes packaged inside spherical viruses is studied. The viral capsid is modeled as a hollowed sphere. The attraction between RNA molecules and the inner viral capsid is assumed to be non-specific and occurs at the inner capsid surface only. For small capsid attraction, it is found that monomer concentration of RNA molecules is maximum at the center of the capsid to maximize their configurational entropy. For stronger capsid attraction, RNA concentration peaks at some distance near the capsid. In the latter case, the competition between the branching of RNA secondary struture and its adsorption to the inner capsid results in the formation of a dense layer of RNA near capsid surface. The layer thickness is a slowly varying (logarithmic) function of the capsid inner radius. Consequently, for immediate strength of RNA-capsid interaction, the amount of RNA packaged inside a virus is proportional to the capsid {\em area} (or the number of proteins) instead of its volume. The numerical profiles describe reasonably well the experimentally observed RNA nucleotide concentration profiles of various viruses.Comment: 5 pages, 2 figures. Abstract, introduction rewritten. Comparison to actual virus profiles added. Submitted to PR

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 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.Comment: 22 pages, 7 figures in J. Phys. Chem. B, 201

Dynamics of Bacteriophage Genome Ejection In Vitro and In Vivo

Bacteriophages, phages for short, are viruses of bacteria. The majority of phages contain a double-stranded DNA genome packaged in a capsid at a density of ~500 mg/ml. This high density requires substantial compression of the normal B form helix, leading to the conjecture that DNA in mature phage virions is under significant pressure, and that pressure is used to eject the DNA during infection. A large number of theoretical, computer simulation and in vitro experimental studies surrounding this conjecture has revealed many --- though often isolated and/or contradictory --- aspects of packaged DNA. This prompts us to present a unified view of the statistical physics and thermodynamics of DNA packaged in phage capsids. We argue that the DNA in a mature phage is in a (meta)stable state, wherein electrostatic self-repulsion is balanced by curvature stress due to confinement in the capsid. We show that in addition to the osmotic pressure associated with the packaged DNA and its counterions, there are four different pressures within the capsid: pressure on the DNA, hydrostatic pressure, the pressure experienced by the capsid, and the pressure associated with the chemical potential of DNA ejection. Significantly, we analyze the mechanism of force transmission in the packaged DNA, and demonstrate that the pressure on DNA is not important for ejection. We derive equations showing a strong hydrostatic pressure difference across the capsid shell. We propose that when a phage is triggered to eject by interaction with its receptor in vitro, the (thermodynamic) incentive of water molecules to enter the phage capsid flushes the DNA out of the capsid. In vivo, the difference between the osmotic pressures in the bacterial cell cytoplasm and the culture medium similarly results in a water flow that drags the DNA out of the capsid and into the bacterial cell.Comment: 28 pages, 6 figures, invited article to Physical Biology, to appea

Density waves theory of the capsid structure of small icosahedral viruses

We apply Landau theory of crystallization to explain and to classify the capsid structures of small viruses with spherical topology and icosahedral symmetry. We develop an explicit method which predicts the positions of centers of mass for the proteins constituting viral capsid shell. Corresponding density distribution function which generates the positions has universal form without any fitting parameter. The theory describes in a uniform way both the structures satisfying the well-known Caspar and Klug geometrical model for capsid construction and those violating it. The quasiequivalence of protein environments in viral capsid and peculiarities of the assembly thermodynamics are also discussed.Comment: 8 pages, 3 figur

Dynamics of polymer ejection from capsid

Polymer ejection from a capsid through a nanoscale pore is an important biological process with relevance to modern biotechnology. Here, we study generic capsid ejection using Langevin dynamics. We show that even when the ejection takes place within the drift-dominated region there is a very high probability for the ejection process not to be completed. Introducing a small aligning force at the pore entrance enhances ejection dramatically. Such a pore asymmetry is a candidate for a mechanism by which a viral ejection is completed. By detailed high-resolution simulations we show that such capsid ejection is an out-of-equilibrium process that shares many common features with the much studied driven polymer translocation through a pore in a wall or a membrane. We find that the escape times scale with polymer length, $\tau \sim N^\alpha$. We show that for the pore without the asymmetry the previous predictions corroborated by Monte Carlo simulations do not hold. For the pore with the asymmetry the scaling exponent varies with the initial monomer density (monomers per capsid volume) $\rho$ inside the capsid. For very low densities $\rho \le 0.002$ the polymer is only weakly confined by the capsid, and we measure $\alpha = 1.33$, which is close to $\alpha = 1.4$ obtained for polymer translocation. At intermediate densities the scaling exponents $\alpha = 1.25$ and $1.21$ for $\rho = 0.01$ and $0.02$, respectively. These scalings are in accord with a crude derivation for the lower limit $\alpha = 1.2$. For the asymmetrical pore precise scaling breaks down, when the density exceeds the value for complete confinement by the capsid, $\rho \gtrapprox 0.25$. The high-resolution data show that the capsid ejection for both pores, analogously to polymer translocation, can be characterized as a multiplicative stochastic process that is dominated by small-scale transitions.Comment: 10 pages, 6 figure

Tolerance to mutations in the foot-and-mouth disease virus integrin-binding RGD region is different in cultured cells and in vivo and depends on the capsid sequence context.

Engineered RNAs carrying substitutions in the integrin receptor-binding Arg-Gly-Asp (RGD) region of foot-and-mouth disease virus (FMDV) were constructed (aa 141-147 of VP1 capsid protein) and their infectivity was assayed in cultured cells and suckling mice. The effect of these changes was studied in the capsid proteins of two FMDVs, C-S8c1, which enters cells through integrins, and 213hs(-), a derivative highly adapted to cell culture whose ability to infect cells using the glycosaminoglycan heparan sulfate (HS) as receptor, acquired by multiple passage on BHK-21 cells, has been abolished. The capsid sequence context determined infectivity in cultured cells and directed the selection of additional replacements in structural proteins. Interestingly, a viral population derived from a C-S8c1/L144A mutant, carrying only three substitutions in the capsid, was able to expand tropism to wild-type (wt) and mutant (mt)glycosaminoglycan-deficient CHO cells. In contrast, the 213hs(-) capsid tolerated all substitutions analysed with no additional mutations, and the viruses recovered maintained the ability of the 213hs(-) parental virus to infect wt and mt CHO cells. Viruses derived from C-S8c1 with atypical RGD regions were virulent and transmissible for mice with no other changes in the capsid. Substitution of Asp143 for Ala in the C-S8c1 capsid eliminated infectivity in cultured cells and mice. Co-inoculation with a neutralizing monoclonal antibody directed against the type C FMDV RGD region abolished infectivity of C-S8c1 virus on suckling mice, suggesting that FMDV can infect mice using integrins. Sequence requirements imposed for viral entry in vitro and in vivo are discussed

Self-assembly of model proteins into virus capsids

We consider self-assembly of proteins into a virus capsid by the methods of molecular dynamics. The capsid corresponds either to SPMV or CCMV and is studied with and without the RNA molecule inside. The proteins are flexible and described by the structure-based coarse-grained model augmented by electrostatic interactions. Previous studies of the capsid self-assembly involved solid objects of a supramolecular scale, e.g. corresponding to capsomeres, with engineered couplings and stochastic movements. In our approach, a single capsid is dissociated by an application of a high temperature for a variable period and then the system is cooled down to allow for self-assembly. The restoration of the capsid proceeds to various extent, depending on the nature of the dissociated state, but is rarely complete because some proteins depart too far unless the process takes place in a confined space.Comment: 13 figure