Viral genome structures are optimal for capsid assembly

Disassembly of Virus-Like Particles and the Stabilizing Role of the Nucleic Acid Cargo

In many simple viruses and virus-like particles, the protein capsid self-assembles around a nucleic-acid genome. Although the assembly process has been studied in detail, relatively little is known about how the capsid disassembles, a potentially important step for infection (in viruses) or cargo delivery (in virus-like particles). We investigate capsid disassembly using a coarse-grained molecular dynamics model of a T = 1 dodecahedral capsid and an RNA-like polymer. We alter the interactions between the subunits of the capsid as well as the ionic strength of the solution to induce partial or complete disassembly of self-assembled particles. We find that disassembly follows nucleation-and-growth kinetics, where the nucleation barrier is related to the interaction strengths as well as to the conformation of the polymer. In particular, we find that polymer segments that interact with adjacent subunits reinforce the subunit-subunit contacts. These results have implications for the design of virus-like particles for applications such as drug delivery. A cargo designed with reinforcement in mind might be used to control the stability of such particles and mediate disassembly.

Resolving viral structural complexity by super-resolution microscopy

In this review, we discuss different super-resolution microscopy (SRM) techniques employed to study viral structures and virus composition with nanometric resolution. We describe the basic principles of the different microscopy methods utilized to break the light diffraction limit, enabling the study of protein composition in viral structures. Finally, we demonstrate for the first time the differential spatial distribution of two structural proteins in an individual baculovirus using single-molecule super-resolution microscopy. We discuss the future of these powerful methods for virology, medicine, and biotechnology applications.

Molecular Dynamics Characterization of the Free and Encapsidated RNA2 of CCMV with the oxRNA Model

The cowpea chlorotic mottle virus (CCMV) has emerged as a model system to assess the balance between electrostatic and topological features of single-stranded RNA viruses, specifically in the context of the viral self-assembly. Yet, despite its biophysical significance, little structural data on the RNA content of the CCMV virion is available. Here, the conformational dynamics of the RNA2 fragment of CCMV was assessed via coarse-grained molecular dynamics simulations, employing the oxRNA2 force field. The behavior of RNA2 was characterized both as a freely-folding molecule and within a mean-field depiction of the capsid. For the former, the role of the salt concentration, the temperature and of ad hoc constraints on the RNA termini was verified on the equilibrium properties of RNA2. For the latter, a multi-scale approach was employed to derive a potential profile of the viral cavity from atomistic structures of the CCMV capsid in solution. The conformational ensembles of the encapsidated RNA2 were significantly altered with respect to the freely-folding counterparts, as shown by the emergence of long-range motifs and pseudoknots. Finally, the role of the N-terminal tails of the CCMV subunits is highlighted as a critical feature in the construction of a proper electrostatic model of the CCMV capsid.

Cryo-EM of human rhinovirus reveals capsid-RNA duplex interactions that provide insights into virus assembly and genome uncoating

The cryo-EM structure of the human rhinovirus B14 determined in this study reveals 13-bp RNA duplexes symmetrically bound to regions around each of the 30 two-fold axes in the icosahedral viral capsid. The RNA duplexes (~12% of the ssRNA genome) define a quasi-dodecahedral cage that line a substantial part of the capsid interior surface. The RNA duplexes establish a complex network of non-covalent interactions with pockets in the capsid inner wall, including coulombic interactions with a cluster of basic amino acid residues that surround each RNA duplex. A direct comparison was made between the cryo-EM structure of RNA-filled virions and that of RNA-free (empty) capsids that resulted from genome release from a small fraction of viruses. The comparison reveals that some specific residues involved in capsid-duplex RNA interactions in the virion undergo remarkable conformational rearrangements upon RNA release from the capsid. RNA release is also associated with the asynchronous opening of channels at the 30 two-fold axes. The results provide further insights into the molecular mechanisms leading to assembly of rhinovirus particles and their genome uncoating during infection. They may also contribute to development of novel antiviral strategies aimed at interfering with viral capsid-genome interactions during the infectious cycle.

Glass-like Relaxation Dynamics during the Disorder-to-Order Transition of Viral Nucleocapsids

Nucleocapsid self-assembly is an essential yet elusive step in virus replication. Using time-resolved small-angle X-ray scattering on a model icosahedral ssRNA virus, we reveal a previously unreported kinetic pathway. Initially, RNA-bound capsid subunits rapidly accumulate beyond the stoichiometry of native virions. This is followed by a disorder-to-order transition characterized by glass-like relaxation dynamics and the release of excess subunits. Our molecular dynamics simulations, employing a coarse-grained elastic model, confirm the physical feasibility of self-ordering accompanied by subunit release. The relaxation can be modeled by an exponential integral decay on the mean squared radius of gyration, with relaxation times varying within the second range depending on RNA type and subunit concentration. A nanogel model suggests that the initially disordered nucleoprotein complexes quickly reach an equilibrium size, while their mass fractal dimension continues to evolve. Understanding virus self-assembly is not only crucial for combating viral infections, but also for designing synthetic virus-inspired nanocages for drug delivery applications.

