All-atom models of the membrane-spanning domain of HIV-1 gp41 from metadynamics

Viral Membrane Fusion and the Transmembrane Domain

Initiation of host cell infection by an enveloped virus requires a viral-to-host cell membrane fusion event. This event is mediated by at least one viral transmembrane glycoprotein, termed the fusion protein, which is a key therapeutic target. Viral fusion proteins have been studied for decades, and numerous critical insights into their function have been elucidated. However, the transmembrane region remains one of the most poorly understood facets of these proteins. In the past ten years, the field has made significant advances in understanding the role of the membrane-spanning region of viral fusion proteins. We summarize developments made in the past decade that have contributed to the understanding of the transmembrane region of viral fusion proteins, highlighting not only their critical role in the membrane fusion process, but further demonstrating their involvement in several aspects of the viral lifecycle.

HIV-1 Env gp41 Transmembrane Domain Dynamics Are Modulated by Lipid, Water, and Ion Interactions

The gp41 transmembrane domain (TMD) of the envelope glycoprotein of the human immunodeficiency virus modulates the conformation of the viral envelope spike, the only druggable target on the surface of the virion. Targeting the envelope glycoprotein with small-molecule and antibody therapies requires an understanding of gp41 TMD dynamics, which is often challenging given the difficulties in describing native membrane properties. Here, atomistic molecular dynamics simulations of a trimeric, prefusion gp41 TMD in a model, asymmetric viral membrane that mimics the native viral envelope were performed. Water and chloride ions were observed to permeate the membrane and interact with the highly conserved arginine bundle, (R696)₃, at the center of the membrane and influenced TMD stability by creating a network of hydrogen bonds and electrostatic interactions. We propose that this (R696)₃ - water - anion network plays an important role in viral fusion with the host cell by modulating protein conformational changes within the membrane. Additionally, R683 and R707 at the exofacial and cytofacial membrane-water interfaces, respectively, are anchored in the lipid headgroup region and serve as a junction point for stabilization of the termini. The membrane thins as a result of the tilting of the gp41 trimer with nearby lipids increasing in volume, leading to an entropic driving force for TMD conformational change. These results provide additional detail and perspective on the influence of certain lipid types on TMD dynamics and a rationale for targeting key residues of the TMD for therapeutic design. These insights into the molecular details of TMD membrane anchoring will build toward a greater understanding of the dynamics that lead to viral fusion with the host cell.

Characterization of Lipid-Protein Interactions and Lipid-Mediated Modulation of Membrane Protein Function through Molecular Simulation

The cellular membrane constitutes one of the most fundamental compartments of a living cell, where key processes such as selective transport of material and exchange of information between the cell and its environment are mediated by proteins that are closely associated with the membrane. The heterogeneity of lipid composition of biological membranes and the effect of lipid molecules on the structure, dynamics, and function of membrane proteins are now widely recognized. Characterization of these functionally important lipid-protein interactions with experimental techniques is however still prohibitively challenging. Molecular dynamics (MD) simulations offer a powerful complementary approach with sufficient temporal and spatial resolutions to gain atomic-level structural information and energetics on lipid-protein interactions. In this review, we aim to provide a broad survey of MD simulations focusing on exploring lipid-protein interactions and characterizing lipid-modulated protein structure and dynamics that have been successful in providing novel insight into the mechanism of membrane protein function.

All-atom virus simulations

The constant threat of viral disease can be combated by the development of novel vaccines and therapeutics designed to disrupt key features of virus structure or infection cycle processes. Such development relies on high-resolution characterization of viruses and their dynamical behaviors, which are often challenging to obtain solely by experiment. In response, all-atom molecular dynamics simulations are widely leveraged to study the structural components of viruses, leading to some of the largest simulation endeavors undertaken to date. The present work reviews exemplary all-atom simulation work on viruses, as well as progress toward simulating entire virions.

All-Atom Structural Models of the Transmembrane Domains of Insulin and Type 1 Insulin-Like Growth Factor Receptors

