Immature HIV-1 lattice assembly dynamics are regulated by scaffolding from nucleic acid and the plasma membrane

Human Immunodeficiency Virus Type 1 Gag Polyprotein Modulates Membrane Physical Properties like a Surfactant: Potential Implications for Virus Assembly

Human immunodeficiency virus (HIV) assembly at an infected cell's plasma membrane requires membrane deformation to organize the near-spherical shape of an immature virus. While the cellular expression of HIV Gag is sufficient to initiate budding of virus-like particles, how Gag generates membrane curvature is not fully understood. Using highly curved lipid nanotubes, we have investigated the physicochemical basis of the membrane activity of recombinant nonmyristoylated Gag-Δp6. Gag protein, upon adsorption onto the membrane, resulted in the shape changes of both charged and uncharged nanotubes. This shape change was more pronounced in the presence of charged lipids, especially phosphatidylinositol bisphosphate (PI(4,5)P2). We found that Gag modified the interfacial tension of phospholipid bilayer membranes, as judged by comparison with the effects of amphipathic peptides and nonionic detergent. Bioinformatic analysis demonstrated that a region of the capsid and SP1 domains junction of Gag is structurally similar to the amphipathic peptide magainin-1. This region accounts for integral changes in the physical properties of the membrane upon Gag adsorption, as we showed with the synthetic CA-SP1 junction peptide. Phenomenologically, membrane-adsorbed Gag could diminish the energetic cost of increasing the membrane area in a way similar to foam formation. We propose that Gag acts as a surface-active substance at the HIV budding site that softens the membrane at the place of Gag adsorption, lowering the energy for membrane bending. Finally, our experimental data and theoretical considerations give a lipid-centric view and common mechanism by which proteins could bend membranes, despite not having intrinsic curvature in their molecular surfaces or assemblies.

Optimizing Coarse-Grained Models for Large-Scale Membrane Protein Simulation

Coarse-grained (CG) models have been developed for studying membrane proteins at physiologically relevant scales. Such methods, including popular CG lipid models, exhibit stability and efficiency at moderate scales, but they can become impractical or even unusable beyond a critical size due to various technical issues. Here, we report that these scale-dependent issues can arise from progressively slower relaxation dynamics and become confounded by unforeseen instabilities observed only at larger scales. To address these issues, we systemically optimized a 4-site solvent-free CG lipid model that is suitable for conducting micron-scale molecular dynamics simulations of membrane proteins under various membrane properties. We applied this lipid model to explore the long-range membrane deformation induced by a large mechanosensitive ion channel, PIEZO. We show that the optimized CG models are powerful in elucidating the structural and dynamic interplay between PIEZO and the membrane. Furthermore, we anticipate that our methodological insights can prove useful for resolving issues stemming from scale-dependent limitations of similar CG methodologies.

Human ESCRT-I and ALIX function as scaffolding helical filaments in vivo

The Endosomal Sorting Complex Required for Transport (ESCRT) is an evolutionarily conserved machinery that performs reverse-topology membrane scission in cells universally required from cytokinesis to budding of enveloped viruses. Upstream acting ESCRT-I and ALIX control these events and link recruitment of viral and cellular partners to late-acting ESCRT-III CHMP4 through incompletely understood mechanisms. Using structure-function analyses combined with super-resolution imaging, we show that ESCRT-I and ALIX function as distinct helical filaments in vivo . Together, they are essential for optimal structural scaffolding of HIV-1 nascent virions, the retention of viral and human genomes through defined functional interfaces, and recruitment of CHMP4 that itself assembles into corkscrew-like filaments intertwined with ESCRT-I or ALIX helices. Disruption of filament assembly or their conformationally clustered RNA binding interfaces in human cells impaired membrane abscission, resulted in major structural instability and leaked nucleic acid from nascent virions and nuclear envelopes. Thus, ESCRT-I and ALIX function as helical filaments in vivo and serve as both nucleic acid-dependent structural scaffolds as well as ESCRT-III assembly templates.

Significance statement

When cellular membranes are dissolved or breached, ESCRT is rapidly deployed to repair membranes to restore the integrity of intracellular compartments. Membrane sealing is ensured by ESCRT-III filaments assembled on the inner face of membrane; a mechanism termed inverse topology membrane scission. This mechanism, initiated by ESCRT-I and ALIX, is universally necessary for cytokinesis, wound repair, budding of enveloped viruses, and more. We show ESCRT-I and ALIX individually oligomerize into helical filaments that cluster newly discovered nucleic acid-binding interfaces and scaffold-in genomes within nascent virions and nuclear envelopes. These oligomers additionally appear to serve as ideal templates for ESCRT-III polymerization, as helical filaments of CHMP4B were found intertwined ESCRT-I or ALIX filaments in vivo . Similarly, corkscrew-like filaments of ALIX are also interwoven with ESCRT-I, supporting a model of inverse topology membrane scission that is synergistically reinforced by inward double filament scaffolding.

Hierarchical Framework for Predicting Entropies in Bottom-Up Coarse-Grained Models

The thermodynamic entropy of coarse-grained (CG) models stands as one of the most important properties for quantifying the missing information during the CG process and for establishing transferable (or extendible) CG interactions. However, performing additional CG simulations on top of model construction often leads to significant additional computational overhead. In this work, we propose a simple hierarchical framework for predicting the thermodynamic entropies of various molecular CG systems. Our approach employs a decomposition of the CG interactions, enabling the estimation of the CG partition function and thermodynamic properties a priori. Starting from the ideal gas description, we leverage classical perturbation theory to systematically incorporate simple yet essential interactions, ranging from the hard sphere model to the generalized van der Waals model. Additionally, we propose an alternative approach based on multiparticle correlation functions, allowing for systematic improvements through higher-order correlations. Numerical applications to molecular liquids validate the high fidelity of our approach, and our computational protocols demonstrate that a reduced model with simple energetics can reasonably estimate the thermodynamic entropy of CG models without performing any CG simulations. Overall, our findings present a systematic framework for estimating not only the entropy but also other thermodynamic properties of CG models, relying solely on information from the reference system.

Cooperative Membrane Binding of HIV-1 Matrix Proteins

The HIV-1 assembly process begins with a newly synthesized Gag polyprotein being targeted to the inner leaflet of the plasma membrane of the infected cells to form immature viral particles. Gag-membrane interactions are mediated through the myristoylated (Myr) N-terminal matrix (MA) domain of Gag, which eventually multimerize on the membrane to form trimers and higher order oligomers. The study of the structure and dynamics of peripheral membrane proteins like MA has been challenging for both experimental and computational studies due to the complex transient dynamics of protein-membrane interactions. Although the roles of anionic phospholipids (PIP2, PS) and the Myr group in the membrane targeting and stable membrane binding of MA are now well-established, the cooperative interactions between the MA monomers and MA-membrane remain elusive in the context of viral assembly and release. Our present study focuses on the membrane binding dynamics of a higher order oligomeric structure of MA protein (a dimer of trimers), which has not been explored before. Employing time-lagged independent component analysis (tICA) to our microsecond-long trajectories, we investigate conformational changes of the matrix protein induced by membrane binding. Interestingly, the Myr switch of an MA monomer correlates with the conformational switch of adjacent monomers in the same trimer. Together, our findings suggest complex protein dynamics during the formation of the immature HIV-1 lattice; while MA trimerization facilitates Myr insertion, MA trimer-trimer interactions in the immature lattice can hinder the same.