Charge Detection Mass Spectrometry Reveals Favored Structures in the Assembly of Virus-Like Particles: Polymorphism in Norovirus GI.1

The main capsid protein (CP) of norovirus, the leading cause of gastroenteritis, is expected to self-assemble into virus-like particles with the same structure as the wild-type virus, a capsid with 180 CPs in a T = 3 icosahedron. Using charge detection mass spectrometry (CD-MS), we find that the norovirus GI.1 variant is structurally promiscuous, forming a wide variety of well-defined structures, some that are icosahedral capsids and others that are not. The structures that are present evolve with time and vary with solution conditions. The presence of icosahedral T = 3 and T = 4 capsids (240 CPs) under some conditions was confirmed by cryo-electron microscopy (cryo-EM). The cryo-EM studies also confirmed the presence of an unexpected prolate geometry based on an elongated T = 4 capsid with 300 CPs. In addition, CD-MS measurements indicate the presence of well-defined peaks with masses corresponding to 420, 480, 600, and 700 CPs. The peak corresponding to 420 CPs is probably due to an icosahedral T = 7 capsid, but this could not be confirmed by cryo-EM. It is possible that the T = 7 particles are too fragile to survive vitrification. There are no mass peaks associated with the T = 9 and T = 12 icosahedra with 540 and 720 CPs. The larger structures with 480, 600, and 700 CPs are not icosahedral; however, their measured charges suggest that they are hollow shells. The use of CD-MS to monitor virus-like particles assembly may have important applications in vaccine development and quality control.

Trends and challenges in electroanalytical biosensing methodologies for infectious viral diseases

Viral pandemic diseases have disruptive global consequences leading to millions of deaths and a severe impact on the global economy. Inadequate preventative protocols have led to an overwhelming demand for intensive care leading to uncontrollable burdens and even breakdown of healthcare sectors across many countries. The rapid detection of viral disease helps in the understanding of the relevant intricacies, helping to tackle infection with improved guidelines. Portable biosensor devices offer promising solutions by facilitating on-site detection of viral pathogens. This review summarizes the latest innovative strategies reported using electroanalytical methods for the screening of viral antigens. The structural components of viruses and their categories are presented followed by the various recognition elements and transduction techniques involved in biosensors. Core sections focus on biosensors reported for viral genomic detection(DNA and RNA) and antigenic capsid protein. Strategies for addressing the challenges of electroanalytical biosensing of viral components are also presented. The advantages, and disadvantages of biorecognition elements and nanozymes for the detection of viral disease are highlighted. Such technical insights will help researchers working in chemistry, and biochemistry as well as clinicians working in medical diagnostics.

Synthesis, insertion, and characterization of SARS-CoV-2 membrane protein within lipid bilayers

Throughout history, coronaviruses have posed challenges to both public health and the global economy; nevertheless, methods to combat them remain rudimentary, primarily due to the absence of experiments to understand the function of various viral components. Among these, membrane (M) proteins are one of the most elusive because of their small size and challenges with expression. Here, we report the development of an expression system to produce tens to hundreds of milligrams of M protein per liter of Escherichia coli culture. These large yields render many previously inaccessible structural and biophysical experiments feasible. Using cryo-electron microscopy and atomic force microscopy, we image and characterize individual membrane-incorporated M protein dimers and discover membrane thinning in the vicinity, which we validated with molecular dynamics simulations. Our results suggest that the resulting line tension, along with predicted induction of local membrane curvature, could ultimately drive viral assembly and budding.

The structure and physical properties of a packaged bacteriophage particle

A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.

Formalizing Coarse-Grained Representations of Anisotropic Interactions at Multimeric Protein Interfaces Using Virtual Sites

Molecular simulations of biomacromolecules that assemble into multimeric complexes remain a challenge due to computationally inaccessible length and time scales. Low-resolution and implicit-solvent coarse-grained modeling approaches using traditional nonbonded interactions (both pairwise and spherically isotropic) have been able to partially address this gap. However, these models may fail to capture the complex anisotropic interactions present at macromolecular interfaces unless higher-order interaction potentials are incorporated at the expense of the computational cost. In this work, we introduce an alternate and systematic approach to represent directional interactions at protein-protein interfaces by using virtual sites restricted to pairwise interactions. We show that virtual site interaction parameters can be optimized within a relative entropy minimization framework by using only information from known statistics between coarse-grained sites. We compare our virtual site models to traditional coarse-grained models using two case studies of multimeric protein assemblies and find that the virtual site models predict pairwise correlations with higher fidelity and, more importantly, assembly behavior that is morphologically consistent with experiments. Our study underscores the importance of anisotropic interaction representations and paves the way for more accurate yet computationally efficient coarse-grained simulations of macromolecular assembly in future research.

