Concerted Interactions between Multiple gp41 Trimers and the Target Cell Lipidome May Be Required for HIV-1 Entry

Molecular dynamics simulations show how antibodies may rescue HIV-1 mutants incapable of infecting host cells

High mutation and replication rates of HIV-1 result in the continuous generation of variants, allowing it to adapt to changing host environments. Mutations often have deleterious effects, but variants carrying them are rapidly purged. Surprisingly, a particular variant incapable of entering host cells was found to be rescued by host antibodies targeting HIV-1. Understanding the molecular mechanism of this rescue is important to develop and improve antibody-based therapies. To unravel the underlying mechanisms, we performed fully atomistic molecular dynamics simulations of the HIV-1 gp41 trimer responsible for viral entry into host cells, its entry-deficient variant, and its complex with the rescuing antibody. We find that the Q563R mutation, which the entry-deficient variant carries, prevents the native conformation of the gp41 6-helix bundle required for entry and stabilizes an alternative conformation instead. This is the consequence of substantial changes in the secondary structure and interactions between the domains of gp41. Binding of the antibody F240 to gp41 reverses these changes and re-establishes the native conformation, resulting in rescue. To test the generality of this mechanism, we performed simulations with the entry-deficient L565A variant and antibody 3D6. We find that 3D6 binding was able to reverse structural and interaction changes introduced by the mutation and restore the native gp41 conformation. Viral variants may not only escape antibodies but be aided by them in their survival, potentially compromising antibody-based therapies, including vaccination and passive immunization. Our simulation framework could serve as a tool to assess the likelihood of such resistance against specific antibodies.

Sequence of the SARS-CoV-2 Spike Transmembrane Domain Encodes Conformational Dynamics

The homotrimeric SARS-CoV-2 spike protein enables viral infection by undergoing a large conformational transition, which facilitates the fusion of the viral envelope with the host cell membrane. The spike protein is anchored to the SARS-CoV-2 envelope by its transmembrane domain (TMD), composed of three TM helices, each contributed by one of the protomers of spike. Although the TMD is known to be important for viral fusion, whether it is a passive anchor of the spike or actively promotes fusion remains unknown. Specifically, it is unclear if the TMD and its dynamics facilitate the prefusion to postfusion conformational transition of the spike. Here, we computationally study the dynamics and self-assembly of the SARS-CoV-2 spike TMD in homogeneous POPC and cholesterol containing membranes. Atomistic simulations of a long TM helix-containing protomer segment show that the membrane-embedded segment bobs, tilts and gains and loses helicity, locally thinning the membrane. Coarse-grained multimerization simulations using representative TM helix structures from the atomistic simulations exhibit diverse trimer populations whose architecture depends on the structure of the TM helix protomer. While a symmetric conformation reflects the symmetry of the resting spike, an asymmetric TMD conformation could promote membrane fusion through the stabilization of a fusion intermediate. Together, our simulations demonstrate that the sequence and length of the SARS-CoV-2 spike TM segment make it inherently dynamic, that trimerization does not abrogate these dynamics and that the various observed TMD conformations may enable viral fusion.

Insertion and Anchoring of the HIV-1 Fusion Peptide into a Complex Membrane Mimicking the Human T-Cell

A fundamental understanding of how the HIV-1 envelope (Env) protein facilitates fusion is still lacking. The HIV-1 fusion peptide, consisting of 15 to 22 residues, is the N-terminus of the gp41 subunit of the Env protein. Further, this peptide, a promising vaccine candidate, initiates viral entry into target cells by inserting and anchoring into human immune cells. The influence of membrane lipid reorganization and the conformational changes of the fusion peptide during the membrane insertion and anchoring processes, which can significantly affect HIV-1 cell entry, remains largely unexplored due to the limitations of experimental measurements. In this work, we investigate the insertion of the fusion peptide into an immune cell membrane mimic through multiscale molecular dynamics simulations. We mimic the native T-cell by constructing a nine-lipid asymmetric membrane, along with geometrical restraints accounting for insertion in the context of gp41. To account for the slow time scale of lipid mixing while enabling conformational changes, we implement a protocol to go back and forth between atomistic and coarse-grained simulations. Our study provides a molecular understanding of the interactions between the HIV-1 fusion peptide and the T-cell membrane, highlighting the importance of the conformational flexibility of fusion peptides and local lipid reorganization in stabilizing the anchoring of gp41 into the targeted host membrane during the early events of HIV-1 cell entry. Importantly, we identify a motif within the fusion peptide critical for fusion that can be further manipulated in future immunological studies.

