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

The tubular cavity of tobacco mosaic virus shields mechanical stress and regulates disassembly

Here we probe Tobacco mosaic virus (TMV) particles immobilized on a solid surface under transversal mechanical stress. We use atomic force microscopy to implement punctual deformation with high force (∼nN) that induces immediate virus rupture (single indentation assay), and continuous cycles of low force (∼100 pN) that generate a gradual disassembly of the virus particle (mechanical fatigue assay). These experiments are interpreted with the help of TMV coarse-grained and finite elements simulations, which indicate that the tubular cavity screens the transmission of mechanical stress from the top to the bottom half of the virion structure. Likewise, mechanical fatigue experiments reveal how TMV disassembles following growing transversal rifts with different dynamics that depend on a combination of the applied force and the tubular geometry of the virus. Our results indicate how the cylindrical cavity of TMV cushions the lower half of the virus structure from mechanical stress and regulates mechanical disassembly. STATEMENT OF SIGNIFICANCE: The inability of plant viruses like tobacco mosaic virus (TMV) to infect mammals makes them ideal for technological applications. While TMV is known for it's durability, it's unclear if this is due solely to its capsid proteins or its tubular structure. Using Atomic Force Microscopy, coarse-grained and finite elements models, we found that the tubular hole screens the transmission of mechanical stress from the top to the bottom half of the virion structure. This characteristic induces a stepwise disassembly process from intact to half virus, finishing in the virion disruption. Since the energies between proteins are comparable to those of other viruses, there is a protective effect of the tubular cavity that transcends the size down to the nanoscale.

Machine Learning of Molecular Dynamics Simulations Provides Insights into the Modulation of Viral Capsid Assembly

An effective approach in the development of novel antivirals is to target the assembly of viral capsids by using capsid assembly modulators (CAMs). CAMs targeting hepatitis B virus (HBV) have two major modes of function: they can either accelerate nucleocapsid assembly, retaining its structure, or misdirect it into noncapsid-like particles. Previous molecular dynamics (MD) simulations of early capsid-assembly intermediates showed differences in protein conformations for the apo and bound states. Here, we have developed and tested several classification machine learning (ML) models to better distinguish between apo-tetramer intermediates and those bound to accelerating or misdirecting CAMs. Models based on tertiary structural properties of the Cp149 tetramers and their interdimer orientation, as well as models based on direct and inverse contact distances between protein residues, were tested. All models distinguished the apo states and the two CAM-bound states with high accuracy. Furthermore, tertiary structure models and residue-distance models highlighted different tetramer regions as being important for classification. Both models can be used to better understand structural transitions that govern the assembly of nucleocapsids and to assist in the development of more potent CAMs. Finally, we demonstrate the utility of classification ML methods in comparing MD trajectories and describe our ML approaches, which can be extended to other systems of interest.

Structural studies of Parvoviridae capsid assembly and evolution: implications for novel AAV vector design

Adeno-associated virus (AAV) vectors have emerged as powerful tools in gene therapy, potentially treating various genetic disorders. Engineering the AAV capsids through computational methods enables the customization of these vectors to enhance their effectiveness and safety. This engineering allows for the development of gene therapies that are not only more efficient but also personalized to unique genetic profiles. When developing, it is essential to understand the structural biology and the vast techniques used to guide vector designs. This review covers the fundamental biology of the Parvoviridae capsids, focusing on modern structural study techniques, including (a) Cryo-electron microscopy and X-ray Crystallography studies and (b) Comparative analysis of capsid structures across different Parvoviridae species. Along with the structure and evolution of the Parvoviridae capsids, computational methods have provided significant insights into the design of novel AAV vector techniques, which include (a) Structure-guided design of AAV capsids with improved properties, (b) Directed Evolution of AAV capsids for specific applications, and (c) Computational prediction of AAV capsid-receptor interactions. Further discussion addressed the ongoing challenges in the AAV vector design and proposed future directions for exploring enhanced computational tools, such as artificial intelligence/machine learning and deep learning.

Induction of hepatitis B core protein aggregation targeting an unconventional binding site

The hepatitis B virus (HBV) infection is a major global health problem, with chronic infection leading to liver complications and high death toll. Current treatments, such as nucleos(t)ide analogs and interferon-α, effectively suppress viral replication but rarely cure the infection. To address this, new antivirals targeting different components of the HBV molecular machinery are being developed. Here we investigated the hepatitis B core protein (HBc) that forms the viral capsids and plays a vital role in the HBV life cycle. We explored two distinct binding pockets on the HBV capsid: the central hydrophobic pocket of HBc-dimers and the pocket at the tips of capsid spikes. We synthesized a geranyl dimer that binds to the central pocket with micromolar affinity, and dimeric peptides that bind the spike-tip pocket with sub-micromolar affinity. Cryo-electron microscopy further confirmed the binding of peptide dimers to the capsid spike tips and their capsid-aggregating properties. Finally, we show that the peptide dimers induce HBc aggregation in vitro and in living cells. Our findings highlight two tractable sites within the HBV capsid and provide an alternative strategy to affect HBV capsids.

Weighted Ensemble Simulations Reveal Novel Conformations and Modulator Effects in Hepatitis B Virus Capsid Assembly

Molecular dynamics (MD) simulations provide a detailed description of biophysical processes allowing mechanistic questions to be addressed at the atomic level. The promise of such approaches is partly hampered by well known sampling issues of typical simulations, where time scales available are significantly shorter than the process of interest. For the system of interest here, the binding of modulators of Hepatitis B virus capsid self-assembly, the binding site is at a flexible protein-protein interface. Characterization of the conformational landscape and how it is altered upon ligand binding is thus a prerequisite for a complete mechanistic description of capsid assembly modulation. However, such a description can be difficult due to the aforementioned sampling issues of standard MD, and enhanced sampling strategies are required. Here we employ the Weighted Ensemble methodology to characterize the free-energy landscape of our earlier determined functionally relevant progress coordinates. It is shown that this approach provides conformations outside those sampled by standard MD, as well as an increased number of structures with correspondingly enlarged binding pockets conducive to ligand binding, illustrating the utility of Weighted Ensemble for computational drug development.