Hierarchical assembly is more robust than egalitarian assembly in synthetic capsids

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

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

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

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.

Kinetic Growth of Multicomponent Microcompartment Shells

An important goal of systems and synthetic biology is to produce high value chemical species in large quantities. Microcompartments, which are protein nanoshells encapsulating catalytic enzyme cargo, could potentially function as tunable nanobioreactors inside and outside cells to generate these high value species. Modifying the morphology of microcompartments through genetic engineering of shell proteins is one viable strategy to tune cofactor and metabolite access to encapsulated enzymes. However, this is a difficult task without understanding how changing interactions between the many different types of shell proteins and enzymes affect microcompartment assembly and shape. Here, we use multiscale molecular dynamics and experimental data to describe assembly pathways available to microcompartments composed of multiple types of shell proteins with varied interactions. As the average interaction between the enzyme cargo and the multiple types of shell proteins is weakened, the shell assembly pathway transitions from (i) nucleating on the enzyme cargo to (ii) nucleating in the bulk and then binding the cargo as it grows to (iii) an empty shell. Atomistic simulations and experiments using the 1,2-propanediol utilization microcompartment system demonstrate that shell protein interactions are highly varied and consistent with our multicomponent, coarse-grained model. Furthermore, our results suggest that intrinsic bending angles control the size of these microcompartments. Overall, our simulations and experiments provide guidance to control microcomparmtent size and assembly by modulating the interactions between shell proteins.

Understanding Virus Structure and Dynamics through Molecular Simulations

Viral outbreaks remain a serious threat to human and animal populations and motivate the continued development of antiviral drugs and vaccines, which in turn benefits from a detailed understanding of both viral structure and dynamics. While great strides have been made in characterizing these systems experimentally, molecular simulations have proven to be an essential, complementary approach. In this work, we review the contributions of molecular simulations to the understanding of viral structure, functional dynamics, and processes related to the viral life cycle. Approaches ranging from coarse-grained to all-atom representations are discussed, including current efforts at modeling complete viral systems. Overall, this review demonstrates that computational virology plays an essential role in understanding these systems.

Self-assembly of dengue virus empty capsid-like particles in solution

Nucleocapsid (NC) assembly is an essential step of the virus replication cycle. It ensures genome protection and transmission among hosts. Flaviviruses are human viruses for which envelope structure is well known, whereas no information on NC organization is available. Here we designed a dengue virus capsid protein (DENVC) mutant in which a highly positive spot conferred by arginine 85 in α4-helix was replaced by a cysteine residue, simultaneously removing the positive charge and restricting the intermolecular motion through the formation of a disulfide cross-link. We showed that the mutant self-assembles into capsid-like particles (CLP) in solution without nucleic acids. Using biophysical techniques, we investigated capsid assembly thermodynamics, showing that an efficient assembly is related to an increased DENVC stability due to α4/α4' motion restriction. To our knowledge, this is the first time that flaviviruses' empty capsid assembly is obtained in solution, revealing the R85C mutant as a powerful tool to understand the NC assembly mechanism.

A binary interaction map between turnip mosaic virus and Arabidopsis thaliana proteomes

Viruses are obligate intracellular parasites that have co-evolved with their hosts to establish an intricate network of protein-protein interactions. Here, we followed a high-throughput yeast two-hybrid screening to identify 378 novel protein-protein interactions between turnip mosaic virus (TuMV) and its natural host Arabidopsis thaliana. We identified the RNA-dependent RNA polymerase NIb as the viral protein with the largest number of contacts, including key salicylic acid-dependent transcription regulators. We verified a subset of 25 interactions in planta by bimolecular fluorescence complementation assays. We then constructed and analyzed a network comprising 399 TuMV-A. thaliana interactions together with intravirus and intrahost connections. In particular, we found that the host proteins targeted by TuMV are enriched in different aspects of plant responses to infections, are more connected and have an increased capacity to spread information throughout the cell proteome, display higher expression levels, and have been subject to stronger purifying selection than expected by chance. The proviral or antiviral role of ten host proteins was validated by characterizing the infection dynamics in the corresponding mutant plants, supporting a proviral role for the transcriptional regulator TGA1. Comparison with similar studies with animal viruses, highlights shared fundamental features in their mode of action.

Bottom-up Coarse-Graining: Principles and Perspectives

