Primary changes of the mechanical properties of Southern Bean Mosaic Virus upon calcium removal

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.

Heat-induced transitions of an empty minute virus of mice capsid in explicit water: all-atom MD simulation

The capsid-like structure of the virus-based protein nanoparticles (NPs) can serve as bionanomaterials, with applications in biomedicines and nanotechnology. Release of packaged material from these nanocontainers is associated with subtle conformational changes of the NP structure, which in vitro, is readily accomplished by heating. Characterizing the structural changes as a function of temperature may provide fresh insights into nanomaterial/antiviral strategies. Here, we have calculated heat induced changes in the properties of an empty minute virus of mice particle using large-scale ≈ 3.0 × 10⁶ all-atom molecular dynamics simulations. We focus on two heat induced structural changes of the NP, namely, dynamical transition (DT) and breathing transition (BT), both characterized by sudden and sharp change of measured parameters at temperatures, T_DT and T_BT, respectively. While DT is assessed by mean-square fluctuation of hydrogen atoms of the NP, BT is monitored through internal volume and permeation rate of water molecules through the NP. Both the transitions, resulting primarily from collective atomistic motion, are found to occur at temperatures widely separated from one another (T_BT>T_DT). The breathing motions, responsible for the translocation events of the packaged materials through the NP to kick off, are further probed by computing atomic resolution stresses from NVE simulations. Distribution of equilibrium atomistic stresses on the NP reveals a largely asymmetric nature and suggests structural breathing may actually represent large dynamic changes in the hotspot regions, far from the NP pores, which is in remarkable resemblance with recently conducted hydrogen-deuterium exchange coupled to mass spectrometry experiment. Communicated by Ramaswamy H. Sarma.

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.

Minimal Cylinder Analysis Reveals the Mechanical Properties of Oncogenic Nucleosomes

Histone variants regulate replication, transcription, DNA damage repair, and chromosome segregation. Though widely accepted as a paradigm, it has not been rigorously demonstrated that histone variants encode unique mechanical properties. Here, we present a new theoretical approach called minimal cylinder analysis that uses strain fluctuations to determine the Young's modulus of nucleosomes from all-atom molecular dynamics simulations. Recently, we validated this computational tool against in vitro single-molecule nanoindentation of histone variant nucleosomes. In this report, we further extend minimal cylinder analysis to study the biophysical properties of hybrid nucleosomes that are known to exist in human cancer cells and contain H3 histone variants CENP-A and H3.3. Here, we report that the heterotypic nucleosome has an intermediate elasticity (8.5 ± 0.5 MPa) compared to CENP-A (6.2 ± 0.4 MPa) and H3 (9.8 ± 0.7 MPa) and that the dynamics of both canonical and CENP-A nucleosomes are preserved and partitioned across the nucleosome pseudodyad. Furthermore, we investigate the mechanism by which the elasticity of these heterotypic nucleosomes augments cryptic binding surfaces. From these analyses, we predict that the heterotypic nucleosome is permissive to the binding of one copy of the kinetochore protein CENP-C while still retaining a closed DNA end configuration required for linker histone H1 to bind. We discuss that the ectopic deposition of CENP-A in cancer by H3.3 chaperones HIRA and DAXX may fortuitously result in hybrid nucleosome formation. Using these results, we propose biological outcomes that might arise when such heterotypic nucleosomes occupy large regions of the genome.

All-Atom Molecular Dynamics Simulations of Whole Viruses

Classical molecular dynamics modeling of whole viruses or their capsids in explicit water is discussed, and known examples from the literature are analyzed. Only works on all-atom modeling in explicit water are included. Physical chemistry of the whole system is the focus, which includes the structure and dynamics of the biomolecules as well as water and ion behavior in and around the virus particle. It was demonstrated that in most investigations molecular phenomena that currently can not be studied experimentally are successfully reproduced and explained by the simulations. These include, for example, transport and distribution of ions inside viruses that ultimately connected to their stability, the hydrodynamic pressure in the capsid related to viruses' elastic properties, the role of metal ions in virus swelling, and others. Current and future tendencies in the development of all-atom virus simulations are outlined.

All-Atom Molecular Dynamics Simulations of Entire Virus Capsid Reveal the Role of Ion Distribution in Capsid's Stability

Present experimental methods do not have sufficient resolution to investigate all processes in virus particles at atomistic details. We report the results of molecular dynamics simulations and analyze the connection between the number of ions inside an empty capsid of PCV2 virus and its stability. We compare the crystallographic structures of the capsids with unresolved N-termini and without them in realistic conditions (room temperature and aqueous solution) and show that the structure is preserved. We find that the chloride ions play a key role in the stability of the capsid. A low number of chloride ions results in loss of the native icosahedral symmetry, while an optimal number of chloride ions create a neutralizing layer next to the positively charged inner surface of the capsid. Understanding the dependence of the capsid stability on the distribution of the ions will help clarify the details of the viral life cycle that is ultimately connected to the role of packaged viral genome inside the capsid.

