Virus capsid dissolution studied by microsecond molecular dynamics simulations

Structure and energetics guide dynamic behaviour in a T = 3 icosahedral virus capsid

Although virus capsids appear as rigid, symmetric particles in experimentally determined structures; biochemical studies suggest a significant degree of structural flexibility in the particles. We carried out all-atom simulations on the icosahedral capsid of an insect virus, Flock House Virus, which show intriguing differences in the degree of flexibility of quasi-equivalent capsid subunits consistent with previously described biological behaviour. The flexibility of all the β and γ subunits of the protein and RNA fragments is analysed and compared. Both γA subunit and RNA fragment exhibit higher flexibility than the γB and γC subunits. The capsid shell is permeable to the bidirectional movement of water molecules, and the movement is heavily influenced by the geometry of the capsid shell along specific symmetry axes. In comparison to the symmetry axes along I5 and I3, the I2 axis exhibits a slightly higher water content. This enriched water environment along I2 could play a pivotal role in facilitating the structural transitions necessary for RNA release, shedding some light on the intricate and dynamic processes underlying the viral life cycle. Our study suggests that the physical characterization of whole virus capsids is the key to identifying biologically relevant transition states in the virus life cycle and understanding the basis of virus infectivity.

Martini on the Rocks: Can a Coarse-Grained Force Field Model Crystals?

Computational chemistry is an important tool in numerous scientific disciplines, including drug discovery and structural biology. Coarse-grained models offer simple representations of molecular systems that enable simulations of large-scale systems. Because there has been an increase in the adoption of such models for simulations of biomolecular systems, critical evaluation is warranted. Here, the stability of the amyloid peptide and organic crystals is evaluated using the Martini 3 coarse-grained force field. The crystals change shape drastically during the simulations. Radial distribution functions show that the distance between backbone beads in β-sheets increases by ∼1 Å, breaking the crystals. The melting points of organic compounds are much too low in the Martini force field. This suggests that Martini 3 lacks the specific interactions needed to accurately simulate peptides or organic crystals without imposing artificial restraints. The problems may be exacerbated by the use of the 12-6 potential, suggesting that a softer potential could improve this model for crystal simulations.

Survival of surface bacteriophages and their hosts in in situ deep-sea environments

IMPORTANCE

The survival of the sinking prokaryotes and viruses in the deep-sea environment is crucial for deep-sea ecosystems and biogeochemical cycles. Through an in situ deep-sea long-term incubation device, our results showed that viral particles and infectivity had still not decayed completely after in situ incubation for 1 year. This suggests that, via infection and lysis, surface viruses with long-term infectious activity in situ deep-sea environments may influence deep-sea microbial populations in terms of activity, function, diversity, and community structure and ultimately affect deep-sea biogeochemical cycles, highlighting the need for additional research in this area.

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.

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.

Optimization of parallel implementation of UNRES package for coarse-grained simulations to treat large proteins

We report major algorithmic improvements of the UNRES package for physics-based coarse-grained simulations of proteins. These include (i) introduction of interaction lists to optimize computations, (ii) transforming the inertia matrix to a pentadiagonal form to reduce computing and memory requirements, (iii) removing explicit angles and dihedral angles from energy expressions and recoding the most time-consuming energy/force terms to minimize the number of operations and to improve numerical stability, (iv) using OpenMP to parallelize those sections of the code for which distributed-memory parallelization involves unfavorable computing/communication time ratio, and (v) careful memory management to minimize simultaneous access of distant memory sections. The new code enables us to run molecular dynamics simulations of protein systems with size exceeding 100,000 amino-acid residues, reaching over 1 ns/day (1 μs/day in all-atom timescale) with 24 cores for proteins of this size. Parallel performance of the code and comparison of its performance with that of AMBER, GROMACS and MARTINI 3 is presented.

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.

One microsecond MD simulations of the SARS-CoV-2 main protease and hydroxychloroquine complex reveal the intricate nature of binding