Machine learning of molecular dynamics simulations provides insights into modulation of viral capsid assembly

An effective approach in the development of novel antivirals is to target the assembly of viral capsids using capsid assembly modulators (CAMs). CAMs targeting hepatitis B virus (HBV) have two major modes of function: they can either accelerate nucleocapsid assembly, retaining its structure, or misdirect it into non-capsid-like particles. Previous molecular dynamics (MD) simulations of early capsid-assembly intermediates showed differences in protein conformations for apo and bound states. Here, we have developed and tested several classification machine learning (ML) models to better distinguish between apo-tetramer intermediates and those bound to accelerating or misdirecting CAMs. Models based on tertiary structural properties of the Cp149 tetramers and their inter-dimer orientation, as well as models based on direct and inverse contact distances between protein residues, were tested. All models distinguished the apo states and the two CAM-bound states with high accuracy. Furthermore, tertiary structure models and residue-distance models highlighted different tetramer regions as important for classification. Both models can be used to better understand structural transitions that govern the assembly of nucleocapsids and to assist the development of more potent CAMs. Finally, we demonstrate the utility of classification ML methods in comparing MD trajectories and describe our ML approaches, which can be extended to other systems of interest.

Enabling Atomistic Modeling and Simulation of Complex Curved Cellular Membranes with xMAS Builder

As more powerful high performance computing resources are becoming available, there is a new opportunity to bring the unique capabilities of molecular dynamics (MD) simulations to cell-scale systems. Membranes are ubiquitous within cells and are responsible for a diverse set of essential biological functions, but building atomistic models of cell-scale membranes for MD simulations is immensely challenging because of their vast sizes, complex geometries, and complex compositions. To meet this challenge, we have developed xMAS Builder (Experimentally-Derived Membranes of Arbitrary Shape Builder), which is designed to take experimental lipidomics and structural (e.g., electron microscopy and tomography) data as input and use them to build MD-ready models of cellular membrane systems. To test xMAS Builder's capabilities, we have used it to build two models (one ~12.0 million atoms and the other ~11.6 million atoms) of a test system with a representative complex lipid composition and geometry. The two models, which differed only in their lipid packing densities, both maintained their membrane integrity during an extended MD simulation (250 ns and 386 ns), but their highly divergent relaxation dynamics indicate that the proper packing density of curved membranes is determined by leaflet volume rather than surface area. These results suggest that xMAS Builder's algorithms produce high quality models and that simulation of these models will provide profound biophysical insights into the behavior of cellular membranes.

From the Cytoplasm into the Nucleus-Hepatitis B Virus Travel and Genome Repair

Hepatitis B virus (HBV) is a major global health concern, affecting millions of people worldwide. HBV is part of the hepadnaviridae family and one of the primary causes of acute and chronic liver infections, leading to conditions such as cirrhosis and hepatocellular carcinoma (HCC). Understanding the intracellular transport and genome repair mechanisms of HBV is crucial for developing new drugs, which-in combination with immune modulators-may contribute to potential cures. This review will explore the current knowledge of HBV intracytoplasmic and nuclear transport, as well as genome repair processes, while drawing comparisons to other viruses with nuclear replication.

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

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

A review on description dynamics and conformational changes of proteins using combination of principal component analysis and molecular dynamics simulation

Understanding how proteins behave dynamically and undergo conformational changes is essential to comprehending their biological roles. This review article examines the potent tool of using Molecular Dynamics simulations in conjunction with Principal Component Analysis (PCA) to explore protein dynamics. Molecular dynamics data can be made easier to read by removing prominent patterns through the use of PCA, a sophisticated dimensionality reduction approach. Researchers can obtain critical insights into the fundamental principles governing protein function by using PCA on MD simulation data. We provide a systematic approach to PCA that includes data collection, input coordinate selection, and result interpretation. Protein collective movements and fundamental dynamics are made visible by PCA, which makes it possible to identify conformational substates that are crucial to function. By means of principal component analysis, scientists are able to observe and measure large-scale movements, like hinge bending and domain motions, as well as pinpoint areas of protein structural stiffness and flexibility. Moreover, PCA allows temporal separation, distinguishing slower global motions from faster local changes. A strong foundation for researching protein dynamics is provided by the combination of PCA and Molecular Dynamics simulations, which have applications in drug development and enhance our comprehension of intricate biological systems.

Effect of Pregenomic RNA on the Mechanical Stability of HBV Capsid by Coarse-Grained Molecular Simulations

Hepatitis B virus (HBV) is a double-stranded DNA virus, but its life cycle involves an intermediate stage, during which pregenomic RNA (pgRNA) is encapsulated in the capsid and then reverse-transcribed into the minus DNA strand. These immature HBV virions are the key target for antiviral drug discovery. In this study, we investigate the flexibility and mechanical stability of the HBV capsid containing pgRNA by employing residue-resolved coarse-grained molecular dynamics simulations. The results showed that the presence of pgRNA tends to decrease the overall flexibility of the capsid. In addition, the symmetrically arranged subunits of the capsid show asymmetry in the dominant modes of the conformational fluctuations with or without the presence of pgRNA. Furthermore, the simulations revealed that the presence of pgRNA enhances the overall mechanical stability of the virion particle. Electrostatic interactions between the disordered CTD of capsid and pgRNA were found to play a crucial role in modulating viral mechanical stability. Decreasing the electrostatic interactions by CTD phosphorylation or high salt concentration significantly reduces the mechanical stability of the HBV capsid containing pgRNA. Finally, the 2-fold symmetric sites have been proposed to be the most vulnerable to rupture during the initial stages of capsid disassembly. These findings could enhance our understanding of the physical basis of viral invasion and provide valuable insights into the development of antiviral drugs.