Stepwise reversible nanomechanical buckling in a viral capsid

Viruses are nanoscale infectious agents constructed of a proteinaceous capsid that protects the packaged genomic material. Nanoindentation experiments using atomic force microscopy have, in recent years, provided unprecedented insight into the elastic properties, structural stability and maturation-dependent mechanical changes in viruses. However, the dynamics of capsid behavior are still unresolved. Here we used high-resolution nanoindentation experiments on mature, DNA-filled T7 bacteriophage particles. The elastic regime of the nanoindentation force trace contained discrete, stepwise transitions that cause buckling of the T7 capsid with magnitudes that are integer multiples of ∼0.6 nm. Remarkably, the transitions are reversible and contribute to the rapid consolidation of the capsid structure against a force during cantilever retraction. The stepwise transitions were present even following the removal of the genomic DNA by heat treatment, indicating that they are related to the structure and dynamics of the capsomeric proteins. Dynamic force spectroscopy experiments revealed that the thermally activated consolidation step is ∼10⁴ times faster than spontaneous buckling, suggesting that the capsid stability is under strong dynamic control. Capsid structural dynamics may play an important role in protecting the genomic material from harsh environmental impacts. The nanomechanics approach employed here may be used to investigate the structural dynamics of other viruses and nanoscale containers as well.

Computational virology: From the inside out

Viruses typically pack their genetic material within a protein capsid. Enveloped viruses also have an outer membrane made up of a lipid bilayer and membrane-spanning glycoproteins. X-ray diffraction and cryoelectron microscopy provide high resolution static views of viral structure. Molecular dynamics (MD) simulations may be used to provide dynamic insights into the structures of viruses and their components. There have been a number of simulations of viral capsids and (in some cases) of the inner core of RNA or DNA packaged within them. These simulations have generally focussed on the structural integrity and stability of the capsid and/or on the influence of the nucleic acid core on capsid stability. More recently there have been a number of simulation studies of enveloped viruses, including HIV-1, influenza A, and dengue virus. These have addressed the dynamic behaviour of the capsid, the matrix, and/or of the outer envelope. Analysis of the dynamics of the lipid bilayer components of the envelopes of influenza A and of dengue virus reveals a degree of biophysical robustness, which may contribute to the stability of virus particles in different environments. Significant computational challenges need to be addressed to aid simulation of complex viruses and their membranes, including the need to integrate structural data from a range of sources to enable us to move towards simulations of intact virions. This article is part of a Special Issue entitled: Membrane Proteins edited by J.C. Gumbart and Sergei Noskov.

Swelling and softening of the cowpea chlorotic mottle virus in response to pH shifts

Cowpea chlorotic mottle virus (CCMV) forms highly elastic icosahedral protein capsids that undergo a characteristic swelling transition when the pH is raised from 5 to 7. Here, we performed nano-indentation experiments using an atomic force microscope to track capsid swelling and measure the shells' Young's modulus at the same time. When we chelated Ca(2+) ions and raised the pH, we observed a gradual swelling of the RNA-filled capsids accompanied by a softening of the shell. Control experiments with empty wild-type virus and a salt-stable mutant revealed that the softening was not strictly coupled to the swelling of the protein shells. Our data suggest that a pH increase and Ca(2+) chelation lead primarily to a loosening of contacts within the protein shell, resulting in a softening of the capsid. This appears to render the shell metastable and make swelling possible when repulsive forces among the capsid proteins become large enough, which is known to be followed by capsid disassembly at even higher pH. Thus, softening and swelling are likely to play a role during inoculation.

Fluctuating Nonlinear Spring Model of Mechanical Deformation of Biological Particles

The mechanical properties of virus capsids correlate with local conformational dynamics in the capsid structure. They also reflect the required stability needed to withstand high internal pressures generated upon genome loading and contribute to the success of important events in viral infectivity, such as capsid maturation, genome uncoating and receptor binding. The mechanical properties of biological nanoparticles are often determined from monitoring their dynamic deformations in Atomic Force Microscopy nanoindentation experiments; but a comprehensive theory describing the full range of observed deformation behaviors has not previously been described. We present a new theory for modeling dynamic deformations of biological nanoparticles, which considers the non-linear Hertzian deformation, resulting from an indenter-particle physical contact, and the bending of curved elements (beams) modeling the particle structure. The beams’ deformation beyond the critical point triggers a dynamic transition of the particle to the collapsed state. This extreme event is accompanied by a catastrophic force drop as observed in the experimental or simulated force (F)-deformation (X) spectra. The theory interprets fine features of the spectra, including the nonlinear components of the FX-curves, in terms of the Young’s moduli for Hertzian and bending deformations, and the structural damage dependent beams’ survival probability, in terms of the maximum strength and the cooperativity parameter. The theory is exemplified by successfully describing the deformation dynamics of natural nanoparticles through comparing theoretical curves with experimental force-deformation spectra for several virus particles. This approach provides a comprehensive description of the dynamic structural transitions in biological and artificial nanoparticles, which is essential for their optimal use in nanotechnology and nanomedicine applications.