Currently, several vaccines and antivirals across the globe are in clinical trials. Hydroxychloroquine (HCQ) was reported to inhibit the SARS-CoV-2 virus in antiviral assays. Here, it raises the curiosity about the molecular target of HCQ inside the cell. It may inhibit some of the viral targets, or some other complex mechanisms must be at disposal towards action mechanisms. In some of the viruses, proteases are experimentally reported to be a potential target of HCQ. However, no in-depth investigations are available in the literature yet. Henceforth, we have carried out extensive, one-microsecond long molecular dynamics simulations of the bound complex of hydroxychloroquine with main protease (Mpro) of SARS-CoV-2. Our analysis found that HCQ binds within the catalytic pocket of Mpro and remains stable upto one-third of simulation time but further causes increased fluctuations in simulation parameters. In the end, the HCQ does not possess any pre-formed hydrogen bond, other non-covalent interactions with Mpro, ultimately showing the unsteadiness in binding at catalytic binding pocket and may suggest that HCQ may not inhibit the Mpro. In the future, this study would require experimental validation on enzyme assays against Mpro, and that may be the final say. Communicated by Ramaswamy H. Sarma.

Applications of Charge Detection Mass Spectrometry in Molecular Biology and Biotechnology

Charge detection mass spectrometry (CDMS) is a single-particle technique where the masses of individual ions are determined from simultaneous measurement of their mass-to-charge ratio (m/z) and charge. Masses are determined for thousands of individual ions, and then the results are binned to give a mass spectrum. Using this approach, accurate mass distributions can be measured for heterogeneous and high-molecular-weight samples that are usually not amenable to analysis by conventional mass spectrometry. Recent applications include heavily glycosylated proteins, protein complexes, protein aggregates such as amyloid fibers, infectious viruses, gene therapies, vaccines, and vesicles such as exosomes.

Analysis of Recombinant Adenovirus Vectors by Ion Trap Charge Detection Mass Spectrometry: Accurate Molecular Weight Measurements beyond 150 MDa

Adenovirus is one of the largest nonenveloped, double-stranded DNA viruses. It is widely used as a gene therapy vector and has recently received a lot of attention as a novel vaccine platform for SARS-CoV-2. Human adenovirus 5 (HAdV5) contains over 2500 protein molecules and has a 36 kbp genome. Adenovirus is well beyond the range of conventional mass spectrometry, and it was unclear how well such a large complex could be desolvated. Here, we report molecular weight (MW) distributions measured for HAdV5 and for 11 recombinant AdV vectors with genomes of varying lengths. The MW distributions were recorded using ion trap charge detection mass spectrometry (CDMS), a single-particle technique where m/z and charge are measured for individual ions. The results show that ions as large as 150 MDa can be effectively desolvated and accurate MW distributions obtained. The MW distribution for HAdV5 contains a narrow peak at 156.1 MDa, assigned to the infectious virus. A smaller peak at 129.6 MDa is attributed to incomplete particles that have not packaged a genome. The ions in the 129.6 MDa peak have a much lower average charge than those in the peak at 156.1 MDa. This is attributed to the empty particles missing some or all of the fibers that decorate the surface of the virion. The MW measured for the mature virus (156.1 MDa) is much larger than that predicted from sequence masses and copy numbers of the constituents (142.5 MDa). Measurements performed for recombinant AdV as a function of genome length show that for every 1 MDa increase in the genome MW, the MW of the mature virus increases by around 2.3 MDa. The additional 1.3 MDa is attributed to core proteins that are copackaged with the DNA. This observation suggests that the discrepancy between the measured and expected MWs for mature HAdV5 is due to an underestimate in the copy numbers of the core proteins.

All-atom molecular dynamics study of hepatitis B virus containing pregenome RNA in solution

Immature hepatitis B virus (HBV) captures nucleotides in its capsid for reverse transcription. The nucleotides and nucleotide analog drugs, which are triphosphorylated and negatively charged in the cell, approach the capsid via diffusion and are absorbed into it. In this study, we performed a long-time molecular dynamics calculation of the entire HBV capsid containing pregenome RNA to investigate the interactions between the capsid and negatively charged substances. Electric field analysis demonstrated that negatively charged substances can approach the HBV capsid by thermal motion, avoiding spikes. The substances then migrate all over the floor of the HBV capsid. Finally, they find pores through which they can pass through the HBV capsid shell. Free energy profiles were calculated along these pores for small ions to understand their permeability through the pores. Anions (Cl^-) showed higher free energy barriers than cations (Na⁺ and K⁺) through all pores, and the permeation rate of Cl^- was eight times slower than that of K⁺ or Na⁺. Furthermore, the ions were more stable in the capsid than in the bulk water. Thus, the HBV capsid exerts ion selectivity for uptake and provides an environment for ions, such as nucleotides and nucleotide analog drugs, to be stabilized within the capsid.