Small Molecule Assembly Agonist Alters the Dynamics of Hepatitis B Virus Core Protein Dimer and Capsid

Chronic hepatitis B virus (HBV) poses a significant public health burden worldwide, encouraging the search for curative antivirals. One approach is capsid assembly modulators (CAMs), which are assembly agonists. CAMs lead to empty and defective capsids, inhibiting the formation of new viruses, and can also lead to defects in the release of the viral genome, inhibiting new infections. In this study, we employed hydrogen-deuterium exchange mass spectrometry (HDX-MS) to assess the impact of one such CAM, HAP18, on HBV dimers, capsids composed of 120 (or 90) capsid protein dimers, and cross-linked capsids (xl-capsids). HDX analysis revealed hydrogen bonding networks within and between the dimers. HAP18 disrupted the hydrogen bonding network of dimers, demonstrating a previously unappreciated impact on the dimer structure. Conversely, HAP18 stabilized both unmodified and cross-linked capsids. Intriguingly, cross-linking the capsid, which was accomplished by forming disulfides between an engineered C-terminal cysteine, increased the overall rate of HDX. Moreover, HAP18 binding induced conformational changes beyond the binding sites. Our findings provide evidence for allosteric communication within and between capsid protein dimers. These results show that CAMs are capable of harnessing this allosteric network to modulate the dimer and capsid dynamics.

Guanidine Hydrochloride-Induced Hepatitis B Virus Capsid Disassembly Hysteresis

Hepatitis B virus (HBV) displays remarkable self-assembly capabilities that interest the scientific community and biotechnological industries as HBV is leading to an annual mortality of up to 1 million people worldwide (especially in Africa and Southeast Asia). When the ionic strength is increased, hepatitis B virus-like particles (VLPs) can assemble from dimers of the first 149 residues of the HBV capsid protein core assembly domain (Cp149). Using solution small-angle X-ray scattering, we investigated the disassembly of the VLPs by titrating guanidine hydrochloride (GuHCl). Measurements were performed with and without 1 M NaCl, added either before or after titrating GuHCl. Fitting the scattering curves to a linear combination of atomic models of Cp149 dimer (the subunit) and T = 3 and T = 4 icosahedral capsids revealed the mass fraction of the dimer in each structure in all the titration points. Based on the mass fractions, the variation in the dimer-dimer association standard free energy was calculated as a function of added GuHCl, showing a linear relation between the interaction strength and GuHCl concentration. Using the data, we estimated the energy barriers for assembly and disassembly and the critical nucleus size for all of the assembly reactions. Extrapolating the standard free energy to [GuHCl] = 0 showed an evident hysteresis in the assembly process, manifested by differences in the dimer-dimer association standard free energy obtained for the disassembly reactions compared with the equivalent assembly reactions. Similar hysteresis was observed in the energy barriers for assembly and disassembly and the critical nucleus size. The results suggest that above 1.5 M, GuHCl disassembled the capsids by attaching to the protein and adding steric repulsion, thereby weakening the hydrophobic attraction.

Fingolimod synergizes and reverses K. pneumoniae resistance to colistin

Klebsiella pneumoniae (K. pneumoniae) infection and the rapid spread of multi-drug resistant (MDR) bacteria pose a serious threat to global healthcare. Polymyxin E (colistin), a group of cationic antimicrobial polypeptides, is currently one of the last resort treatment options against carbapenem-resistant Gram-negative pathogens. The effectiveness of colistin has been compromised due to its intensive use. This study found that fingolimod (FLD), a natural product derivative, exhibited a significant synergistic bactericidal effect on K. pneumoniae when combined with colistin, both in vitro and in vivo. The checkerboard method was employed to assess the in vitro synergistic effect of FLD with colistin. FLD enhanced the susceptibility of bacteria to colistin and lowered effectively minimum inhibitory concentrations (MIC) when compared to colistin MIC, and the fractional inhibitory concentrations (FIC) value was less than 0.3. The time-kill curve demonstrated that the combination treatment of FLD and colistin had significant bactericidal efficacy. The in vitro concurrent administration of colistin and FLD resulted in heightening membrane permeability, compromising cell integrity, diminishing membrane fluidity, and perturbing membrane homeostasis. They also induced alterations in membrane potential, levels of reactive oxygen species, and adenosine triphosphate synthesis, ultimately culminating in bacterial death. Moreover, the combination of FLD with colistin significantly influenced fatty acid metabolism. In the mouse infection model, the survival rate of mice injected with K. pneumoniae was significantly improved to 67% and pathological damage was significantly relieved with combination treatment of FLD and colistin when compared with colistin treatment. This study highlights the potential of FLD in combining with colistin for treating infections caused by MDR isolates of K. pneumoniae.

Biophysics-Guided Lead Discovery of HBV Capsid Assembly Modifiers

Hepatitis B virus (HBV) is the leading cause of chronic liver pathologies worldwide. HBV nucleocapsid, a key structural component, is formed through the self-assembly of the capsid protein units. Therefore, interfering with the self-assembly process is a promising approach for the development of novel antiviral agents. Applied to HBV, this approach has led to several classes of capsid assembly modulators (CAMs). Here, we report structurally novel CAMs with moderate activity and low toxicity, discovered through a biophysics-guided approach combining docking, molecular dynamics simulations, and a series of assays with a particular emphasis on biophysical experiments. Several of the identified compounds induce the formation of aberrant capsids and inhibit HBV DNA replication in vitro, suggesting that they possess modest capsid assembly modulation effects. The synergistic computational and experimental approaches provided key insights that facilitated the identification of compounds with promising activities. The discovery of preclinical CAMs presents opportunities for subsequent optimization efforts, thereby opening new avenues for HBV inhibition.

