The impact of viral RNA on the association rates of capsid protein assembly: bacteriophage MS2 as a case study

Stability and conformational memory of electrosprayed and rehydrated bacteriophage MS2 virus coat proteins

Proteins are innately dynamic, which is important for their functions, but which also poses significant challenges when studying their structures. Gas-phase techniques can utilise separation and a range of sample manipulations to transcend some of the limitations of conventional techniques for structural biology in crystalline or solution phase, and isolate different states for separate interrogation. However, the transfer from solution to the gas phase risks affecting the structures, and it is unclear to what extent different conformations remain distinct in the gas phase, and if resolution in silico can recover the native conformations and their differences. Here, we use extensive molecular dynamics simulations to study the two distinct conformations of dimeric capsid protein of the MS2 bacteriophage. The protein undergoes notable restructuring of its peripheral parts in the gas phase, but subsequent simulation in solvent largely recovers the native structure. Our results suggest that despite some structural loss due to the experimental conditions, gas-phase structural biology techniques provide meaningful data that inform not only about the structures but also conformational dynamics of proteins.

•

Presented extensive molecular dynamics (MD) simulation data investigating protein vacuum exposure and rehydration dynamics.

•

Demonstrated that the majority of the protein structure recovers their initial solution conformation after vacuum exposure.

•

Explored the potential gain for structural biology of using MD simulation to refine gas-phase determined protein structures.

The diversity of protein-protein interaction interfaces within T=3 icosahedral viral capsids

Some non-enveloped virus capsids assemble from multiple copies of a single type of coat-protein (CP). The comparative energetics of the diverse CP-CP interfaces present in such capsids likely govern virus assembly-disassembly mechanisms. The T = 3 icosahedral capsids comprise 180 CP copies arranged about two-, three-, five- and six-fold axes of (quasi-)rotation symmetry. Structurally diverse CPs can assemble into T = 3 capsids. Specifically, the Leviviridae CPs are structurally distinct from the Bromoviridae, Tombusviridae and Tymoviridae CPs which fold into the classic “jelly-roll” fold. However, capsids from across the four families are known to disassemble into dimers. To understand whether the overall symmetry of the capsid or the structural details of the CP determine virus assembly-disassembly mechanisms, we analyze the different CP-CP interfaces that occur in the four virus families. Previous work studied protein homodimer interfaces using interface size (relative to the monomer) and hydrophobicity. Here, we analyze all CP-CP interfaces using these two parameters and find that the dimerization interface (present between two CPs congruent through a two-fold axis of rotation) has a larger relative size in the Leviviridae than in the other viruses. The relative sizes of the other Leviviridae interfaces and all the jelly-roll interfaces are similar. However, the dimerization interfaces across families have slightly higher hydrophobicity, potentially making them stronger than other interfaces. Finally, although the CP-monomers of the jelly-roll viruses are structurally similar, differences in their dimerization interfaces leads to varied dimer flexibility. Overall, differences in CP-structures may induce different modes of swelling and assembly-disassembly in the T = 3 viruses.

A Multiscale Model for the Self-Assembly of Coat Proteins in Bacteriophage MS2

The self-assembly of viral capsids is an essential step to the formation of infectious viruses. Elucidating the kinetic mechanisms of how a capsid or virus-like particle assembles could advance our knowledge about the viral lifecycle, as well as the general principles in self-assembly of biomaterials. However, current understanding of capsid assembly remains incomplete for many viruses due to the fact that the transient intermediates along the assembling pathways are experimentally difficult to be detected. In this paper, we constructed a new multiscale computational framework to simulate the self-assembly of virus-like particles. We applied our method to the coat proteins of bacteriophage MS2 as a specific model system. This virus-like particle of bacteriophage MS2 has a unique feature that its 90 sequence-identical dimers can be classified into two structurally various groups: one is the symmetric CC dimer, and the other is the asymmetric AB dimer. The homotypic interactions between AB dimers result in a 5-fold symmetric contact, while the heterotypic interactions between AB and CC dimers result in 6-fold symmetric contact. We found that the assembly can be described as a physical process of phase transition that is regulated by various factors such as concentration and specific stoichiometry between AB and CC dimers. Our simulations also demonstrate that heterotypic and homotypic interfaces play distinctive roles in modulating the assembling kinetics. The interaction between AB and CC dimers is much more dynamic than that between two AB dimers. We therefore suggest that the alternate growth of viral capsid through the heterotypic dimer interactions dominates the assembling pathways. This is, to the best of our knowledge, the first multiscale model to simulate the assembling process of coat proteins in bacteriophage MS2. The generality of this approach opens the door to its further applications in assembly of other viral capsids, virus-like particles, and novel drug delivery systems.