Roles of RNA scaffolding in nanoscale Gag multimerization and selective protein sorting at HIV membranes

HIV-1 Gag proteins can multimerize upon the viral genomic RNA or multiple random cellular messenger RNAs to form a virus particle or a virus-like particle, respectively. To date, whether the two types of particles form via the same Gag multimerization process has remained unclarified. Using photoactivated localization microscopy to illuminate Gag organizations and dynamics at the nanoscale, here, we showed that genomic RNA mediates Gag multimerization in a more cluster-centric, cooperative, and spatiotemporally coordinated fashion, with the ability to drive dense Gag clustering dependent on its ability to act as a long-stranded scaffold not easily attainable by cellular messenger RNAs. These differences in Gag multimerization were further shown to affect downstream selective protein sorting into HIV membranes, indicating that the choice of RNA for packaging can modulate viral membrane compositions. These findings should advance the understanding of HIV assembly and further benefit the development of virus-like particle-based therapeutics.

Formalizing Coarse-Grained Representations of Anisotropic Interactions at Multimeric Protein Interfaces Using Virtual Sites

Molecular simulations of biomacromolecules that assemble into multimeric complexes remain a challenge due to computationally inaccessible length and time scales. Low-resolution and implicit-solvent coarse-grained modeling approaches using traditional nonbonded interactions (both pairwise and spherically isotropic) have been able to partially address this gap. However, these models may fail to capture the complex anisotropic interactions present at macromolecular interfaces unless higher-order interaction potentials are incorporated at the expense of the computational cost. In this work, we introduce an alternate and systematic approach to represent directional interactions at protein-protein interfaces by using virtual sites restricted to pairwise interactions. We show that virtual site interaction parameters can be optimized within a relative entropy minimization framework by using only information from known statistics between coarse-grained sites. We compare our virtual site models to traditional coarse-grained models using two case studies of multimeric protein assemblies and find that the virtual site models predict pairwise correlations with higher fidelity and, more importantly, assembly behavior that is morphologically consistent with experiments. Our study underscores the importance of anisotropic interaction representations and paves the way for more accurate yet computationally efficient coarse-grained simulations of macromolecular assembly in future research.

Molecular dynamics simulations of HIV-1 matrix-membrane interactions at different stages of viral maturation

Although the structural rearrangement of the membrane-bound matrix (MA) protein trimers upon HIV-1 maturation has been reported, the consequences of MA maturation on the MA-lipid interactions are not well understood. Long-timescale molecular dynamics simulations of the MA multimeric assemblies of immature and mature virus particles with our realistic asymmetric membrane model have explored MA-lipid interactions and lateral organization of lipids around MA complexes. The number of stable MA-phosphatidylserine and MA-phosphatidylinositol 4,5-bisphosphate (PIP2) interactions at the trimeric interface of the mature MA complex is observed to be greater compared to that of the immature MA complex. Our simulations identified an alternative PIP2-binding site in the immature MA complex where the multivalent headgroup of a PIP2 lipid with a greater negative charge binds to multiple basic amino acid residues such as ARG3 residues of both the MA monomers at the trimeric interface and highly basic region (HBR) residues (LYS29, LYS31) of one of the MA monomers. Our enhanced sampling simulations have explored the conformational space of phospholipids at different binding sites of the trimer-trimer interface of MA complexes that are not accessible by conventional unbiased molecular dynamics. Unlike the immature MA complex, the 2' acyl tail of two PIP2 lipids at the trimeric interface of the mature MA complex is observed to sample stable binding pockets of MA consisting of helix-4 residues. Together, our results provide molecular-level insights into the interactions of MA trimeric complexes with membrane and different lipid conformations at the specific binding sites of MA protein before and after viral maturation.

HIV-1 capsid shape, orientation, and entropic elasticity regulate translocation into the nuclear pore complex

Nuclear import and uncoating of the viral capsid are critical steps in the HIV-1 life cycle that serve to transport and release genomic material into the nucleus. Viral core import involves translocating the HIV-1 capsid at the nuclear pore complex (NPC). Notably, the central channel of the NPC appears to often accommodate and allow passage of intact HIV-1 capsid, though mechanistic details of the process remain to be fully understood. Here, we investigate the molecular interactions that operate in concert between the HIV-1 capsid and the NPC that regulate capsid translocation through the central channel. To this end, we develop a "bottom-up" coarse-grained (CG) model of the human NPC from recently released cryo-electron tomography structure and then construct composite membrane-embedded CG NPC models. We find that successful translocation from the cytoplasmic side to the NPC central channel is contingent on the compatibility of the capsid morphology and channel dimension and the proper orientation of the capsid approach to the channel from the cytoplasmic side. The translocation dynamics is driven by maximizing the contacts between phenylalanine-glycine nucleoporins at the central channel and the capsid. For the docked intact capsids, structural analysis reveals correlated striated patterns of lattice disorder likely related to the intrinsic capsid elasticity. Uncondensed genomic material inside the docked capsid augments the overall lattice disorder of the capsid. Our results suggest that the intrinsic "elasticity" can also aid the capsid to adapt to the stress and remain structurally intact during translocation.

Computational Prediction of Coiled-Coil Protein Gelation Dynamics and Structure

Protein hydrogels represent an important and growing biomaterial for a multitude of applications, including diagnostics and drug delivery. We have previously explored the ability to engineer the thermoresponsive supramolecular assembly of coiled-coil proteins into hydrogels with varying gelation properties, where we have defined important parameters in the coiled-coil hydrogel design. Using Rosetta energy scores and Poisson-Boltzmann electrostatic energies, we iterate a computational design strategy to predict the gelation of coiled-coil proteins while simultaneously exploring five new coiled-coil protein hydrogel sequences. Provided this library, we explore the impact of in silico energies on structure and gelation kinetics, where we also reveal a range of blue autofluorescence that enables hydrogel disassembly and recovery. As a result of this library, we identify the new coiled-coil hydrogel sequence, Q5, capable of gelation within 24 h at 4 °C, a more than 2-fold increase over that of our previous iteration Q2. The fast gelation time of Q5 enables the assessment of structural transition in real time using small-angle X-ray scattering (SAXS) that is correlated to coarse-grained and atomistic molecular dynamics simulations revealing the supramolecular assembling behavior of coiled-coils toward nanofiber assembly and gelation. This work represents the first system of hydrogels with predictable self-assembly, autofluorescent capability, and a molecular model of coiled-coil fiber formation.

Conformational transitions of the HIV-1 Gag polyprotein upon multimerization and gRNA binding

During the HIV-1 assembly process, the Gag polyprotein multimerizes at the producer cell plasma membrane, resulting in the formation of spherical immature virus particles. Gag-genomic RNA (gRNA) interactions play a crucial role in the multimerization process, which is yet to be fully understood. We performed large-scale all-atom molecular dynamics simulations of membrane-bound full-length Gag dimer, hexamer, and 18-mer. The inter-domain dynamic correlation of Gag, quantified by the heterogeneous elastic network model applied to the simulated trajectories, is observed to be altered by implicit gRNA binding, as well as the multimerization state of the Gag. The lateral dynamics of our simulated membrane-bound Gag proteins, with and without gRNA binding, agree with prior experimental data and help to validate our simulation models and methods. The gRNA binding is observed to affect mainly the SP1 domain of the 18-mer and the matrix-capsid linker domain of the hexamer. In the absence of gRNA binding, the independent dynamical motion of the nucleocapsid domain results in a collapsed state of the dimeric Gag. Unlike stable SP1 helices in the six-helix bundle, without IP6 binding, the SP1 domain undergoes a spontaneous helix-to-coil transition in the dimeric Gag. Together, our findings reveal conformational switches of Gag at different stages of the multimerization process and predict that the gRNA binding reinforces an efficient binding surface of Gag for multimerization, and also regulates the dynamic organization of the local membrane region itself.