Impact of DNA on interactions between core proteins of Hepatitis B virus-like particles comprising different C-terminals

Hepatitis B virus (HBV) virus-like particles (VLPs) are promising therapeutic agents derived from HBV core proteins (Cp). This study investigates the assembly dynamics of HBV VLPs, which is crucial for their potential as drug carriers or gene delivery systems. Coarse-grained molecular dynamics simulations explore the impact of C-terminal domain length (in the Cp ranging from Cp149 to wild-type Cp183) on Cp assembly and stability, particularly in the presence of DNA. Our findings reveal that the C-terminal nucleic acid binding region significantly influences Cp assembly and stability of trimers comprising Cp dimers. Shorter C-terminal domains (Cp164, Cp167) enhance stability and protein-protein interactions, while interactions between naturally occurring Cp183 are destabilized in the absence of DNA. Interestingly, DNA addition further stabilizes Cp assemblies, and this effect is influenced by the length of the nucleic acid binding region. Shorter C-terminal domains show less dependency on DNA content. This stabilization is attributed to electrostatic forces between positively charged C-terminal chains and negatively charged nucleic acids. Our study sheds light on the molecular mechanisms governing protein-protein and protein-DNA interactions in HBV VLP assembly, providing insights into Cp processability and informing the development of efficient gene therapy carriers using VLP technology.

Role of Protonation States in the Stability of Molecular Dynamics Simulations of High-Resolution Membrane Protein Structures

Classical molecular dynamics (MD) simulations provide unmatched spatial and time resolution of protein structure and function. However, the accuracy of MD simulations often depends on the quality of force field parameters and the time scale of sampling. Another limitation of conventional MD simulations is that the protonation states of titratable amino acid residues remain fixed during simulations, even though protonation state changes coupled to conformational dynamics are central to protein function. Due to the uncertainty in selecting protonation states, classical MD simulations are sometimes performed with all amino acids modeled in their standard charged states at pH 7. Here, we performed and analyzed classical MD simulations on high-resolution cryo-EM structures of two large membrane proteins that transfer protons by catalyzing protonation/deprotonation reactions. In simulations performed with titratable amino acids modeled in their standard protonation (charged) states, the structure diverges far from its starting conformation. In comparison, MD simulations performed with predetermined protonation states of amino acid residues reproduce the structural conformation, protein hydration, and protein-water and protein-protein interactions of the structure much better. The results support the notion that it is crucial to perform basic protonation state calculations, especially on structures where protonation changes play an important functional role, prior to the launch of any conventional MD simulations. Furthermore, the combined approach of fast protonation state prediction and MD simulations can provide valuable information about the charge states of amino acids in the cryo-EM sample. Even though accurate prediction of protonation states in proteinaceous environments currently remains a challenge, we introduce an approach of combining pKa prediction with cryo-EM density map analysis that helps in improving not only the protonation state predictions but also the atomic modeling of density data.

The structure and physical properties of a packaged bacteriophage particle

A string of nucleotides confined within a protein capsid contains all the instructions necessary to make a functional virus particle, a virion. Although the structure of the protein capsid is known for many virus species1,2, the three-dimensional organization of viral genomes has mostly eluded experimental probes3,4. Here we report all-atom structural models of an HK97 virion5, including its entire 39,732 base pair genome, obtained through multiresolution simulations. Mimicking the action of a packaging motor6, the genome was gradually loaded into the capsid. The structure of the packaged capsid was then refined through simulations of increasing resolution, which produced a 26 million atom model of the complete virion, including water and ions confined within the capsid. DNA packaging occurs through a loop extrusion mechanism7 that produces globally different configurations of the packaged genome and gives each viral particle individual traits. Multiple microsecond-long all-atom simulations characterized the effect of the packaged genome on capsid structure, internal pressure, electrostatics and diffusion of water, ions and DNA, and revealed the structural imprints of the capsid onto the genome. Our approach can be generalized to obtain complete all-atom structural models of other virus species, thereby potentially revealing new drug targets at the genome-capsid interface.

pyCapsid: identifying dominant dynamics and quasi-rigid mechanical units in protein shells

SUMMARY

pyCapsid is a Python package developed to facilitate the characterization of the dynamics and quasi-rigid mechanical units of protein shells and other protein complexes. The package was developed in response to the rapid increase of high-resolution structures, particularly capsids of viruses, requiring multiscale biophysical analyses. Given a protein shell, pyCapsid generates the collective vibrations of its amino-acid residues, identifies quasi-rigid mechanical regions associated with the disassembly of the structure, and maps the results back to the input proteins for interpretation. pyCapsid summarizes the main results in a report that includes publication-quality figures.

AVAILABILITY AND IMPLEMENTATION

pyCapsid's source code is available under MIT License on GitHub. It is compatible with Python 3.8-3.10 and has been deployed in two leading Python package-management systems, PIP and Conda. Installation instructions and tutorials are available in the online documentation and in the pyCapsid's YouTube playlist. In addition, a cloud-based implementation of pyCapsid is available as a Google Colab notebook. pyCapsid Colab does not require installation and generates the same report and outputs as the installable version. Users can post issues regarding pyCapsid in the repository's issues section.

Structure of the Hepatitis B virus capsid quasi-6-fold with a trapped C-terminal domain reveals capsid movements associated with domain exit