Dynamic force experiments, which have become available to explore the physical properties of biological assemblies, oftentimes reveal results that are difficult to understand without theoretical framework. We employed a multiscale modeling approach—a combination of Molecular Dynamics simulations of atomic structures with Langevin simulations of coarse-grained models of virus shells—to characterize the degrees of freedom defining the deformation and structural collapse of biological particles tested mechanically. This enabled us to develop an analytical model that provides meaningful interpretation of force-deformation spectra available from single-particle nanoindentation experiments. The Fluctuating Nonlinear Spring (FNS) model of uniaxial particle’s deformation captures essential features of the force-deformation spectra as observed in nanomanipulations in vitro and in silico: initial non-linearity, then a subsequent force decrease transition due to structural collapse. Our theory uniquely combines the elements of continuum mechanics with the statistics of extremes, enabling one to gather mechanical and statistical characteristics of nanoparticles, which determine the Hertzian deformation of the particle’s protein layer, and bending deformation and structural damage to the particle structure. We have demonstrated how the FNS theory can accurately model the deformation of several viral shells, showing promising model applications for describing a variety of natural and synthetic nanoparticles.

A hybrid molecular dynamics/fluctuating hydrodynamics method for modelling liquids at multiple scales in space and time

A new 3D implementation of a hybrid model based on the analogy with two-phase hydrodynamics has been developed for the simulation of liquids at microscale. The idea of the method is to smoothly combine the atomistic description in the molecular dynamics zone with the Landau-Lifshitz fluctuating hydrodynamics representation in the rest of the system in the framework of macroscopic conservation laws through the use of a single "zoom-in" user-defined function s that has the meaning of a partial concentration in the two-phase analogy model. In comparison with our previous works, the implementation has been extended to full 3D simulations for a range of atomistic models in GROMACS from argon to water in equilibrium conditions with a constant or a spatially variable function s. Preliminary results of simulating the diffusion of a small peptide in water are also reported.

Overview on Sobemoviruses and a Proposal for the Creation of the Family Sobemoviridae

The genus Sobemovirus, unassigned to any family, consists of viruses with single-stranded plus-oriented single-component RNA genomes and small icosahedral particles. Currently, 14 species within the genus have been recognized by the International Committee on Taxonomy of Viruses (ICTV) but several new species are to be recognized in the near future. Sobemovirus genomes are compact with a conserved structure of open reading frames and with short untranslated regions. Several sobemoviruses are important pathogens. Moreover, over the last decade sobemoviruses have become important model systems to study plant virus evolution. In the current review we give an overview of the structure and expression of sobemovirus genomes, processing and functions of individual proteins, particle structure, pathology and phylogenesis of sobemoviruses as well as of satellite RNAs present together with these viruses. Based on a phylogenetic analysis we propose that a new family Sobemoviridae should be recognized including the genera Sobemovirus and Polemovirus. Finally, we outline the future perspectives and needs for the research focusing on sobemoviruses.

Nothing to sneeze at: a dynamic and integrative computational model of an influenza A virion

The influenza virus is surrounded by an envelope composed of a lipid bilayer and integral membrane proteins. Understanding the structural dynamics of the membrane envelope provides biophysical insights into aspects of viral function, such as the wide-ranging survival times of the virion in different environments. We have combined experimental data from X-ray crystallography, nuclear magnetic resonance spectroscopy, cryo-electron microscopy, and lipidomics to build a model of the intact influenza A virion. This is the basis of microsecond-scale coarse-grained molecular dynamics simulations of the virion, providing simulations at different temperatures and with varying lipid compositions. The presence of the Forssman glycolipid alters a number of biophysical properties of the virion, resulting in reduced mobility of bilayer lipid and protein species. Reduced mobility in the virion membrane may confer physical robustness to changes in environmental conditions. Our simulations indicate that viral spike proteins do not aggregate and thus are competent for multivalent immunoglobulin G interactions.

Elastic properties and heterogeneous stiffness of the phi29 motor connector channel

