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.

SlyB encapsulates outer membrane proteins in stress-induced lipid nanodomains.

The outer membrane in Gram-negative bacteria consists of an asymmetric phospholipid-lipopolysaccharide bilayer that is densely packed with outer-membrane ß-barrel proteins (OMPs) and lipoproteins1. The architecture and composition of this bilayer is closely monitored and is essential to cell integrity and survival2-4. Here we find that SlyB, a lipoprotein in the PhoPQ stress regulon, forms stable stress-induced complexes with the outer-membrane proteome. SlyB comprises a 10 kDa periplasmic ß-sandwich domain and a glycine zipper domain that forms a transmembrane α-helical hairpin with discrete phospholipid- and lipopolysaccharide-binding sites. After loss in lipid asymmetry, SlyB oligomerizes into ring-shaped transmembrane complexes that encapsulate ß-barrel proteins into lipid nanodomains of variable size. We find that the formation of SlyB nanodomains is essential during lipopolysaccharide destabilization by antimicrobial peptides or acute cation shortage, conditions that result in a loss of OMPs and compromised outer-membrane barrier function in the absence of a functional SlyB. Our data reveal that SlyB is a compartmentalizing transmembrane guard protein that is involved in cell-envelope proteostasis and integrity, and suggest that SlyB represents a larger family of broadly conserved lipoproteins with 2TM glycine zipper domains with the ability to form lipid nanodomains.

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.

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.

Critical mechanistic features of HIV-1 viral capsid assembly.

The maturation of HIV-1 capsid protein (CA) into a cone-shaped lattice capsid is critical for viral infectivity. CA can self-assemble into a range of capsid morphologies made of ~175 to 250 hexamers and 12 pentamers. The cellular polyanion inositol hexakisphosphate (IP6) has recently been demonstrated to facilitate conical capsid formation by coordinating a ring of arginine residues within the central cavity of capsid hexamers and pentamers. However, the kinetic interplay of events during IP6 and CA coassembly is unclear. In this work, we use coarse-grained molecular dynamics simulations to elucidate the molecular mechanism of capsid formation, including the role played by IP6. We show that IP6, in small quantities at first, promotes curvature generation by trapping pentameric defects in the growing lattice and shifts assembly behavior toward kinetically favored outcomes. Our analysis also suggests that IP6 can stabilize metastable capsid intermediates and can induce structural pleomorphism in mature capsids.

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.

Seipin transmembrane segments critically function in triglyceride nucleation and lipid droplet budding from the membrane.

Lipid droplets (LDs) are organelles formed in the endoplasmic reticulum (ER) to store triacylglycerol (TG) and sterol esters. The ER protein seipin is key for LD biogenesis. Seipin forms a cage-like structure, with each seipin monomer containing a conserved hydrophobic helix and two transmembrane (TM) segments. How the different parts of seipin function in TG nucleation and LD budding is poorly understood. Here, we utilized molecular dynamics simulations of human seipin, along with cell-based experiments, to study seipin's functions in protein-lipid interactions, lipid diffusion, and LD maturation. An all-atom simulation indicates that seipin TM segment residues and hydrophobic helices residues located in the phospholipid tail region of the bilayer attract TG. Simulating larger, growing LDs with coarse-grained models, we find that the seipin TM segments form a constricted neck structure to facilitate conversion of a flat oil lens into a budding LD. Using cell experiments and simulations, we also show that conserved, positively charged residues at the end of seipin's TM segments affect LD maturation. We propose a model in which seipin TM segments critically function in TG nucleation and LD growth.

Cooperative multivalent receptor binding promotes exposure of the SARS-CoV-2 fusion machinery core.

The molecular events that permit the spike glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to bind and enter cells are important to understand for both fundamental and therapeutic reasons. Spike proteins consist of S1 and S2 domains, which recognize angiotensin-converting enzyme 2 (ACE2) receptors and contain the viral fusion machinery, respectively. Ostensibly, the binding of spike trimers to ACE2 receptors promotes dissociation of the S1 domains and exposure of the fusion machinery, although the molecular details of this process have yet to be observed. We report the development of bottom-up coarse-grained (CG) models consistent with cryo-electron tomography data, and the use of CG molecular dynamics simulations to investigate viral binding and S2 core exposure. We show that spike trimers cooperatively bind to multiple ACE2 dimers at virion-cell interfaces in a manner distinct from binding between soluble proteins, which processively induces S1 dissociation. We also simulate possible variant behavior using perturbed CG models, and find that ACE2-induced S1 dissociation is primarily sensitive to conformational state populations and the extent of S1/S2 cleavage, rather than ACE2 binding affinity. These simulations reveal an important concerted interaction between spike trimers and ACE2 dimers that primes the virus for membrane fusion and entry.

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.

Cooperative multivalent receptor binding promotes exposure of the SARS-CoV-2 fusion machinery core.

The molecular events that permit the spike glycoprotein of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) to bind, fuse, and enter cells are important to understand for both fundamental and therapeutic reasons. Spike proteins consist of S1 and S2 domains, which recognize angiotensin-converting enzyme 2 (ACE2) receptors and contain the viral fusion machinery, respectively. Ostensibly, the binding of spike trimers to ACE2 receptors promotes the preparation of the fusion machinery by dissociation of the S1 domains. We report the development of bottom-up coarse-grained (CG) models validated with cryo-electron tomography (cryo-ET) data, and the use of CG molecular dynamics simulations to investigate the dynamical mechanisms involved in viral binding and exposure of the S2 trimeric core. We show that spike trimers cooperatively bind to multiple ACE2 dimers at virion-cell interfaces. The multivalent interaction cyclically and processively induces S1 dissociation, thereby exposing the S2 core containing the fusion machinery. Our simulations thus reveal an important concerted interaction between spike trimers and ACE2 dimers that primes the virus for membrane fusion and entry.