Encapsulation of Negatively Charged Cargo in MS2 Viral Capsids

Encapsulation into virus-like particles is an efficient way of loading cargo of interest for delivery applications. Here, we describe the encapsulation of proteins with tags comprising anionic amino acids or DNA and gold nanoparticles with negative surface charges inside MS2 bacteriophage capsids to obtain homogeneous nanoparticles with a diameter of 27 nm.

The conformational feasibility for the formation of reaching dimer in ASV and HIV integrase: a molecular dynamics study

Retroviral integrases are reported to form alternate dimer assemblies like the core-core dimer and reaching dimer. The core-core dimer is stabilized predominantly by an extensive interface between two catalytic core domains. The reaching dimer is stabilized by N-terminal domains that reach to form intermolecular interfaces with the other subunit's core and C-terminal domains (CTD), as well as CTD-CTD interactions. In this study, molecular dynamics (MD), Brownian dynamics (BD) simulations, and free energy analyses, were performed to elucidate determinants for the stability of the reaching dimer forms of full-length Avian Sarcoma Virus (ASV) and Human Immunodeficiency Virus (HIV) IN, and to examine the role of the C-tails (the last ~16-18 residues at the C-termini) in their structural dynamics. The dynamics of an HIV reaching dimer derived from small angle X-ray scattering and protein crosslinking data, was compared with the dynamics of a core-core dimer model derived from combining the crystal structures of two-domain fragments. The results showed that the core domains in the ASV reaching dimer express free dynamics, whereas those in the HIV reaching dimer are highly stable. BD simulations suggest a higher rate of association for the HIV core-core dimer than the reaching dimer. The predicted stability of these dimers was therefore ranked in the following order: ASV reaching dimer < HIV reaching dimer < composite core-core dimer. Analyses of MD trajectories have suggested residues that are critical for intermolecular contacts in each reaching dimer. Tests of these predictions and insights gained from these analyses could reveal a potential pathway for the association and dissociation of full-length IN multimers.

The allosteric switching mechanism in bacteriophage MS2

We use all-atom simulations to elucidate the mechanisms underlying conformational switching and allostery within the coat protein of the bacteriophage MS2. Assembly of most icosahedral virus capsids requires that the capsid protein adopts different conformations at precise locations within the capsid. It has been shown that a 19 nucleotide stem loop (TR) from the MS2 genome acts as an allosteric effector, guiding conformational switching of the coat protein during capsid assembly. Since the principal conformational changes occur far from the TR binding site, it is important to understand the molecular mechanism underlying this allosteric communication. To this end, we use all-atom simulations with explicit water combined with a path sampling technique to sample the MS2 coat protein conformational transition, in the presence and absence of TR-binding. The calculations find that TR binding strongly alters the transition free energy profile, leading to a switch in the favored conformation. We discuss changes in molecular interactions responsible for this shift. We then identify networks of amino acids with correlated motions to reveal the mechanism by which effects of TR binding span the protein. We find that TR binding strongly affects residues located at the 5-fold and quasi-sixfold interfaces in the assembled capsid, suggesting a mechanism by which the TR binding could direct formation of the native capsid geometry. The analysis predicts amino acids whose substitution by mutagenesis could alter populations of the conformational substates or their transition rates.

SDA 7: A modular and parallel implementation of the simulation of diffusional association software

The simulation of diffusional association (SDA) Brownian dynamics software package has been widely used in the study of biomacromolecular association. Initially developed to calculate bimolecular protein-protein association rate constants, it has since been extended to study electron transfer rates, to predict the structures of biomacromolecular complexes, to investigate the adsorption of proteins to inorganic surfaces, and to simulate the dynamics of large systems containing many biomacromolecular solutes, allowing the study of concentration-dependent effects. These extensions have led to a number of divergent versions of the software. In this article, we report the development of the latest version of the software (SDA 7). This release was developed to consolidate the existing codes into a single framework, while improving the parallelization of the code to better exploit modern multicore shared memory computer architectures. It is built using a modular object-oriented programming scheme, to allow for easy maintenance and extension of the software, and includes new features, such as adding flexible solute representations. We discuss a number of application examples, which describe some of the methods available in the release, and provide benchmarking data to demonstrate the parallel performance.

webSDA: a web server to simulate macromolecular diffusional association