The DNA packaging motor of the bacteriophage ϕ29, comprising head-tail connector, ATPase, and pRNA, transports the viral DNA inside the procapsid against pressure differences of up to ∼60 atm during replication. Several models for the DNA packaging mechanism have been proposed, which attribute different roles to the connector, and require specific mechanical properties of the connector. To characterize these properties at the atomic level, and to understand how the connector withstands this large pressure, we have carried out molecular dynamics simulations of the whole connector both in equilibrium and under mechanical stress. The simulations revealed a quite heterogeneous distribution of stiff and soft regions, resembling that of typical composite materials that are also optimized to resist mechanical stress. In particular, the conserved middle α-helical region is found to be remarkably stiff, similar only to structural proteins forming viral shell, silk, or collagen. In contrast, large parts of the peripheral interface to the ϕ29 procapsid turned out to be rather soft. Force probe and umbrella sampling simulations showed that large connector deformations are remarkably reversible, and served to calculate the free energies required for these deformations. In particular, for an untwisting deformation by 12°, as postulated by the untwist-twist model, more than four times' larger energy is required than is available from hydrolysis of one ATP molecule. Combined with previous experiments, this result is incompatible with the untwist-twist model. In contrast, our simulations support the recently proposed one-way revolution model and suggest in structural terms how the connector blocks DNA leakage. In particular, conserved loops at the rim of the central channel, which are in direct contact with the DNA, are found to be rather flexible and tightly anchored to the rigid central region. These findings suggest a check-valve mechanism, with the flexible loops obstructing the channel by interacting with the viral DNA.

Universal relaxation governs the nonequilibrium elasticity of biomolecules

Experimental and computational dynamic force spectroscopy is widely used to determine the mechanical properties of single biomolecules. Whereas so far the focus has mainly been on rupture or unfolding forces, recent force-probe molecular dynamics simulations have revealed a strong loading rate dependence of biomolecular elasticities, which cannot be explained by the established one-dimensional transition-state treatments. We show that this nonequilibrium behavior can be explained by a theory that includes relaxation effects. For three structurally and mechanically quite diverse systems, a single relaxation mode suffices to quantitatively describe their loading-rate-dependent elastic behavior. Atomistic simulations of these systems revealed the microscopic nature of the respective relaxation modes. This result suggests a new type of "elasticity spectroscopy" experiment, which should render nonequilibrium properties of structured macromolecules accessible to single-molecule force spectroscopy.

Virus capsid dissolution studied by microsecond molecular dynamics simulations

Dissolution of many plant viruses is thought to start with swelling of the capsid caused by calcium removal following infection, but no high-resolution structures of swollen capsids exist. Here we have used microsecond all-atom molecular simulations to describe the dynamics of the capsid of satellite tobacco necrosis virus with and without the 92 structural calcium ions. The capsid expanded 2.5% upon removal of the calcium, in good agreement with experimental estimates. The water permeability of the native capsid was similar to that of a phospholipid membrane, but the permeability increased 10-fold after removing the calcium, predominantly between the 2-fold and 3-fold related subunits. The two calcium binding sites close to the icosahedral 3-fold symmetry axis were pivotal in the expansion and capsid-opening process, while the binding site on the 5-fold axis changed little structurally. These findings suggest that the dissociation of the capsid is initiated at the 3-fold axis.

Microscopic Symmetry Imposed by Rotational Symmetry Boundary Conditions in Molecular Dynamics Simulation

A large number of viral capsids, as well as other macromolecular assemblies, have icosahedral structure or structures with other rotational symmetries. This symmetry can be exploited during molecular dynamics (MD) to model in effect the full viral capsid using only a subset of primary atoms plus copies of image atoms generated from rotational symmetry boundary conditions (RSBC). A pure rotational symmetry operation results in both primary and image atoms at short range, and within nonbonded interaction distance of each other, so that nonbonded interactions can not be specified by the minimum image convention and explicit treatment of image atoms is required. As such an unavoidable consequence of RSBC is that the enumeration of nonbonded interactions in regions surrounding certain rotational axes must include both a primary atom and its copied image atom, thereby imposing microscopic symmetry for some forces. We examined the possibility of artifacts arising from this imposed microscopic symmetry of RSBC using two simulation systems: a water shell and human rhinovirus 14 (HRV14) capsid with explicit water. The primary unit was a pentamer of the icosahedron, which has the advantage of direct comparison of icosahedrally equivalent spatial regions, for example regions near a 2-fold symmetry axis with imposed symmetry and a 2-fold axis without imposed symmetry. Analysis of structural and dynamic properties of water molecules and protein atoms found similar behavior near symmetry axes with imposed symmetry and where the minimum image convention fails compared with that in other regions in the simulation system, even though an excluded volume effect was detected for water molecules near the axes with imposed symmetry. These results validate the use of RSBC for icosahedral viral capsids or other rotationally symmetric systems.