Outcomes of the EMDataResource Cryo-EM Ligand Modeling Challenge

The EMDataResource Ligand Model Challenge aimed to assess the reliability and reproducibility of modeling ligands bound to protein and protein/nucleic-acid complexes in cryogenic electron microscopy (cryo-EM) maps determined at near-atomic (1.9-2.5 Å) resolution. Three published maps were selected as targets: E. coli beta-galactosidase with inhibitor, SARS-CoV-2 RNA-dependent RNA polymerase with covalently bound nucleotide analog, and SARS-CoV-2 ion channel ORF3a with bound lipid. Sixty-one models were submitted from 17 independent research groups, each with supporting workflow details. We found that (1) the quality of submitted ligand models and surrounding atoms varied, as judged by visual inspection and quantification of local map quality, model-to-map fit, geometry, energetics, and contact scores, and (2) a composite rather than a single score was needed to assess macromolecule+ligand model quality. These observations lead us to recommend best practices for assessing cryo-EM structures of liganded macromolecules reported at near-atomic resolution.

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.

Improved AlphaFold modeling with implicit experimental information

Machine-learning prediction algorithms such as AlphaFold and RoseTTAFold can create remarkably accurate protein models, but these models usually have some regions that are predicted with low confidence or poor accuracy. We hypothesized that by implicitly including new experimental information such as a density map, a greater portion of a model could be predicted accurately, and that this might synergistically improve parts of the model that were not fully addressed by either machine learning or experiment alone. An iterative procedure was developed in which AlphaFold models are automatically rebuilt on the basis of experimental density maps and the rebuilt models are used as templates in new AlphaFold predictions. We show that including experimental information improves prediction beyond the improvement obtained with simple rebuilding guided by the experimental data. This procedure for AlphaFold modeling with density has been incorporated into an automated procedure for interpretation of crystallographic and electron cryo-microscopy maps.

Putting AlphaFold models to work with phenix.process_predicted_model and ISOLDE

AlphaFold has recently become an important tool in providing models for experimental structure determination by X-ray crystallography and cryo-EM. Large parts of the predicted models typically approach the accuracy of experimentally determined structures, although there are frequently local errors and errors in the relative orientations of domains. Importantly, residues in the model of a protein predicted by AlphaFold are tagged with a predicted local distance difference test score, informing users about which regions of the structure are predicted with less confidence. AlphaFold also produces a predicted aligned error matrix indicating its confidence in the relative positions of each pair of residues in the predicted model. The phenix.process_predicted_model tool downweights or removes low-confidence residues and can break a model into confidently predicted domains in preparation for molecular replacement or cryo-EM docking. These confidence metrics are further used in ISOLDE to weight torsion and atom-atom distance restraints, allowing the complete AlphaFold model to be interactively rearranged to match the docked fragments and reducing the need for the rebuilding of connecting regions.

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.

Evolution of Intermediates during Capsid Assembly of Hepatitis B Virus with Phenylpropenamide-Based Antivirals

Self-assembly of virus capsids is a potential target for antivirals due to its importance in the virus lifecycle. Here, we investigate the effect of phenylpropenamide derivatives B-21 and AT-130 on the assembly of hepatitis B virus (HBV) core protein. Phenylpropenamides are widely believed to yield assembly of spherical particles resembling native, empty HBV capsids. Because the details of assembly can be overlooked with ensemble measurements, we performed resistive-pulse sensing on nanofluidic devices with four pores in series to characterize the size distributions of the products in real time. With its single particle sensitivity and compatibility with typical assembly buffers, resistive-pulse sensing is well-suited for analyzing virus assembly in vitro. We observed that assembly with B-21 and AT-130 produced a large fraction of partially complete virus particles that may be on-path, off-path, or trapped. For both B-21 and AT-130, capsid assembly was more sensitive to disruption under conditions where the interprotein association energy was low at lower salt concentrations. Dilution of the reaction solutions led to the rearrangement of the incomplete particles and demonstrated that these large intermediates may be on-path, but are labile, and exist in a frustrated dynamic equilibrium. During capsid assembly, phenylpropenamide molecules modestly increase the association energy of dimers, prevent intermediates from dissociating, and lead to kinetic trapping where the formation of too many capsids has been initiated, which results in both empty and incomplete particles.

Competition between Normative and Drug-Induced Virus Self-Assembly Observed with Single-Particle Methods

Disruption of virus capsid assembly has compelling antiviral potential that has been applied to hepatitis B virus (HBV). HBV core protein assembly can be modulated by heteroaryldihydropyrimidines (HAPs), and such molecules are collectively termed core protein allosteric modulators (CpAMs). Although the antiviral effects of CpAMs are acknowledged, the mechanism of action remains an open question. Challenging aspects of characterizing misdirected assembly are the large size and nonuniform nature of the final particles. In this study of HBV assembly, we observed a competition between normative and CpAM-induced aberrant assembly with electron microscopy and resistive-pulse sensing on nanofluidic devices. This competition was a function of the strength of the association energy between individual core proteins, which is proportional to ionic strength. At strong association energy, assembly reactions primarily yielded morphologically normal HBV capsids, despite the presence of HAP-TAMRA. At weak association energy, HAP-TAMRA led to increased assembly product size and disrupted morphology. The smallest particles were T = 4 icosahedra, whereas the larger particles were defective spheres, ellipsoids, and bacilliform cylinders, with regions of T = 4 geometry interspersed with flat regions. Deviation from spherical geometry progressively increased with particle size, which is consistent with the interpretation of a competition between two alternative assembly pathways.

An Assembly-Activating Site in the Hepatitis B Virus Capsid Protein Can Also Trigger Disassembly

The Hepatitis B Virus (HBV) core protein homodimers self-assemble to form an icosahedral capsid that packages the viral genome. Disassembly occurs in the nuclear basket to release the mature genome to the nucleus. Small molecules have been developed that bind to a pocket at the interdimer interface to accelerate assembly and strengthen interactions between subunits; these are under development as antiviral agents. Here, we explore the role of the dimer-dimer interface by mutating sites in the drug-binding pocket to cysteine and examining the effect of covalently linking small molecules to them. We find that ligands bound to the pocket may trigger capsid disassembly in a dose-dependent manner. This result indicates that, at least transiently, the pocket adopts a destabilizing conformation. We speculate that this pocket also plays a role in virus disassembly and genome release by binding ligands that are incompatible with virus stability, "unwanted guests." Investigating protein-protein interactions, especially large protein polymers, offers new and unique challenges. By using an engineered addressable thiol, we provide a means to examine the effects of modifying an interface without requiring drug-like properties for the ligand.