Theory and Practice of Coarse-Grained Molecular Dynamics of Biologically Important Systems

Molecular dynamics with coarse-grained models is nowadays extensively used to simulate biomolecular systems at large time and size scales, compared to those accessible to all-atom molecular dynamics. In this review article, we describe the physical basis of coarse-grained molecular dynamics, the coarse-grained force fields, the equations of motion and the respective numerical integration algorithms, and selected practical applications of coarse-grained molecular dynamics. We demonstrate that the motion of coarse-grained sites is governed by the potential of mean force and the friction and stochastic forces, resulting from integrating out the secondary degrees of freedom. Consequently, Langevin dynamics is a natural means of describing the motion of a system at the coarse-grained level and the potential of mean force is the physical basis of the coarse-grained force fields. Moreover, the choice of coarse-grained variables and the fact that coarse-grained sites often do not have spherical symmetry implies a non-diagonal inertia tensor. We describe selected coarse-grained models used in molecular dynamics simulations, including the most popular MARTINI model developed by Marrink's group and the UNICORN model of biological macromolecules developed in our laboratory. We conclude by discussing examples of the application of coarse-grained molecular dynamics to study biologically important processes.

Molecular mechanism of capsid disassembly in hepatitis B virus

The disassembly of a viral capsid leading to the release of its genetic material into the host cell is a fundamental step in viral infection. In hepatitis B virus (HBV), the capsid consists of identical protein monomers that dimerize and then arrange themselves into pentamers or hexamers on the capsid surface. By applying atomistic molecular dynamics simulation to an entire solvated HBV capsid subjected to a uniform mechanical stress protocol, we monitor the capsid-disassembly process and analyze the process down to the level of individual amino acids in 20 independent simulation replicas. The strain of an isotropic external force, combined with structural fluctuations, causes structurally heterogeneous cracks to appear in the HBV capsid. Analysis of the monomer-monomer interfaces reveals that, in contrast to the expectation from purely mechanical considerations, the cracks mainly occur within hexameric sites, whereas pentameric sites remain largely intact. Only a small subset of the capsid protein monomers, different in each simulation, are engaged in each instance of disassembly. We identify specific residues whose interactions are most readily lost during disassembly; R127, I139, Y132, N136, A137, and V149 are among the hot spots at the interfaces between dimers that lie within hexamers, leading to disassembly. The majority of these hot-spot residues are conserved by evolution, hinting to their importance for disassembly by avoiding overstabilization of capsids.

Modeling protein structures with the coarse-grained UNRES force field in the CASP14 experiment

The UNited RESidue (UNRES) force field was tested in the 14th Community Wide Experiment on the Critical Assessment of Techniques for Protein Structure Prediction (CASP14), in which larger oligomeric and multimeric targets were present compared to previous editions. Three prediction modes were tested (i) ab initio (the UNRES group), (ii) contact-assisted (the UNRES-contact group), and (iii) template-assisted (the UNRES-template group). For most of the targets, the contact restraints were derived from the server models top-ranked by the DeepQA method, while the DNCON2 method was used for 11 targets. Our consensus-fragment procedure was used to run template-assisted predictions. Each group also processed the Nuclear Magnetic Resonance (NMR)- and Small Angle X-Ray Scattering (SAXS)-data assisted targets. The average Global Distance Test Total Score (GDT_TS) of the 'Model 1' predictions were 29.17, 39.32, and 56.37 for the UNRES, UNRES-contact, and UNRES-template predictions, respectively, increasing by 0.53, 2.24, and 3.76, respectively, compared to CASP13. It was also found that the GDT_TS of the UNRES models obtained in ab initio mode and in the contact-assisted mode decreases with the square root of chain length, while the exponent in this relationship is 0.20 for the UNRES-template group models and 0.11 for the best performing AlphaFold2 models, which suggests that incorporation of database information, which stems from protein evolution, brings in long-range correlations, thus enabling the correction of force-field inaccuracies.

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.

Structural and dynamic asymmetry in icosahedrally symmetric virus capsids