Large-scale computational molecular models provide scientists a means to investigate the effect of microscopic details on emergent mesoscopic behavior. Elucidating the relationship between variations on the molecular scale and macroscopic observable properties facilitates an understanding of the molecular interactions driving the properties of real world materials and complex systems (e.g., those found in biology, chemistry, and materials science). As a result, discovering an explicit, systematic connection between microscopic nature and emergent mesoscopic behavior is a fundamental goal for this type of investigation. The molecular forces critical to driving the behavior of complex heterogeneous systems are often unclear. More problematically, simulations of representative model systems are often prohibitively expensive from both spatial and temporal perspectives, impeding straightforward investigations over possible hypotheses characterizing molecular behavior. While the reduction in resolution of a study, such as moving from an atomistic simulation to that of the resolution of large coarse-grained (CG) groups of atoms, can partially ameliorate the cost of individual simulations, the relationship between the proposed microscopic details and this intermediate resolution is nontrivial and presents new obstacles to study. Small portions of these complex systems can be realistically simulated. Alone, these smaller simulations likely do not provide insight into collectively emergent behavior. However, by proposing that the driving forces in both smaller and larger systems (containing many related copies of the smaller system) have an explicit connection, systematic bottom-up CG techniques can be used to transfer CG hypotheses discovered using a smaller scale system to a larger system of primary interest. The proposed connection between different CG systems is prescribed by (i) the CG representation (mapping) and (ii) the functional form and parameters used to represent the CG energetics, which approximate potentials of mean force (PMFs). As a result, the design of CG methods that facilitate a variety of physically relevant representations, approximations, and force fields is critical to moving the frontier of systematic CG forward. Crucially, the proposed connection between the system used for parametrization and the system of interest is orthogonal to the optimization used to approximate the potential of mean force present in all systematic CG methods. The empirical efficacy of machine learning techniques on a variety of tasks provides strong motivation to consider these approaches for approximating the PMF and analyzing these approximations.

Biophysical Modeling of SARS-CoV-2 Assembly: Genome Condensation and Budding

The COVID-19 pandemic caused by the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) has spurred unprecedented and concerted worldwide research to curtail and eradicate this pathogen. SARS-CoV-2 has four structural proteins: Envelope (E), Membrane (M), Nucleocapsid (N), and Spike (S), which self-assemble along with its RNA into the infectious virus by budding from intracellular lipid membranes. In this paper, we develop a model to explore the mechanisms of RNA condensation by structural proteins, protein oligomerization and cellular membrane-protein interactions that control the budding process and the ultimate virus structure. Using molecular dynamics simulations, we have deciphered how the positively charged N proteins interact and condense the very long genomic RNA resulting in its packaging by a lipid envelope decorated with structural proteins inside a host cell. Furthermore, considering the length of RNA and the size of the virus, we find that the intrinsic curvature of M proteins is essential for virus budding. While most current research has focused on the S protein, which is responsible for viral entry, and it has been motivated by the need to develop efficacious vaccines, the development of resistance through mutations in this crucial protein makes it essential to elucidate the details of the viral life cycle to identify other drug targets for future therapy. Our simulations will provide insight into the viral life cycle through the assembly of viral particles de novo and potentially identify therapeutic targets for future drug development.

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.

Transient RNA Interactions Leave a Covalent Imprint on a Viral Capsid Protein

The hepatitis B virus (HBV) is the leading cause of persistent liver infections. Its DNA-based genome is synthesized through reverse transcription of an RNA template inside the assembled capsid shell. In addition to the structured assembly domain, the capsid protein harbors a C-terminal extension that mediates both the enclosure of RNA during capsid assembly and the nuclear entry of the capsid during infection. The arginine-rich motifs within this extension, though common to many viruses, have largely escaped atomic-scale investigation. Here, we leverage solution and solid-state nuclear magnetic resonance spectroscopy at ambient and cryogenic temperatures, under dynamic nuclear polarization signal enhancement, to investigate the organization of the genome within the capsid. Transient interactions with phosphate groups of the RNA backbone confine the arginine-rich motifs to the interior capsid space. While no secondary structure is induced in the C-terminal extension, interactions with RNA counteract the formation of a disulfide bond, which covalently tethers this peptide arm onto the inner capsid surface. Electrostatic and covalent contributions thus compete in the spatial regulation of capsid architecture. This disulfide switch represents a coupling mechanism between the structured assembly domain of the capsid and the enclosed nucleic acids. In particular, it enables the redox-dependent regulation of the exposure of the C-terminal extension on the capsid surface, which is required for nuclear uptake of the capsid. Phylogenetic analysis of capsid proteins from hepadnaviruses points toward a function of this switch in the persistence of HBV infections.

Cargo Release from Nonenveloped Viruses and Virus-like Nanoparticles: Capsid Rupture or Pore Formation