Many viruses undergo transient conformational change to surveil their environments for receptors and host factors. In Hepatitis B virus (HBV) infection, after the virus enters the cell, it is transported to the nucleus by interaction of the HBV capsid with an importin α/β complex. The interaction between virus and importins is mediated by nuclear localization signals on the capsid protein's C-terminal domain (CTD). However, CTDs are located inside the capsid. In this study, we asked where does a CTD exit the capsid, are all quasi-equivalent CTDs created equal, and does the capsid structure deform to facilitate CTD egress from the capsid? Here, we used Impβ as a tool to trap transiently exposed CTDs and examined this complex by cryo-electron microscopy. We examined an asymmetric reconstruction of a T = 4 icosahedral capsid and a focused reconstruction of a quasi-6-fold vertex (3.8 and 4.0 Å resolution, respectively). Both approaches showed that a subset of CTDs extended through a pore in the center of the quasi-6-fold complex. CTD egress was accompanied by enlargement of the pore and subtle changes in quaternary and tertiary structure of the quasi-6-fold. When compared to molecular dynamics simulations, structural changes were within the normal range of capsid flexibility. Although pore diameter was enlarged in the Impβ-bound reconstruction, simulations indicate that CTD egress does not exclusively depend on enlarged pores. In summary, we find that HBV surveillance of its environment by transient exposure of its CTD requires only modest conformational change of the capsid.

Biology of the hepatitis B virus (HBV) core and capsid assembly modulators (CAMs) for chronic hepatitis B (CHB) cure

Hepatitis B virus (HBV) is a hepadnavirus, a small DNA virus that infects liver tissue, with some unusual replication steps that share similarities to retroviruses. HBV infection can lead to chronic hepatitis B (CHB), a life-long infection associated with significant risks of liver disease, especially if untreated. HBV is a significant global health problem, with hundreds of millions currently living with CHB. Currently approved strategies to prevent or inhibit HBV are highly effective, however, a cure for CHB has remained elusive. To achieve a cure, elimination of the functionally integrated HBV covalently closed chromosomal DNA (cccDNA) genome is required. The capsid core is an essential component of HBV replication, serving roles when establishing infection and in creating new virions. Over the last two and a half decades, significant efforts have been made to find and characterize antivirals that target the capsid, specifically the HBV core protein (Cp). The antivirals that interfere with the kinetics and morphology of the capsid, termed capsid assembly modulators (CAMs), are extremely potent, and clinical investigations indicate they are well tolerated and highly effective. Several CAMs offer the potential to cure CHB by decreasing the cccDNA pools. Here, we review the biology of the HBV capsid, focused on Cp, and the development of inhibitors that target it.

Kinetic Growth of Multicomponent Microcompartment Shells

An important goal of systems and synthetic biology is to produce high value chemical species in large quantities. Microcompartments, which are protein nanoshells encapsulating catalytic enzyme cargo, could potentially function as tunable nanobioreactors inside and outside cells to generate these high value species. Modifying the morphology of microcompartments through genetic engineering of shell proteins is one viable strategy to tune cofactor and metabolite access to encapsulated enzymes. However, this is a difficult task without understanding how changing interactions between the many different types of shell proteins and enzymes affect microcompartment assembly and shape. Here, we use multiscale molecular dynamics and experimental data to describe assembly pathways available to microcompartments composed of multiple types of shell proteins with varied interactions. As the average interaction between the enzyme cargo and the multiple types of shell proteins is weakened, the shell assembly pathway transitions from (i) nucleating on the enzyme cargo to (ii) nucleating in the bulk and then binding the cargo as it grows to (iii) an empty shell. Atomistic simulations and experiments using the 1,2-propanediol utilization microcompartment system demonstrate that shell protein interactions are highly varied and consistent with our multicomponent, coarse-grained model. Furthermore, our results suggest that intrinsic bending angles control the size of these microcompartments. Overall, our simulations and experiments provide guidance to control microcomparmtent size and assembly by modulating the interactions between shell proteins.

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.

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.

Diffusion and molecular partitioning in hierarchically complex virus-like particles

Viruses are diverse infectious agents found in virtually every type of natural environment. Due to the range of conditions in which viruses have evolved, they exhibit a wide range of structure and function which has been exploited for biotechnology. The self-assembly process of virus-like particles (VLPs), derived from structural virus components, allows for the assembly of a hierarchy of materials. Because VLPs are robust in both their assembly and the final product, functionality can be incorporated through design of their building blocks or chemical modification after their synthesis and assembly. In particular, encapsulation of active enzymes inside VLP results in macromolecular concentration approximating that of cells, introducing excluded volume effects on encapsulated cargo which are not present in traditional experiments done on dilute proteins. This work reviews the hierarchical assembly of VLPs, experiments investigating diffusion in VLP systems, and methods for partitioning of chemical species in VLPs as functional biomaterials.

Molecular architecture and dynamics of SARS-CoV-2 envelope by integrative modeling

Despite tremendous efforts, the exact structure of SARS-CoV-2 and related betacoronaviruses remains elusive. SARS-CoV-2 envelope is a key structural component of the virion that encapsulates viral RNA. It is composed of three structural proteins, spike, membrane (M), and envelope, which interact with each other and with the lipids acquired from the host membranes. Here, we developed and applied an integrative multi-scale computational approach to model the envelope structure of SARS-CoV-2 with near atomistic detail, focusing on studying the dynamic nature and molecular interactions of its most abundant, but largely understudied, M protein. The molecular dynamics simulations allowed us to test the envelope stability under different configurations and revealed that the M dimers agglomerated into large, filament-like, macromolecular assemblies with distinct molecular patterns. These results are in good agreement with current experimental data, demonstrating a generic and versatile approach to model the structure of a virus de novo.

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

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

Molecular elucidation of drug-induced abnormal assemblies of the hepatitis B virus capsid protein by solid-state NMR