Molecular Sensing and Manipulation of Protein Oligomerization in Membrane Nanotubes with Bolaamphiphilic Foldamers

Adaptive and reversible self-assembly of supramolecular protein structures is a fundamental characteristic of dynamic living matter. However, the quantitative detection and assessment of the emergence of mesoscale protein complexes from small and dynamic oligomeric precursors remains highly challenging. Here, we present a novel approach utilizing a short membrane nanotube (sNT) pulled from a planar membrane reservoir as nanotemplates for molecular reconstruction, manipulation, and sensing of protein oligomerization and self-assembly at the mesoscale. The sNT reports changes in membrane shape and rigidity caused by membrane-bound proteins as variations of the ionic conductivity of the sNT lumen. To confine oligomerization to the sNT, we have designed and synthesized rigid oligoamide foldamer tapes (ROFTs). Charged ROFTs incorporate into the planar and sNT membranes, mediate protein binding to the membranes, and, driven by the luminal electric field, shuttle the bound proteins between the sNT and planar membranes. Using Annexin-V (AnV) as a prototype, we show that the sNT detects AnV oligomers shuttled into the nanotube by ROFTs. Accumulation of AnV on the sNT induces its self-assembly into a curved lattice, restricting the sNT geometry and inhibiting the material uptake from the reservoir during the sNT extension, leading to the sNT fission. By comparing the spontaneous and ROFT-mediated entry of AnV into the sNT, we reveal how intricate membrane curvature sensing by small AnV oligomers controls the lattice self-assembly. These results establish sNT-ROFT as a powerful tool for molecular reconstruction and functional analyses of protein oligomerization and self-assembly, with broad application to various membrane processes.

Gaussian representation of coarse-grained interactions of liquids: Theory, parametrization, and transferability

Coarse-grained (CG) interactions determined via bottom-up methodologies can faithfully reproduce the structural correlations observed in fine-grained (atomistic resolution) systems, yet they can suffer from limited extensibility due to complex many-body correlations. As part of an ongoing effort to understand and improve the applicability of bottom-up CG models, we propose an alternative approach to address both accuracy and transferability. Our main idea draws from classical perturbation theory to partition the hard sphere repulsive term from effective CG interactions. We then introduce Gaussian basis functions corresponding to the system's characteristic length by linking these Gaussian sub-interactions to the local particle densities at each coordination shell. The remaining perturbative long-range interaction can be treated as a collective solvation interaction, which we show exhibits a Gaussian form derived from integral equation theories. By applying this numerical parametrization protocol to CG liquid systems, our microscopic theory elucidates the emergence of Gaussian interactions in common phenomenological CG models. To facilitate transferability for these reduced descriptions, we further infer equations of state to determine the sub-interaction parameter as a function of the system variables. The reduced models exhibit excellent transferability across the thermodynamic state points. Furthermore, we propose a new strategy to design the cross-interactions between distinct CG sites in liquid mixtures. This involves combining each Gaussian in the proper radial domain, yielding accurate CG potentials of mean force and structural correlations for multi-component systems. Overall, our findings establish a solid foundation for constructing transferable bottom-up CG models of liquids with enhanced extensibility.

Effect of ionic strength on the assembly of simian vacuolating virus capsid protein around poly(styrene sulfonate)

Virus-like particles (VLPs) are noninfectious nanocapsules that can be used for drug delivery or vaccine applications. VLPs can be assembled from virus capsid proteins around a condensing agent, such as RNA, DNA, or a charged polymer. Electrostatic interactions play an important role in the assembly reaction. VLPs assemble from many copies of capsid protein, with a combinatorial number of intermediates. Hence, the mechanism of the reaction is poorly understood. In this paper, we combined solution small-angle X-ray scattering (SAXS), cryo-transmission electron microscopy (TEM), and computational modeling to determine the effect of ionic strength on the assembly of Simian Vacuolating Virus 40 (SV40)-like particles. We mixed poly(styrene sulfonate) with SV40 capsid protein pentamers at different ionic strengths. We then characterized the assembly product by SAXS and cryo-TEM. To analyze the data, we performed Langevin dynamics simulations using a coarse-grained model that revealed incomplete, asymmetric VLP structures consistent with the experimental data. We found that close to physiological ionic strength, [Formula: see text] VLPs coexisted with VP1 pentamers. At lower or higher ionic strengths, incomplete particles coexisted with pentamers and [Formula: see text] particles. Including the simulated structures was essential to explain the SAXS data in a manner that is consistent with the cryo-TEM images.

OpenMSCG: A Software Tool for Bottom-Up Coarse-Graining

The "bottom-up" approach to coarse-graining, for building accurate and efficient computational models to simulate large-scale and complex phenomena and processes, is an important approach in computational chemistry, biophysics, and materials science. As one example, the Multiscale Coarse-Graining (MS-CG) approach to developing CG models can be rigorously derived using statistical mechanics applied to fine-grained, i.e., all-atom simulation data for a given system. Under a number of circumstances, a systematic procedure, such as MS-CG modeling, is particularly valuable. Here, we present the development of the OpenMSCG software, a modularized open-source software that provides a collection of successful and widely applied bottom-up CG methods, including Boltzmann Inversion (BI), Force-Matching (FM), Ultra-Coarse-Graining (UCG), Relative Entropy Minimization (REM), Essential Dynamics Coarse-Graining (EDCG), and Heterogeneous Elastic Network Modeling (HeteroENM). OpenMSCG is a high-performance and comprehensive toolset that can be used to derive CG models from large-scale fine-grained simulation data in file formats from common molecular dynamics (MD) software packages, such as GROMACS, LAMMPS, and NAMD. OpenMSCG is modularized in the Python programming framework, which allows users to create and customize modeling "recipes" for reproducible results, thus greatly improving the reliability, reproducibility, and sharing of bottom-up CG models and their applications.

Cooperative Membrane Binding of HIV-1 Matrix Proteins

The HIV-1 assembly process begins with a newly synthesized Gag polyprotein being targeted to the inner leaflet of the plasma membrane of the infected cells to form immature viral particles. Gag-membrane interactions are mediated through the myristoylated(Myr) N-terminal matrix (MA) domain of Gag which eventually multimerize on the membrane to form trimers and higher-order oligomers. The study of the structure and dynamics of peripheral membrane proteins like MA has been challenging for both experimental and computational studies due to the complex dynamics of protein-membrane interactions. Although the roles of anionic phospholipids (PIP2, PS) and the Myr group in the membrane targeting and stable membrane binding of MA are now well-established, the cooperative interactions between MA monomers and MA-membrane still remain elusive. Our present study focuses on the membrane binding dynamics of a higher-order oligomeric structure of MA protein (a dimer of trimers), which has not been explored before. Employing time-lagged independent component analysis (tICA) to our microsecond-long trajectories, we investigate conformational changes of the matrix protein induced by membrane binding. Interestingly, the Myr switch of a MA monomer correlates with the conformational switch of adjacent monomers in the same trimer. Together, our findings suggest that MA trimerization facilitates Myr insertion, but MA trimer-trimer interactions in the lattice of immature HIV-1 particles can hinder the same. Additionally, local lipid density patterns of different lipid species provide a signature of the initial stage of lipid-domain formation upon membrane binding of the protein complex.