Virus-like nanoparticles are protein shells similar to wild-type viruses, and both aim to deliver their content into a cell. Unfortunately, the release mechanism of their cargo/genome remains elusive. Pores on the symmetry axes were proposed to enable the slow release of the viral genome. In contrast, cryo-EM images showed that capsids of nonenveloped RNA viruses can crack open and rapidly release the genome. We combined in vitro cryo-EM observations of the genome release of three viruses with coarse-grained simulations of generic virus-like nanoparticles to investigate the cargo/genome release pathways. Simulations provided details on both slow and rapid release pathways, including the success rates of individual releases. Moreover, the simulated structures from the rapid release pathway were in agreement with the experiment. Slow release occurred when interactions between capsid subunits were long-ranged, and the cargo/genome was noncompact. In contrast, rapid release was preferred when the interaction range was short and/or the cargo/genome was compact. These findings indicate a design strategy of virus-like nanoparticles for drug delivery.

Phase transitions in virology

Viruses have established relationships with almost every other living organism on Earth and at all levels of biological organization: from other viruses up to entire ecosystems. In most cases, they peacefully coexist with their hosts, but in most relevant cases, they parasitize them and induce diseases and pandemics, such as the AIDS and the most recent avian influenza and COVID-19 pandemic events, causing a huge impact on health, society, and economy. Viruses play an essential role in shaping the eco-evolutionary dynamics of their hosts, and have been also involved in some of the major evolutionary innovations either by working as vectors of genetic information or by being themselves coopted by the host into their genomes. Viruses can be studied at different levels of biological organization, from the molecular mechanisms of genome replication, gene expression and encapsidation, to global pandemics. All these levels are different and yet connected through the presence of threshold conditions allowing for the formation of a capsid, the loss of genetic information or epidemic spreading. These thresholds, as occurs with temperature separating phases in a liquid, define sharp qualitative types of behaviour. Thesephase transitionsare very well known in physics. They have been studied by means of simple, but powerful models able to capture their essential properties, allowing us to better understand them. Can the physics of phase transitions be an inspiration for our understanding of viral dynamics at different scales? Here we review well-known mathematical models of transition phenomena in virology. We suggest that the advantages of abstract, simplified pictures used in physics are also the key to properly understanding the origins and evolution of complexity in viruses. By means of several examples, we explore this multilevel landscape and how minimal models provide deep insights into a diverse array of problems. The relevance of these transitions in connecting dynamical patterns across scales and their evolutionary and clinical implications are outlined.

SARS-CoV-2 nucleocapsid protein forms condensates with viral genomic RNA

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) infection causes Coronavirus Disease 2019 (COVID-19), a pandemic that seriously threatens global health. SARS-CoV-2 propagates by packaging its RNA genome into membrane enclosures in host cells. The packaging of the viral genome into the nascent virion is mediated by the nucleocapsid (N) protein, but the underlying mechanism remains unclear. Here, we show that the N protein forms biomolecular condensates with viral genomic RNA both in vitro and in mammalian cells. While the N protein forms spherical assemblies with homopolymeric RNA substrates that do not form base pairing interactions, it forms asymmetric condensates with viral RNA strands. Cross-linking mass spectrometry (CLMS) identified a region that drives interactions between N proteins in condensates, and deletion of this region disrupts phase separation. We also identified small molecules that alter the size and shape of N protein condensates and inhibit the proliferation of SARS-CoV-2 in infected cells. These results suggest that the N protein may utilize biomolecular condensation to package the SARS-CoV-2 RNA genome into a viral particle.

RNA Structures and Their Role in Selective Genome Packaging

To generate infectious viral particles, viruses must specifically select their genomic RNA from milieu that contains a complex mixture of cellular or non-genomic viral RNAs. In this review, we focus on the role of viral encoded RNA structures in genome packaging. We first discuss how packaging signals are constructed from local and long-range base pairings within viral genomes, as well as inter-molecular interactions between viral and host RNAs. Then, how genome packaging is regulated by the biophysical properties of RNA. Finally, we examine the impact of RNA packaging signals on viral evolution.

Molecular determinants of Ebola nucleocapsid stability from molecular dynamics simulations