Spontaneous Dimerization and Distinct Packing Modes of Transmembrane Domains in Receptor Tyrosine Kinases

The insulin receptor (IR) and the insulin-like growth factor-1 receptor (IGF1R) are homodimeric transmembrane glycoproteins that transduce signals across the membrane on binding of extracellular peptide ligands. The structures of IR/IGF1R fragments in apo and liganded states have revealed that the extracellular subunits of these receptors adopt Λ-shaped configurations to which are connected the intracellular tyrosine kinase (TK) domains. The binding of peptide ligands induces structural transitions in the extracellular subunits leading to potential dimerization of transmembrane domains (TMDs) and autophosphorylation in TKs. However, the activation mechanisms of IR/IGF1R, especially the role of TMDs in coordinating signal-inducing structural transitions, remain poorly understood, in part due to the lack of structures of full-length receptors in apo or liganded states. While atomistic simulations of IR/IGF1R TMDs showed that these domains can dimerize in single component membranes, spontaneous unbiased dimerization in a plasma membrane having a physiologically representative lipid composition has not been observed. We address this limitation by employing coarse-grained (CG) molecular dynamics simulations to probe the dimerization propensity of IR/IGF1R TMDs. We observed that TMDs in both receptors spontaneously dimerized independent of their initial orientations in their dissociated states, signifying their natural propensity for dimerization. In the dimeric state, IR TMDs predominantly adopted X-shaped configurations with asymmetric helical packing and significant tilt relative to the membrane normal, while IGF1R TMDs adopted symmetric V-shaped or parallel configurations with either no tilt or a small tilt relative to the membrane normal. Our results suggest that IR/IGF1R TMDs spontaneously dimerize and adopt distinct dimerized configurations.

Membrane lateral organization from potential energy disconnectivity graph

Understanding the thermodynamic and kinetic properties of biomolecules requires elucidation of their complex energy landscape. A disconnectivity graph analysis of the energy landscape provides a framework for mapping the multi-dimensional landscape onto a two-dimensional representation while preserving the key features of the energy landscape. Several studies show that the structure or shape of the disconnectity graph is directly associated with the function of protein and nucleic acid molecules. In this review, we discuss how disconnectivity analysis of the potential energy surface can be extended to lipid molecules to glean important information about membrane organization. The shape of the disconnectivity graphs can be used to predict the lateral organization of multi-component lipid bilayer. We hope that this review encourages the use of disconnectivity graphs routinely by membrane biophysicists to predict the lateral organization of lipids.

Insertion and Anchoring of HIV-1 Fusion Peptide into Complex Membrane Mimicking Human T-cell

A fundamental understanding of how HIV-1 envelope (Env) protein facilitates fusion is still lacking. The HIV-1 fusion peptide, consisting of 15 to 22 residues, is the N-terminus of the gp41 subunit of the Env protein. Further, this peptide, a promising vaccine candidate, initiates viral entry into target cells by inserting and anchoring into human immune cells. The influence of membrane lipid reorganization and the conformational changes of the fusion peptide during the membrane insertion and anchoring processes, which can significantly affect HIV-1 cell entry, remains largely unexplored due to the limitations of experimental measurements. In this work, we investigate the insertion of the fusion peptide into an immune cell membrane mimic through multiscale molecular dynamics simulations. We mimic the native T-cell by constructing a 9-lipid asymmetric membrane, along with geometrical restraints accounting for insertion in the context of gp41. To account for the slow timescale of lipid mixing while enabling conformational changes, we implement a protocol to go back and forth between atomistic and coarse-grained simulations. Our study provides a molecular understanding of the interactions between the HIV-1 fusion peptide and the T-cell membrane, highlighting the importance of conformational flexibility of fusion peptides and local lipid reorganization in stabilizing the anchoring of gp41 into the targeted host membrane during the early events of HIV-1 cell entry. Importantly, we identify a motif within the fusion peptide critical for fusion that can be further manipulated in future immunological studies.

Mathematical model for force and energy of virion-cell interactions during full engulfment in HIV: Impact of virion maturation and host cell morphology

