Mature HIV-1 capsid structure by cryo-electron microscopy and all-atom molecular dynamics

Molecular size dependence on achievable resolution from XFEL single-particle 3D reconstruction

Single-particle analysis using x-ray free-electron lasers (XFELs) is a novel method for obtaining structural information of samples in a state close to nature. In particular, it is suitable for observing the inner structure of large biomolecules by taking advantage of the high transmittance of x-rays. However, systematic studies on the resolution achievable for large molecules are lacking. In this study, the molecular size dependence of the resolution of a three-dimensional (3D) structure resulting from XFEL single-particle reconstruction is evaluated using synthetic data. Evidently, 3D structures of larger molecules can be restored with higher detail (defined relative to the molecular sizes) than smaller ones; however, reconstruction with high absolute resolution (defined in nm^-1) is challenging. Our results provide useful information for the experimental design of 3D structure reconstruction using coherent x-ray diffraction patterns of single-particles.

Understanding Retroviral Life Cycle and its Genomic RNA Packaging

Members of the family Retroviridae are important animal and human pathogens. Being obligate parasites, their replication involves a series of steps during which the virus hijacks the cellular machinery. Additionally, many of the steps of retrovirus replication are unique among viruses, including reverse transcription, integration, and specific packaging of their genomic RNA (gRNA) as a dimer. Progress in retrovirology has helped identify several molecular mechanisms involved in each of these steps, but many are still unknown or remain controversial. This review summarizes our present understanding of the molecular mechanisms involved in various stages of retrovirus replication. Furthermore, it provides a comprehensive analysis of our current understanding of how different retroviruses package their gRNA into the assembling virions. RNA packaging in retroviruses holds a special interest because of the uniqueness of packaging a dimeric genome. Dimerization and packaging are highly regulated and interlinked events, critical for the virus to decide whether its unspliced RNA will be packaged as a "genome" or translated into proteins. Finally, some of the outstanding areas of exploration in the field of RNA packaging are highlighted, such as the role of epitranscriptomics, heterogeneity of transcript start sites, and the necessity of functional polyA sequences. An in-depth knowledge of mechanisms that interplay between viral and cellular factors during virus replication is critical in understanding not only the virus life cycle, but also its pathogenesis, and development of new antiretroviral compounds, vaccines, as well as retroviral-based vectors for human gene therapy.

Investigational drugs for HIV: trends, opportunities and key players

INTRODUCTION

Since the first antiretroviral drug was described, the field of HIV treatment and prevention has undergone two drug-based revolutions: the first one, enabled by the virtually concomitant discovery of non-nucleoside reverse transcriptase and protease inhibitors, was the inception of combined antiretroviral therapy. The second followed the creation of integrase strand-transfer inhibitors with improved safety, potency, and resistance profiles. Long-acting antiretroviral drugs, including broadly neutralizing antibodies, now offer the opportunity for a third transformational change in HIV management.

AREAS COVERED

Our review focused on HIV treatment and prevention with investigational drugs that offer the potential for infrequent dosing, including drugs not yet approved for clinical use. We also discussed approved drugs for which administration modalities or formulations are being optimized. We performed a literature search in published manuscripts, conference communications, and registered clinical trials.

EXPERT OPINION

While the field focuses on extending dosing intervals, we identify drug tissue penetration as an understudied opportunity to improve HIV care. We repeat that self-administration remains an essential milestone to reach the full potential of long-acting drugs. Treatments and prevention strategies based on broadly neutralizing antibodies require a deeper understanding of their antiretroviral properties.

Molecular Modeling Insights into the Structure and Behavior of Integrins: A Review

Integrins are heterodimeric glycoproteins crucial to the physiology and pathology of many biological functions. As adhesion molecules, they mediate immune cell trafficking, migration, and immunological synapse formation during inflammation and cancer. The recognition of the vital roles of integrins in various diseases revealed their therapeutic potential. Despite the great effort in the last thirty years, up to now, only seven integrin-based drugs have entered the market. Recent progress in deciphering integrin functions, signaling, and interactions with ligands, along with advancement in rational drug design strategies, provide an opportunity to exploit their therapeutic potential and discover novel agents. This review will discuss the molecular modeling methods used in determining integrins' dynamic properties and in providing information toward understanding their properties and function at the atomic level. Then, we will survey the relevant contributions and the current understanding of integrin structure, activation, the binding of essential ligands, and the role of molecular modeling methods in the rational design of antagonists. We will emphasize the role played by molecular modeling methods in progress in these areas and the designing of integrin antagonists.

Computing the Volume, Surface Area, Mean, and Gaussian Curvatures of Molecules and Their Derivatives

Geometry is crucial in our efforts to comprehend the structures and dynamics of biomolecules. For example, volume, surface area, and integrated mean and Gaussian curvature of the union of balls representing a molecule are used to quantify its interactions with the water surrounding it in the morphometric implicit solvent models. The Alpha Shape theory provides an accurate and reliable method for computing these geometric measures. In this paper, we derive homogeneous formulas for the expressions of these measures and their derivatives with respect to the atomic coordinates, and we provide algorithms that implement them into a new software package, AlphaMol. The only variables in these formulas are the interatomic distances, making them insensitive to translations and rotations. AlphaMol includes a sequential algorithm and a parallel algorithm. In the parallel version, we partition the atoms of the molecule of interest into 3D rectangular blocks, using a kd-tree algorithm. We then apply the sequential algorithm of AlphaMol to each block, augmented by a buffer zone to account for atoms whose ball representations may partially cover the block. The current parallel version of AlphaMol leads to a 20-fold speed-up compared to an independent serial implementation when using 32 processors. For instance, it takes 31 s to compute the geometric measures and derivatives of each atom in a viral capsid with more than 26 million atoms on 32 Intel processors running at 2.7 GHz. The presence of the buffer zones, however, leads to redundant computations, which ultimately limit the impact of using multiple processors. AlphaMol is available as an OpenSource software.

Determining Spatial Variability of Elastic Properties for Biological Samples Using AFM