Ebola virus (EBOV) is a human pathogen with the ability to cause hemorrhagic fever and bleeding diathesis in hosts. The life cycle of EBOV depends on its nucleocapsid. The Ebola nucleocapsid consists of a helical assembly of nucleoproteins (NPs) encapsidating single-stranded viral RNA (ssRNA). Knowledge of the molecular determinants of Ebola nucleocapsid stability is essential for the development of therapeutics against EBOV. However, large degrees of freedom associated with the Ebola nucleocapsid helical assembly pose a computational challenge, thereby limiting the previous simulation studies to the level of monomers. In the present work, we have performed all atom molecular dynamics (MD) simulations of the helical assembly of EBOV nucleoproteins in the absence and presence of ssRNA. We found that ssRNA is essential for maintaining structural integrity of the nucleocapsid. Other molecular determinants observed to stabilize the nucleocapsid include NP-RNA and NP-NP interactions and ion distributions. Additionally, the structural and dynamical behavior of the nucleocapsid monomer depends on its position in the helical assembly. NP monomers present on the longitudinal edges of the helical tube are more exposed, flexible, and have weaker NP-NP interactions than those residing in the center. This work provides key structural features stabilizing the nucleocapsid that may serve as therapeutic targets.

Scaling relation between genome length and particle size of viruses provides insights into viral life history

In terms of genome and particle sizes, viruses exhibit great diversity. With the discovery of several nucleocytoplasmic large DNA viruses (NCLDVs) and jumbo phages, the relationship between particle and genome sizes has emerged as an important criterion for understanding virus evolution. We use allometric scaling of capsid volume with the genome length of different groups of viruses to shed light on its relationship with virus life history. The allometric exponents for icosahedral dsDNA bacteriophages and NCDLVs were found to be 1 and 2, respectively, indicating that with increasing capsid size DNA packaging density remains the same in bacteriophages but decreases for NCLDVs. We argue that the exponents are largely shaped by their entry mechanism and capsid mechanical stability. We further show that these allometric size parameters are also intricately linked to the relative energy costs of translation and replication in viruses and can have further implications on viral life history.

•

Capsid and genome size allometric exponent gives insights into viral life history

•

The allometric exponent of NCLDVs is almost twice that of bacteriophages

•

The exponent is largely shaped by the viral entry mechanism and capsid stability

•

The relaxed genome size constraint allows large viruses to evolve greater autonomy

Atomistic dynamics of a viral infection process: Release of membrane lytic peptides from a non-enveloped virus

Molecular simulations have played an instrumental role in uncovering the structural dynamics and physical properties of virus capsids. In this work, we move beyond equilibrium physicochemical characterization of a virus system to study a stage of the infection process that is required for viral proliferation. Despite many biochemical and functional studies, the molecular mechanism of host cell entry by non-enveloped viruses remains largely unresolved. Flock House virus (FHV) is a model system for non-enveloped viruses and is the subject of the current study. FHV infects through the acid-dependent endocytic pathway, where low pH triggers externalization of membrane-disrupting (γ) peptides from the capsid interior. Using all-atom equilibrium and enhanced sampling simulations, the mechanism and energetics of γ peptide liberation and the effect of pH on this process are investigated. Our computations agree with experimental findings and reveal nanoscopic details regarding the pH control mechanism, which are not readily accessible in experiments.

SARS-CoV-2 nucleocapsid protein forms condensates with viral genomic RNA

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection causes COVID-19, a pandemic that seriously threatens global health. SARS-CoV-2 propagates by packaging its RNA genome into membrane enclosures in host cells. The packaging of the viral genome into the nascent virion is mediated by the nucleocapsid (N) protein, but the underlying mechanism remains unclear. Here, we show that the N protein forms biomolecular condensates with viral genomic RNA both in vitro and in mammalian cells. Phase separation is driven, in part, by hydrophobic and electrostatic interactions. While the N protein forms spherical assemblies with unstructured RNA, it forms asymmetric condensates with viral RNA strands that contain secondary structure elements. Cross-linking mass spectrometry identified a region that forms interactions between N proteins in condensates, and truncation of this region disrupts phase separation. We also identified small molecules that alter the formation of N protein condensates. These results suggest that the N protein may utilize biomolecular condensation to package the SARS-CoV-2 RNA genome into a viral particle.

An Efficient Brome mosaic virus-Based Gene Silencing Protocol for Hexaploid Wheat (Triticum aestivum L.)