The receptor tyrosine kinase superfamily comprises many cell-surface receptors including the insulin receptor (IR) and type 1 insulin-like growth factor receptor (IGF1R) that are constitutively homodimeric transmembrane glycoproteins. Therefore, these receptors require ligand-triggered domain rearrangements rather than receptor dimerization for activation. Specifically, binding of peptide ligands to receptor ectodomains transduces signals across the transmembrane domains for trans-autophosphorylation in cytoplasmic kinase domains. The molecular details of these processes are poorly understood in part due to the absence of structures of full-length receptors. Using MD simulations and enhanced conformational sampling algorithms, we present all-atom structural models of peptides containing 51 residues from the transmembrane and juxtamembrane regions of IR and IGF1R. In our models, the transmembrane regions of both receptors adopt helical conformations with kinks at Pro961 (IR) and Pro941 (IGF1R), but the C-terminal residues corresponding to the juxtamembrane region of each receptor adopt unfolded and flexible conformations in IR as opposed to a helix in IGF1R. We also observe that the N-terminal residues in IR form a kinked-helix sitting at the membrane-solvent interface, while homologous residues in IGF1R are unfolded and flexible. These conformational differences result in a larger tilt-angle of the membrane-embedded helix in IGF1R in comparison to IR to compensate for interactions with water molecules at the membrane-solvent interfaces. Our metastable/stable states for the transmembrane domain of IR, observed in a lipid bilayer, are consistent with a known NMR structure of this domain determined in detergent micelles, and similar states in IGF1R are consistent with a previously reported model of the dimerized transmembrane domains of IGF1R. Our all-atom structural models suggest potentially unique structural organization of kinase domains in each receptor.

The Atomic Structure of the HIV-1 gp41 Transmembrane Domain and Its Connection to the Immunogenic Membrane-proximal External Region

The membrane-proximal external region (MPER) C-terminal segment and the transmembrane domain (TMD) of gp41 are involved in HIV-1 envelope glycoprotein-mediated fusion and modulation of immune responses during viral infection. However, the atomic structure of this functional region remains unsolved. Here, based on the high resolution NMR data obtained for peptides spanning the C-terminal segment of MPER and the TMD, we report two main findings: (i) the conformational variability of the TMD helix at a membrane-buried position; and (ii) the existence of an uninterrupted α-helix spanning MPER and the N-terminal region of the TMD. Thus, our structural data provide evidence for the bipartite organization of TMD predicted by previous molecular dynamics simulations and functional studies, but they do not support the breaking of the helix at Lys-683, as was suggested by some models to mark the initiation of the TMD anchor. Antibody binding energetics examined with isothermal titration calorimetry and humoral responses elicited in rabbits by peptide-based vaccines further support the relevance of a continuous MPER-TMD helix for immune recognition. We conclude that the transmembrane anchor of HIV-1 envelope is composed of two distinct subdomains: 1) an immunogenic helix at the N terminus also involved in promoting membrane fusion; and 2) an immunosuppressive helix at the C terminus, which might also contribute to the late stages of the fusion process. The unprecedented high resolution structural data reported here may guide future vaccine and inhibitor developments.

Fusion-competent state induced by a C-terminal HIV-1 fusion peptide in cholesterol-rich membranes

The replicative cycle of the human immunodeficiency virus type-1 begins after fusion of the viral and target-cell membranes. The envelope glycoprotein gp41 transmembrane subunit contains conserved hydrophobic domains that engage and perturb the merging lipid bilayers. In this work, we have characterized the fusion-committed state generated in vesicles by CpreTM, a synthetic peptide derived from the sequence connecting the membrane-proximal external region (MPER) and the transmembrane domain (TMD) of gp41. Pre-loading cholesterol-rich vesicles with CpreTM rendered them competent for subsequent lipid-mixing with fluorescently-labeled target vesicles. Highlighting the physiological relevance of the lasting fusion-competent state, the broadly neutralizing antibody 4E10 bound to the CpreTM-primed vesicles and inhibited lipid-mixing. Heterotypic fusion assays disclosed dependence on the lipid composition of the vesicles that acted either as virus or cell membrane surrogates. Lipid-mixing exhibited above all a critical dependence on the cholesterol content in those experiments. We infer that the fusion-competent state described herein resembles bona-fide perturbations generated by the pre-hairpin MPER-TMD connection within the viral membrane.

Dynamics of lipids, cholesterol, and transmembrane α-helices from microsecond molecular dynamics simulations

