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.

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.

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.

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.

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.

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.

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

The packaging and budding of Gag polyprotein and viral RNA is a critical step in the HIV-1 life cycle. High-resolution structures of the Gag polyprotein have revealed that the capsid (CA) and spacer peptide 1 (SP1) domains contain important interfaces for Gag self-assembly. However, the molecular details of the multimerization process, especially in the presence of RNA and the cell membrane, have remained unclear. In this work, we investigate the mechanisms that work in concert between the polyproteins, RNA, and membrane to promote immature lattice growth. We develop a coarse-grained (CG) computational model that is derived from subnanometer resolution structural data. Our simulations recapitulate contiguous and hexameric lattice assembly driven only by weak anisotropic attractions at the helical CA-SP1 junction. Importantly, analysis from CG and single-particle tracking photoactivated localization (spt-PALM) trajectories indicates that viral RNA and the membrane are critical constituents that actively promote Gag multimerization through scaffolding, while overexpression of short competitor RNA can suppress assembly. We also find that the CA amino-terminal domain imparts intrinsic curvature to the Gag lattice. As a consequence, immature lattice growth appears to be coupled to the dynamics of spontaneous membrane deformation. Our findings elucidate a simple network of interactions that regulate the early stages of HIV-1 assembly and budding.

Unusual Organization of I-BAR Proteins on Tubular and Vesicular Membranes.

Protein-mediated membrane remodeling is a ubiquitous and critical process for proper cellular function. Inverse Bin/Amphiphysin/Rvs (I-BAR) domains drive local membrane deformation as a precursor to large-scale membrane remodeling. We employ a multiscale approach to provide the molecular mechanism of unusual I-BAR domain-driven membrane remodeling at a low protein surface concentration with near-atomistic detail. We generate a bottom-up coarse-grained model that demonstrates similar membrane-bound I-BAR domain aggregation behavior as our recent Mesoscopic Membrane with Explicit Proteins model. Together, these models bridge several length scales and reveal an aggregation behavior of I-BAR domains. We find that at low surface coverage (i.e., low bound protein density), I-BAR domains form transient, tip-to-tip strings on periodic flat membrane sheets. Inside of lipid bilayer tubules, we find linear aggregates parallel to the axis of the tubule. Finally, we find that I-BAR domains form tip-to-tip aggregates around the edges of membrane domes. These results are supported by in vitro experiments showing low curvature bulges surrounded by I-BAR domains on giant unilamellar vesicles. Overall, our models reveal new I-BAR domain aggregation behavior in membrane tubules and on the surface of vesicles at low surface concentration that add insight into how I-BAR domain proteins may contribute to certain aspects of membrane remodeling in cells.

Off-Pathway Assembly: A Broad-Spectrum Mechanism of Action for Drugs That Undermine Controlled HIV-1 Viral Capsid Formation.

The early and late stages of human immunodeficiency virus (HIV) replication are orchestrated by the capsid (CA) protein, which self-assembles into a conical protein shell during viral maturation. Small molecule drugs known as capsid inhibitors (CIs) impede the highly regulated activity of CA. Intriguingly, a few CIs, such as PF-3450074 (PF74) and GS-CA1, exhibit effects at multiple stages of the viral lifecycle at effective concentrations in the pM to nM regimes, while the majority of CIs target a single stage of the viral lifecycle and are effective at nM to µM concentrations. In this work, we use coarse-grained molecular dynamics simulations to elucidate the molecular mechanisms that enable CIs to have such curious broad-spectrum activity. Our quantitatively analyzed findings show that CIs can have a profound impact on the hierarchical self-assembly of CA by perturbing populations of small CA oligomers. The self-assembly process is accelerated by the emergence of alternative assembly pathways that favor the rapid incorporation of CA pentamers, and leads to increased structural pleomorphism in mature capsids. Two relevant phenotypes are observed: (1) eccentric capsid formation that may fail to encase the viral genome and (2) rapid disassembly of the capsid, which express at late and early stages of infection, respectively. Finally, our study emphasizes the importance of adopting a dynamical perspective on inhibitory mechanisms and provides a basis for the design of future therapeutics that are effective at low stoichiometric ratios of drug to protein.

The Bright Future for Electrode Materials of Energy Devices: Highly Conductive Porous Na-Embedded Carbon.