Virus-induced gene silencing (VIGS) is a rapid and powerful method to evaluate gene function, especially for species like hexaploid wheat that have large, redundant genomes and are difficult and time-consuming to transform. The Brome mosaic virus (BMV)-based VIGS vector is widely used in monocotyledonous species but not wheat. Here we report the establishment of a simple and effective VIGS procedure in bread wheat using BMVCP5, the most recently improved BMV silencing vector, and wheat genes PHYTOENE DESATURASE (TaPDS) and PHOSPHATE2 (TaPHO2) as targets. Time-course experiments revealed that smaller inserts (~100 nucleotides, nt) were more stable in BMVCP5 and conferred higher silencing efficiency and longer silencing duration, compared with larger inserts. When using a 100-nt insert and a novel coleoptile inoculation method, BMVCP5 induced extensive silencing of TaPDS transcript and a visible bleaching phenotype in the 2nd to 5th systemically-infected leaves from nine to at least 28 days post inoculation (dpi). For TaPHO2, the ability of BMVCP5 to simultaneously silence all three homoeologs was demonstrated. To investigate the feasibility of BMV VIGS in wheat roots, ectopically expressed enhanced GREEN FLUORESCENT PROTEIN (eGFP) in a transgenic wheat line was targeted for silencing. Silencing of eGFP fluorescence was observed in both the maturation and elongation zones of roots. BMVCP5 mediated significant silencing of eGFP and TaPHO2 mRNA expression in roots at 14 and 21 dpi, and TaPHO2 silencing led to the doubling of inorganic phosphate concentration in the 2nd through 4th systemic leaves. All 54 wheat cultivars screened were susceptible to BMV infection. BMVCP5-mediated TaPDS silencing resulted in the expected bleaching phenotype in all eight cultivars examined, and decreased TaPDS transcript was detected in all three cultivars examined. This BMVCP5 VIGS technology may serve as a rapid and effective functional genomics tool for high-throughput gene function studies in aerial and root tissues and in many wheat cultivars.

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.

Assembly and Stability of Simian Virus 40 Polymorphs

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

Revealing in real-time a multistep assembly mechanism for SV40 virus-like particles

Many viruses use their genome as template for self-assembly into an infectious particle. However, this reaction remains elusive because of the transient nature of intermediate structures. To elucidate this process, optical tweezers and acoustic force spectroscopy are used to follow viral assembly in real time. Using Simian virus 40 (SV40) virus-like particles as model system, we reveal a multistep assembly mechanism. Initially, binding of VP1 pentamers to DNA leads to a significantly decreased persistence length. Moreover, the pentamers seem able to stabilize DNA loops. Next, formation of interpentamer interactions results in intermediate structures with reduced contour length. These structures stabilize into objects that permanently decrease the contour length to a degree consistent with DNA compaction in wild-type SV40. These data indicate that a multistep mechanism leads to fully assembled cross-linked SV40 particles. SV40 is studied as drug delivery system. Our insights can help optimize packaging of therapeutic agents in these particles.

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.

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.

RNA Homopolymers Form Higher-Curvature Virus-like Particles Than Do Normal-Composition RNAs

Unlike double-stranded DNA, single-stranded RNA can be spontaneously packaged into spherical capsids by viral capsid protein (CP) because it is a more compact and flexible polymer. Many systematic investigations of this self-assembly process have been carried out using CP from cowpea chlorotic mottle virus, with a wide range of sequences and lengths of single-stranded RNA. Among these studies are measurements of the relative packaging efficiencies of these RNAs into spherical capsids. In this work, we address a fundamental issue that has received very little attention, namely the question of the preferred curvature of the capsid formed around different RNA molecules. We show in particular that homopolymers of RNA-polyribouridylic acid and polyriboadenylic acid-form exclusively T = 2-sized (∼22-nm diameter) virus-like particles (VLPs) when mixed with cowpea chlorotic mottle virus CP, independent of their length, ranging from 500 to more than 4000 nucleotides. This is in contrast to "normal-composition" RNAs (i.e., molecules with comparable numbers of each of the four nucleotides and hence capable of developing a large amount of secondary structure because of intramolecular complementarity/basepairing); a curvature corresponding to T = 3-size (∼28 nm in diameter) is preferred for the VLPs formed with such RNAs. Our work is consistent with the preferred curvature of VLPs being a consequence of interaction of CP with RNA-in particular, the presence or absence of short RNA duplexes-and suggests that the equilibrium size of the capsid results from a trade-off between this optimum size and the cost of confinement.

Real-Time Assembly of Viruslike Nucleocapsids Elucidated at the Single-Particle Level

While the structure of a multitude of viral particles has been resolved to atomistic detail, their assembly pathways remain largely elusive. Key unresolved issues are particle nucleation, particle growth, and the mode of genome compaction. These issues are difficult to address in bulk approaches and are effectively only accessible by the real-time tracking of assembly dynamics of individual particles. This we do here by studying the assembly into rod-shaped viruslike particles (VLPs) of artificial capsid polypeptides. Using fluorescence optical tweezers, we establish that small oligomers perform one-dimensional diffusion along the DNA. Larger oligomers are immobile and nucleate VLP growth. A multiplexed acoustic force spectroscopy approach reveals that DNA is compacted in regular steps, suggesting packaging via helical wrapping into a nucleocapsid. By reporting how real-time assembly tracking elucidates viral nucleation and growth principles, our work opens the door to a fundamental understanding of the complex assembly pathways of both VLPs and naturally evolved viruses.

Minimal coarse-grained models for molecular self-organisation in biology