A common characteristic of virus capsids is icosahedral symmetry, yet these highly symmetric structures can display asymmetric features within their virions and undergo asymmetric dynamics. The fields of structural and computational biology have entered a new realm in the investigation of virus infection mechanisms, with the ability to observe symmetry-breaking features. This review will cover important studies on icosahedral virus structure and dynamics, covering both symmetric and asymmetric conformational changes. However, the main emphasis of the review will be towards recent studies employing cryo-electron microscopy or molecular dynamics simulations, which can uncover asymmetric aspects of these systems relevant to understanding viral physical-chemical properties and their biological impact.

Multiscale modelling and simulation of viruses

In recent years, advances in structural biology, integrative modelling, and simulation approaches have allowed us to gain unprecedented insights into viral structure and dynamics. In this article we survey recent studies utilizing this wealth of structural information to build computational models of partial or complete viruses and to elucidate mechanisms of viral function. Additionally, the close interplay of viral pathogens with host factors - such as cellular and intracellular membranes, receptors, antibodies, and other host proteins - makes accurate models of viral interactions and dynamics essential. As viruses continue to pose severe challenges in prevention and treatment, enhancing our mechanistic understanding of viral infection is vital to enable the development of novel therapeutic strategies.

Rotational and Translational Diffusion of Proteins as a Function of Concentration

Atomistic simulations of three different proteins at different concentrations are performed to obtain insight into protein mobility as a function of protein concentration. We report on simulations of proteins from diluted to the physiological water concentration (about 70% of the mass). First, the viscosity was computed and found to increase by a factor of 7-9 going from pure water to the highest protein concentration, in excellent agreement with in vivo nuclear magnetic resonance results. At a physiological concentration of proteins, the translational diffusion is found to be slowed down to about 30% of the in vitro values. The slow-down of diffusion found here using atomistic models is slightly more than that of a hard sphere model that neglects the electrostatic interactions. Interestingly, rotational diffusion of proteins is slowed down somewhat more (by about 80-95% compared to in vitro values) than translational diffusion, in line with experimental findings and consistent with the increased viscosity. The finding that rotation is retarded more than translation is attributed to solvent-separated clustering. No direct interactions between the proteins are found, and the clustering can likely be attributed to dispersion interactions that are stronger between proteins than between protein and water. Based on these simulations, we can also conclude that the internal dynamics of the proteins in our study are affected only marginally under crowding conditions, and the proteins become somewhat more stable at higher concentrations. Simulations were performed using a force field that was tuned for dealing with crowding conditions by strengthening the protein-water interactions. This force field seems to lead to a reproducible partial unfolding of an α-helix in one of the proteins, an effect that was not observed in the unmodified force field.

High-throughput computational pipeline for 3-D structure preparation and in silico protein surface property screening: A case study on HBcAg dimer structures

Knowledge-based experimental design can aid biopharmaceutical high-throughput screening (HTS) experiments needed to identify critical manufacturability parameters. Prior knowledge can be obtained via computational methods such as protein property extraction from 3-D protein structures. This study presents a high-throughput 3-D structure preparation and refinement pipeline that supports structure screenings with an automated and data-dependent workflow. As a case study, three chimeric virus-like particle (VLP) building blocks, hepatitis B core antigen (HBcAg) dimers, were constructed. Molecular dynamics (MD) refinement quality, speed, stability, and correlation to zeta potential data was evaluated using different MD simulation settings. Settings included 2 force fields (YASARA2 and AMBER03) and 2 pKa computation methods (YASARA and H++). MD simulations contained a data-dependent termination via identification of a 2 ns Window of Stability, which was also used for robust descriptor extraction. MD simulation with YASARA2, independent of pKa computation method, was found to be most stable and computationally efficient. These settings resulted in a fast refinement (6.6-37.5 h), a good structure quality (-1.17--1.13) and a strong linear dependence between dimer surface charge and complete chimeric HBcAg VLP zeta potential. These results indicate the computational pipeline's applicability for early-stage candidate assessment and design optimization of HTS manufacturability or formulability experiments.

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.

Single-molecule studies of flavivirus envelope dynamics: Experiment and computation