Extensive all-atom molecular dynamics (∼24 μs total) allowed exploration of configurational space and calculation of lateral diffusion coefficients of the components of a protein-embedded, cholesterol-containing model bilayer. The three model membranes are composed of an ∼50/50 (by mole) dipalmitoylphosphatidylcholine (DPPC)/cholesterol bilayer and contained an α-helical transmembrane protein (HIV-1 gp41 TM). Despite the high concentration of cholesterol, normal Brownian motion was observed and the calculated diffusion coefficients (on the order of 10(-9) cm(2)/s) are consistent with experiments. Diffusion is sensitive to a variety of parameters, and a temperature difference of ∼4 K from thermostat artifacts resulted in 2-10-fold differences in diffusion coefficients and significant differences in lipid order, membrane thickness, and unit cell area. Also, the specific peptide sequence likely underlies the consistently observed faster diffusion in one leaflet. Although the simulations here present molecular dynamics (MD) an order of magnitude longer than those from previous studies, the three systems did not approach ergodicity. The distributions of cholesterol and DPPC around the peptides changed on the microsecond time scale, but not significantly enough to thoroughly explore configurational space. These simulations support conclusions of other recent microsecond MD in that even longer time scales are needed for equilibration of model membranes and simulations of more realistic cellular or viral bilayers.

Cholesterol-dependent membrane fusion induced by the gp41 membrane-proximal external region-transmembrane domain connection suggests a mechanism for broad HIV-1 neutralization

The HIV-1 glycoprotein 41 promotes fusion of the viral membrane with that of the target cell. Structural, biochemical, and biophysical studies suggest that its membrane-proximal external region (MPER) may interact with the HIV-1 membrane and induce its disruption and/or deformation during the process. However, the high cholesterol content of the envelope (ca. 40 to 50 mol%) imparts high rigidity, thereby acting against lipid bilayer restructuring. Here, based on the outcome of vesicle stability assays, all-atom molecular dynamics simulations, and atomic force microscopy observations, we propose that the conserved sequence connecting the MPER with the N-terminal residues of the transmembrane domain (TMD) is involved in HIV-1 fusion. This junction would function by inducing phospholipid protrusion and acyl-chain splay in the cholesterol-enriched rigid envelope. Supporting the functional relevance of such a mechanism, membrane fusion was inhibited by the broadly neutralizing 4E10 antibody but not by a nonneutralizing variant with the CDR-H3 loop deleted. We conclude that the MPER-TMD junction embodies an envelope-disrupting C-terminal fusion peptide that can be targeted by broadly neutralizing antibodies.

IMPORTANCE

Fusion of the cholesterol-enriched viral envelope with the cell membrane marks the beginning of the infectious HIV-1 replicative cycle. Consequently, the Env glycoprotein-mediated fusion function constitutes an important clinical target for inhibitors and preventive vaccines. Antibodies 4E10 and 10E8 bind to one Env vulnerability site located at the gp41 membrane-proximal external region (MPER)-transmembrane domain (TMD) junction and block infection. These antibodies display broad viral neutralization, which underscores the conservation and functionality of the MPER-TMD region. In this work, we combined biochemical assays with molecular dynamics simulations and microscopy observations to characterize the unprecedented fusogenic activity of the MPER-TMD junction. The fact that such activity is dependent on cholesterol and inhibited by the broadly neutralizing 4E10 antibody emphasizes its physiological relevance. Discovery of this functional element adds to our understanding of the mechanisms underlying HIV-1 infection and its blocking by antibodies.

Characterization of the water defect at the HIV-1 gp41 membrane spanning domain in bilayers with and without cholesterol using molecular simulations

The membrane spanning domain (MSD) of human immunodeficiency virus 1 (HIV-1) envelope glycoprotein gp41 is important for fusion and infection. We used molecular dynamics (MD) simulations (3.4 μs total) to relate membrane and peptide properties that lead to water solvation of the α-helical gp41 MSD's midspan arginine in pure dipalmitoylphosphatidylcholine (DPPC) and in 50/50 DPPC/cholesterol membranes. We find that the midspan arginine is solvated by water that penetrates the inner leaflet, leading to a so-called water defect. The water defect is surprisingly robust across initial conditions and membrane compositions, but the presence of cholesterol modulates its behavior in several key ways. In the cholesterol-containing membranes, fluctuations in membrane thickness and water penetration depth are localized near the midspan arginine, and the MSD helices display a tightly regulated tilt angle. In the cholesterol-free membranes, thickness fluctuations are not as strongly correlated to the peptide position and tilt angles vary significantly depending on protein position relative to boundaries between domains of differing thickness. Cholesterol in an HIV-1 viral membrane is required for infection. Therefore, this work suggests that the colocalized water defect and membrane thickness fluctuations in cholesterol-containing viral membranes play an important role in fusion by bringing the membrane closer to a stability limit that must be crossed for fusion to occur.