Immature HIV-1 assembles from Gag dimers leaving partial hexamers at lattice edges as potential substrates for proteolytic 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.

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.

A new one-site coarse-grained model for water: Bottom-up many-body projected water (BUMPer). I. General theory and model.

Water is undoubtedly one of the most important molecules for a variety of chemical and physical systems, and constructing precise yet effective coarse-grained (CG) water models has been a high priority for computer simulations. To recapitulate important local correlations in the CG water model, explicit higher-order interactions are often included. However, the advantages of coarse-graining may then be offset by the larger computational cost in the model parameterization and simulation execution. To leverage both the computational efficiency of the CG simulation and the inclusion of higher-order interactions, we propose a new statistical mechanical theory that effectively projects many-body interactions onto pairwise basis sets. The many-body projection theory presented in this work shares similar physics from liquid state theory, providing an efficient approach to account for higher-order interactions within the reduced model. We apply this theory to project the widely used Stillinger-Weber three-body interaction onto a pairwise (two-body) interaction for water. Based on the projected interaction with the correct long-range behavior, we denote the new CG water model as the Bottom-Up Many-Body Projected Water (BUMPer) model, where the resultant CG interaction corresponds to a prior model, the iteratively force-matched model. Unlike other pairwise CG models, BUMPer provides high-fidelity recapitulation of pair correlation functions and three-body distributions, as well as N-body correlation functions. BUMPer extensively improves upon the existing bottom-up CG water models by extending the accuracy and applicability of such models while maintaining a reduced computational cost.

Lipid-Composition-Mediated Forces Can Stabilize Tubular Assemblies of I-BAR Proteins.

Collective action by inverse-Bin/Amphiphysin/Rvs (I-BAR) domains drive micron-scale membrane remodeling. The macroscopic curvature sensing and generation behavior of I-BAR domains is well characterized, and computational models have suggested various mechanisms on simplified membrane systems, but there remain missing connections between the complex environment of the cell and the models proposed thus far. Here, we show a connection between the role of protein curvature and lipid clustering in the relaxation of large membrane deformations. When we include phosphatidylinositol 4,5-bisphosphate-like lipids that preferentially interact with the charged ends of an I-BAR domain, we find clustering of phosphatidylinositol 4,5-bisphosphate-like lipids that induce a directional membrane-mediated interaction between membrane-bound I-BAR domains. Lipid clusters mediate I-BAR domain interactions and cause I-BAR domain aggregates that would not arise through membrane fluctuation-based or curvature-based interactions. Inside of membrane protrusions, lipid cluster-mediated interaction draws long side-by-side aggregates together, resulting in more cylindrical protrusions as opposed to bulbous, irregularly shaped protrusions.

A multiscale coarse-grained model of the SARS-CoV-2 virion.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and ongoing development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data become publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion.

A Multiscale Coarse-grained Model of the SARS-CoV-2 Virion.

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the causative agent of the COVID-19 pandemic. Computer simulations of complete viral particles can provide theoretical insights into large-scale viral processes including assembly, budding, egress, entry, and fusion. Detailed atomistic simulations, however, are constrained to shorter timescales and require billion-atom simulations for these processes. Here, we report the current status and on-going development of a largely "bottom-up" coarse-grained (CG) model of the SARS-CoV-2 virion. Structural data from a combination of cryo-electron microscopy (cryo-EM), x-ray crystallography, and computational predictions were used to build molecular models of structural SARS-CoV-2 proteins, which were then assembled into a complete virion model. We describe how CG molecular interactions can be derived from all-atom simulations, how viral behavior difficult to capture in atomistic simulations can be incorporated into the CG models, and how the CG models can be iteratively improved as new data becomes publicly available. Our initial CG model and the detailed methods presented are intended to serve as a resource for researchers working on COVID-19 who are interested in performing multiscale simulations of the SARS-CoV-2 virion. SIGNIFICANCE STATEMENT: This study reports the construction of a molecular model for the SARS-CoV-2 virion and details our multiscale approach towards model refinement. The resulting model and methods can be applied to and enable the simulation of SARS-CoV-2 virions.

A helical assembly of human ESCRT-I scaffolds reverse-topology membrane scission.

The ESCRT complexes drive membrane scission in HIV-1 release, autophagosome closure, multivesicular body biogenesis, cytokinesis, and other cell processes. ESCRT-I is the most upstream complex and bridges the system to HIV-1 Gag in virus release. The crystal structure of the headpiece of human ESCRT-I comprising TSG101-VPS28-VPS37B-MVB12A was determined, revealing an ESCRT-I helical assembly with a 12-molecule repeat. Electron microscopy confirmed that ESCRT-I subcomplexes form helical filaments in solution. Mutation of VPS28 helical interface residues blocks filament formation in vitro and autophagosome closure and HIV-1 release in human cells. Coarse-grained (CG) simulations of ESCRT assembly at HIV-1 budding sites suggest that formation of a 12-membered ring of ESCRT-I molecules is a geometry-dependent checkpoint during late stages of Gag assembly and HIV-1 budding and templates ESCRT-III assembly for membrane scission. These data show that ESCRT-I is not merely a bridging adaptor; it has an essential scaffolding and mechanical role in its own right.

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α.