All-atom virus 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.

Integrative approaches for characterizing protein dynamics: NMR, CryoEM, and computer simulations

Proteins are inherently dynamic and their internal motions are essential for biological function. Protein motions cover a broad range of timescales: 10-14-10 s, spanning from sub-picosecond vibrational motions of atoms via microsecond loop conformational rearrangements to millisecond large amplitude domain reorientations. Observing protein dynamics over all timescales and connecting motions and structure to biological mechanisms requires integration of multiple experimental and computational techniques. This review reports on state-of-the-art approaches for assessing dynamics in biological systems using recent examples of virus assemblies, enzymes, and molecular machines. By integrating NMR spectroscopy in solution and the solid state, cryo electron microscopy, and molecular dynamics simulations, atomistic pictures of protein motions are obtained, not accessible from any single method in isolation. This information provides fundamental insights into protein behavior that can guide the development of future therapeutics.

Hexagonal Lattices of HIV Capsid Proteins Explored by Simulations Based on a Thermodynamically Consistent Model

HIV capsid proteins (CAs) may self-assemble into a variety of shapes under in vivo and in vitro conditions. Here, we employed simulations based on a residue-level coarse-grained (CG) model with full conformational flexibility to investigate hexagonal lattices, which are the underlying structural pattern for CA aggregations. Facilitated by enhanced sampling simulations to rigorously calculate CA dimerization and polymerization affinities, we calibrated our model to reproduce the experimentally measured affinities. Using the calibrated model, we performed unbiased simulations on several large systems consisting of 1512 CA subunits, allowing reversible binding and unbinding of the CAs in a thermodynamically consistent manner. In one simulation, a preassembled hexagonal CA sheet developed spontaneous curvatures reminiscent of those observed in experiments, and the edges of the sheet exhibited local curvatures larger than those of the interior. In other simulations starting with randomly distributed CAs at different concentrations, existing CA assemblies grew by binding free capsomeres to the edges and by merging with other assemblies. At high CA concentrations, rapid establishment of predominant aggregates was followed by much slower adjustments toward more regular hexagonal lattices, with increasing numbers of intact CA hexamers and pentamers being formed. Our approach of adapting a general CG model to specific systems by using experimental binding data represents a practical and effective strategy for simulating and elucidating intricate protein aggregations.

The SIRAH force field: A suite for simulations of complex biological systems at the coarse-grained and multiscale levels

The different combinations of molecular dynamics simulations with coarse-grained representations have acquired considerable popularity among the scientific community. Especially in biocomputing, the significant speedup granted by simplified molecular models opened the possibility of increasing the diversity and complexity of macromolecular systems, providing realistic insights on large assemblies for more extended time windows. However, a holistic view of biological ensembles' structural and dynamic features requires a self-consistent force field, namely, a set of equations and parameters that describe the intra and intermolecular interactions among moieties of diverse chemical nature (i.e., nucleic and amino acids, lipids, solvent, ions, etc.). Nevertheless, examples of such force fields are scarce in the literature at the fully atomistic and coarse-grained levels. Moreover, the number of force fields capable of handling simultaneously different scales is restricted to a handful. Among those, the SIRAH force field, developed in our group, furnishes a set of topologies and tools that facilitate the setting up and running of molecular dynamics simulations at the coarse-grained and multiscale levels. SIRAH uses the same classical pairwise Hamiltonian function implemented in the most popular molecular dynamics software. In particular, it runs natively in AMBER and Gromacs engines, and porting it to other simulation packages is straightforward. This review describes the underlying philosophy behind the development of SIRAH over the years and across families of biological molecules, discussing current limitations and future implementations.

Molecular Dynamics Simulations of Deformable Viral Capsomers

Most coarse-grained models of individual capsomers associated with viruses employ rigid building blocks that do not exhibit shape adaptation during self-assembly. We develop a coarse-grained general model of viral capsomers that incorporates their stretching and bending energies while retaining many features of the rigid-body models, including an overall trapezoidal shape with attractive interaction sites embedded in the lateral walls to favor icosahedral capsid assembly. Molecular dynamics simulations of deformable capsomers reproduce the rich self-assembly behavior associated with a general T=1 icosahedral virus system in the absence of a genome. Transitions from non-assembled configurations to icosahedral capsids to kinetically-trapped malformed structures are observed as the steric attraction between capsomers is increased. An assembly diagram in the space of capsomer-capsomer steric attraction and capsomer deformability reveals that assembling capsomers of higher deformability into capsids requires increasingly large steric attraction between capsomers. Increasing capsomer deformability can reverse incorrect capsomer-capsomer binding, facilitating transitions from malformed structures to symmetric capsids; however, making capsomers too soft inhibits assembly and yields fluid-like structures.