Flaviviruses are simple enveloped viruses exhibiting complex structural and functional heterogeneities. Decades of research have provided crucial basic insights, antiviral medication and moderately successful gene therapy trials. The most infectious particle is, however, not always the most abundant one in a population, questioning the utility of classic ensemble-averaging virology approaches. Indeed, viral replication is often not particularly efficient, prone to errors or containing parallel routes. Here, we review different single-molecule sensitive fluorescence methods that are employed to investigate flaviviruses. In particular, we review how (i) time-resolved Förster resonance energy transfer (trFRET) was applied to probe dengue envelope conformations; (ii) FRET-fluorescence correlation spectroscopy to investigate dengue envelope intrinsic dynamics and (iii) single particle tracking to follow the path of dengue viruses in cells. We also discuss how such methods may be supported by molecular dynamics (MD) simulations over a range of spatio-temporal scales, to provide complementary data on the structure and dynamics of flaviviral systems. We describe recent improvements in multiscale MD approaches that allowed the simulation of dengue particle envelopes in near-atomic resolution. We hope this review is an incentive for setting up and applying similar single-molecule studies and combine them with MD simulations to investigate structural dynamics of entire flavivirus particles over the nanosecond-to-millisecond time-scale and follow viruses during infection in cells over milliseconds to minutes.

Scaffolds as Structural Tools for Bone-Targeted Drug Delivery

Although bone has a high potential to regenerate itself after damage and injury, the efficacious repair of large bone defects resulting from resection, trauma or non-union fractures still requires the implantation of bone grafts. Materials science, in conjunction with biotechnology, can satisfy these needs by developing artificial bones, synthetic substitutes and organ implants. In particular, recent advances in materials science have provided several innovations, underlying the increasing importance of biomaterials in this field. To address the increasing need for improved bone substitutes, tissue engineering seeks to create synthetic, three-dimensional scaffolds made from organic or inorganic materials, incorporating drugs and growth factors, to induce new bone tissue formation. This review emphasizes recent progress in materials science that allows reliable scaffolds to be synthesized for targeted drug delivery in bone regeneration, also with respect to past directions no longer considered promising. A general overview concerning modeling approaches suitable for the discussed systems is also provided.

All-atom virus simulations

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

All-atom molecular dynamics of the HBV capsid reveals insights into biological function and cryo-EM resolution limits

The hepatitis B virus capsid represents a promising therapeutic target. Experiments suggest the capsid must be flexible to function; however, capsid structure and dynamics have not been thoroughly characterized in the absence of icosahedral symmetry constraints. Here, all-atom molecular dynamics simulations are leveraged to investigate the capsid without symmetry bias, enabling study of capsid flexibility and its implications for biological function and cryo-EM resolution limits. Simulation results confirm flexibility and reveal a propensity for asymmetric distortion. The capsid’s influence on ionic species suggests a mechanism for modulating the display of cellular signals and implicates the capsid’s triangular pores as the location of signal exposure. A theoretical image reconstruction performed using simulated conformations indicates how capsid flexibility may limit the resolution of cryo-EM. Overall, the present work provides functional insight beyond what is accessible to experimental methods and raises important considerations regarding asymmetry in structural studies of icosahedral virus capsids.

How Computational Models Enable Mechanistic Insights into Virus Infection

An implicit aim in cellular infection biology is to understand the mechanisms how viruses, microbes, eukaryotic parasites, and fungi usurp the functions of host cells and cause disease. Mechanistic insight is a deep understanding of the biophysical and biochemical processes that give rise to an observable phenomenon. It is typically subject to falsification, that is, it is accessible to experimentation and empirical data acquisition. This is different from logic and mathematics, which are not empirical, but built on systems of inherently consistent axioms. Here, we argue that modeling and computer simulation, combined with mechanistic insights, yields unprecedented deep understanding of phenomena in biology and especially in virus infections by providing a way of showing sufficiency of a hypothetical mechanism. This ideally complements the necessity statements accessible to empirical falsification by additional positive evidence. We discuss how computational implementations of mathematical models can assist and enhance the quantitative measurements of infection dynamics of enveloped and non-enveloped viruses and thereby help generating causal insights into virus infection biology.

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.

Ammonia as an In Situ Sanitizer: Influence of Virus Genome Type on Inactivation