TOC

Conformational transitions of the HIV-1 Gag polyprotein upon multimerization and gRNA binding

During the HIV-1 assembly process, the Gag polyprotein multimerizes at the producer cell plasma membrane, resulting in the formation of spherical immature virus particles. Gag-gRNA interactions play a crucial role in the multimerization process, which is yet to be fully understood. We have performed large-scale all-atom molecular dynamics simulations of membrane-bound full-length Gag dimer, hexamer, and 18-mer. The inter-domain dynamic correlation of Gag, quantified by the heterogeneous elastic network model (hENM) applied to the simulated trajectories, is observed to be altered by implicit gRNA binding, as well as the multimerization state of the Gag. The lateral dynamics of our simulated membrane-bound Gag proteins, with and without gRNA binding, agree with prior experimental data and help to validate our simulation models and methods. The gRNA binding is observed to impact mainly the SP1 domain of the 18-mer and the MA-CA linker domain of the hexamer. In the absence of gRNA binding, the independent dynamical motion of the NC domain results in a collapsed state of the dimeric Gag. Unlike stable SP1 helices in the six-helix bundle, without IP6 binding, the SP1 domain undergoes a spontaneous helix-to-coil transition in the dimeric Gag. Together, our findings reveal conformational switches of Gag at different stages of the multimerization process and predict that the gRNA binding reinforces an efficient binding surface of Gag for multimerization, as well as regulates the dynamic organization of the local membrane region itself.

Significance

Gag(Pr 55 Gag ) polyprotein orchestrates many essential events in HIV-1 assembly, including packaging of the genomic RNA (gRNA) in the immature virion. Although various experimental techniques, such as cryo-ET, X-ray, and NMR, have revealed structural properties of individual domains in the immature Gag clusters, structural and biophysical characterization of a full-length Gag molecule remains a challenge for existing experimental techniques. Using atomistic molecular dynamics simulations of the different model systems of Gag polyprotein, we present here a detailed structural characterization of Gag molecules in different multimerization states and interrogate the synergy between Gag-Gag, Gag-membrane, and Gag-gRNA interactions during the viral assembly process.

Temporal control by cofactors prevents kinetic trapping in retroviral Gag lattice assembly

For retroviruses like HIV to proliferate, they must form virions shaped by the self-assembly of Gag polyproteins into a rigid lattice. This immature Gag lattice has been structurally characterized and reconstituted in vitro, revealing the sensitivity of lattice assembly to multiple cofactors. Due to this sensitivity, the energetic criterion for forming stable lattices is unknown, as are their corresponding rates. Here, we use a reaction-diffusion model designed from the cryo-ET structure of the immature Gag lattice to map a phase diagram of assembly outcomes controlled by experimentally constrained rates and free energies, over experimentally relevant timescales. We find that productive assembly of complete lattices in bulk solution is extraordinarily difficult due to the large size of this ∼3700 monomer complex. Multiple Gag lattices nucleate before growth can complete, resulting in loss of free monomers and frequent kinetic trapping. We therefore derive a time-dependent protocol to titrate or "activate" the Gag monomers slowly within the solution volume, mimicking the biological roles of cofactors. This general strategy works remarkably well, yielding productive growth of self-assembled lattices for multiple interaction strengths and binding rates. By comparing to the in vitro assembly kinetics, we can estimate bounds on rates of Gag binding to Gag and the cellular cofactor IP6. Our results show that Gag binding to IP6 can provide the additional time delay necessary to support smooth growth of the immature lattice with relatively fast assembly kinetics, mostly avoiding kinetic traps. Our work provides a foundation for predicting and disrupting formation of the immature Gag lattice via targeting specific protein-protein binding interactions.

Structure of the HIV immature lattice allows for essential lattice remodeling within budded virions

For HIV virions to become infectious, the immature lattice of Gag polyproteins attached to the virion membrane must be cleaved. Cleavage cannot initiate without the protease formed by the homo-dimerization of domains linked to Gag. However, only 5% of the Gag polyproteins, termed Gag-Pol, carry this protease domain, and they are embedded within the structured lattice. The mechanism of Gag-Pol dimerization is unknown. Here, we use spatial stochastic computer simulations of the immature Gag lattice as derived from experimental structures, showing that dynamics of the lattice on the membrane is unavoidable due to the missing 1/3 of the spherical protein coat. These dynamics allow for Gag-Pol molecules carrying the protease domains to detach and reattach at new places within the lattice. Surprisingly, dimerization timescales of minutes or less are achievable for realistic binding energies and rates despite retaining most of the large-scale lattice structure. We derive a formula allowing extrapolation of timescales as a function of interaction free energy and binding rate, thus predicting how additional stabilization of the lattice would impact dimerization times. We further show that during assembly, dimerization of Gag-Pol is highly likely and therefore must be actively suppressed to prevent early activation. By direct comparison to recent biochemical measurements within budded virions, we find that only moderately stable hexamer contacts (-12kBT<∆G<-8kBT) retain both the dynamics and lattice structures that are consistent with experiment. These dynamics are likely essential for proper maturation, and our models quantify and predict lattice dynamics and protease dimerization timescales that define a key step in understanding formation of infectious viruses.

Dynamics of upstream ESCRT organization at the HIV-1 budding site

In the late stages of the HIV-1 life cycle, membrane localization and self-assembly of Gag polyproteins induce membrane deformation and budding. Release of the virion requires direct interaction between immature Gag lattice and upstream ESCRT machinery at the viral budding site, followed by assembly of downstream ESCRT-III factors, culminating in membrane scission. However, molecular details of upstream ESCRT assembly dynamics at the viral budding site remain unclear. In this work, using coarse-grained (CG) molecular dynamics (MD) simulations, we investigated the interactions between Gag, ESCRT-I, ESCRT-II, and membrane to delineate the dynamical mechanisms by which upstream ESCRTs assemble templated by late-stage immature Gag lattice. We first systematically derived "bottom-up" CG molecular models and interactions of upstream ESCRT proteins from experimental structural data and extensive all-atom MD simulations. Using these molecular models, we performed CG MD simulations of ESCRT-I oligomerization and ESCRT-I/II supercomplex formation at the neck of the budding virion. Our simulations demonstrate that ESCRT-I can effectively oligomerize to higher-order complexes templated by the immature Gag lattice both in the absence of ESCRT-II and when multiple copies of ESCRT-II are localized at the bud neck. The ESCRT-I/II supercomplexes formed in our simulations exhibit predominantly columnar structures, which has important implications for the nucleation pathway of downstream ESCRT-III polymers. Importantly, ESCRT-I/II supercomplexes bound to Gag initiate membrane neck constriction by pulling the inner edge of the bud neck closer to the ESCRT-I headpiece ring. Our findings serve to elucidate a network of interactions between upstream ESCRT machinery, immature Gag lattice, and membrane neck that regulate protein assembly dynamics at the HIV-1 budding site.

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.

Understanding dynamics in coarse-grained models. I. Universal excess entropy scaling relationship