Measuring the mechanical properties (i.e., elasticity in terms of Young's modulus) of biological samples using Atomic Force Microscopy (AFM) indentation at the nanoscale has opened new horizons in studying and detecting various pathological conditions at early stages, including cancer and osteoarthritis. It is expected that AFM techniques will play a key role in the future in disease diagnosis and modeling using rigorous mathematical criteria (i.e., automated user-independent diagnosis). In this review, AFM techniques and mathematical models for determining the spatial variability of elastic properties of biological materials at the nanoscale are presented and discussed. Significant issues concerning the rationality of the elastic half-space assumption, the possibility of monitoring the depth-dependent mechanical properties, and the construction of 3D Young's modulus maps are also presented.

Low-molecular-weight anti-HIV-1 agents targeting HIV-1 capsid proteins

The HIV-1 capsid is a shell that encapsulates viral RNA, and forms a conical structure by assembling oligomers of capsid (CA) proteins. Since the CA proteins are highly conserved among many strains of HIV-1, the inhibition of the CA function could be an appropriate goal for suppression of HIV-1 replication, but to date, no drug targeting CA has been developed. Hydrophobic interactions between two CA molecules through Trp184 and Met185 in the protein are known to be indispensable for conformational stabilization of the CA multimer. In our previous study, a small molecule designed by in silico screening as a dipeptide mimic of Trp184 and Met185 in the interaction site was synthesized and found to have significant anti-HIV-1 activity. In the present study, molecules with different scaffolds based on a dipeptide mimic of Trp184 and Met185 have been designed and synthesized. Their significant anti-HIV activity and their advantages compared to the previous compounds were examined. The present results should be useful in the design of novel CA-targeting anti-HIV agents.

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.

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.

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.

Using physical property surrogate models to perform accelerated multi-fidelity optimization of force field parameters†

Accurate representations of van der Waals dispersion–repulsion interactions play an important role in high-quality molecular dynamics simulations. Training the force field parameters used in the Lennard Jones (LJ) potential typically used to represent these interactions is challenging, generally requiring adjustment based on simulations of macroscopic physical properties. The large computational expense of these simulations, especially when many parameters must be trained simultaneously, limits the size of training data set and number of optimization steps that can be taken, often requiring modelers to perform optimizations within a local parameter region. To allow for more global LJ parameter optimization against large training sets, we introduce a multi-fidelity optimization technique which uses Gaussian process surrogate modeling to build inexpensive models of physical properties as a function of LJ parameters. This approach allows for fast evaluation of approximate objective functions, greatly accelerating searches over parameter space and enabling the use of optimization algorithms capable of searching more globally. In this study, we use an iterative framework which performs global optimization with differential evolution at the surrogate level, followed by validation at the simulation level and surrogate refinement. Using this technique on two previously studied training sets, containing up to 195 physical property targets, we refit a subset of the LJ parameters for the OpenFF 1.0.0 (Parsley) force field. We demonstrate that this multi-fidelity technique can find improved parameter sets compared to a purely simulation-based optimization by searching more broadly and escaping local minima. Additionally, this technique often finds significantly different parameter minima that have comparably accurate performance. In most cases, these parameter sets are transferable to other similar molecules in a test set. Our multi-fidelity technique provides a platform for rapid, more global optimization of molecular models against physical properties, as well as a number of opportunities for further refinement of the technique.

We present a multi-fidelity method for optimizing nonbonded force field parameters against physical property data. Leveraging fast surrogate models, we accelerate the parameter search and find novel solutions that improve force field performance.

In Situ Imaging of Virus-Infected Cells by Cryo-Electron Tomography: An Overview

Cryo-electron tomography (cryo-ET) has emerged as a powerful tool in structural biology to study viruses and is undergoing a resolution revolution. Enveloped viruses comprise several RNA and DNA pleomorphic viruses that are pathogens of clinical importance to humans and animals. Considerable efforts in cryogenic correlative light and electron microscopy (cryo-CLEM), cryogenic focused ion beam milling (cryo-FIB), and integrative structural techniques are helping to identify virus structures within cells leading to a rise of in situ discoveries shedding light on how viruses interact with their hosts during different stages of infection. This chapter reviews recent advances in the application of cryo-ET in imaging enveloped viruses and the structural and mechanistic insights revealed studying the viral infection cycle within their eukaryotic cellular hosts, with particular attention to viral entry, replication, assembly, and egress during infection.

Poxvirus under the eyes of electron microscope

Zoonotic poxvirus infections pose significant threat to human health as we have witnessed recent spread of monkeypox. Therefore, insights into molecular mechanism behind poxvirus replication cycle are needed for the development of efficient antiviral strategies. Virion assembly is one of the key steps that determine the fate of replicating poxviruses. However, in-depth understanding of poxvirus assembly is challenging due to the complex nature of multi-step morphogenesis and heterogeneous virion structures. Despite these challenges, decades of research have revealed virion morphologies at various maturation stages, critical protein components and interactions with host cell compartments. Transmission electron microscopy has been employed as an indispensable tool for the examination of virion morphology, and more recently for the structure determination of protein complexes. In this review, we describe some of the major findings in poxvirus morphogenesis and the contributions of continuously advancing electron microscopy techniques.

A click chemistry amplified nanopore assay for ultrasensitive quantification of HIV-1 p24 antigen in clinical samples

Despite major advances in HIV testing, ultrasensitive detection of early infection remains challenging, especially for the viral capsid protein p24, which is an early virological biomarker of HIV-1 infection. Here, To improve p24 detection in patients missed by immunological tests that dominate the diagnostics market, we show a click chemistry amplified nanopore (CAN) assay for ultrasensitive quantitative detection. This strategy achieves a 20.8 fM (0.5 pg/ml) limit of detection for HIV-1 p24 antigen in human serum, demonstrating 20~100-fold higher analytical sensitivity than nanocluster-based immunoassays and clinically used enzyme-linked immunosorbent assay, respectively. Clinical validation of the CAN assay in a pilot cohort shows p24 quantification at ultra-low concentration range and correlation with CD4 count and viral load. We believe that this strategy can improve the utility of p24 antigen in detecting early infection and monitoring HIV progression and treatment efficacy, and also can be readily modified to detect other infectious diseases.

Bridging length scales from molecules to the whole organism by cryoCLEM and cryoET

Resolving atomic structures of isolated proteins has uncovered mechanisms and fundamental processes in biology. However, many functions can only be tested in the context of intact cells and tissues that are many orders of magnitude larger than the macromolecules on which they depend. Therefore, methods that interrogate macromolecular structure in situ provide a means of directly relating structure to function across length scales. Here, we developed several workflows using cryogenic correlated light and electron microscopy (cryoCLEM) and electron tomography (cryoET) that can bridge this gap to reveal the molecular infrastructure that underlies higher order functions within cells and tissues. We also describe experimental design considerations, including cryoCLEM labelling, sample preparation, and quality control, for determining the in situ molecular architectures within native, hydrated cells and tissues.

Human Immunodeficiency Virus Type 2 Capsid Protein Mutagenesis Reveals Amino Acid Residues Important for Virus Particle Assembly

Human immunodeficiency virus (HIV) Gag drives virus particle assembly. The capsid (CA) domain is critical for Gag multimerization mediated by protein-protein interactions. The Gag protein interaction network defines critical aspects of the retroviral lifecycle at steps such as particle assembly and maturation. Previous studies have demonstrated that the immature particle morphology of HIV-2 is intriguingly distinct relative to that of HIV-1. Based upon this observation, we sought to determine the amino acid residues important for virus assembly that might help explain the differences between HIV-1 and HIV-2. To do this, we conducted site-directed mutagenesis of targeted locations in the HIV-2 CA domain of Gag and analyzed various aspects of virus particle assembly. A panel of 31 site-directed mutants of residues that reside at the HIV-2 CA inter-hexamer interface, intra-hexamer interface and CA inter-domain linker were created and analyzed for their effects on the efficiency of particle production, particle morphology, particle infectivity, Gag subcellular distribution and in vitro protein assembly. Seven conserved residues between HIV-1 and HIV-2 (L19, A41, I152, K153, K157, N194, D196) and two non-conserved residues (G38, N127) were found to significantly impact Gag multimerization and particle assembly. Taken together, these observations complement structural analyses of immature HIV-2 particle morphology and Gag lattice organization as well as provide important comparative insights into the key amino acid residues that can help explain the observed differences between HIV immature particle morphology and its association with virus replication and particle infectivity.

Prion-like low complexity regions enable avid virus-host interactions during HIV-1 infection

Cellular proteins CPSF6, NUP153 and SEC24C play crucial roles in HIV-1 infection. While weak interactions of short phenylalanine-glycine (FG) containing peptides with isolated capsid hexamers have been characterized, how these cellular factors functionally engage with biologically relevant mature HIV-1 capsid lattices is unknown. Here we show that prion-like low complexity regions (LCRs) enable avid CPSF6, NUP153 and SEC24C binding to capsid lattices. Structural studies revealed that multivalent CPSF6 assembly is mediated by LCR-LCR interactions, which are templated by binding of CPSF6 FG peptides to a subset of hydrophobic capsid pockets positioned along adjoining hexamers. In infected cells, avid CPSF6 LCR-mediated binding to HIV-1 cores is essential for functional virus-host interactions. The investigational drug lenacapavir accesses unoccupied hydrophobic pockets in the complex to potently impair HIV-1 inside the nucleus without displacing the tightly bound cellular cofactor from virus cores. These results establish previously undescribed mechanisms of virus-host interactions and antiviral action.

Analysis and comparison of electron radiation damage assessments in Cryo-EM by single particle analysis and micro-crystal electron diffraction

Electron radiation damage to macromolecules is an inevitable resolution limit factor in all major structural determination applications using cryo-electron microscopy (cryo-EM). Single particle analysis (SPA) and micro-crystal electron diffraction (MicroED) have been employed to assess radiation damage with a variety of protein complexes. Although radiation induced sidechain density loss and resolution decay were observed by both methods, the minimum dose of electron irradiation reducing high-resolution limit reported by SPA is more than ten folds higher than measured by MicroED using the conventional dose concept, and there is a gap between the attained resolutions assessed by these two methods. We compared and analyzed these two approaches side-by-side in detail from several aspects to identify some crucial determinants and to explain this discrepancy. Probability of a high energy electron being inelastically scattered by a macromolecule is proportional to number of layers of the molecules in its transmission path. As a result, the same electron dose could induce much more site-specific damage to macromolecules in 3D protein crystal than single particle samples. Major differences in data collection and processing scheme are the key factors to different levels of sensitivity to radiation damage at high resolution between the two methods. High resolution electron diffraction in MicroED dataset is very sensitive to global damage to 3D protein crystals with low dose accumulation, and its intensity attenuation rates at atomic resolution shell could be applied for estimating ratio of damaged and total selected single particles for SPA. More in-depth systematically radiation damage assessments using SPA and MicroED will benefit all applications of cryo-EM, especially cellular structure analysis by tomography.

HIV-1 Preintegration Complex Preferentially Integrates the Viral DNA into Nucleosomes Containing Trimethylated Histone 3-Lysine 36 Modification and Flanking Linker DNA

HIV-1 DNA is preferentially integrated into chromosomal hot spots by the preintegration complex (PIC). To understand the mechanism, we measured the DNA integration activity of PICs-extracted from infected cells-and intasomes, biochemically assembled PIC substructures using a number of relevant target substrates. We observed that PIC-mediated integration into human chromatin is preferred compared to genomic DNA. Surprisingly, nucleosomes lacking histone modifications were not preferred integration compared to the analogous naked DNA. Nucleosomes containing the trimethylated histone 3 lysine 36 (H3K36me3), an epigenetic mark linked to active transcription, significantly stimulated integration, but the levels remained lower than the naked DNA. Notably, H3K36me3-modified nucleosomes with linker DNA optimally supported integration mediated by the PIC but not by the intasome. Interestingly, optimal intasome-mediated integration required the cellular cofactor LEDGF. Unexpectedly, LEDGF minimally affected PIC-mediated integration into naked DNA but blocked integration into nucleosomes. The block for the PIC-mediated integration was significantly relieved by H3K36me3 modification. Mapping the integration sites in the preferred substrates revealed that specific features of the nucleosome-bound DNA are preferred for integration, whereas integration into naked DNA was random. Finally, biochemical and genetic studies demonstrate that DNA condensation by the H1 protein dramatically reduces integration, providing further evidence that features inherent to the open chromatin are preferred for HIV-1 integration. Collectively, these results identify the optimal target substrate for HIV-1 integration, report a mechanistic link between H3K36me3 and integration preference, and importantly, reveal distinct mechanisms utilized by the PIC for integration compared to the intasomes. IMPORTANCE HIV-1 infection is dependent on integration of the viral DNA into the host chromosomes. The preintegration complex (PIC) containing the viral DNA, the virally encoded integrase (IN) enzyme, and other viral/host factors carries out HIV-1 integration. HIV-1 integration is not dependent on the target DNA sequence, and yet the viral DNA is selectively inserted into specific "hot spots" of human chromosomes. A growing body of literature indicates that structural features of the human chromatin are important for integration targeting. However, the mechanisms that guide the PIC and enable insertion of the PIC-associated viral DNA into specific hot spots of the human chromosomes are not fully understood. In this study, we describe a biochemical mechanism for the preference of the HIV-1 DNA integration into open chromatin. Furthermore, our study defines a direct role for the histone epigenetic mark H3K36me3 in HIV-1 integration preference and identify an optimal substrate for HIV-1 PIC-mediated viral DNA integration.

Multiscale Modeling of Hepatitis B Virus Capsid Assembly and Its Dimorphism

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

A self-complementary macrocycle by a dual interaction system

Self-complementary assembly is one of the most promising phenomena for the formation of discrete assemblies, e.g., proteins and capsids. However, self-complementary assembly based on multiple host-guest systems has been scarcely reported due to the difficulty in controlling each assembly. Herein, we report a dual interaction system in which the key assembly direction is well regulated by both π-π stacking and hydrogen bonding to construct a self-complementary macrocycle. Continuous host-guest behavior of anthracene-based molecular tweezers during crystallization leads to successful construction of a cyclic hexamer, which is reminiscent of Kekulé’s monkey model. Furthermore, the cyclic hexamer in a tight and triple-layered fashion shows hierarchical assembly into cuboctahedron and rhombohedral assemblies in the presence of trifluoroacetic acid. Our findings would be potentially one of metal-free strategies for constructing anthracene-based supramolecular assemblies with higher-order structure.

In nature, HIV capsid consists of single class of protein unit by self-complementarity. Here, the authors find that a molecular tweezer forms a cyclic hexamer by its continuous host-guest behavior, and constructs a large cuboctahedron by hierarchical assembly.

The Discovery of GSK3640254, a Next-Generation Inhibitor of HIV-1 Maturation

GSK3640254 is an HIV-1 maturation inhibitor (MI) that exhibits significantly improved antiviral activity toward a range of clinically relevant polymorphic variants with reduced sensitivity toward the second-generation MI GSK3532795 (BMS-955176). The key structural difference between GSK3640254 and its predecessor is the replacement of the para-substituted benzoic acid moiety attached at the C-3 position of the triterpenoid core with a cyclohex-3-ene-1-carboxylic acid substituted with a CH2F moiety at the carbon atom α- to the pharmacophoric carboxylic acid. This structural element provided a new vector with which to explore structure-activity relationships (SARs) and led to compounds with improved polymorphic coverage while preserving pharmacokinetic (PK) properties. The approach to the design of GSK3640254, the development of a synthetic route and its preclinical profile are discussed. GSK3640254 is currently in phase IIb clinical trials after demonstrating a dose-related reduction in HIV-1 viral load over 7-10 days of dosing to HIV-1-infected subjects.

Design, Synthesis and Structure-Activity Relationships of Phenylalanine-Containing Peptidomimetics as Novel HIV-1 Capsid Binders Based on Ugi Four-Component Reaction

As a key structural protein, HIV capsid (CA) protein plays multiple roles in the HIV life cycle, and is considered a promising target for anti-HIV treatment. Based on the structural information of CA modulator PF-74 bound to HIV-1 CA hexamer, 18 novel phenylalanine derivatives were synthesized via the Ugi four-component reaction. In vitro anti-HIV activity assays showed that most compounds exhibited low-micromolar-inhibitory potency against HIV. Among them, compound I-19 exhibited the best anti-HIV-1 activity (EC50 = 2.53 ± 0.84 μM, CC50 = 107.61 ± 27.43 μM). In addition, I-14 displayed excellent HIV-2 inhibitory activity (EC50 = 2.30 ± 0.11 μM, CC50 > 189.32 μM) with relatively low cytotoxicity, being more potent than that of the approved drug nevirapine (EC50 > 15.02 μM, CC50 > 15.2 μM). Additionally, surface plasmon resonance (SPR) binding assays demonstrated direct binding to the HIV CA protein. Moreover, molecular docking and molecular dynamics simulations provided additional information on the binding mode of I-19 to HIV-1 CA. In summary, we further explored the structure—activity relationships (SARs) and selectivity of anti-HIV-1/HIV-2 of PF-74 derivatives, which is conducive to discovering efficient anti-HIV drugs.

Dielectric dispersion characteristics of the phospholipid bilayer with subnanometer resolution from terahertz to mid-infrared

There is growing interest in whether the myelinated nerve fiber acts as a dielectric waveguide to propagate terahertz to mid-infrared electromagnetic waves, which are presumed stable signal carrier for neurotransmission. The myelin sheath is formed as a multilamellar biomembrane structure, hence insights into the dielectric properties of the phospholipid bilayer is essential for a complete understanding of the myelinated fiber functioning. In this work, by means of atomistic molecular dynamics simulations of the dimyristoylphosphatidylcholine (DMPC) bilayer in water and numerical calculations of carefully layered molecules along with calibration of optical dielectric constants, we for the first time demonstrate the spatially resolved (in sub-nm) dielectric spectrum of the phospholipid bilayer in a remarkably wide range from terahertz to mid-infrared. More specifically, the membrane head regions exhibit both larger real and imaginary permittivities than that of the tail counterparts in the majority of the 1–100 THz band. In addition, the spatial variation of dielectric properties suggests advantageous propagation characteristics of the phospholipid bilayer in a relatively wide band of 55–85 THz, where the electromagnetic waves are well confined within the head regions.

Intrahost evolution of the HIV-2 capsid correlates with progression to AIDS

HIV-2 infection will progress to AIDS in most patients without treatment, albeit at approximately half the rate of HIV-1 infection. HIV-2 capsid (p26) amino acid polymorphisms are associated with lower viral loads and enhanced processing of T cell epitopes, which may lead to protective Gag-specific T cell responses common in slower progressors. Lower virus evolutionary rates, and positive selection on conserved residues in HIV-2 env have been associated with slower progression to AIDS. In this study we analysed 369 heterochronous HIV-2 p26 sequences from 12 participants with a median age of 30 years at enrolment. CD4% change over time was used to stratify participants into relative faster and slower progressor groups. We analysed p26 sequence diversity evolution, measured site-specific selection pressures and evolutionary rates, and determined if these evolutionary parameters were associated with progression status. Faster progressors had lower CD4% and faster CD4% decline rates. Median pairwise sequence diversity was higher in faster progressors (5.7x10-3 versus 1.4x10-3 base substitutions per site, P<0.001). p26 evolved under negative selection in both groups (dN/dS=0.12). Median virus evolutionary rates were higher in faster than slower progressors - synonymous rates: 4.6x10-3 vs. 2.3x10-3; and nonsynonymous rates: 6.9x10-4 vs. 2.7x10-4 substitutions/site/year, respectively. Virus evolutionary rates correlated negatively with CD4% change rates (ρ = -0.8, P=0.02), but not CD4% level. The signature amino acid at p26 positions 6, 12 and 119 differed between faster (6A, 12I, 119A) and slower (6G, 12V, 119P) progressors. These amino acid positions clustered near to the TRIM5α/p26 hexamer interface surface. p26 evolutionary rates were associated with progression to AIDS and were mostly driven by synonymous substitutions. Nonsynonymous evolutionary rates were an order of magnitude lower than synonymous rates, with limited amino acid sequence evolution over time within hosts. These results indicate HIV-2 p26 may be an attractive therapeutic target.

The potential role of exosomal miRNAs and membrane proteins in acute HIV-infected people

Exosomes play an important role during human immunodeficiency virus (HIV) acute infection. Yet, information regarding its cargo and its association with HIV rapid progressors (RPs) and typical progressors (TPs) remain largely unknown. In this study, exosomal miRNAs sequencing and mass cytometry were used to identify differential exosomal miRNAs and membrane proteins that participate in the pathogenesis of TPs and RPs. We discovered that miR-144-5p, miR-1180-3p, miR-451a, miR-362-5p, and miR-625-5p are associated with the TPs and miR-362-5p with the RPs. Decreased autophagy, amino acid metabolism, immune response, and IL-6 are closely related to RPs. In addition, SP1 was selected as the most significant transcription factor (TF) associated with disease progression. CD49D, CD5, CCR5, CD40, CD14, and CD86 were selected as the differential exosomal membrane proteins between TPs and RPs. This study provides valuable information for clarifying the mechanism in people with acute HIV infection.

The HIV-1 Gag Protein Displays Extensive Functional and Structural Roles in Virus Replication and Infectivity

Once merely thought of as the protein responsible for the overall physical nature of the human immunodeficiency virus type 1 (HIV-1), the Gag polyprotein has since been elucidated to have several roles in viral replication and functionality. Over the years, extensive research into the polyproteins’ structure has revealed that Gag can mediate its own trafficking to the plasma membrane, it can interact with several host factors and can even aid in viral genome packaging. Not surprisingly, Gag has also been associated with HIV-1 drug resistance and even treatment failure. Therefore, this review provides an extensive overview of the structural and functional roles of the HIV-1 Gag domains in virion integrity, functionality and infectivity.

Electrostatic Map of the SARS-CoV-2 Virion Specifies Binding Sites of the Antiviral Cationic Photosensitizer

Electrostatics is an important part of virus life. Understanding the detailed distribution of charges over the surface of a virus is important to predict its interactions with host cells, antibodies, drugs, and different materials. Using a coarse-grained model of the entire viral envelope developed by D. Korkin and S.-J. Marrink’s scientific groups, we created an electrostatic map of the external surface of SARS-CoV-2 and found a highly heterogeneous distribution of the electrostatic potential field of the viral envelope. Numerous negative patches originate mainly from negatively charged lipid domains in the viral membrane and negatively charged areas on the “stalks” of the spike (S) proteins. Membrane (M) and envelope (E) proteins with the total positive charge tend to colocalize with the negatively charged lipids. In the E protein pentamer exposed to the outer surface, negatively charged glutamate residues and surrounding lipids form a negative electrostatic potential ring around the channel entrance. We simulated the interaction of the antiviral octacationic photosensitizer octakis(cholinyl)zinc phthalocyanine with the surface structures of the entire model virion using the Brownian dynamics computational method implemented in ProKSim software (version r661). All mentioned negatively charged envelope components attracted the photosensitizer molecules and are thus potential targets for reactive oxygen generated in photosensitized reactions.

CP-MAS and solution NMR studies of allosteric communication in CA-assemblies of HIV-1

Solution and solid-state NMR spectroscopy are highly complementary techniques for studying structure and dynamics in very high molecular weight systems. Here we have analysed the dynamics of HIV-1 capsid (CA) assemblies in presence of the cofactors IP6 and ATPγS and the host-factor CPSF6 using a combination of solution state and cross polarisation magic angle spinning (CP-MAS) solid-state NMR. In particular, dynamical effects on ns to µs and µs to ms timescales are observed revealing diverse motions in assembled CA. Using CP-MAS NMR, we exploited the sensitivity of the amide/Cα-Cβ backbone chemical shifts in DARR and NCA spectra to observe the plasticity of the HIV-1 CA tubular assemblies and also map the binding of cofactors and the dynamics of cofactor-CA complexes. In solution, we measured how the addition of host- and co-factors to CA -hexamers perturbed the chemical shifts and relaxation properties of CA-Ile and -Met methyl groups using transverse-relaxation-optimized NMR spectroscopy to exploit the sensitivity of methyl groups as probes in high-molecular weight proteins. These data show how dynamics of the CA protein assembly over a range of spatial and temporal scales play a critical role in CA function. Moreover, we show that binding of IP6, ATPγS and CPSF6 results in local chemical shift as well as dynamic changes for a significant, contiguous portion of CA, highlighting how allosteric pathways communicate ligand interactions between adjacent CA protomers.

PDBx/mmCIF Ecosystem: Foundational Semantic Tools for Structural Biology

PDBx/mmCIF, Protein Data Bank Exchange (PDBx) macromolecular Crystallographic Information Framework (mmCIF), has become the data standard for structural biology. With its early roots in the domain of small-molecule crystallography, PDBx/mmCIF provides an extensible data representation that is used for deposition, archiving, remediation, and public dissemination of experimentally determined three-dimensional (3D) structures of biological macromolecules by the Worldwide Protein Data Bank (wwPDB, wwpdb.org). Extensions of PDBx/mmCIF are similarly used for computed structure models by ModelArchive (modelarchive.org), integrative/hybrid structures by PDB-Dev (pdb-dev.wwpdb.org), small angle scattering data by Small Angle Scattering Biological Data Bank SASBDB (sasbdb.org), and for models computed generated with the AlphaFold 2.0 deep learning software suite (alphafold.ebi.ac.uk). Community-driven development of PDBx/mmCIF spans three decades, involving contributions from researchers, software and methods developers in structural sciences, data repository providers, scientific publishers, and professional societies. Having a semantically rich and extensible data framework for representing a wide range of structural biology experimental and computational results, combined with expertly curated 3D biostructure data sets in public repositories, accelerates the pace of scientific discovery. Herein, we describe the architecture of the PDBx/mmCIF data standard, tools used to maintain representations of the data standard, governance, and processes by which data content standards are extended, plus community tools/software libraries available for processing and checking the integrity of PDBx/mmCIF data. Use cases exemplify how the members of the Worldwide Protein Data Bank have used PDBx/mmCIF as the foundation for its pipeline for delivering Findable, Accessible, Interoperable, and Reusable (FAIR) data to many millions of users worldwide.

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.

Pre-exascale HPC approaches for molecular dynamics simulations. Covid-19 research: A use case

Exascale computing has been a dream for ages and is close to becoming a reality that will impact how molecular simulations are being performed, as well as the quantity and quality of the information derived for them. We review how the biomolecular simulations field is anticipating these new architectures, making emphasis on recent work from groups in the BioExcel Center of Excellence for High Performance Computing. We exemplified the power of these simulation strategies with the work done by the HPC simulation community to fight Covid-19 pandemics. This article is categorized under:Data Science > Computer Algorithms and ProgrammingData Science > Databases and Expert SystemsMolecular and Statistical Mechanics > Molecular Dynamics and Monte-Carlo Methods.

HIV Capsid Protein Genetic Diversity Across HIV-1 Variants and Impact on New Capsid-Inhibitor Lenacapavir

The HIV p24 capsid protein has an essential, structural, and functional role in the viral replication cycle, being an interesting target for vaccine design, diagnostic tests, and new antiretroviral drugs (ARVs). The HIV-1 variability poses a challenge for the accuracy and efficiency of diagnostic and treatment tools. This study analyzes p24 diversity among HIV-1 variants and within its secondary structure in HIV-1 M, O, P, and N groups. All available HIV-1 p24 nucleotide sequences were downloaded from the Los Alamos HIV Sequence Database, selecting 23,671 sequences belonging to groups O, N, P, and M (9 subtypes, 7 sub-sub types, and 109 circulating recombinant forms or CRFs). Using a bioinformatics tool developed in our laboratory (EpiMolBio program), we analyzed the amino acid conservation compared to the HXB2 subtype B reference sequence and the V-markers, or amino acid changes that were specific for each variant with at least 10 available sequences. We inferred the p24 consensus sequence for HIV-1 and for each group to analyze the overall conservation in p24 main structural regions, reporting the percentage of substitutions per variant affecting the capsid assembly and molecule-binding, including those associated with resistance to the new capsid-inhibitor lenacapavir, and the key residues involved in lenacapavir-p24 interaction, according to the bibliography. Although the overall structure of p24 was highly conserved, the conservation in the secondary structure varied between HIV-1 variants and the type of secondary structure. All HIV-1 variants presented >80% amino acid conservation vs. HXB2 reference sequence, except for group M sub-subtype F1 (69.27%). Mutants affecting the capsid assembly or lenacapavir capsid-binding were found in <1% of the p24 consensus sequence. Our study reports the HIV-1 variants carrying 14 unique single V-markers in 9/38 group M variants and the level of p24 conservation in each secondary structure region among the 4 HIV-1 groups and group M variants, revealing no natural resistance to lenacapavir in any HIV-1 variant. We present a thorough analysis of p24 variability among all HIV-1 variants circulating to date. Since p24 genetic variability can impact the viral replication cycle and the efficacy of new p24-based diagnostic, therapeutic, and vaccine strategies, conservation studies must consider all HIV-1 variants circulating worldwide.

HIV-1 mutants that escape the cytotoxic T-lymphocytes are defective in viral DNA integration

HIV-1 replication is durably controlled without antiretroviral therapy (ART) in certain infected individuals called elite controllers (ECs). These individuals express specific human leukocyte antigens (HLA) that tag HIV-infected cells for elimination by presenting viral epitopes to CD8+ cytotoxic T-lymphocytes (CTL). In HIV-infected individuals expressing HLA-B27, CTLs primarily target the viral capsid protein (CA)-derived KK10 epitope. While selection of CA mutation R264K helps HIV-1 escape this potent CTL response, the accompanying fitness cost severely diminishes virus infectivity. Interestingly, selection of a compensatory CA mutation S173A restores HIV-1 replication. However, the molecular mechanism(s) underlying HIV-1 escape from this ART-free virus control by CTLs is not fully understood. Here, we report that the R264K mutation-associated infectivity defect arises primarily from impaired HIV-1 DNA integration, which is restored by the S173A mutation. Unexpectedly, the integration defect of the R264K variant was also restored upon depletion of the host cyclophilin A. These findings reveal a nuclear crosstalk between CA and HIV-1 integration as well as identify a previously unknown role of cyclophilin A in viral DNA integration. Finally, our study identifies a novel immune escape mechanism of an HIV-1 variant escaping a CA-directed CTL response.

Multimodal Functionalities of HIV-1 Integrase

Integrase is the retroviral protein responsible for integrating reverse transcripts into cellular genomes. Co-packaged with viral RNA and reverse transcriptase into capsid-encased viral cores, human immunodeficiency virus 1 (HIV-1) integrase has long been implicated in reverse transcription and virion maturation. However, the underlying mechanisms of integrase in these non-catalytic-related viral replication steps have remained elusive. Recent results have shown that integrase binds genomic RNA in virions, and that mutational or pharmacological disruption of integrase-RNA binding yields eccentric virion particles with ribonucleoprotein complexes situated outside of the capsid shell. Such viruses are defective for reverse transcription due to preferential loss of integrase and viral RNA from infected target cells. Parallel research has revealed defective integrase-RNA binding and eccentric particle formation as common features of class II integrase mutant viruses, a phenotypic grouping of viruses that display defects at steps beyond integration. In light of these new findings, we propose three new subclasses of class II mutant viruses (a, b, and c), all of which are defective for integrase-RNA binding and particle morphogenesis, but differ based on distinct underlying mechanisms exhibited by the associated integrase mutant proteins. We also assess how these findings inform the role of integrase in HIV-1 particle maturation.

Molecular dynamics simulations reveal the selectivity mechanism of structurally similar agonists to TLR7 and TLR8

TLR7 and TLR8 are key members of the Toll-like receptor family, playing crucial roles in the signaling pathways of innate immunity, and thus become attractive therapeutic targets of many diseases including infections and cancer. Although TLR7 and TLR8 show a high degree of sequence homology, their biological response to small molecule binding is very different. Aiming to understand the mechanism of selective profiles of small molecule modulators against TLR7 and TLR8, we carried out molecular dynamic simulations on three imidazoquinoline derivatives bound to the receptors separately. They are Resiquimod (R), Hybrid-2 (H), and Gardiquimod (G), selective agonists of TLR7 and TLR8. Our MD trajectories indicated that in the complex of TLR7-R and TLR7-G, the two chains forming the TLR7 dimer tended to remain “open” conformation, while the rest systems maintained in the closed format. The agonists R, H, and G developed conformational deviation mainly on the aliphatic tail. Furthermore, we attempted to quantify the selectivity between TLR7 and TLR8 by binding free energies via MM-GBSA method. It showed that the three selected modulators were more favorable for TLR7 than TLR8, and the ranking from the strongest to the weakest was H, R and G, aligning well with experimental data. In the TLR7, the flexible and hydrophobic aliphatic side chain of H has stronger van der Waals interactions with V381 and F351 but only pick up interaction with one amino acid residue i.e. Y353 of TLR8. Unsurprisingly, the positively charged side chain of G has less favorable interaction with I585 of TLR7 and V573 of TLR8 explaining G is weak agonist of both TLR7 and TLR8. All three imidazoquinoline derivatives can form stable hydrogen bonds with D555 of TLR7 and the corresponding D543 of TLR8. In brief, the set of total 400ns MD studies sheds light on the potential selectivity mechanisms of agonists towards TLR7 and TLR8, indicating the van der Waals interaction as the driving force for the agonists binding, thus provides us insights for designing more potent and selective modulators to cooperate with the hydrophobic nature of the binding pocket.

Graph neural networks accelerated molecular dynamics

Molecular Dynamics (MD) simulation is a powerful tool for understanding the dynamics and structure of matter. Since the resolution of MD is atomic-scale, achieving long timescale simulations with femtosecond integration is very expensive. In each MD step, numerous iterative computations are performed to calculate energy based on different types of interaction and their corresponding spatial gradients. These repetitive computations can be learned and surrogated by a deep learning model, such as a Graph Neural Network (GNN). In this work, we developed a GNN Accelerated MD (GAMD) model that directly predicts forces, given the state of the system (atom positions, atom types), bypassing the evaluation of potential energy. By training the GNN on a variety of data sources (simulation data derived from classical MD and density functional theory), we show that GAMD can predict the dynamics of two typical molecular systems, Lennard-Jones system and water system, in the NVT ensemble with velocities regulated by a thermostat. We further show that GAMD’s learning and inference are agnostic to the scale, where it can scale to much larger systems at test time. We also perform a comprehensive benchmark test comparing our implementation of GAMD to production-level MD software, showing GAMD’s competitive performance on the large-scale simulation.

The use of vector formalism in the analysis of hydrophobic and electric driving forces in biological assemblies

Hydrophobic forces are known to have a crucial part not only in the conformation of the three-dimensional structure of proteins, but also in the build-up of DNA–protein complexes. Electric forces also play an important role both in the tertiary as well in the quaternary structure of macromolecular associations. Sometimes both hydrophobic and electric interactions add up their strengths to accomplish these structures but in most cases they act in opposite directions. This fact, together with being overall interactions with different ranges, provides a nuanced equilibrium also modulated by the need to comply with steric hindrances and geometric frustration effects. This review focuses on the utility of using the hydrophobic and electrical dipole moment vectors to describe the interactions that give rise to the structures of biological macromolecules. Although different definitions of both electric dipole and hydrophobic moments have been described in the literature, results obtained in biological assemblies demonstrate the principle of the biological membrane model. According to this model, postulated by our group, biological macromolecules tend to associate by aligning their hydrophobic moments in a similar manner to phospholipids in a membrane. Examples of both closed and open structures are used to assess the predictability of our model. We seek agreement between our results with those described in the current literature. The review ends with possible future projections using this formalism.

Applications of Atomic Force Microscopy in HIV-1 Research

Obtaining an understanding of the mechanism underlying the interrelations between the structure and function of HIV-1 is of pivotal importance. In previous decades, this mechanism was addressed extensively in a variety of studies using conventional approaches. More recently, atomic force microscopy, which is a relatively new technique with unique capabilities, has been utilized to study HIV-1 biology. Atomic force microscopy can generate high-resolution images at the nanometer-scale and analyze the mechanical properties of individual HIV-1 virions, virus components (e.g., capsids), and infected live cells under near-physiological environments. This review describes the working principles and various imaging and analysis modes of atomic force microscopy, and elaborates on its distinctive contributions to HIV-1 research in areas such as mechanobiology and the physics of infection.

Modeling the Dynamics of Protein–Protein Interfaces, How and Why?

Protein–protein assemblies act as a key component in numerous cellular processes. Their accurate modeling at the atomic level remains a challenge for structural biology. To address this challenge, several docking and a handful of deep learning methodologies focus on modeling protein–protein interfaces. Although the outcome of these methods has been assessed using static reference structures, more and more data point to the fact that the interaction stability and specificity is encoded in the dynamics of these interfaces. Therefore, this dynamics information must be taken into account when modeling and assessing protein interactions at the atomistic scale. Expanding on this, our review initially focuses on the recent computational strategies aiming at investigating protein–protein interfaces in a dynamic fashion using enhanced sampling, multi-scale modeling, and experimental data integration. Then, we discuss how interface dynamics report on the function of protein assemblies in globular complexes, in fuzzy complexes containing intrinsically disordered proteins, as well as in active complexes, where chemical reactions take place across the protein–protein interface.

Structure-based virtual screening workflow to identify antivirals targeting HIV-1 capsid

We have identified novel HIV-1 capsid inhibitors targeting the PF74 binding site. Acting as the building block of the HIV-1 capsid core, the HIV-1 capsid protein plays an important role in the viral life cycle and is an attractive target for antiviral development. A structure-based virtual screening workflow for hit identification was employed, which includes docking 1.6 million commercially-available drug-like compounds from the ZINC database to the capsid dimer, followed by applying two absolute binding free energy (ABFE) filters on the 500 top-ranked molecules from docking. The first employs the Binding Energy Distribution Analysis Method (BEDAM) in implicit solvent. The top-ranked compounds are then refined using the Double Decoupling method in explicit solvent. Both docking and BEDAM refinement were carried out on the IBM World Community Grid as part of the FightAIDS@Home project. Using this virtual screening workflow, we identified 24 molecules with calculated binding free energies between − 6 and − 12 kcal/mol. We performed thermal shift assays on these molecules to examine their potential effects on the stability of HIV-1 capsid hexamer and found that two compounds, ZINC520357473 and ZINC4119064 increased the melting point of the latter by 14.8 °C and 33 °C, respectively. These results support the conclusion that the two ZINC compounds are primary hits targeting the capsid dimer interface. Our simulations also suggest that the two hit molecules may bind at the capsid dimer interface by occupying a new sub-pocket that has not been exploited by existing CA inhibitors. The possible causes for why other top-scored compounds suggested by ABFE filters failed to show measurable activity are discussed.

Natural Food Colorants and Preservatives: A Review, a Demand, and a Challenge

The looming urgency of feeding the growing world population along with the increasing consumers' awareness and expectations have driven the evolution of food production systems and the processes and products applied in the food industry. Although substantial progress has been made on food additives, the controversy in which some of them are still shrouded has encouraged research on safer and healthier next generations. These additives can come from natural sources and confer numerous benefits for health, beyond serving the purpose of coloring or preserving, among others. As limiting factors, these additives are often related to stability, sustainability, and cost-effectiveness issues, which justify the need for innovative solutions. In this context, and with the advances witnessed in computers and computational methodologies for in silico experimental aid, the development of new safer and more efficient natural additives with dual functionality (colorant and preservative), for instance by the copigmentation phenomena, may be achieved more efficiently, circumventing the current difficulties.

Strain and rupture of HIV-1 capsids during uncoating

SignificanceThe mature capsids of HIV-1 are transiently stable complexes that self-assemble around the viral genome during maturation, and uncoat to release preintegration complexes that archive a double-stranded DNA copy of the virus in the host cell genome. However, a detailed view of how HIV cores rupture remains lacking. Here, we elucidate the physical properties involved in capsid rupture using a combination of large-scale all-atom molecular dynamics simulations and cryo-electron tomography. We find that intrinsic strain on the capsid forms highly correlated patterns along the capsid surface, along which cracks propagate. Capsid rigidity also increases with high strain. Our findings provide fundamental insight into viral capsid uncoating.

Three-dimensional insights into human enveloped viruses in vitro and in situ

Viruses can be enveloped or non-enveloped, and require a host cell to replicate and package their genomes into new virions to infect new cells. To accomplish this task, viruses hijack the host-cell machinery to facilitate their replication by subverting and manipulating normal host cell function. Enveloped viruses can have severe consequences for human health, causing various diseases such as acquired immunodeficiency syndrome (AIDS), seasonal influenza, COVID-19, and Ebola virus disease. The complex arrangement and pleomorphic architecture of many enveloped viruses pose a challenge for the more widely used structural biology techniques, such as X-ray crystallography. Cryo-electron tomography (cryo-ET), however, is a particularly well-suited tool for overcoming the limitations associated with visualizing the irregular shapes and morphology enveloped viruses possess at macromolecular resolution. The purpose of this review is to explore the latest structural insights that cryo-ET has revealed about enveloped viruses, with particular attention given to their architectures, mechanisms of entry, replication, assembly, maturation and egress during infection. Cryo-ET is unique in its ability to visualize cellular landscapes at 3-5 nanometer resolution. Therefore, it is the most suited technique to study asymmetric elements and structural rearrangements of enveloped viruses during infection in their native cellular context.

Advances in HIV-1 Assembly

The assembly of HIV-1 particles is a concerted and dynamic process that takes place on the plasma membrane of infected cells. An abundance of recent discoveries has advanced our understanding of the complex sequence of events leading to HIV-1 particle assembly, budding, and release. Structural studies have illuminated key features of assembly and maturation, including the dramatic structural transition that occurs between the immature Gag lattice and the formation of the mature viral capsid core. The critical role of inositol hexakisphosphate (IP6) in the assembly of both the immature and mature Gag lattice has been elucidated. The structural basis for selective packaging of genomic RNA into virions has been revealed. This review will provide an overview of the HIV-1 assembly process, with a focus on recent advances in the field, and will point out areas where questions remain that can benefit from future investigation.

Exploring cryo-electron microscopy with molecular dynamics

Single particle analysis cryo-electron microscopy (EM) and molecular dynamics (MD) have been complimentary methods since cryo-EM was first applied to the field of structural biology. The relationship started by biasing structural models to fit low-resolution cryo-EM maps of large macromolecular complexes not amenable to crystallization. The connection between cryo-EM and MD evolved as cryo-EM maps improved in resolution, allowing advanced sampling algorithms to simultaneously refine backbone and sidechains. Moving beyond a single static snapshot, modern inferencing approaches integrate cryo-EM and MD to generate structural ensembles from cryo-EM map data or directly from the particle images themselves. We summarize the recent history of MD innovations in the area of cryo-EM modeling. The merits for the myriad of MD based cryo-EM modeling methods are discussed, as well as, the discoveries that were made possible by the integration of molecular modeling with cryo-EM. Lastly, current challenges and potential opportunities are reviewed.

Assembly and Analysis of Cell-Scale Membrane Envelopes

The march toward exascale computing will enable routine molecular simulation of larger and more complex systems, for example, simulation of entire viral particles, on the scale of approximately billions of atoms─a simulation size commensurate with a small bacterial cell. Anticipating the future hardware capabilities that will enable this type of research and paralleling advances in experimental structural biology, efforts are currently underway to develop software tools, procedures, and workflows for constructing cell-scale structures. Herein, we describe our efforts in developing and implementing an efficient and robust workflow for construction of cell-scale membrane envelopes and embedding membrane proteins into them. A new approach for construction of massive membrane structures that are stable during the simulations is built on implementing a subtractive assembly technique coupled with the development of a structure concatenation tool (fastmerge), which eliminates overlapping elements based on volumetric criteria rather than adding successive molecules to the simulation system. Using this approach, we have constructed two "protocells" consisting of MARTINI coarse-grained beads to represent cellular membranes, one the size of a cellular organelle and another the size of a small bacterial cell. The membrane envelopes constructed here remain whole during the molecular dynamics simulations performed and exhibit water flux only through specific proteins, demonstrating the success of our methodology in creating tight cell-like membrane compartments. Extended simulations of these cell-scale structures highlight the propensity for nonspecific interactions between adjacent membrane proteins leading to the formation of protein microclusters on the cell surface, an insight uniquely enabled by the scale of the simulations. We anticipate that the experiences and best practices presented here will form the basis for the next generation of cell-scale models, which will begin to address the addition of soluble proteins, nucleic acids, and small molecules essential to the function of a cell.

Molecular basis of the potential interaction of SARS-CoV-2 spike protein to CD147 in COVID-19 associated-lymphopenia

Lymphopenia is considered one of the most characteristic clinical features of the coronavirus disease 2019 (COVID-19). SARS-CoV-2 infects host cells via the interaction of its spike protein with the human angiotensin-converting enzyme 2 (hACE2) receptor. Since T lymphocytes display a very low expression level of hACE2, a novel receptor might be involved in the entry of SARS-CoV-2 into T cells. The transmembrane glycoprotein CD147 is highly expressed by activated T lymphocytes, and was recently proposed as a probable route for SARS-CoV-2 invasion. To understand the molecular basis of the potential interaction of SARS-CoV-2 to CD147, we have investigated the binding of the viral spike protein to this receptor in-silico. The results showed that this binding is dominated by electrostatic interactions involving residues Arg403, Asn481, and the backbone of Gly502. The overall binding arrangement shows the CD147 C-terminal domain interacting with the spike external subdomain in the grove between the short antiparallel β strands, β1' and β2', and the small helix α1'. This proposed interaction was further confirmed using MD simulation and binding free energy calculation. These data contribute to a better understanding of the mechanism of infection of SARS-CoV-2 to T lymphocytes and could provide valuable insights for the rational design of adjuvant treatment for COVID-19. Communicated by Ramaswamy H. Sarma.