Macromolecular interactions play a crucial role in biological systems. Simulation of diffusional association (SDA) is a software for carrying out Brownian dynamics simulations that can be used to study the interactions between two or more biological macromolecules. webSDA allows users to run Brownian dynamics simulations with SDA to study bimolecular association and encounter complex formation, to compute association rate constants, and to investigate macromolecular crowding using atomically detailed macromolecular structures. webSDA facilitates and automates the use of the SDA software, and offers user-friendly visualization of results. webSDA currently has three modules: 'SDA docking' to generate structures of the diffusional encounter complexes of two macromolecules, 'SDA association' to calculate bimolecular diffusional association rate constants, and 'SDA multiple molecules' to simulate the diffusive motion of hundreds of macromolecules. webSDA is freely available to all users and there is no login requirement. webSDA is available at http://mcm.h-its.org/webSDA/.

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.

Applying molecular crowding models to simulations of virus capsid assembly in vitro

Virus capsid assembly has been widely studied as a biophysical system, both for its biological and medical significance and as an important model for complex self-assembly processes. No current technology can monitor assembly in detail and what information we have on assembly kinetics comes exclusively from in vitro studies. There are many differences between the intracellular environment and that of an in vitro assembly assay, however, that might be expected to alter assembly pathways. Here, we explore one specific feature characteristic of the intracellular environment and known to have large effects on macromolecular assembly processes: molecular crowding. We combine prior particle simulation methods for estimating crowding effects with coarse-grained stochastic models of capsid assembly, using the crowding models to adjust kinetics of capsid simulations to examine possible effects of crowding on assembly pathways. Simulations suggest a striking difference depending on whether or not a system uses nucleation-limited assembly, with crowding tending to promote off-pathway growth in a nonnucleation-limited model but often enhancing assembly efficiency at high crowding levels even while impeding it at lower crowding levels in a nucleation-limited model. These models may help us understand how complicated assembly systems may have evolved to function with high efficiency and fidelity in the densely crowded environment of the cell.

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.

A new paradigm for the roles of the genome in ssRNA viruses

Recent work with RNA phages and an ssRNA plant satellite virus challenges the widely held view that the sequences and structures of genomic RNAs are unimportant for virion assembly. In the T=3 phages, RNA–coat protein interactions occur throughout the genome, defining the quasiconformers of their protein shells. In the plant virus, there are multiple packaging signals dispersed throughout the genome that overcome electrostatic barriers to protein self-assembly. Both viral coat proteins cause the solution structures of their cognate genomes to collapse into a form that is readily encapsidated in a two-stage assembly process. Such similar behavior in two structurally unrelated viral protein folds implies that this might be a conserved feature of many viral assembly reactions. These results suggest a highly defined structure for the RNA in the virions, consistent with recent structural studies. They also have implications both for subsequent genome release during infection and for the evolution of viral sequences.

A two-stage mechanism of viral RNA compaction revealed by single molecule fluorescence

Long RNAs often exist as multiple conformers in equilibrium. For the genomes of single-stranded RNA viruses, one of these conformers must include a compacted state allowing the RNA to be confined within the virion. We have used single molecule fluorescence correlation spectroscopy to monitor the conformations of viral genomes and sub-fragments in the absence and presence of coat proteins. Cognate RNA-coat protein interactions in two model viruses cause a rapid collapse in the hydrodynamic radii of their respective RNAs. This is caused by protein binding at multiple sites on the RNA that facilitate additional protein-protein contacts. The collapsed species recruit further coat proteins to complete capsid assembly with great efficiency and fidelity. The specificity in RNA-coat protein interactions seen at single-molecule concentrations reflects the packaging selectivity seen for such viruses in vivo. This contrasts with many in vitro reassembly measurements performed at much higher concentrations. RNA compaction by coat protein or polycation binding are distinct processes, implying that defined RNA-coat protein contacts are required for assembly.

Osmolyte-mediated encapsulation of proteins inside MS2 viral capsids

The encapsulation of enzymes in nanometer-sized compartments has the potential to enhance and control enzymatic activity, both in vivo and in vitro. Despite this potential, there is little quantitative data on the effect of encapsulation in a well-defined compartment under varying conditions. To gain more insight into these effects, we have characterized two improved methods for the encapsulation of heterologous molecules inside bacteriophage MS2 viral capsids. First, attaching DNA oligomers to a molecule of interest and incubating it with MS2 coat protein dimers yielded reassembled capsids that packaged the tagged molecules. The addition of a protein-stabilizing osmolyte, trimethylamine-N-oxide, significantly increased the yields of reassembly. Second, we found that expressed proteins with genetically encoded negatively charged peptide tags could also induce capsid reassembly, resulting in high yields of reassembled capsids containing the protein. This second method was used to encapsulate alkaline phosphatase tagged with a 16 amino acid peptide. The purified encapsulated enzyme was found to have the same K(m) value and a slightly lower k(cat) value than the free enzyme, indicating that this method of encapsulation had a minimal effect on enzyme kinetics. This method provides a practical and potentially scalable way of studying the complex effects of encapsulating enzymes in protein-based compartments.