Coarse-grained (CG) models facilitate an efficient exploration of complex systems by reducing the unnecessary degrees of freedom of the fine-grained (FG) system while recapitulating major structural correlations. Unlike structural properties, assessing dynamic properties in CG modeling is often unfeasible due to the accelerated dynamics of the CG models, which allows for more efficient structural sampling. Therefore, the ultimate goal of the present series of articles is to establish a better correspondence between the FG and CG dynamics. To assess and compare dynamical properties in the FG and the corresponding CG models, we utilize the excess entropy scaling relationship. For Paper I of this series, we provide evidence that the FG and the corresponding CG counterpart follow the same universal scaling relationship. By carefully reviewing and examining the literature, we develop a new theory to calculate excess entropies for the FG and CG systems while accounting for entropy representability. We demonstrate that the excess entropy scaling idea can be readily applied to liquid water and methanol systems at both the FG and CG resolutions. For both liquids, we reveal that the scaling exponents remain unchanged from the coarse-graining process, indicating that the scaling behavior is universal for the same underlying molecular systems. Combining this finding with the concept of mapping entropy in CG models, we show that the missing entropy plays an important role in accelerating the CG dynamics.

The dawn of interoperating spatial models in cell biology

Spatial simulations are becoming an increasingly ubiquitous component in the cycle of discovery, experimentation, and communication across the sciences. In cell biology, many researchers share a vision of developing multiscale models that recapitulate observable behaviors spanning from atoms to cells to tissues. For this dream to become a reality, however, simulation technologies must provide a means for integration and interoperability as they advance. Already, the field has developed numerous methods that span scales of length, time, and complexity to create an extensive body of effective simulation approaches, and although these approaches rarely interoperate, they collectively cover a large spectrum of knowledge that future models may handle in a more unified manner. Here, we discuss the importance of making the data, workflows, and outputs of spatial simulations shareable and interoperable; and how democratization could encourage diverse biologists to participate more easily in developing models to advance our understanding of biological systems.

Bottom-up Coarse-Graining: Principles and Perspectives

Large-scale computational molecular models provide scientists a means to investigate the effect of microscopic details on emergent mesoscopic behavior. Elucidating the relationship between variations on the molecular scale and macroscopic observable properties facilitates an understanding of the molecular interactions driving the properties of real world materials and complex systems (e.g., those found in biology, chemistry, and materials science). As a result, discovering an explicit, systematic connection between microscopic nature and emergent mesoscopic behavior is a fundamental goal for this type of investigation. The molecular forces critical to driving the behavior of complex heterogeneous systems are often unclear. More problematically, simulations of representative model systems are often prohibitively expensive from both spatial and temporal perspectives, impeding straightforward investigations over possible hypotheses characterizing molecular behavior. While the reduction in resolution of a study, such as moving from an atomistic simulation to that of the resolution of large coarse-grained (CG) groups of atoms, can partially ameliorate the cost of individual simulations, the relationship between the proposed microscopic details and this intermediate resolution is nontrivial and presents new obstacles to study. Small portions of these complex systems can be realistically simulated. Alone, these smaller simulations likely do not provide insight into collectively emergent behavior. However, by proposing that the driving forces in both smaller and larger systems (containing many related copies of the smaller system) have an explicit connection, systematic bottom-up CG techniques can be used to transfer CG hypotheses discovered using a smaller scale system to a larger system of primary interest. The proposed connection between different CG systems is prescribed by (i) the CG representation (mapping) and (ii) the functional form and parameters used to represent the CG energetics, which approximate potentials of mean force (PMFs). As a result, the design of CG methods that facilitate a variety of physically relevant representations, approximations, and force fields is critical to moving the frontier of systematic CG forward. Crucially, the proposed connection between the system used for parametrization and the system of interest is orthogonal to the optimization used to approximate the potential of mean force present in all systematic CG methods. The empirical efficacy of machine learning techniques on a variety of tasks provides strong motivation to consider these approaches for approximating the PMF and analyzing these approximations.

Inositol Hexakisphosphate (IP6) Accelerates Immature HIV-1 Gag Protein Assembly toward Kinetically Trapped Morphologies

During the late stages of the HIV-1 lifecycle, immature virions are produced by the concerted activity of Gag polyproteins, primarily mediated by the capsid (CA) and spacer peptide 1 (SP1) domains, which assemble into a spherical lattice, package viral genomic RNA, and deform the plasma membrane. Recently, inositol hexakisphosphate (IP6) has been identified as an essential assembly cofactor that efficiently produces both immature virions in vivo and immature virus-like particles in vitro. To date, however, several distinct mechanistic roles for IP6 have been proposed on the basis of independent functional, structural, and kinetic studies. In this work, we investigate the molecular influence of IP6 on the structural outcomes and dynamics of CA/SP1 assembly using coarse-grained (CG) molecular dynamics (MD) simulations and free energy calculations. Here, we derive a bottom-up, low-resolution, and implicit-solvent CG model of CA/SP1 and IP6, and simulate their assembly under conditions that emulate both in vitro and in vivo systems. Our analysis identifies IP6 as an assembly accelerant that promotes curvature generation and fissure-like defects throughout the lattice. Our findings suggest that IP6 induces kinetically trapped immature morphologies, which may be physiologically important for later stages of viral morphogenesis and potentially useful for virus-like particle technologies.

Molecular Shape Solution for Mesoscopic Remodeling of Cellular Membranes

Cellular membranes self-assemble from and interact with various molecular species. Each molecule locally shapes the lipid bilayer, the soft elastic core of cellular membranes. The dynamic architecture of intracellular membrane systems is based on elastic transformations and lateral redistribution of these elementary shapes, driven by chemical and curvature stress gradients. The minimization of the total elastic stress by such redistribution composes the most basic, primordial mechanism of membrane curvature-composition coupling (CCC). Although CCC is generally considered in the context of dynamic compositional heterogeneity of cellular membrane systems, in this article we discuss a broader involvement of CCC in controlling membrane deformations. We focus specifically on the mesoscale membrane transformations in open, reservoir-governed systems, such as membrane budding, tubulation, and the emergence of highly curved sites of membrane fusion and fission. We reveal that the reshuffling of molecular shapes constitutes an independent deformation mode with complex rheological properties.This mode controls effective elasticity of local deformations as well as stationary elastic stress, thus emerging as a major regulator of intracellular membrane remodeling.

Relationship between HIV-1 Gag Multimerization and Membrane Binding

HIV-1 viral particle assembly occurs specifically at the plasma membrane and is driven primarily by the viral polyprotein Gag. Selective association of Gag with the plasma membrane is a key step in the viral assembly pathway, which is traditionally attributed to the MA domain. MA regulates specific plasma membrane binding through two primary mechanisms including: (1) specific interaction of the MA highly basic region (HBR) with the plasma membrane phospholipid phosphatidylinositol (4,5) bisphosphate [PI(4,5)P2], and (2) tRNA binding to the MA HBR, which prevents Gag association with non-PI(4,5)P2 containing membranes. Gag multimerization, driven by both CA–CA inter-protein interactions and NC-RNA binding, also plays an essential role in viral particle assembly, mediating the establishment and growth of the immature Gag lattice on the plasma membrane. In addition to these functions, the multimerization of HIV-1 Gag has also been demonstrated to enhance its membrane binding activity through the MA domain. This review provides an overview of the mechanisms regulating Gag membrane binding through the MA domain and multimerization through the CA and NC domains, and examines how these two functions are intertwined, allowing for multimerization mediated enhancement of Gag membrane binding.