Mesoscale simulations: An indispensable approach to understand biomembranes

Computer simulation techniques form a versatile tool, a computational microscope, for exploring biological processes. This tool has been particularly effective in exploring different features of biological membranes. In recent years, thanks to elegant multiscale simulation schemes, some fundamental limitations of investigations by distinct simulation techniques have been resolved. As a result, we are now capable of exploring processes spanning multiple scales beyond the capacity of any single technique. In this perspective, we argue that mesoscale simulations require more attention and must be further developed to fill evident gaps in a quest toward simulating and modeling living cell membranes.

The stressed life of a lipid in the Zika virus membrane

Protein-lipid interactions modulate a plethora of physiopathologic processes and have been the subject of countless studies. However, these kinds of interactions in the context of viral envelopes have remained relatively unexplored, partially because the intrinsically small dimensions of the molecular systems escape to the current resolution of experimental techniques. However, coarse-grained and multiscale simulations may fill that niche, providing nearly atomistic resolution at an affordable computational price. Here we use multiscale simulations to characterize the lipid-protein interactions in the envelope of the Zika Virus, a prominent member of the Flavivirus genus. Comparisons between the viral envelope and simpler molecular systems indicate that the viral membrane is under extreme pressures and asymmetric forces. Furthermore, the dense net of protein-protein contacts established by the envelope proteins creates poorly solvated regions that destabilize the external leaflet leading to a decoupled dynamics between both membrane layers. These findings lead to the idea that the Flaviviral membrane may store a significant amount of elastic energy, playing an active role in the membrane fusion process.

Molecular dynamics of the viral life cycle: progress and prospects

Molecular dynamics (MD) simulations across spatiotemporal resolutions are widely applied to study viruses and represent the central technique uniting the field of computational virology. We discuss the progress of MD in elucidating the dynamics of the viral life cycle, including the status of modeling intact extracellular virions and leveraging advanced simulations to mimic active life cycle processes. We further remark on the prospects of MD for continued contributions to the basic science characterization of viruses, especially given the increasing availability of high-quality experimental data and supercomputing power. Overall, integrative computational methods that are closely guided by experiments are unmatched in the level of detail they provide, enabling-now and in the future-new discoveries relevant to thwarting viral infection.

Challenging the Existing Model of the Hexameric HIV-1 Gag Lattice and MA Shell Superstructure: Implications for Viral Entry

Despite type 1 human immunodeficiency virus (HIV-1) being discovered in the early 1980s, significant knowledge gaps remain in our understanding of the superstructure of the HIV-1 matrix (MA) shell. Current viral assembly models assume that the MA shell originates via recruitment of group-specific antigen (Gag) polyproteins into a hexagonal lattice but fails to resolve and explain lattice overlapping that occurs when the membrane is folded into a spherical/ellipsoidal shape. It further fails to address how the shell recruits, interacts with and encompasses the viral spike envelope (Env) glycoproteins. These Env glycoproteins are crucial as they facilitate viral entry by interacting with receptors and coreceptors located on T-cells. In our previous publication, we proposed a six-lune hosohedral structure, snowflake-like model for the MA shell of HIV-1. In this article, we improve upon the six-lune hosohedral structure by incorporating into our algorithm the recruitment of complete Env glycoproteins. We generated the Env glycoprotein assembly using a combination of predetermined Env glycoprotein domains from X-ray crystallography, nuclear magnetic resonance (NMR), cryoelectron tomography, and three-dimensional prediction tools. Our novel MA shell model comprises 1028 MA trimers and 14 Env glycoproteins. Our model demonstrates the movement of Env glycoproteins in the interlunar spaces, with effective clustering at the fusion hub, where multiple Env complexes bind to T-cell receptors during the process of viral entry. Elucidating the HIV-1 MA shell structure and its interaction with the Env glycoproteins is a key step toward understanding the mechanism of HIV-1 entry.

Algorithm to minimize MPI communications in the parallelized fast multipole method combined with molecular dynamics calculations

In the era of exascale supercomputers, large-scale, and long-time molecular dynamics (MD) calculations are expected to make breakthroughs in various fields of science and technology. Here, we propose a new algorithm to improve the parallelization performance of message passing interface (MPI)-communication in the MPI-parallelized fast multipole method (FMM) combined with MD calculations under three-dimensional periodic boundary conditions. Our approach enables a drastic reduction in the amount of communication data, including the atomic coordinates and multipole coefficients, both of which are required to calculate the electrostatic interaction by using the FMM. In communications of multipole coefficients, the reduction rate of communication data in the new algorithm relative to the amount of data in the conventional one increases as both the number of FMM levels and the number of MPI processes increase. The aforementioned rate increase could exceed 50% as the number of MPI processes becomes larger for very large systems. The proposed algorithm, named the minimum-transferred data (MTD) method, should enable large-scale and long-time MD calculations to be calculated efficiently, under the condition of massive MPI-parallelization on exascale supercomputers.