C-terminal tail of human immunodeficiency virus gp41: functionally rich and structurally enigmatic

The human immunodeficiency virus (HIV) and acquired immunodeficiency syndrome (AIDS) pandemic is amongst the most important current worldwide public health threats. While much research has been focused on AIDS vaccines that target the surface viral envelope (Env) protein, including gp120 and the gp41 ectodomain, the C-terminal tail (CTT) of gp41 has received relatively little attention. Despite early studies highlighting the immunogenicity of a particular CTT sequence, the CTT has been classically portrayed as a type I membrane protein limited to functioning in Env trafficking and virion incorporation. Recent studies demonstrate, however, that the Env CTT has other important functions. The CTT has been shown to additionally modulate Env ectodomain structure on the cell and virion surface, affect Env reactivity and viral sensitivity to conformation-dependent neutralizing antibodies, and alter cell–cell and virus–cell fusogenicity of Env. This review provides an overview of the Env structure and function with a particular emphasis on the CTT and recent studies that highlight its functionally rich nature.

Selected amino acid changes in HIV-1 subtype-C gp41 are associated with specific gp120(V3) signatures in the regulation of co-receptor usage

The majority of studies have characterized the tropism of HIV-1 subtype-B isolates, but little is known about the determinants of tropism in other subtypes. So, the goal of the present study was to genetically characterize the envelope of viral proteins in terms of co-receptor usage by analyzing 356 full-length env sequences derived from HIV-1 subtype-C infected individuals. The co-receptor usage of V3 sequences was inferred by using the Geno2Pheno and PSSM algorithms, and also analyzed to the "11/25 rule". All reported env sequences were also analyzed with regard to N-linked glycosylation sites, net charge and hydrophilicity, as well as the binomial correlation phi coefficient to assess covariation among gp120(V3) and gp41 signatures and the average linkage hierarchical agglomerative clustering were also performed. Among env sequences present in Los Alamos Database, 255 and 101 sequences predicted as CCR5 and CXCR4 were selected, respectively. The classical V3 signatures at positions 11 and 25, and other specific V3 and gp41 amino acid changes were found statistically associated with different co-receptor usage. Furthermore, several statistically significant associations between V3 and gp41 signatures were also observed. The dendrogram topology showed a cluster associated with CCR5-usage composed by five gp41 mutated positions, A22V, R133M, E136G, N140L, and N166Q that clustered with T2V(V3) and G24T(V3) (bootstrap=1). Conversely, a heterogeneous cluster with CXCR4-usage, involving S11GR(V3), 13-14insIG/LG(V3), P16RQ(V3), Q18KR(V3), F20ILV(V3), D25KRQ(V3), Q32KR(V3) along with A30T(gp41), S107N(gp41), D148E(gp41), A189S(gp41) was identified (bootstrap=0.86). Our results show that as observed for HIV-1 subtype-B, also in subtype-C specific and different gp41 and gp120V3 amino acid changes are associated individually or together with CXCR4 and/or CCR5 usage. These findings strengthen previous observations that determinants of tropism may also reside in the gp41 protein.

Evidence against extracellular exposure of a highly immunogenic region in the C-terminal domain of the simian immunodeficiency virus gp41 transmembrane protein

The generally accepted model for human immunodeficiency virus type 1 (HIV-1) envelope glycoprotein topology includes a single membrane-spanning domain. An alternate model has been proposed which features multiple membrane-spanning domains. Consistent with the alternate model, a high percentage of HIV-1-infected individuals produce unusually robust antibody responses to a region of envelope, the so-called "Kennedy epitope," that in the conventional model should be in the cytoplasm. Here we show analogous, robust antibody responses in simian immunodeficiency virus SIVmac239-infected rhesus macaques to a region of SIVmac239 envelope located in the C-terminal domain, which in the conventional model should be inside the cell. Sera from SIV-infected rhesus macaques consistently reacted with overlapping oligopeptides corresponding to a region located within the cytoplasmic domain of gp41 by the generally accepted model, at intensities comparable to those observed for immunodominant areas of the surface component gp120. Rabbit serum raised against this highly immunogenic region (HIR) reacted with SIV envelope in cell surface-staining experiments, as did monoclonal anti-HIR antibodies isolated from an SIVmac239-infected rhesus macaque. However, control experiments demonstrated that this surface staining could be explained in whole or in part by the release of envelope protein from expressing cells into the supernatant and the subsequent attachment to the surfaces of cells in the culture. Serum and monoclonal antibodies directed against the HIR failed to neutralize even the highly neutralization-sensitive strain SIVmac316. Furthermore, a potential N-linked glycosylation site located close to the HIR and postulated to be outside the cell in the alternate model was not glycosylated. An artificially introduced glycosylation site within the HIR was also not utilized for glycosylation. Together, these data support the conventional model of SIV envelope as a type Ia transmembrane protein with a single membrane-spanning domain and without any extracellular loops.