Hepatitis B virus (HBV) capsid assembly modulators (CAMs) represent a recent class of anti-HBV antivirals. CAMs disturb proper nucleocapsid assembly, by inducing formation of either aberrant assemblies (CAM-A) or of apparently normal but genome-less empty capsids (CAM-E). Classical structural approaches have revealed the CAM binding sites on the capsid protein (Cp), but conformational information on the CAM-induced off-path aberrant assemblies is lacking. Here we show that solid-state NMR can provide such information, including for wild-type full-length Cp183, and we find that in these assemblies, the asymmetric unit comprises a single Cp molecule rather than the four quasi-equivalent conformers typical for the icosahedral T = 4 symmetry of the normal HBV capsids. Furthermore, while in contrast to truncated Cp149, full-length Cp183 assemblies appear, on the mesoscopic level, unaffected by CAM-A, NMR reveals that on the molecular level, Cp183 assemblies are equally aberrant. Finally, we use a eukaryotic cell-free system to reveal how CAMs modulate capsid-RNA interactions and capsid phosphorylation. Our results establish a structural view on assembly modulation of the HBV capsid, and they provide a rationale for recently observed differences between in-cell versus in vitro capsid assembly modulation.

Breathing and Tilting: Mesoscale Simulations Illuminate Influenza Glycoprotein Vulnerabilities

Influenza virus has resurfaced recently from inactivity during the early stages of the COVID-19 pandemic, raising serious concerns about the nature and magnitude of future epidemics. The main antigenic targets of influenza virus are two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Whereas the structural and dynamical properties of both glycoproteins have been studied previously, the understanding of their plasticity in the whole-virion context is fragmented. Here, we investigate the dynamics of influenza glycoproteins in a crowded protein environment through mesoscale all-atom molecular dynamics simulations of two evolutionary-linked glycosylated influenza A whole-virion models. Our simulations reveal and kinetically characterize three main molecular motions of influenza glycoproteins: NA head tilting, HA ectodomain tilting, and HA head breathing. The flexibility of HA and NA highlights antigenically relevant conformational states, as well as facilitates the characterization of a novel monoclonal antibody, derived from convalescent human donor, that binds to the underside of the NA head. Our work provides previously unappreciated views on the dynamics of HA and NA, advancing the understanding of their interplay and suggesting possible strategies for the design of future vaccines and antivirals against influenza.

Multiscale Modeling of Hepatitis B Virus Capsid Assembly and Its Dimorphism

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

Breathing and tilting: mesoscale simulations illuminate influenza glycoprotein vulnerabilities

Influenza virus has resurfaced recently from inactivity during the early stages of the COVID-19 pandemic, raising serious concerns about the nature and magnitude of future epidemics. The main antigenic targets of influenza virus are two surface glycoproteins, hemagglutinin (HA) and neuraminidase (NA). Whereas the structural and dynamical properties of both glycoproteins have been studied previously, the understanding of their plasticity in the whole-virion context is fragmented. Here, we investigate the dynamics of influenza glycoproteins in a crowded protein environment through mesoscale all-atom molecular dynamics simulations of two evolutionary-linked glycosylated influenza A whole-virion models. Our simulations reveal and kinetically characterize three main molecular motions of influenza glycoproteins: NA head tilting, HA ectodomain tilting, and HA head breathing. The flexibility of HA and NA highlights antigenically relevant conformational states, as well as facilitates the characterization of a novel monoclonal antibody, derived from human convalescent plasma, that binds to the underside of the NA head. Our work provides previously unappreciated views on the dynamics of HA and NA, advancing the understanding of their interplay and suggesting possible strategies for the design of future vaccines and antivirals against influenza.

How Pore Architecture Regulates the Function of Nanoscale Protein Compartments

Self-assembling proteins can form porous compartments that adopt well-defined architectures at the nanoscale. In nature, protein compartments act as semipermeable barriers to enable spatial separation and organization of complex biochemical processes. The compartment pores play a key role in their overall function by selectively controlling the influx and efflux of important biomolecular species. By engineering the pores, the functionality of compartments can be tuned to facilitate non-native applications, such as artificial nanoreactors for catalysis. In this review, we analyze how protein structure determines the porosity and impacts the function of both native and engineered compartments, highlighting the wealth of structural data recently obtained by cryo-EM and X-ray crystallography. Through this analysis, we offer perspectives on how current structural insights can inform future studies into the design of artificial protein compartments as nanoreactors with tunable porosity and function.

Protein Cages: From Fundamentals to Advanced Applications

Proteins that self-assemble into polyhedral shell-like structures are useful molecular containers both in nature and in the laboratory. Here we review efforts to repurpose diverse protein cages, including viral capsids, ferritins, bacterial microcompartments, and designed capsules, as vaccines, drug delivery vehicles, targeted imaging agents, nanoreactors, templates for controlled materials synthesis, building blocks for higher-order architectures, and more. A deep understanding of the principles underlying the construction, function, and evolution of natural systems has been key to tailoring selective cargo encapsulation and interactions with both biological systems and synthetic materials through protein engineering and directed evolution. The ability to adapt and design increasingly sophisticated capsid structures and functions stands to benefit the fields of catalysis, materials science, and medicine.

Hysteresis in Hepatitis B Virus (HBV) Requires Assembly of Near-Perfect Capsids

The hepatitis B virus (HBV) must release its contents to initiate infection, making capsid disassembly critical to the viral life cycle. Capsid assembly proceeds through a cascade of weak interactions between copies of capsid protein (Cp) to yield uniform particles. However, there is a hysteresis to capsid dissociation that allows capsids to persist under conditions where they could not assemble. In this study, we have sought to define the basis of hysteresis by examining urea-induced dissociation of in vitro-assembled HBV capsids. In general, capsid samples show a mixture of two pools, differentiated by stability. Labile capsid dissociation corresponds to an ∼5 μM pseudocritical concentration of assembly (pcc), the same as that observed in assembly reactions. Dissociation of the stable pool corresponds to a subfemtomolar pcc, indicative of hysteresis. The fraction of stable capsids in an assembly reaction increases with the integrity of the Cp preparation and when association is performed at a higher ionic strength, which modifies the Cp conformation. Labile complexes are more prevalent when assembly conditions yield many kinetically trapped (incomplete and overgrown) products. Cp isolated from stable capsids reassembles into a mixture of stable and labile capsids. These results suggest that hysteresis arises from an ideal capsid lattice, even when some of the substituents in that lattice have defects. Consistent with structural studies that show a subtle difference between Cp dimers and Cp in capsid, we propose that hysteresis arises when HBV capsids undergo a lattice-dependent structural transition.