Treatment of human excreta and animal manure (HEAM) is key in controlling the spread of persistent enteric pathogens, such as viruses. The extent of virus inactivation during HEAM storage and treatment appears to vary with virus genome type, although the reasons for this variability are not clear. Here, we investigated the inactivation of viruses of different genome types under conditions representative of HEAM storage or mesophilic digestion. The goals were to characterize the influence of HEAM solution conditions on inactivation and to determine the potential mechanisms involved. Specifically, eight viruses representing the four viral genome types (single-stranded RNA [ssRNA], double-stranded RNA [dsRNA], single-stranded DNA [ssDNA], and double-stranded DNA [dsDNA]) were exposed to synthetic solutions with well-controlled temperature (20 to 35°C), pH (8 to 9), and ammonia (NH3) concentrations (0 to 40 mmol liter(-1)). DNA and dsRNA viruses were considerably more resistant than ssRNA viruses, resulting in up to 1,000-fold-longer treatment times to reach a 4-log inactivation. The apparently slower inactivation of DNA viruses was rationalized by the higher stability of DNA than that of ssRNA in HEAM. Pushing the system toward harsher pH (>9) and temperature (>35°C) conditions, such as those encountered in thermophilic digestion and alkaline treatments, led to more consistent inactivation kinetics among ssRNA and other viruses. This suggests that the dependence of inactivation on genome type disappeared in favor of protein-mediated inactivation mechanisms common to all viruses. Finally, we recommend the use of MS2 as a conservative indicator to assess the inactivation of ssRNA viruses and the stable ΦX174 or dsDNA phages as indicators for persistent viruses.

IMPORTANCE

Viruses are among the most environmentally persistent pathogens. They can be present in high concentrations in human excreta and animal manure (HEAM). Therefore, appropriate treatment of HEAM is important prior to its reuse or discharge into the environment. Here, we investigated the factors that determine the persistence of viruses in HEAM, and we determined the main mechanisms that lead to their inactivation. Unlike other organisms, viruses can have four different genome types (double- or single-stranded RNA or DNA), and the viruses studied herein represent all four types. Genome type appeared to be the major determinant for persistence. Single-stranded RNA viruses are the most labile, because this genome type is susceptible to degradation in HEAM. In contrast, the other genome types are more stable; therefore, inactivation is slower and mainly driven by the degradation of viral proteins. Overall, this study allows us to better understand the behavior of viruses in HEAM.

Physical Ingredients Controlling Stability and Structural Selection of Empty Viral Capsids

One of the crucial steps in the viral replication cycle is the self-assembly of its protein shell. Typically, each native virus adopts a unique architecture, but the coat proteins of many viruses have the capability to self-assemble in vitro into different structures by changing the assembly conditions. However, the mechanisms determining which of the possible capsid shapes and structures is selected by a virus are still not well-known. We present a coarse-grained model to analyze and understand the physical mechanisms controlling the size and structure selection in the assembly of empty viral capsids. Using this model and Monte Carlo simulations, we have characterized the phase diagram and stability of T = 1,3,4,7 and snub cube shells. In addition, we have studied the tolerance of different shells to changes in physical parameters related to ambient conditions, identifying possible strategies to induce misassembly or failure. Finally, we discuss the factors that select the shape of a capsid as spherical, faceted, elongated, or decapsidated. Our model sheds important light on the ingredients that control the assembly and stability of viral shells. This knowledge is essential to get capsids with well-defined size and structure that could be used for promising applications in medicine or bionanotechnology.

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.

Principles and Overview of Sampling Methods for Modeling Macromolecular Structure and Dynamics

Investigation of macromolecular structure and dynamics is fundamental to understanding how macromolecules carry out their functions in the cell. Significant advances have been made toward this end in silico, with a growing number of computational methods proposed yearly to study and simulate various aspects of macromolecular structure and dynamics. This review aims to provide an overview of recent advances, focusing primarily on methods proposed for exploring the structure space of macromolecules in isolation and in assemblies for the purpose of characterizing equilibrium structure and dynamics. In addition to surveying recent applications that showcase current capabilities of computational methods, this review highlights state-of-the-art algorithmic techniques proposed to overcome challenges posed in silico by the disparate spatial and time scales accessed by dynamic macromolecules. This review is not meant to be exhaustive, as such an endeavor is impossible, but rather aims to balance breadth and depth of strategies for modeling macromolecular structure and dynamics for a broad audience of novices and experts.

Picosecond infrared laser-induced all-atom nonequilibrium molecular dynamics simulation of dissociation of viruses

Laser-induced all-atom nonequilibrium molecular dynamics simulation of virus dissociation.

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.

Prediction of stability changes upon mutation in an icosahedral capsid