Approach to map nanotopography of cell surface receptors

Cells communicate with their environment via surface receptors, but nanoscopic receptor organization with respect to complex cell surface morphology remains unclear. This is mainly due to a lack of accessible, robust and high-resolution methods. Here, we present an approach for mapping the topography of receptors at the cell surface with nanometer precision. The method involves coating glass coverslips with glycine, which preserves the fine membrane morphology while allowing immobilized cells to be positioned close to the optical surface. We developed an advanced and simplified algorithm for the analysis of single-molecule localization data acquired in a biplane detection scheme. These advancements enable direct and quantitative mapping of protein distribution on ruffled plasma membranes with near isotropic 3D nanometer resolution. As demonstrated successfully for CD4 and CD45 receptors, the described workflow is a straightforward quantitative technique to study molecules and their interactions at the complex surface nanomorphology of differentiated metazoan cells.

A super-resolution localisation-based method is shown to map receptor topography in immune cells, which allows quantitative study of molecules and their interactions at the complex surface nanomorphology of cells.

Molecular Determinants in tRNA D-arm Required for Inhibition of HIV-1 Gag Membrane Binding

Plasma-membrane-specific localization of Gag, an essential step in HIV-1 particle assembly, is regulated by the interaction of the Gag MA domain with PI(4,5)P₂ and tRNA-mediated inhibition of non-specific or premature membrane binding. Different tRNAs inhibit PI(4,5)P₂-independent membrane binding to varying degrees in vitro; however, the structural determinants for this difference remain unknown. Here we demonstrate that membrane binding of full-length Gag synthesized in vitro using reticulocyte lysates is inhibited when RNAs that contain the anticodon arm of tRNA^Pro, but not that of tRNA^Lys3, are added exogenously. In contrast, in the context of a liposome binding assay in which the effects of tRNAs on purified MA were tested, full-length tRNA^Lys3 showed greater inhibition of MA membrane binding than full-length tRNA^Pro. While transplantation of the D loop sequence of tRNA^Lys3 into tRNA^Pro resulted in a modest increase in the inhibitory effect relative to WT tRNA^Pro, replacing the entire D arm sequence with that of tRNA^Lys3 was necessary to confer the full inhibitory effects upon tRNA^Pro. Together, these results demonstrate that the D arm of tRNA^Lys3 is a major determinant of strong inhibition of MA membrane binding and that this inhibitory effect requires not only the D loop, which was recently reported to contact the MA highly basic region, but the loop sequence in the context of the D arm structure.

Rational design of self-assembled RNA nanostructures for HIV-1 virus assembly blockade

Many pathological processes are driven by RNA-protein interactions, making such interactions promising targets for molecular interventions. HIV-1 assembly is one such process, in which the viral genomic RNA interacts with the viral Gag protein and serves as a scaffold to drive Gag multimerization that ultimately leads to formation of a virus particle. Here, we develop self-assembled RNA nanostructures that can inhibit HIV-1 virus assembly, achieved through hybridization of multiple artificial small RNAs with a stem-loop structure (STL) that we identify as a prominent ligand of Gag that can inhibit virus particle production via STL-Gag interactions. The resulting STL-decorated nanostructures (double and triple stem-loop structures denoted as Dumbbell and Tribell, respectively) can elicit more pronounced viral blockade than their building blocks, with the inhibition arising as a result of nanostructures interfering with Gag multimerization. These findings could open up new avenues for RNA-based therapy.

Preservation of HIV-1 Gag Helical Bundle Symmetry by Bevirimat Is Central to Maturation Inhibition

The assembly and maturation of human immunodeficiency virus type 1 (HIV-1) require proteolytic cleavage of the Gag polyprotein. The rate-limiting step resides at the junction between the capsid protein CA and spacer peptide 1, which assembles as a six-helix bundle (6HB). Bevirimat (BVM), the first-in-class maturation inhibitor drug, targets the 6HB and impedes proteolytic cleavage, yet the molecular mechanisms of its activity, and relatedly, the escape mechanisms of mutant viruses, remain unclear. Here, we employed extensive molecular dynamics (MD) simulations and free energy calculations to quantitatively investigate molecular structure-activity relationships, comparing wild-type and mutant viruses in the presence and absence of BVM and inositol hexakisphosphate (IP6), an assembly cofactor. Our analysis shows that the efficacy of BVM is directly correlated with preservation of 6-fold symmetry in the 6HB, which exists as an ensemble of structural states. We identified two primary escape mechanisms, and both lead to loss of symmetry, thereby facilitating helix uncoiling to aid access of protease. Our findings also highlight specific interactions that can be targeted for improved inhibitor activity and support the use of MD simulations for future inhibitor design.

Betacoronavirus Assembly: Clues and Perspectives for Elucidating SARS-CoV-2 Particle Formation and Egress

In 2019, a new pandemic virus belonging to the betacoronavirus family emerged, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). This new coronavirus appeared in Wuhan, China, and is responsible for severe respiratory pneumonia in humans, namely, coronavirus disease 2019 (COVID-19). Having infected almost 200 million people worldwide and caused more than 4.1 million deaths as of today, this new disease has raised a significant number of questions about its molecular mechanism of replication and, in particular, how infectious viral particles are produced. Although viral entry is well characterized, the full assembly steps of SARS-CoV-2 have still not been fully described. Coronaviruses, including SARS-CoV-2, have four main structural proteins, namely, the spike glycoprotein (S), the membrane glycoprotein (M), the envelope protein (E), and the nucleocapsid protein (N). All these proteins have key roles in the process of coronavirus assembly and budding. In this review, we gathered the current knowledge about betacoronavirus structural proteins involved in viral particle assembly, membrane curvature and scission, and then egress in order to suggest and question a coherent model for SARS-CoV-2 particle production and release.

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.

A Background Assessable and Correctable Bimolecular Fluorescence Complementation System for Nanoscopic Single-Molecule Imaging of Intracellular Protein-Protein Interactions

Bimolecular Fluorescence Complementation (BiFC) is a versatile approach for intracellular analysis of protein-protein interactions (PPIs), but the tendency of the split fluorescent protein (FP) fragments to self-assemble when brought into close proximity of each other by random collision can lead to generation of false-positive signals that hamper high-definition imaging of PPIs occurring on the nanoscopic level. While it is thought that expressing the fusion proteins at a low level can remove false positives without impacting specific signals, there has been no effective strategy to test this possibility. Here, we present a system capable of assessing and removing BiFC false positives, termed Background Assessable and Correctable-BiFC (BAC-BiFC), in which one of the split FP fragments is fused with an optically distinct FP that serves as a reference marker, and the single-cell fluorescence ratio of the BiFC signal to the reference signal is used to gauge an optimal transfection condition. We showed that when BAC-BiFC is designed to image PPIs regulating Human Immunodeficiency Virus type 1 (HIV-1) assembly, the fluorescence ratio could decrease with decreasing probe quantity, and ratios approaching the limit of detection could allow physiologically relevant characterization of the assembly process on the nanoscale by single-molecule localization microscopy (SMLM). With much improved clarity, previously undescribed features of HIV-1 assembly were revealed.

When in Need of an ESCRT: The Nature of Virus Assembly Sites Suggests Mechanistic Parallels between Nuclear Virus Egress and Retroviral Budding