The Mechanism of Action of Hepatitis B Virus Capsid Assembly Modulators Can Be Predicted from Binding to Early Assembly Intermediates

Interfering with the self-assembly of virus nucleocapsids is a promising approach for the development of novel antiviral agents. Applied to hepatitis B virus (HBV), this approach has led to several classes of capsid assembly modulators (CAMs) that target the virus by either accelerating nucleocapsid assembly or misdirecting it into noncapsid-like particles, thereby inhibiting the HBV replication cycle. Here, we have assessed the structures of early nucleocapsid assembly intermediates, bound with and without CAMs, using molecular dynamics simulations. We find that distinct conformations of the intermediates are induced depending on whether the bound CAM accelerates or misdirects assembly. Specifically, the assembly intermediates with bound misdirecting CAMs appear to be flattened relative to those with bound accelerators. Finally, the potency of CAMs within the same class was studied. We find that an increased number of contacts with the capsid protein and favorable binding energies inferred from free energy perturbation calculations are indicative of increased potency.

The ugly, bad, and good stories of large-scale biomolecular simulations

Molecular modeling of large biomolecular assemblies exemplifies a disruptive area holding both promises and contentions. Propelled by peta and exascale computing, several simulation methodologies have now matured into user-friendly tools that are successfully employed for modeling viruses, membranous nano-constructs, and key pieces of the genetic machinery. We present three unifying biophysical themes that emanate from some of the most recent multi-million atom simulation endeavors. Despite connecting molecular changes with phenotypic outcomes, the quality measures of these simulations remain questionable. We discuss the existing and upcoming strategies for constructing representative ensembles of large systems, how new computing technologies will boost this area, and make a point that integrative modeling guided by experimental data is the future of biomolecular computations.

Experimental Characterization of the Hepatitis B Virus Capsid Dynamics by Solid-State NMR

Protein plasticity and dynamics are important aspects of their function. Here we use solid-state NMR to experimentally characterize the dynamics of the 3.5 MDa hepatitis B virus (HBV) capsid, assembled from  240 copies of the Cp149 core protein. We measure both T₁ and T_1ρ relaxation times, which we use to establish detectors on the nanosecond and microsecond timescale. We compare our results to those from a 1 microsecond all-atom Molecular Dynamics (MD) simulation trajectory for the capsid. We show that, for the constituent residues, nanosecond dynamics are faithfully captured by the MD simulation. The calculated values can be used in good approximation for the NMR-non-detected residues, as well as to extrapolate into the range between the nanosecond and microsecond dynamics, where NMR has a blind spot at the current state of technology. Slower motions on the microsecond timescale are difficult to characterize by all-atom MD simulations owing to computational expense, but are readily accessed by NMR. The two methods are, thus, complementary, and a combination thereof can reliably characterize motions covering correlation times up to a few microseconds.

Disassembly of Single Virus Capsids Monitored in Real Time with Multicycle Resistive-Pulse Sensing

Virus assembly and disassembly are critical steps in the virus lifecycle; however, virus disassembly is much less well understood than assembly. For hepatitis B virus (HBV) capsids, disassembly of the virus capsid in the presence of guanidine hydrochloride (GuHCl) exhibits strong hysteresis that requires additional chemical energy to initiate disassembly and disrupt the capsid structure. To study disassembly of HBV capsids, we mixed T = 4 HBV capsids with 1.0-3.0 M GuHCl, monitored the reaction over time by randomly selecting particles, and measured their size with resistive-pulse sensing. Particles were cycled forward and backward multiple times to increase the observation time and likelihood of observing a disassembly event. The four-pore device used for resistive-pulse sensing produces four current pulses for each particle during translocation that improves tracking and identification of single particles and increases the precision of particle-size measurements when pulses are averaged. We studied disassembly at GuHCl concentrations below and above denaturing conditions of the dimer, the fundamental unit of HBV capsid assembly. As expected, capsids showed little disassembly at low GuHCl concentrations (e.g., 1.0 M GuHCl), whereas at higher GuHCl concentrations (≥1.5 M), capsids exhibited disassembly, sometimes as a complex series of events. In all cases, disassembly was an accelerating process, where capsids catastrophically disassembled within a few 100 ms of reaching critical stability; disassembly rates reached tens of dimers per second just before capsids fell apart. Some disassembly events exhibited metastable intermediates that appeared to lose one or more trimers of dimers in a stepwise fashion.

Initiating Coarse-Grained MD Simulations for Membrane-Bound Proteins

Molecular dynamics (MD) simulations have become a widely used tool in the scientific community for understanding molecular scale phenomena that are challenging to address with wet-lab techniques. Coarse-grained simulations, in which multiple atoms are combined into single beads, allow for larger systems and longer time scales to be explored than atomistic techniques. Here, we describe the procedures and equipment required to set up coarse-grained simulations of membrane-bound proteins in a lipid bilayer that can mimic many membrane environments.

Current Progress in the Development of Hepatitis B Virus Capsid Assembly Modulators: Chemical Structure, Mode-of-Action and Efficacy