Structural dynamics of the S4 voltage-sensor helix in lipid bilayers lacking phosphate groups

Voltage-dependent K(+) (Kv) channels require lipid phosphates for functioning. The S4 helix, which carries the gating charges in the voltage-sensing domain (VSD), inserts into membranes while being stabilized by a protein-lipid interface in which lipid phosphates play an essential role. To examine the physical basis of the protein-lipid interface in the absence of lipid phosphates, we performed molecular dynamics (MD) simulations of a KvAP S4 variant (S4mut) in bilayers with and without lipid phosphates. We find that, in dioleoyltrimethylammoniumpropane (DOTAP) bilayers lacking lipid phosphates, the gating charges are solvated by anionic counterions and, hence, lack the bilayer support provided by phosphate-containing palmitoyloleoylglycerophosphocholine (POPC) bilayers. The result is a water-permeable bilayer with significantly smaller deformations around the peptide. Together, these results provide an explanation for the nonfunctionality of VSDs in terms of a destabilizing protein-lipid interface.

Conserved arginine residue in the membrane-spanning domain of HIV-1 gp41 is required for efficient membrane fusion

Despite the high mutation rate of HIV-1, the amino acid sequences of the membrane-spanning domain (MSD) of HIV-1 gp41 are well conserved. Arginine residues are rarely found in single membrane-spanning domains, yet an arginine residue, R(696) (the numbering is based on that of HXB2), is highly conserved in HIV-1 gp41. To examine the role of R(696), it was mutated to K, A, I, L, D, E, N, and Q. Most of these substitutions did not affect the expression, processing or surface distribution of the envelope protein (Env). However, a syncytia formation assay showed that the substitution of R(696) with amino acid residues other than K, a naturally observed mutation in the gp41 MSD, decreased fusion activity. Substitution with hydrophobic amino acid residues (A, I, and L) resulted in a modest decrease, while substitution with D or E, potentially negatively-charged residues, almost abolished the syncytia formation. All the fusion-defective mutants showed slower kinetics with the cell-based dual split protein (DSP) assay that scores the degree of membrane fusion based on pore formation between fusing cells. Interestingly, the D and E substitutions did show some fusion activity in the DSP assays, suggesting that proteins containing D or E substitutions retained some fusion pore-forming capability. However, nascent pores failed to develop, due probably to impaired activity in the pore enlargement process. Our data show the importance of this conserved arginine residue for efficient membrane fusion.

Highly conserved structural properties of the C-terminal tail of HIV-1 gp41 protein despite substantial sequence variation among diverse clades: implications for functions in viral replication

Although the HIV-1 Env gp120 and gp41 ectodomain have been extensively characterized in terms of structure and function, similar characterizations of the C-terminal tail (CTT) of HIV gp41 remain relatively limited and contradictory. The current study was designed to examine in detail CTT sequence conservation relative to gp120 and the gp41 ectodomain and to examine the conservation of predicted physicochemical and structural properties across a number of divergent HIV clades and groups. Results demonstrate that CTT sequences display intermediate levels of sequence evolution and diversity in comparison to the more diverse gp120 and the more conserved gp41 ectodomain. Despite the relatively high level of CTT sequence variation, the physicochemical properties of the lentivirus lytic peptide domains (LLPs) within the CTT are evidently highly conserved across clades/groups. Additionally, predictions using PEP-FOLD indicate a high level of structural similarity in the LLP regions that was confirmed by circular dichroism measurements of secondary structure of LLP peptides from clades B, C, and group O. Results demonstrate that LLP peptides adopt helical structure in the presence of SDS or trifluoroethanol but are predominantly unstructured in aqueous buffer. Thus, these data for the first time demonstrate strong conservations of characteristic CTT physicochemical and structural properties despite substantial sequence diversity, apparently indicating a delicate balance between evolutionary pressures and the conservation of CTT structure and associated functional roles in virus replication.