Viral endocytosis involves elastic cell deformation, driven by chemical adhesion energy, and depends on physical interactions between the virion and cell membrane. These interactions are not easy to quantify experimentally. Hence, this study aimed to develop a mathematical model of the interactions of HIV particles with host cells and explore the effects of mechanical and morphological parameters during full virion engulfment. The invagination force and engulfment energy were described as viscoelastic and linear-elastic functions of radius and elastic modulus of virion and cell, ligand-receptor energy density and engulfment depth. The influence of changes in the virion-cell contact geometry representing different immune cells and ultrastructural membrane features and the decrease in virion radius and shedding of gp120 proteins during maturation on invagination force and engulfment energy was investigated. A low invagination force and high ligand-receptor energy are associated with high virion entry ability. The required invagination force was the same for immune cells of different sizes but lower for a local convex geometry of the cell membrane at the virion length scale. This suggests that localized membrane features of immune cells play a role in viral entry ability. The available engulfment energy decreased during virion maturation, indicating the involvement of additional biological or biochemical changes in viral entry. The developed mathematical model offers potential for the mechanobiological assessment of the invagination of enveloped viruses towards improving the prevention and treatment of viral infections.

Enhanced Expansion and Reduced Kiss-and-Run Events in Fusion Pores Steered by Synaptotagmin-1 C2B Domains

The fusion pore controls the release of exocytotic vesicle contents through a precise orchestration of lipids from the fusing membranes and proteins. There is a major lipid reorganization during the different stages in life of the fusion pore (membrane fusion, nucleation, and expansion) that can be scrutinized thermodynamically. In this work, using umbrella sampling simulations we describe the expansion of the fusion pore. We have calculated free energy profiles to drive a nascent, just nucleated, fusion pore to its expanded configuration. We have quantified the effects on the free energy of one and two Synaptotagmin-1 C2B domains in the cytosolic space. We show that C2B domains cumulatively reduce the cost for expansion, favoring the system to evolve toward full fusion. Finally, by conducting thousands of unbiased molecular dynamics simulations, we show that C2B domains significantly decrease the probability of kiss-and-run events.

Role of Protein-Lipid Interactions in Viral Entry

The viral entry consists of several sequential events that ensure the attachment of the virus to the host cell and the introduction of its genetic material for the continuation of the replication cycle. Both cellular and viral lipids have gained a wider focus in recent years in the field of viral entry, as they are found to play key roles in different steps of the process. The specific role is summarized that lipids and lipid membrane nanostructures play in viral attachment, fusion, and immune evasion and how they can be targeted with antiviral therapies. Finally, some of the limitations of techniques commonly used for protein-lipid interactions studies are discussed, and new emerging tools are reviewed that can be applied to this field.

Cholesterol-Mediated Clustering of the HIV Fusion Protein gp41 in Lipid Bilayers

The envelope glycoprotein (Env) of the human immunodeficient virus (HIV-1) is known to cluster on the viral membrane surface to attach to target cells and cause membrane fusion for HIV-1 infection. However, the molecular structural mechanisms that drive Env clustering remain opaque. Here, we use solid-state NMR spectroscopy and molecular dynamics (MD) simulations to investigate nanometer-scale clustering of the membrane-proximal external region (MPER) and transmembrane domain (TMD) of gp41, the fusion protein component of Env. Using ¹⁹F solid-state NMR experiments of mixed fluorinated peptides, we show that MPER-TMD trimers form clusters with interdigitated MPER helices in cholesterol-containing membranes. Inter-trimer ¹⁹F-¹⁹F cross peaks, which are indicative of spatial contacts within ∼2 nm, are observed in cholesterol-rich virus-mimetic membranes but are suppressed in cholesterol-free model membranes. Water-peptide and lipid-peptide cross peaks in 2D ¹H-¹⁹F correlation spectra indicate that the MPER is well embedded in model phosphocholine membranes but is more exposed to the surface of the virus-mimetic membrane. These experimental results are reproduced in coarse-grained and atomistic molecular dynamics simulations, which suggest that the effects of cholesterol on gp41 clustering is likely via indirect modulation of the MPER orientation. Cholesterol binding to the helix-turn-helix region of the MPER-TMD causes a parallel orientation of the MPER with the membrane surface, thus allowing MPERs of neighboring trimers to interact with each other to cause clustering. These solid-state NMR data and molecular dynamics simulations suggest that MPER and cholesterol cooperatively govern the clustering of gp41 trimers during virus-cell membrane fusion.