Hepatitis B virus (HBV) is a major causative agent of human hepatitis. Its viral genome comprises partially double-stranded DNA, which is complexed with viral polymerase within an icosahedral capsid consisting of a dimeric core protein. Here, we describe the effects of capsid assembly modulators (CAMs) on the geometric or kinetic disruption of capsid construction and the virus life cycle. We highlight classical, early-generation CAMs such as heteroaryldihydropyrimidines, phenylpropenamides or sulfamoylbenzamides, and focus on the chemical structure and antiviral efficacy of recently identified non-classical CAMs, which consist of carboxamides, aryl ureas, bithiazoles, hydrazones, benzylpyridazinones, pyrimidines, quinolines, dyes, and antimicrobial compounds. We summarize the therapeutic efficacy of four representative classical compounds with data from clinical phase 1 studies in chronic HBV patients. Most of these compounds are in phase 2 trials, either as monotherapy or in combination with approved nucleos(t)ides drugs or other immunostimulatory molecules. As followers of the early CAMs, the therapeutic efficacy of several non-classical CAMs has been evaluated in humanized mouse models of HBV infection. It is expected that these next-generation HBV CAMs will be promising candidates for a series of extended human clinical trials.

Viral Phrenology

We introduce Viral Phrenology, a new scheme for understanding the genomic composition of spherical viruses based on the locations of their structural protrusions. We used icosahedral point arrays to classify 135 distinct viral capsids collected from over 600 capsids available in the VIPERdb. Using gauge points of point arrays, we found 149 unique structural protrusions. We then show how to use the locations of these protrusions to determine the genetic composition of the virus. We then show that ssDNA, dsDNA, dsRNA and ssRNA viruses use different arrangements for distributing their protrusions. We also found that Triangulation number is also partially dependent on the structural protrusions. This analysis begins to tie together Baltimore Classification and Triangulation number using point arrays.

The Hepatitis B Virus Nucleocapsid-Dynamic Compartment for Infectious Virus Production and New Antiviral Target

Hepatitis B virus (HBV) is a small enveloped DNA virus which replicates its tiny 3.2 kb genome by reverse transcription inside an icosahedral nucleocapsid, formed by a single ~180 amino acid capsid, or core, protein (Cp). HBV causes chronic hepatitis B (CHB), a severe liver disease responsible for nearly a million deaths each year. Most of HBV's only seven primary gene products are multifunctional. Though less obvious than for the multi-domain polymerase, P protein, this is equally crucial for Cp with its multiple roles in the viral life-cycle. Cp provides a stable genome container during extracellular phases, allows for directed intracellular genome transport and timely release from the capsid, and subsequent assembly of new nucleocapsids around P protein and the pregenomic (pg) RNA, forming a distinct compartment for reverse transcription. These opposing features are enabled by dynamic post-transcriptional modifications of Cp which result in dynamic structural alterations. Their perturbation by capsid assembly modulators (CAMs) is a promising new antiviral concept. CAMs inappropriately accelerate assembly and/or distort the capsid shell. We summarize the functional, biochemical, and structural dynamics of Cp, and discuss the therapeutic potential of CAMs based on clinical data. Presently, CAMs appear as a valuable addition but not a substitute for existing therapies. However, as part of rational combination therapies CAMs may bring the ambitious goal of a cure for CHB closer to reality.

Binding of a Pocket Factor to Hepatitis B Virus Capsids Changes the Rotamer Conformation of Phenylalanine 97

(1) Background: During maturation of the Hepatitis B virus, a viral polymerase inside the capsid transcribes a pre-genomic RNA into a partly double stranded DNA-genome. This is followed by envelopment with surface proteins inserted into a membrane. Envelopment is hypothetically regulated by a structural signal that reports the maturation state of the genome. NMR data suggest that such a signal can be mimicked by the binding of the detergent Triton X 100 to hydrophobic pockets in the capsid spikes. (2) Methods: We have used electron cryo-microscopy and image processing to elucidate the structural changes that are concomitant with the binding of Triton X 100. (3) Results: Our maps show that Triton X 100 binds with its hydrophobic head group inside the pocket. The hydrophilic tail delineates the outside of the spike and is coordinated via Lys-96. The binding of Triton X 100 changes the rotamer conformation of Phe-97 in helix 4, which enables a π-stacking interaction with Trp-62 in helix 3. Similar changes occur in mutants with low secretion phenotypes (P5T and L60V) and in a mutant with a pre-mature secretion phenotype (F97L). (4) Conclusion: Binding of Triton X 100 is unlikely to mimic structural maturation because mutants with different secretion phenotypes show similar structural responses.

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.

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.

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.

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.

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.

Conformational Plasticity of Hepatitis B Core Protein Spikes Promotes Peptide Binding Independent of the Secretion Phenotype

Hepatitis B virus is a major human pathogen, which forms enveloped virus particles. During viral maturation, membrane-bound hepatitis B surface proteins package hepatitis B core protein capsids. This process is intercepted by certain peptides with an “LLGRMKG” motif that binds to the capsids at the tips of dimeric spikes. With microcalorimetry, electron cryo microscopy and peptide microarray-based screens, we have characterized the structural and thermodynamic properties of peptide binding to hepatitis B core protein capsids with different secretion phenotypes. The peptide “GSLLGRMKGA” binds weakly to hepatitis B core protein capsids and mutant capsids with a premature (F97L) or low-secretion phenotype (L60V and P5T). With electron cryo microscopy, we provide novel structures for L60V and P5T and demonstrate that binding occurs at the tips of the spikes at the dimer interface, splaying the helices apart independent of the secretion phenotype. Peptide array screening identifies “SLLGRM” as the core binding motif. This shortened motif binds only to one of the two spikes in the asymmetric unit of the capsid and induces a much smaller conformational change. Altogether, these comprehensive studies suggest that the tips of the spikes act as an autonomous binding platform that is unaffected by mutations that affect secretion phenotypes.