The proper assembly and dissemination of progeny virions is a fundamental step in virus replication. As a whole, viruses have evolved a myriad of strategies to exploit cellular compartments and mechanisms to ensure a successful round of infection. For enveloped viruses such as retroviruses and herpesviruses, acquisition and incorporation of cellular membrane is an essential process during the formation of infectious viral particles. To do this, these viruses have evolved to hijack the host Endosomal Sorting Complexes Required for Transport (ESCRT-I, -II, and -III) to coordinate the sculpting of cellular membrane at virus assembly and dissemination sites, in seemingly different, yet fundamentally similar ways. For instance, at the plasma membrane, ESCRT-I recruitment is essential for HIV-1 assembly and budding, while it is dispensable for the release of HSV-1. Further, HSV-1 was shown to recruit ESCRT-III for nuclear particle assembly and egress, a process not used by retroviruses during replication. Although the cooption of ESCRTs occurs in two separate subcellular compartments and at two distinct steps for these viral lifecycles, the role fulfilled by ESCRTs at these sites appears to be conserved. This review discusses recent findings that shed some light on the potential parallels between retroviral budding and nuclear egress and proposes a model where HSV-1 nuclear egress may occur through an ESCRT-dependent mechanism.

Full assembly of HIV-1 particles requires assistance of the membrane curvature factor IRSp53

During HIV-1 particle formation, the requisite plasma membrane curvature is thought to be solely driven by the retroviral Gag protein. Here, we reveal that the cellular I-BAR protein IRSp53 is required for the progression of HIV-1 membrane curvature to complete particle assembly. siRNA-mediated knockdown of IRSp53 gene expression induces a decrease in viral particle production and a viral bud arrest at half completion. Single-molecule localization microscopy at the cell plasma membrane shows a preferential localization of IRSp53 around HIV-1 Gag assembly sites. In addition, we observe the presence of IRSp53 in purified HIV-1 particles. Finally, HIV-1 Gag protein preferentially localizes to curved membranes induced by IRSp53 I-BAR domain on giant unilamellar vesicles. Overall, our data reveal a strong interplay between IRSp53 I-BAR and Gag at membranes during virus assembly. This highlights IRSp53 as a crucial host factor in HIV-1 membrane curvature and its requirement for full HIV-1 particle assembly.

Immature HIV-1 assembles from Gag dimers leaving partial hexamers at lattice edges as potential substrates for proteolytic maturation

HIV-1 particle assembly is driven by the viral Gag protein, which oligomerizes into a hexameric array on the inner surface of the viral envelope, forming a truncated spherical lattice containing large and small gaps. Gag is then cut by the viral protease, disassembles, and rearranges to form the mature, infectious virus. Here, we present structures and molecular dynamics simulations of the edges of the immature Gag lattice. Our analysis shows that Gag dimers are the basic assembly unit of the HIV-1 particle, lattice edges are partial hexamers, and partial hexamers are prone to structural changes allowing protease to cut Gag. These findings provide insights into assembly of the immature virus, its structure, and how it disassembles during maturation.

The CA (capsid) domain of immature HIV-1 Gag and the adjacent spacer peptide 1 (SP1) play a key role in viral assembly by forming a lattice of CA hexamers, which adapts to viral envelope curvature by incorporating small lattice defects and a large gap at the site of budding. This lattice is stabilized by intrahexameric and interhexameric CA-CA interactions, which are important in regulating viral assembly and maturation. We applied subtomogram averaging and classification to determine the oligomerization state of CA at lattice edges and found that CA forms partial hexamers. These structures reveal the network of interactions formed by CA-SP1 at the lattice edge. We also performed atomistic molecular dynamics simulations of CA-CA interactions stabilizing the immature lattice and partial CA-SP1 helical bundles. Free energy calculations reveal increased propensity for helix-to-coil transitions in partial hexamers compared to complete six-helix bundles. Taken together, these results suggest that the CA dimer is the basic unit of lattice assembly, partial hexamers exist at lattice edges, these are in a helix-coil dynamic equilibrium, and partial helical bundles are more likely to unfold, representing potential sites for HIV-1 maturation initiation.

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.

A new one-site coarse-grained model for water: Bottom-up many-body projected water (BUMPer). II. Temperature transferability and structural properties at low temperature

A number of studies have constructed coarse-grained (CG) models of water to understand its anomalous properties. Most of these properties emerge at low temperatures, and an accurate CG model needs to be applicable to these low-temperature ranges. However, direct use of CG models parameterized from other temperatures, e.g., room temperature, encounters a problem known as transferability, as the CG potential essentially follows the form of the many-body CG free energy function. Therefore, temperature-dependent changes to CG interactions must be accounted for. The collective behavior of water at low temperature is generally a many-body process, which often motivates the use of expensive many-body terms in the CG interactions. To surmount the aforementioned problems, we apply the Bottom-Up Many-Body Projected Water (BUMPer) CG model constructed from Paper I to study the low-temperature behavior of water. We report for the first time that the embedded three-body interaction enables BUMPer, despite its pairwise form, to capture the growth of ice at the ice/water interface with corroborating many-body correlations during the crystal growth. Furthermore, we propose temperature transferable BUMPer models that are indirectly constructed from the free energy decomposition scheme. Changes in CG interactions and corresponding structures are faithfully recapitulated by this framework. We further extend BUMPer to examine its ability to predict the structure, density, and diffusion anomalies by employing an alternative analysis based on structural correlations and pairwise potential forms to predict such anomalies. The presented analysis highlights the existence of these anomalies in the low-temperature regime and overcomes potential transferability problems.

Structural characterization of HIV-1 matrix mutants implicated in envelope incorporation

During the late phase of HIV-1 infection, viral Gag polyproteins are targeted to the plasma membrane (PM) for assembly. Gag localization at the PM is a prerequisite for the incorporation of the envelope protein (Env) into budding particles. Gag assembly and Env incorporation are mediated by the N-terminal myristoylated matrix (MA) domain of Gag. Nonconservative mutations in the trimer interface of MA (A45E, T70R, and L75G) were found to impair Env incorporation and infectivity, leading to the hypothesis that MA trimerization is an obligatory step for Env incorporation. Conversely, Env incorporation can be rescued by a compensatory mutation in the MA trimer interface (Q63R). The impact of these MA mutations on the structure and trimerization properties of MA is not known. In this study, we employed NMR spectroscopy, X-ray crystallography, and sedimentation techniques to characterize the structure and trimerization properties of HIV-1 MA A45E, Q63R, T70R, and L75G mutant proteins. NMR data revealed that these point mutations did not alter the overall structure and folding of MA but caused minor structural perturbations in the trimer interface. Analytical ultracentrifugation data indicated that mutations had a minimal effect on the MA monomer–trimer equilibrium. The high-resolution X-ray structure of the unmyristoylated MA Q63R protein revealed hydrogen bonding between the side chains of adjacent Arg-63 and Ser-67 on neighboring MA molecules, providing the first structural evidence for an additional intermolecular interaction in the trimer interface. These findings advance our knowledge of the interplay of MA trimerization and Env incorporation into HIV-1 particles.

Single-molecule localization microscopy