Molecular exclusion limits for diffusion across a porous capsid

Molecular communication across physical barriers requires pores to connect the environments on either side and discriminate between the diffusants. Here we use porous virus-like particles (VLPs) derived from bacteriophage P22 to investigate the range of molecule sizes able to gain access to its interior. Although there are cryo-EM models of the VLP, they may not accurately depict the parameters of the molecules able to pass across the pores due to the dynamic nature of the P22 particles in the solution. After encapsulating the enzyme AdhD within the P22 VLPs, we use a redox reaction involving PAMAM dendrimer modified NADH/NAD+ to examine the size and charge limitations of molecules entering P22. Utilizing the three different accessible morphologies of the P22 particles, we determine the effective pore sizes of each and demonstrate that negatively charged substrates diffuse across more readily when compared to those that are neutral, despite the negatively charge exterior of the particles.

Molecular determinants of Ebola nucleocapsid stability from molecular dynamics simulations

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

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.

All-Atom MD Simulations of the HBV Capsid Complexed with AT130 Reveal Secondary and Tertiary Structural Changes and Mechanisms of Allostery

The hepatitis B virus (HBV) capsid is an attractive drug target, relevant to combating viral hepatitis as a major public health concern. Among small molecules known to interfere with capsid assembly, the phenylpropenamides, including AT130, represent an important antiviral paradigm based on disrupting the timing of genome packaging. Here, all-atom molecular dynamics simulations of an intact AT130-bound HBV capsid reveal that the compound increases spike flexibility and improves recovery of helical secondary structure in the spike tips. Regions of the capsid-incorporated dimer that undergo correlated motion correspond to established sub-domains that pivot around the central chassis. AT130 alters patterns of correlated motion and other essential dynamics. A new conformational state of the dimer is identified, which can lead to dramatic opening of the intradimer interface and disruption of communication within the spike tip. A novel salt bridge is also discovered, which can mediate contact between the spike tip and fulcrum even in closed conformations, revealing a mechanism of direct communication across these sub-domains. Altogether, results describe a dynamical connection between the intra- and interdimer interfaces and enable mapping of allostery traversing the entire core protein dimer.

Physics of viral dynamics

Viral capsids are often regarded as inert structural units, but in actuality they display fascinating dynamics during different stages of their life cycle. With the advent of single-particle approaches and high-resolution techniques, it is now possible to scrutinize viral dynamics during and after their assembly and during the subsequent development pathway into infectious viruses. In this Review, the focus is on the dynamical properties of viruses, the different physical virology techniques that are being used to study them, and the physical concepts that have been developed to describe viral dynamics.

Permeability of the HIV-1 capsid to metabolites modulates viral DNA synthesis

Reverse transcription, an essential event in the HIV-1 life cycle, requires deoxynucleotide triphosphates (dNTPs) to fuel DNA synthesis, thus requiring penetration of dNTPs into the viral capsid. The central cavity of the capsid protein (CA) hexamer reveals itself as a plausible channel that allows the passage of dNTPs into assembled capsids. Nevertheless, the molecular mechanism of nucleotide import into the capsid remains unknown. Employing all-atom molecular dynamics (MD) simulations, we established that cooperative binding between nucleotides inside a CA hexamer cavity results in energetically favorable conditions for passive translocation of dNTPs into the HIV-1 capsid. Furthermore, binding of the host cell metabolite inositol hexakisphosphate (IP₆) enhances dNTP import, while binding of synthesized molecules like benzenehexacarboxylic acid (BHC) inhibits it. The enhancing effect on reverse transcription by IP₆ and the consequences of interactions between CA and nucleotides were corroborated using atomic force microscopy, transmission electron microscopy, and virological assays. Collectively, our results provide an atomistic description of the permeability of the HIV-1 capsid to small molecules and reveal a novel mechanism for the involvement of metabolites in HIV-1 capsid stabilization, nucleotide import, and reverse transcription.

This study shows that the HIV-1 capsid protein, in addition to its structural role, regulates reverse transcription, an essential metabolic process of the virus, by mediating the import of nucleotides. In addition, host cell metabolites such as inositol phosphates are recruited by the capsid to regulate viral DNA synthesis.

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.

Scalable Analysis of Authentic Viral Envelopes on FRONTERA