Identifying the contributions to thermodynamic stability of capsids is of fundamental and practical importance. Here we use simulation to assess how mutations affect the stability of lumazine synthase from the hyperthermophile Aquifex aeolicus, a T = 1 icosahedral capsid; in the simulations the icosahedral symmetry of the capsid is preserved by simulating a single pentamer and imposing crystal symmetry, in effect simulating an infinite cubic lattice of icosahedral capsids. The stability is assessed by estimating the free energy of association using an empirical method previously proposed to identify biological units in crystal structures. We investigate the effect on capsid formation of seven mutations, for which it has been experimentally assessed whether they disrupt capsid formation or not. With one exception, our approach predicts the effect of the mutations on the capsid stability. The method allows the identification of interaction networks, which drive capsid assembly, and highlights the plasticity of the interfaces between subunits in the capsid. Proteins 2015; 83:1733–1741. © 2015 The Authors. Proteins: Structure, Function, and Bioinformatics Published by Wiley Periodicals, Inc

SPECTRUS: A Dimensionality Reduction Approach for Identifying Dynamical Domains in Protein Complexes from Limited Structural Datasets

Identifying dynamical, quasi-rigid domains in proteins provides a powerful means for characterizing functionally oriented structural changes via a parsimonious set of degrees of freedom. In fact, the relative displacements of few dynamical domains usually suffice to rationalize the mechanics underpinning biological functionality in proteins and can even be exploited for structure determination or refinement purposes. Here we present SPECTRUS, a general scheme that, by solely using amino acid distance fluctuations, can pinpoint the innate quasi-rigid domains of single proteins or large complexes in a robust way. Consistent domains are usually obtained by using either a pair of representative structures or thousands of conformers. The functional insights offered by the approach are illustrated for biomolecular systems of very different size and complexity such as kinases, ion channels, and viral capsids. The decomposition tool is available as a software package and web server at spectrus.sissa.it.

Qgui: A high-throughput interface for automated setup and analysis of free energy calculations and empirical valence bond simulations in biological systems

Structural information and activity data has increased rapidly for many protein targets during the last decades. In this paper, we present a high-throughput interface (Qgui) for automated free energy and empirical valence bond (EVB) calculations that use molecular dynamics (MD) simulations for conformational sampling. Applications to ligand binding using both the linear interaction energy (LIE) method and the free energy perturbation (FEP) technique are given using the estrogen receptor (ERα) as a model system. Examples of free energy profiles obtained using the EVB method for the rate-limiting step of the enzymatic reaction catalyzed by trypsin are also shown. In addition, we present calculation of high-precision Arrhenius plots to obtain the thermodynamic activation enthalpy and entropy with Qgui from running a large number of EVB simulations.

Molecular dynamics simulations of large macromolecular complexes

Connecting dynamics to structural data from diverse experimental sources, molecular dynamics simulations permit the exploration of biological phenomena in unparalleled detail. Advances in simulations are moving the atomic resolution descriptions of biological systems into the million-to-billion atom regime, in which numerous cell functions reside. In this opinion, we review the progress, driven by large-scale molecular dynamics simulations, in the study of viruses, ribosomes, bioenergetic systems, and other diverse applications. These examples highlight the utility of molecular dynamics simulations in the critical task of relating atomic detail to the function of supramolecular complexes, a task that cannot be achieved by smaller-scale simulations or existing experimental approaches alone.

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.

Cholesterol modulates the dimer interface of the β₂-adrenergic receptor via cholesterol occupancy sites

The β2-adrenergic receptor is an important member of the G-protein-coupled receptor (GPCR) superfamily, whose stability and function are modulated by membrane cholesterol. The recent high-resolution crystal structure of the β2-adrenergic receptor revealed the presence of possible cholesterol-binding sites in the receptor. However, the functional relevance of cholesterol binding to the receptor remains unexplored. We used MARTINI coarse-grained molecular-dynamics simulations to explore dimerization of the β2-adrenergic receptor in lipid bilayers containing cholesterol. A novel (to our knowledge) aspect of our results is that receptor dimerization is modulated by membrane cholesterol. We show that cholesterol binds to transmembrane helix IV, and cholesterol occupancy at this site restricts its involvement at the dimer interface. With increasing cholesterol concentration, an increased presence of transmembrane helices I and II, but a reduced presence of transmembrane helix IV, is observed at the dimer interface. To our knowledge, this study is one of the first to explore the correlation between cholesterol occupancy and GPCR organization. Our results indicate that dimer plasticity is relevant not just as an organizational principle but also as a subtle regulatory principle for GPCR function. We believe these results constitute an important step toward designing better drugs for GPCR dimer targets.