Single-molecule localization microscopy (SMLM) describes a family of powerful imaging techniques that dramatically improve spatial resolution over standard, diffraction-limited microscopy techniques and can image biological structures at the molecular scale. In SMLM, individual fluorescent molecules are computationally localized from diffraction-limited image sequences and the localizations are used to generate a super-resolution image or a time course of super-resolution images, or to define molecular trajectories. In this Primer, we introduce the basic principles of SMLM techniques before describing the main experimental considerations when performing SMLM, including fluorescent labelling, sample preparation, hardware requirements and image acquisition in fixed and live cells. We then explain how low-resolution image sequences are computationally processed to reconstruct super-resolution images and/or extract quantitative information, and highlight a selection of biological discoveries enabled by SMLM and closely related methods. We discuss some of the main limitations and potential artefacts of SMLM, as well as ways to alleviate them. Finally, we present an outlook on advanced techniques and promising new developments in the fast-evolving field of SMLM. We hope that this Primer will be a useful reference for both newcomers and practitioners of SMLM.

HIV-1 Gag release from yeast reveals ESCRT interaction with the Gag N-terminal protein region

The HIV-1 protein Gag assembles at the plasma membrane and drives virion budding, assisted by the cellular endosomal complex required for transport (ESCRT) proteins. Two ESCRT proteins, TSG101 and ALIX, bind to the Gag C-terminal p6 peptide. TSG101 binding is important for efficient HIV-1 release, but how ESCRTs contribute to the budding process and how their activity is coordinated with Gag assembly is poorly understood. Yeast, allowing genetic manipulation that is not easily available in human cells, has been used to characterize the cellular ESCRT function. Previous work reported Gag budding from yeast spheroplasts, but Gag release was ESCRT-independent. We developed a yeast model for ESCRT-dependent Gag release. We combined yeast genetics and Gag mutational analysis with Gag-ESCRT binding studies and the characterization of Gag-plasma membrane binding and Gag release. With our system, we identified a previously unknown interaction between ESCRT proteins and the Gag N-terminal protein region. Mutations in the Gag-plasma membrane-binding matrix domain that reduced Gag-ESCRT binding increased Gag-plasma membrane binding and Gag release. ESCRT knockout mutants showed that the release enhancement was an ESCRT-dependent effect. Similarly, matrix mutation enhanced Gag release from human HEK293 cells. Release enhancement partly depended on ALIX binding to p6, although binding site mutation did not impair WT Gag release. Accordingly, the relative affinity for matrix compared with p6 in GST-pulldown experiments was higher for ALIX than for TSG101. We suggest that a transient matrix-ESCRT interaction is replaced when Gag binds to the plasma membrane. This step may activate ESCRT proteins and thereby coordinate ESCRT function with virion assembly.

The impact of Gag non-cleavage site mutations on HIV-1 viral fitness from integrative modelling and simulations

The high mutation rate in retroviruses is one of the leading causes of drug resistance. In human immunodeficiency virus type-1 (HIV-1), synergistic mutations in its protease and the protease substrate - the Group-specific antigen (Gag) polyprotein - work together to confer drug resistance against protease inhibitors and compensate the mutations affecting viral fitness. Some Gag mutations can restore Gag-protease binding, yet most Gag-protease correlated mutations occur outside of the Gag cleavage site. To investigate the molecular basis for this, we now report multiscale modelling approaches to investigate various sequentially cleaved Gag products in the context of clinically relevant mutations that occur outside of the cleavage sites, including simulations of the largest Gag proteolytic product in its viral membrane-bound state. We found that some mutations, such as G123E and H219Q, involve direct interaction with cleavage site residues to influence their local environment, while certain mutations in the matrix domain lead to the enrichment of lipids important for Gag targeting and assembly. Collectively, our results reveal why non-cleavage site mutations have far-reaching implications outside of Gag proteolysis, with important consequences for drugging Gag maturation intermediates and tackling protease inhibitor resistance.

Protein crowding mediates membrane remodeling in upstream ESCRT-induced formation of intraluminal vesicles

As part of the lysosomal degradation pathway, the endosomal sorting complexes required for transport (ESCRT-0 to -III/VPS4) sequester receptors at the endosome and simultaneously deform the membrane to generate intraluminal vesicles (ILVs). Whereas ESCRT-III/VPS4 have an established function in ILV formation, the role of upstream ESCRTs (0 to II) in membrane shape remodeling is not understood. Combining experimental measurements and electron microscopy analysis of ESCRT-III-depleted cells with a mathematical model, we show that upstream ESCRT-induced alteration of the Gaussian bending rigidity and their crowding in concert with the transmembrane cargo on the membrane induce membrane deformation and facilitate ILV formation: Upstream ESCRT-driven budding does not require ATP consumption as only a small energy barrier needs to be overcome. Our model predicts that ESCRTs do not become part of the ILV, but localize with a high density at the membrane neck, where the steep decline in the Gaussian curvature likely triggers ESCRT-III/VPS4 assembly to enable neck constriction and scission.

Binding mechanism of the matrix domain of HIV-1 gag on lipid membranes

Specific protein-lipid interactions are critical for viral assembly. We present a molecular dynamics simulation study on the binding mechanism of the membrane targeting domain of HIV-1 Gag protein. The matrix (MA) domain drives Gag onto the plasma membrane through electrostatic interactions at its highly-basic-region (HBR), located near the myristoylated (Myr) N-terminus of the protein. Our study suggests Myr insertion is involved in the sorting of membrane lipids around the protein-binding site to prepare it for viral assembly. Our realistic membrane models confirm interactions with PIP₂ and PS lipids are highly favored around the HBR and are strong enough to keep the protein bound even without Myr insertion. We characterized Myr insertion events from microsecond trajectories and examined the membrane response upon initial membrane targeting by MA. Insertion events only occur with one of the membrane models, showing a combination of surface charge and internal membrane structure modulate this process.

Simple, General Criterion for Onset of Disclination Disorder on Curved Surfaces

Determining the positions of lattice defects on bounded elastic surfaces with Gaussian curvature is a nontrivial task of mechanical energy optimization. We introduce a simple way to predict the onset of disclination disorder from the shape of the surface. The criterion fixes the value of a weighted integral Gaussian curvature to a universal constant and proves accurate across a great variety of shapes. It provides improved understanding of the limitations to crystalline order in many natural and engineering contexts, such as the assembly of viral capsids.

TRIM5α self-assembly and compartmentalization of the HIV-1 viral capsid

The tripartite-motif protein, TRIM5α, is an innate immune sensor that potently restricts retrovirus infection by binding to human immunodeficiency virus capsids. Higher-ordered oligomerization of this protein forms hexagonally patterned structures that wrap around the viral capsid, despite an anomalously low affinity for the capsid protein (CA). Several studies suggest TRIM5α oligomerizes into a lattice with a symmetry and spacing that matches the underlying capsid, to compensate for the weak affinity, yet little is known about how these lattices form. Using a combination of computational simulations and electron cryo-tomography imaging, we reveal the dynamical mechanisms by which these lattices self-assemble. Constrained diffusion allows the lattice to reorganize, whereas defects form on highly curved capsid surfaces to alleviate strain and lattice symmetry mismatches. Statistical analysis localizes the TRIM5α binding interface at or near the CypA binding loop of CA. These simulations elucidate the molecular-scale mechanisms of viral capsid cellular compartmentalization by TRIM5α.

Tripartite-motif containing (TRIM) proteins modulate cellular responses to viral infection. Here the authors use molecular dynamics simulations to demonstrate that TRIM5α uses a two-dimensional lattice hopping mechanism to aggregate on the HIV capsid surface and initiate lattice growth.

Multiscale modelling and simulation of viruses

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