Enveloped viruses, such as SARS-CoV-2, infect cells via fusion of their envelope with the host membrane. By employing molecular simulations to characterize viral envelopes, researchers can gain insights into key determinants of infection. Here, the Frontera supercomputer is leveraged for large-scale modeling and analysis of authentic viral envelopes, whose lipid compositions are complex and realistic. Visual Molecular Dynamics (VMD) with support for MPI is employed, overcoming previous computational limitations and enabling investigation into virus biology at an unprecedented scale. The techniques applied here to an authentic HIV-1 envelope at two levels of spatial resolution (29 million particles and 280 million atoms) are broadly applicable to the study of other viruses. The authors are actively employing these techniques to develop and characterize an authentic SARS-CoV-2 envelope. A general framework for carrying out scalable analysis of simulation trajectories on Frontera is presented, expanding the utility of the machine in humanity’s ongoing fight against infectious diseases.

Coronavirus through Delaware's Computational Microscope

The Perilla/Hadden-Perilla research team at the University of Delaware presents an overview of computational structural biology, their efforts to model the SARS-CoV-2 viral particle, and their perspective on how their work and training endeavors can contribute to public health.

Mesoscale All-Atom Influenza Virus Simulations Suggest New Substrate Binding Mechanism

Influenza virus circulates in human, avian, and swine hosts, causing seasonal epidemic and occasional pandemic outbreaks. Influenza neuraminidase, a viral surface glycoprotein, has two sialic acid binding sites. The catalytic (primary) site, which also binds inhibitors such as oseltamivir carboxylate, is responsible for cleaving the sialic acid linkages that bind viral progeny to the host cell. In contrast, the functional annotation of the secondary site remains unclear. Here, we better characterize these two sites through the development of an all-atom, explicitly solvated, and experimentally based integrative model of the pandemic influenza A H1N1 2009 viral envelope, containing ∼160 million atoms and spanning ∼115 nm in diameter. Molecular dynamics simulations of this crowded subcellular environment, coupled with Markov state model theory, provide a novel framework for studying realistic molecular systems at the mesoscale and allow us to quantify the kinetics of the neuraminidase 150-loop transition between the open and closed states. An analysis of chloride ion occupancy along the neuraminidase surface implies a potential new role for the neuraminidase secondary site, wherein the terminal sialic acid residues of the linkages may bind before transfer to the primary site where enzymatic cleavage occurs. Altogether, our work breaks new ground for molecular simulation in terms of size, complexity, and methodological analyses of the components. It also provides fundamental insights into the understanding of substrate recognition processes for this vital influenza drug target, suggesting a new strategy for the development of anti-influenza therapeutics.

High-Performance Analysis of Biomolecular Containers to Measure Small-Molecule Transport, Transbilayer Lipid Diffusion, and Protein Cavities

Compartmentalization is a central theme in biology. Cells are composed of numerous membrane-enclosed structures, evolved to facilitate specific biochemical processes; viruses act as containers of genetic material, optimized to drive infection. Molecular dynamics simulations provide a mechanism to study biomolecular containers and the influence they exert on their environments; however, trajectory analysis software generally lacks knowledge of container interior versus exterior. Further, many relevant container analyses involve large-scale particle tracking endeavors, which may become computationally prohibitive with increasing system size. Here, a novel method based on 3-D ray casting is presented, which rapidly classifies the space surrounding biomolecular containers of arbitrary shape, enabling fast determination of the identities and counts of particles (e.g., solvent molecules) found inside and outside. The method is broadly applicable to the study of containers and enables high-performance characterization of properties such as solvent density, small-molecule transport, transbilayer lipid diffusion, and topology of protein cavities. The method is implemented in VMD, a widely used simulation analysis tool that supports personal computers, clouds, and parallel supercomputers, including ORNL's Summit and Titan and NCSA's Blue Waters, where the method can be employed to efficiently analyze trajectories encompassing millions of particles. The ability to rapidly characterize the spatial relationships of particles relative to a biomolecular container over many trajectory frames, irrespective of large particle counts, enables analysis of containers on a scale that was previously unfeasible, at a level of accuracy that was previously unattainable.

MS2 bacteriophage capsid studied using all-atom molecular dynamics

The all-atom model of an MS2 bacteriophage particle without its genome (the capsid) was built using high-resolution cryo-electron microscopy (EM) measurements for initial conformation. The structural characteristics of the capsid and the dynamics of the surrounding solution were examined using molecular dynamics simulation. The model demonstrates the overall preservation of the cryo-EM structure of the capsid at physiological conditions (room temperature and ions composition). The formation of a dense anion layer near the inner surface and a diffuse cation layer near the outer surface of the capsid was detected. The flow of water molecules and ions across the capsid through its pores were quantified, which was considerable for water and substantial for ions.