Structural and dynamic asymmetry in icosahedrally symmetric virus capsids

Modeling Viral Capsid Assembly: A Review of Computational Strategies and Applications

Viral capsid assembly is a complex and critical process, essential for understanding viral behavior, evolution, and the development of antiviral treatments, vaccines, and nanotechnology. Significant progress in studying viral capsid assembly has been achieved through various computational approaches, including molecular dynamics (MD) simulations, stochastic dynamics simulations, coarse-grained (CG) models, electrostatic analyses, lattice models, hybrid techniques, machine learning methods, and kinetic models. Each of these techniques offers unique advantages, and by integrating these diverse computational strategies, researchers can more accurately model the dynamic behaviors and structural features of viral capsids, deepening our understanding of the assembly process. This review provides a comprehensive overview of studies on viral capsid assembly, emphasizing their critical role in advancing our knowledge. It examines the contributions, strengths, and limitations of different computational methods, presents key computational works in the field, and analyzes milestone studies that have shaped current research.

Theoretical Studies on Assembly, Physical Stability, and Dynamics of Viruses

All matter must obey the general laws of physics and living matter is not an exception. Viruses have not only learnt how to cope with them but have managed to use them for their own survival. In this chapter, we will review some of the exciting physics that are behind viruses and discuss simple physical models that can shed some light on different aspects of the viral life cycle and viral properties. In particular, we will focus on how the structure and shape of the viral capsid, its assembly and stability, and the entry and exit of viral particles and their genomes can be explained using fundamental physics theories.

Non-equilibrium virus particle dynamics: Microsecond MD simulations of the complete Flock House virus capsid under different conditions

Flock House virus (FHV) is an animal virus and considered a model system for non-enveloped viruses. It has a small, icosahedral capsid (T=3) and a bipartite positive-sense RNA genome. We present an extensive study of the FHV capsid dynamics from all-atom molecular dynamics simulations of the complete capsid. The simulations explore different biologically relevant conditions (neutral/low pH, with/without RNA in the capsid) using the CHARMM force field. The results show that low pH destabilizes the capsid, causing radial expansion, and RNA stabilizes the capsid. The finding of low pH destabilization is biologically relevant because the capsid is exposed to low pH in the endosome, where conformational changes occur leading to genome release. We also observe structural changes at the fivefold and twofold symmetry axes that likely relate to the externalization of membrane active γ peptides through the fivefold vertex and extrusion of RNA at the twofold axis. Simulations using the Amber force field at neutral pH are also performed and display similar characteristics to the CHARMM simulations.

The balance between fitness advantages and costs drives adaptation of bacteriophage Qβ to changes in host density at different temperatures

Introduction

Host density is one of the main factors affecting the infective capacity of viruses. When host density is low, it is more difficult for the virus to find a susceptible cell, which increases its probability of being damaged by the physicochemical agents of the environment. Nevertheless, viruses can adapt to variations in host density through different strategies that depend on the particular characteristics of the life cycle of each virus. In a previous work, using the bacteriophage Qβ as an experimental model, we found that when bacterial density was lower than optimal the virus increased its capacity to penetrate into the bacteria through a mutation in the minor capsid protein (A1) that is not described to interact with the cell receptor.

Results

Here we show that the adaptive pathway followed by Qβ in the face of similar variations in host density depends on environmental temperature. When the value for this parameter is lower than optimal (30°C), the mutation selected is the same as at the optimal temperature (37°C). However, when temperature increases to 43°C, the mutation selected is located in a different protein (A2), which is involved both in the interaction with the cell receptor and in the process of viral progeny release. The new mutation increases the entry of the phage into the bacteria at the three temperatures assayed. However, it also considerably increases the latent period at 30 and 37°C, which is probably the reason why it is not selected at these temperatures.

Conclusion

The conclusion is that the adaptive strategies followed by bacteriophage Qβ, and probably other viruses, in the face of variations in host density depend not only on their advantages at this selective pressure, but also on the fitness costs that particular mutations may present in function of the rest of environmental parameters that influence viral replication and stability.

Long-time scale simulations of virus-like particles from three human-norovirus strains

The dynamics of the virus like particles (VLPs) corresponding to the GII.4 Houston, GII.2 SMV, and GI.1 Norwalk strains of human noroviruses (HuNoV) that cause gastroenteritis was investigated by means of long-time (about 30 μs in the laboratory timescale) molecular dynamics simulations with the coarse-grained UNRES force field. The main motion of VLP units turned out to be the bending at the junction between the P1 subdomain (that sits in the VLP shell) and the P2 subdomain (that protrudes outside) of the major VP1 protein, this resulting in a correlated wagging motion of the P2 subdomains with respect to the VLP surface. The fluctuations of the P2 subdomain were found to be more pronounced and the P2 domain made a greater angle with the normal to the VLP surface for the GII.2 strain, which could explain the inability of this strain to bind the histo-blood group antigens (HBGAs).

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.

Multiscale Modeling of Hepatitis B Virus Capsid Assembly and Its Dimorphism

Hepatitis B virus (HBV) is an endemic, chronic virus that leads to 800000 deaths per year. Central to the HBV lifecycle, the viral core has a protein capsid assembled from many copies of a single protein. The capsid protein adopts different (quasi-equivalent) conformations to form icosahedral capsids containing 180 or 240 proteins: T = 3 or T = 4, respectively, in Caspar-Klug nomenclature. HBV capsid assembly has become an important target for recently developed antivirals; nonetheless, the assembly pathways and mechanisms that control HBV dimorphism remain unclear. We describe computer simulations of the HBV assembly, using a coarse-grained model that has parameters learned from all-atom molecular dynamics simulations of a complete HBV capsid and yet is computationally tractable. Dynamical simulations with the resulting model reproduce experimental observations of HBV assembly pathways and products. By constructing Markov state models and employing transition path theory, we identify pathways leading to T = 3, T = 4, and other experimentally observed capsid morphologies. The analysis shows that capsid polymorphism is promoted by the low HBV capsid bending modulus, where the key factors controlling polymorphism are the conformational energy landscape and protein-protein binding affinities.

Fighting viruses with computers, right now

The synergistic conjunction of various technological revolutions with the accumulated knowledge and workflows is rapidly transforming several scientific fields. Particularly, Virology can now feed from accurate physical models, polished computational tools, and massive computational power to readily integrate high-resolution structures into biological representations of unprecedented detail. That preparedness allows for the first time to get crucial information for vaccine and drug design from in-silico experiments against emerging pathogens of worldwide concern at relevant action windows. The present work reviews some of the main milestones leading to these breakthroughs in Computational Virology, providing an outlook for future developments in capacity building and accessibility to computational resources.

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.