High electrical conductivity and large accessible surface area, which are required for ideal electrode materials of energy conversion and storage devices, are opposed to each other in current materials. It is a long-term goal to solve this issue. Herein, we report highly conductive porous Na-embedded carbon (Na@C) nanowalls with large surface areas, which have been synthesized by an invented reaction of CO with liquid Na. Their electrical conductivities are 2 orders of magnitude larger than highly conductive 3D graphene. Furthermore, almost all their surface areas are accessible for electrolyte ions. These unique properties make them ideal electrode materials for energy devices, which significantly surpass expensive Pt. Consequently, the dye-sensitized solar cells (DSSCs) with the Na@C counter electrode has reached a high power conversion efficiency of 11.03%. The Na@C also exhibited excellent performance for supercapacitors, leading to high capacitance of 145 F g-1 at current density of 1 A g-1.

What is the thermal conductivity limit of silicon germanium alloys?

The lowest possible thermal conductivity of silicon-germanium (SiGe) bulk alloys achievable through alloy scattering, or the so-called alloy limit, is important to identify for thermoelectric applications. However, this limit remains a subject of contention as both experimentally-reported and theoretically-predicted values tend to be widely scattered and inconclusive. In this work, we present a possible explanation for these discrepancies by demonstrating that the thermal conductivity can vary significantly depending on the degree of randomness in the spatial arrangement of the constituent atoms. Our study suggests that the available experimental data, obtained from alloy samples synthesized using ball-milling techniques, and previous first-principles calculations, restricted by small supercell sizes, may not have accessed the alloy limit. We find that low-frequency anharmonic phonon modes can persist unless the spatial distribution of Si and Ge atoms is completely random at the atomic scale, in which case the lowest possible thermal conductivity may be achieved. Our theoretical analysis predicts that the alloy limit of SiGe could be around 1-2 W m(-1) K(-1) with an optimal composition around 25 at% Ge, which is substantially lower than previously reported values from experiments and first-principles calculations.

Thickness-Dependent Dielectric Constant of Few-Layer In2Se3 Nanoflakes.

The dielectric constant or relative permittivity (ε(r)) of a dielectric material, which describes how the net electric field in the medium is reduced with respect to the external field, is a parameter of critical importance for charging and screening in electronic devices. Such a fundamental material property is intimately related to not only the polarizability of individual atoms but also the specific atomic arrangement in the crystal lattice. In this Letter, we present both experimental and theoretical investigations on the dielectric constant of few-layer In2Se3 nanoflakes grown on mica substrates by van der Waals epitaxy. A nondestructive microwave impedance microscope is employed to simultaneously quantify the number of layers and local electrical properties. The measured ε(r) increases monotonically as a function of the thickness and saturates to the bulk value at around 6-8 quintuple layers. The same trend of layer-dependent dielectric constant is also revealed by first-principles calculations. Our results of the dielectric response, being ubiquitously applicable to layered 2D semiconductors, are expected to be significant for this vibrant research field.

On the influence of polarization effects in predicting the interfacial structure and capacitance of graphene-like electrodes in ionic liquids.

The electric double layer (CD) and electrode quantum (CQ) capacitances of graphene-based supercapacitors are investigated using a combined molecular dynamics and density functional theory approach. In particular, we compare an approach that includes electronic polarization to one that is polarization-free by evaluating both CD and CQ using [EMIM][BF4] ionic liquid as a model electrolyte. Our results indicate that the inclusion of polarization effects can yield higher CD values-in this study by up to 40% around ±2 V-which we attribute primarily to the presence of charge smearing at the electrode-electrolyte interface. On the other hand, we find that the polarization-induced distortion of the electronic structure of graphene does not noticeably alter the predicted CQ. Our analysis suggests that an accurate description of the spatial charge distribution at the graphene interface due to polarization is necessary to improve our predictive capabilities, though more notably for CD. However, the conventional polarization-free approximation can serve as an efficient tool to study trends associated with both the CQ and CD at the interface of various graphene-like materials.

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.

Relative contributions of quantum and double layer capacitance to the supercapacitor performance of carbon nanotubes in an ionic liquid.

Motivated by promising demonstrations of carbon nanotube (CNT) electrodes in supercapacitors, we evaluate the capacitive performance of a (6,6) CNT in [BMIM][PF6] ionic liquid (IL), with particular attention to the relative contributions of the electric double layer (EDL) capacitance (CD) at the CNT/IL interface and the quantum capacitance (CQ) of the CNT. Our classical molecular dynamics simulations reveal that the use of the CNT improves CD when compared to planar graphene, which we discuss in terms of how the electrode curvature affects both the electric field strength and IL packing density. In addition, according to density functional theory calculations, the CQ of the CNT is constant and significantly larger than that of graphene near the Fermi level, which is a consequence of the larger number of available electron states in the CNT. Our study also shows that the relative performance of the CNT- and graphene-based electrodes can be a strong function of applied voltage, which we attribute to the shifting contributions of CQ and CD.

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.