Mechanism of membrane curvature sensing by amphipathic helix containing proteins

There are several examples of membrane-associated protein domains that target curved membranes. This behavior is believed to have functional significance in a number of essential pathways, such as clathrin-mediated endocytosis, which involve dramatic membrane remodeling and require the recruitment of various cofactors at different stages of the process. This work is motivated in part by recent experiments that demonstrated that the amphipathic N-terminal helix of endophilin (H0) targets curved membranes by binding to hydrophobic lipid bilayer packing defects which increase in number with increasing membrane curvature. Here we use state-of-the-art atomistic simulation to explore the packing defect structure of curved membranes, and the effect of this structure on the folding of H0. We find that not only are packing defects increased in number with increasing membrane curvature, but also that their size distribution depends nontrivially on the curvature, falling off exponentially with a decay constant that depends on the curvature, and crucially that even on highly curved membranes defects large enough to accommodate the hydrophobic face of H0 are never observed. We furthermore find that a percolation model for the defects explains the defect size distribution, which implies that larger defects are formed by coalescence of noninteracting smaller defects. We also use the recently developed metadynamics algorithm to study in detail the effect of such defects on H0 folding. It is found that the comparatively larger defects found on a convex membrane promote H0 folding by several kcal/mol, while the smaller defects found on flat and concave membrane surfaces inhibit folding by kinetically trapping the peptide. Together, these observations suggest H0 folding is a cooperative process in which the folding peptide changes the defect structure relative to an unperturbed membrane.

Selected amino acid mutations in HIV-1 B subtype gp41 are associated with specific gp120v₃ signatures in the regulation of co-receptor usage

Background

The third variable loop (V3) of the HIV-1 gp120 surface protein is a major determinant of cellular co-receptor binding. However, HIV-1 can also modulate its tropism through other regions in gp120, such as V1, V2 and C4 regions, as well as in the gp41 protein. Moreover, specific changes in gp41 are likely to be responsible for of damage in gp120-CCR5 interactions, resulting in potential resistance to CCR5 inhibitors.

In order to genetically characterize the two envelope viral proteins in terms of co-receptor usage, we have analyzed 526 full-length env sequences derived from HIV-1 subtype-B infected individuals, from our and public (Los Alamos) databases. The co-receptor usage was predicted by the analysis of V3 sequences using Geno2Pheno (G2P) algorithm. The binomial correlation phi coefficient was used to assess covariation among gp120_V3 and gp41 mutations; subsequently the average linkage hierarchical agglomerative clustering was performed.

Results

According to G2P false positive rate (FPR) values, among 526 env-sequences analyzed, we further characterized 196 sequences: 105 with FPR <5% and 91 with FPR >70%, for X4-using and R5-using viruses, respectively.

Beyond the classical signatures at 11/25 V3 positions (S11S and E25D, R5-tropic viruses; S11KR and E25KRQ, X4-tropic viruses), other specific V3 and gp41 mutations were found statistically associated with the co-receptor usage. Almost all of these specific gp41 positions are exposed on the surface of the glycoprotein. By the covariation analysis, we found several statistically significant associations between V3 and gp41 mutations, especially in the context of CXCR4 viruses. The topology of the dendrogram showed the existence of a cluster associated with R5-usage involving E25D_V3, S11S_V3, T22A_V3, S129DQ_gp41 and A96N_gp41 signatures (bootstrap = 0.88). Conversely, a large cluster was found associated with X4-usage involving T8I_V3, S11KR_V3, F20IVY_V3, G24EKR_V3, E25KR_V3, Q32KR_V3, A30T_gp41, A189S_gp41, N195K_gp41 and L210P_gp41 mutations (bootstrap = 0.84).

Conclusions

Our results show that gp120_V3 and several specific amino acid changes in gp41 are associated together with CXCR4 and/or CCR5 usage. These findings implement previous observations that determinants of tropism may reside outside the V3-loop, even in the gp41. Further studies will be needed to confirm the degree to which these gp41 mutations contribute directly to co-receptor use.