The molecular machinery of life is largely created via self-organisation of individual molecules into functional assemblies. Minimal coarse-grained models, in which a whole macromolecule is represented by a small number of particles, can be of great value in identifying the main driving forces behind self-organisation in cell biology. Such models can incorporate data from both molecular and continuum scales, and their results can be directly compared to experiments. Here we review the state of the art of models for studying the formation and biological function of macromolecular assemblies in living organisms. We outline the key ingredients of each model and their main findings. We illustrate the contribution of this class of simulations to identifying the physical mechanisms behind life and diseases, and discuss their future developments.

Packaging of DNA origami in viral capsids

Here we show the encapsulation of 35 nm diameter, nearly-spherical, DNA origami by self-assembly of SV40-like (simian virus 40) particles. The self-assembly of this new type of nanoparticles is highly reproducible and efficient. The structure of these particles was determined by cryo-EM. The capsid forms a regular SV40 lattice of T = 7d icosahedral symmetry and the structural features of encapsulated DNA origami are fully visible. These particles are a promising biomaterial for use in various medical applications.

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

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

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 equilibrium structure of self-assembled protein nano-cages

Understanding how highly symmetric, robust, monodisperse protein nano-cages self-assemble can have major applications in various areas of bio-nanotechnology, such as drug delivery, biomedical imaging and gene therapy. We develop a model to investigate the assembly of protein subunits into the structures with different sizes and symmetries. Using Monte Carlo simulation, we obtain global minimum energy structures. Our results suggest that the physical properties including the spontaneous curvature, flexibility and bending rigidity of coat proteins are sufficient to predict the size, symmetry and shape selectivity of the assembly products. Further, on a thermodynamic basis, we discuss the polymorphism of nano-cages observed in assembly experiments.

Genome packaging in multi-segmented dsRNA viruses: distinct mechanisms with similar outcomes

Segmented double-stranded (ds)RNA viruses share remarkable similarities in their replication strategy and capsid structure. During virus replication, positive-sense single-stranded (+)RNAs are packaged into procapsids, where they serve as templates for dsRNA synthesis, forming progeny particles containing a complete equimolar set of genome segments. How the +RNAs are recognized and stoichiometrically packaged remains uncertain. Whereas bacteriophages of the Cystoviridae family rely on specific RNA-protein interactions to select appropriate +RNAs for packaging, viruses of the Reoviridae instead rely on specific inter-molecular interactions between +RNAs that guide multi-segmented genome assembly. While these families use distinct mechanisms to direct +RNA packaging, both yield progeny particles with a complete set of genomic dsRNAs.

Why large icosahedral viruses need scaffolding proteins

While small single-stranded viral shells encapsidate their genome spontaneously, many large viruses, such as the herpes simplex virus or infectious bursal disease virus (IBDV), typically require a template, consisting of either scaffolding proteins or an inner core. Despite the proliferation of large viruses in nature, the mechanisms by which hundreds or thousands of proteins assemble to form structures with icosahedral order (IO) is completely unknown. Using continuum elasticity theory, we study the growth of large viral shells (capsids) and show that a nonspecific template not only selects the radius of the capsid, but also leads to the error-free assembly of protein subunits into capsids with universal IO. We prove that as a spherical cap grows, there is a deep potential well at the locations of disclinations that later in the assembly process will become the vertices of an icosahedron. Furthermore, we introduce a minimal model and simulate the assembly of a viral shell around a template under nonequilibrium conditions and find a perfect match between the results of continuum elasticity theory and the numerical simulations. Besides explaining available experimental results, we provide a number of predictions. Implications for other problems in spherical crystals are also discussed.

Experimental and Theoretical Determination of the pH inside the Confinement of a Virus-Like Particle

In biology, a variety of highly ordered nanometer-size protein cages is found. Such structures find increasing application in, for example, vaccination, drug delivery, and catalysis. Understanding the physiochemical properties, particularly inside the confinement of a protein cage, helps to predict the behavior and properties of new materials based on such particles. Here, the relation between the bulk solution pH and the local pH inside a model protein cage, based on virus-like particles (VLPs) built from the coat proteins of the cowpea chlorotic mottle virus, is investigated. The pH is a crucial parameter in a variety of processes and is potentially significantly influenced by the high concentration of charges residing on the interior of the VLPs. The data show a systematic more acidic pH of 0.5 unit inside the VLP compared to that of the bulk solution for pH values above pH 6, which is explained using a theoretical model based on a Donnan equilibrium. The model agrees with the experimental data over almost two orders of magnitude, while below pH 6 the experimental data point to a buffering capacity of the VLP. These results are a first step in a better understanding of the physiochemical conditions inside a protein cage.

The role of the encapsulated cargo in microcompartment assembly

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