Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Next Article in Journal
Letea Virus: Comparative Genomics and Phylogenetic Analysis of a Novel Reassortant Orbivirus Discovered in Grass Snakes (Natrix natrix)
Previous Article in Journal
Measles Vaccines Designed for Enhanced CD8+ T Cell Activation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Identification of HIV-1 Envelope Mutations that Enhance Entry Using Macaque CD4 and CCR5

1
Divisions of Human Biology, Fred Hutchinson Cancer Research Center, Seattle, WA 98109, USA
2
Basic Sciences, Fred Hutchinson Cancer Research Center, Seattle, WA 98109, USA
*
Author to whom correspondence should be addressed.
Present Address: Berkeley Brewing Science, Inc., Berkeley, CA 94704, USA.
Viruses 2020, 12(2), 241; https://doi.org/10.3390/v12020241
Submission received: 27 December 2019 / Revised: 13 February 2020 / Accepted: 19 February 2020 / Published: 21 February 2020
(This article belongs to the Section Animal Viruses)

Abstract

:
Although Rhesus macaques are an important animal model for HIV-1 vaccine development research, most transmitted HIV-1 strains replicate poorly in macaque cells. A major genetic determinant of this species-specific restriction is a non-synonymous mutation in macaque CD4 that results in reduced HIV-1 Envelope (Env)-mediated viral entry compared to human CD4. Recent research efforts employing either laboratory evolution or structure-guided design strategies have uncovered several mutations in Env’s gp120 subunit that enhance binding of macaque CD4 by transmitted/founder HIV-1 viruses. In order to identify additional Env mutations that promote infection of macaque cells, we utilized deep mutational scanning to screen thousands of Env point mutants for those that enhance HIV-1 entry via macaque receptors. We identified many uncharacterized amino acid mutations in the N-terminal heptad repeat (NHR) and C-terminal heptad repeat (CHR) regions of gp41 that increased entry into cells bearing macaque receptors up to 9-fold. Many of these mutations also modestly increased infection of cells bearing human CD4 and CCR5 (up to 1.5-fold). NHR/CHR mutations identified by deep mutational scanning that enhanced entry also increased sensitivity to neutralizing antibodies targeting the MPER epitope, and to inactivation by cold-incubation, suggesting that they promote sampling of an intermediate trimer conformation between closed and receptor bound states. Identification of this set of mutations can inform future macaque model studies, and also further our understanding of the relationship between Env structure and function.

1. Introduction

Entry receptors are a major determinant of the host cell specificity of viruses [1]. Within an infected individual, the distribution of receptors on distinct cell types defines the population of cells susceptible to viral infection. In the case of HIV-1, human T lymphocytes, monocytes and macrophages are the principal host cells due to their expression of the HIV-1 receptor, CD4, and preferred co-receptor, CCR5. Similarly, variation in receptor orthologs among species determines the range of organisms susceptible to viral infection [2]. For HIV-1, although the human CD4 receptor permits efficient infection, CD4 receptors from other nonhuman primates (NHPs) are suboptimal, and only weakly support viral entry [3,4,5,6].
The inability of HIV-1 to efficiently utilize CD4 receptors from many NHPs has limited the use of NHPs as models for studying HIV-1 pathogenesis. The Rhesus macaque has been the favored NHP used to date, with SIV/HIV-1 hybrid viruses (SHIVs), which encode HIV-1 Env in place of SIV Env, used to model HIV-1 infection. As SHIVs encoding Envs from circulating HIV-1 strains replicate poorly in macaque cells [7], most studies employing the SHIV/macaque model have used HIV-1 Envs that were lab adapted to enhance macaque infection, and/or that were derived from HIV-1 variants isolated from the chronic stages of infection [7]. These evolved and chronic-stage Envs are highly pathogenic and transmissible in macaques [8]. While these characteristics enabled robust infection of macaques, evolved and chronic-stage Envs display altered structural conformations and increased susceptibility to neutralizing antibodies compared to Envs from circulating HIV-1 strains [8,9]. Use of these Envs in SHIVs therefore compromises the relevance of the SHIV/macaque model for informing HIV-1 infection in humans [7].
In 2012, it was discovered that a single amino acid mutation in macaque CD4 was responsible for the poor ability of HIV-1 to utilize macaque receptors for entry [3]. In vitro evolution experiments also revealed that single amino acid substitutions at sites 204 or 312 in Env’s gp120 subunit allow HIV-1 to more effectively use macaque CD4 as an entry receptor [4]. These mutations disrupted quaternary contacts between Env trimer subunits, and similar to Envs encoded by highly pathogenic SHIVs, increased sensitivity to neutralizing antibodies [10].
Identification of these mutants prompted further efforts to identify additional Env mutations that would enhance usage of macaque CD4. In 2016, Li et al. used a structural modeling approach to discover that multiple substitutions at Env site 375 permitted robust infection of macaque cells [11]. Later, in 2018, Del Prete et al. utilized an in vivo evolution strategy to identify a mutation at site 281, located within the CD4 binding site, that enhanced the ability of HIV-1 to use macaque CD4 [12]. Mutations at both 375 and 281 had minimal impact on the trimer structure and antibody neutralization profile.
The discovery of mutations that allow HIV-1 to use the macaque CD4 receptor without altering neutralization susceptibility has opened the possibility of improving the relevance of the SHIV/macaque model. As these previous studies each utilized a different experimental approach and discovered distinct adaptive mutations, it is likely that additional mutations also exist that would enhance Env usage of macaque CD4. Identification of these mutations could serve to further improve the relevance of the SHIV/macaque model, while also advancing our understanding of the relationship between Env sequence and function.
Deep mutational scanning (DMS) is a powerful experimental tool that allows researchers to simultaneously assay the effects that thousands of individual mutations have on a protein during selection for a desired function [13]. DMS has previously been used to determine the effects of all single amino acid substitutions in HIV-1 Env in the context of HIV-1 replication in cell culture [14,15], neutralization by antibodies [16,17], and inhibition by fusion inhibitor drugs [18]. Here, we utilized DMS to identify Env mutations that enhance infection of cells expressing macaque CD4 and CCR5 receptors. Importantly, the Env strain background used for these studies was cloned directly from a recently infected individual and is thus representative of transmitted/founder (T/F), circulating HIV-1 strains. We identify a suite of novel mutations in Env’s gp41 subunit that increase HIV-1 infection of cells expressing macaque CD4 and CCR5 receptors.

2. Materials and Methods

2.1. Cells

HEK293T cells (referred to as 293Ts) were used for production of all replication-competent viruses and pseudoviruses. The 293Ts were maintained in D10 media (Dulbecco’s modified Eagle medium supplemented with 10% heat inactivated fetal calf serum and 2 mM L-glutamine).

2.2. Generation of Stable 293Trhm and 293Thu Cells Lines

Plasmids encoding Rhesus macaque CD4 (rhmCD4) and human CD4 (huCD4) were obtained from [4]. Plasmids encoding Rhesus macaque CCR5 (rhmCCR5) and human CCR5 (huCCR5) were obtained from the NIH-ARP as pBABE.Rh-CCR5 (catalog #3599, GenBank Accession U73739), and pBABE.CCR5 (catalog #3331), respectively. Receptor genes encoded by these four plasmids were each cloned into the pHAGE2 lentiviral vector under a CMV promoter. Transfection of 293Ts with these lentiviral vectors and Fugene-6 (Promega, Madison, WI; E2691) was used to generate virus-like particles (VLPs) carrying receptor gene coding sequences. VLPs carrying rhmCCR5 or huCCR5 genes were used to transduce 293T cells, followed by sorting of cells to isolate populations expressing intermediate levels of either rhmCCR5 (293TrhmCCR5) or huCCR5 (293ThuCCR5). The 293TrhmCCR5 cells were then transduced with VLPs carrying the rhmCD4 gene and sorted to isolate 293Ts expressing both rhmCD4 and rhmCCR5 (293Trhm). The 293Thu cells were analogously generated by transduction of 293ThuCCR5 with VLPs carrying the huCD4 gene, followed by sorting. Antibodies used for sorting and quantifying receptor expression were CD4: BD#551980, and CCR5: BD#556042. PCR primers used to amplify and clone receptor sequences are available upon request.

2.3. TZM-bl Titering Assays

Serial dilutions of viral supernatants were added to 2 × 104 TZM-bl cells in the presence of 10 μg/mL DEAE-dextran. After 48 h, cells were fixed and stained for beta-galactosidase, and blue foci were counted to quantify titer.

2.4. Preparation of Green Fluorescent Protein (GFP) Reporter Pseudovirus

The 293T cells were plated at 5 × 105 cells/well in a 6-well dish approximately 24 h prior to transfection. Cells were co-transfected using 6 ul of Fugene-6, 1.2 ug of plasmid Q23ΔEnv-GFP [4], and 8 ug of the Env expression plasmid of interest. Env clones BF520.W14.C2 [19] and BG505.W6.C2.T332N [20] were used for generation of BF520 and BG505 Env pseudoviruses, respectively. Two days after transfection, supernatants were collected and filtered through a 0.2 μM filter, then aliquoted and stored at −80 °C. Supernatants were titered by measuring p24gag levels with an HIV-1 p24 antigen capture assay kit (ABI, Foster City, CA, USA; 5447), and also by TZM-bl assay.
For generation of Env point mutant pseudoviruses, Env expression plasmids were mutagenized by PCR-mediated site-directed mutagenesis, and Sanger sequencing used to confirm the desired mutation. Primers used for site-directed mutagenesis are available upon request.

2.5. Preparation of Replication-Competent BF20.W14M.C2 Viruses

The 293T cells were plated at 5 × 105 cells/well in a 6 well dish approximately 24 h prior to transfection. Cells were co-transfected using 6 μL of Fugene-6, and 2 μg of plasmid Q23.BsmBI.BF520.C2, described in [17]. Two days after transfection, supernatants were collected and filtered through a 0.2 μM filter, then aliquoted and stored at −80 °C. Supernatants were titered by TZM-bl assay.

2.6. Preparation of Mutant Env Libraries

The 293T cells were plated at 5 × 105 cells/well in a 6-well dish approximately 24 h prior to transfection. Cells were transfected using 6 μL of Fugene-6 and 2 μg of BF520.C2 mutant plasmid library, described in [17]. Two days after transfection, supernatants from 16 replicate wells were pooled, filtered through a 2 μM filter, and DNase treated (Roche, Basel, Switzerland; 4716728001) as previously described [14]. Supernatants were titered by TZM-bl assay.
To establish a genotype–phenotype link for viruses in each replicate mutant library, 1.16 × 106 infectious units of transfection supernatant were used to infect a 3:1 mixture of 293Thu and 293Trhm cells at an MOI of 0.01 in 5 layer flasks (Falcon, Franklin Lakes, NJ, USA; 353144) supplemented with 10 μg/mL DEAE-dextran. Twenty-four hours after infection, the media was replaced with fresh D10 media supplemented with 10 μg/mL DEAE-dextran. Four days after infection, media from all flasks was pooled and filtered through a 0.2 μM filter. Viruses were then concentrated 50-fold by ultracentrifugation over a sucrose cushion at 23,000 RPM for 1 hour at 4°C using a SW28 rotor (Beckman Coulter, Brea California; 342207). Concentrated virus was re-suspended in D10 media, and titered by TZM-bl assay.

2.7. Infection of 293Trhm and 293Thu Cells with Wild-Type BF520 Env

T25 flasks were seeded with 3.3 × 106 of either 293T, 293Trhm, 293Thu, or a 3:1 mixture of 293Thu:293Trhm cells. All flasks were infected with Q23.BsmBI.BF520.C2 virus at an MOI of 0.01. Twenty-four hours after infection, cells were split 1:4. At time points two, three and four days after infection, 60 μL samples were taken for titering on TZM-bl cells.

2.8. 293Trhm and 293Thu Selection of Mutant Env Libraries and Deep Sequencing

Wells of a 12-well dish were seeded with 3 × 105 of either 293Trhm or 293Thu cells. The following day, both 293Trhm and 293Thu cells were infected with 1.44 × 106 (replicate 1) or 2.06 × 106 (replicate 2) infectious units of genotype–phenotype linked mutant virus library in the presence of 100 μg/mL DEAE-dextran. We additionally infected one well of 293Thu cells with 1.72 × 106 infectious units of wild-type BF20.W14M.C2 to use as a control during DMS analysis. Three hours after infection, the media was removed and replaced with fresh D10 media. Twelve hours after infection, cells were harvested, pelleted, washed once with PBS, and non-integrated viral DNA was isolated by miniprep. Deep sequencing of viral DNA to quantify the functional effects of Env mutations was performed as described in [17].

2.9. Env Sequence Numbering

BF520 and BG505 Env residues are numbered according to the reference strain numbering system [21].

2.10. Infection of 293Trhm and 293Thu Cells by GFP Reporter Pseudovirus

Wells of a 96-well dish were seeded with 1.2 × 104 of either 293Trhm, 293Thu, or 293TrhmCCR5 cells. The following day, 8 × 103 pg of GFP reporter pseudovirus was added to wells of the plate (corresponding to an MOI of 5) in the presence of 10 μg/mL DEAE-dextran. Two days after infection, cells were harvested, washed once in PBS and re-suspended in 1% paraformaldehyde. Cells were analyzed for GFP expression by flow cytometry.

2.11. Neutralization Assays

Neutralization assays using TZM-bl cells were performed as described in [22]. Reported values reflect the mean of two independent replicates and are representative of at least two neutralization assays.

2.12. Cold Inactivation Assays

Recombinant pseudoviruses were diluted to 1.2 × 105 infectious units/mL, split into 50 μL aliquots, and frozen at −80 °C overnight. At different time points, aliquots were thawed in a 37 °C water bath for two minutes, then placed on ice. Following incubation on ice for different periods of time, all pseudovirus aliquots were incubated in a 37 °C water bath for two minutes, and 6 × 102 IPs were added to wells of a 48-well plate containing 2 × 104 TZM-bl cells in 10 μg/mL DEAE-dextran. Two days later, infection of TZM-bl cells was quantified by a firefly luciferase expression assay [23].

2.13. Structural Analysis

All structural analyses used the X-ray diffraction model of the BG505 SOSIP trimer in complex with a PGT121 precursor Fab [24]. Molecular graphics and analyses were performed using UCSF Chimera, developed by the Resource for Biocomputing, Visualization, and Informatics at the University of California, San Francisco, with support from NIH P41-GM103311 [25].

2.14. DMS Data Analysis and Source Code

DMS sequencing data analysis was performed using the dms_tools2 software package version 2.4.1 [26] (https://jbloomlab.github.io/dms_tools2/). A jupyter notebook that describes this analysis and generates Figures 2,3 and 5 is available at: https://github.com/jbloomlab/HIV_Env_macaque_receptor_DMS. The deep sequencing data are available on the Sequence Read Archive under accession numbers SRA: SRX5199658–SRX5199665.

3. Results

3.1. Deep Mutational Scanning Identifies Env Mutations That Improve Usage of Macaque Receptors

The Env used in our DMS, subtype A BF520.W14M.C2 env (hereafter referred to as BF520), is a (T/F) strain cloned directly from a Kenyan infant at the time point when HIV was first detected following mother-to-child transmission [19]. Like Envs from other circulating HIV-1 strains, the wild type BF520 Env supported infection of cells engineered to express human CD4 and CCR5 (293Thu) but was severely restricted in its ability to infect cells expressing macaque receptors (293Trhm) (Figure S1) [3].
In our DMS, we screened a library of BF520 Env variants to identify mutations that increased entry into cells expressing macaque receptors. The DMS strategy we employed is outlined in Figure 1a, with all steps performed in biological duplicate. Viruses bearing mutant env genes were generated by transfection of two independently constructed mutant plasmid libraries into 293T cells. Each plasmid in these replicate plasmid libraries encodes a codon level mutational variant of BF520 env, in the context of the proviral genome from the subtype A strain, Q23, which was isolated during the first year of infection [17]. Prior deep sequencing demonstrated that these libraries collectively encode over 95% of the 12,559 possible single amino acid mutations to Env’s ectodomain (residues 32 to 702), with an average of 1.1 codon mutations per each env variant [17]. Mutant viruses were then passaged for four days at a low MOI (0.01) in a 3:1 mixture of 293Thu and 293Trhm cells in order to establish a genotype–phenotype link between mutant env genes and Env proteins on each virus. A mixture of 293Thu and 293Trhm cells was used for this passage so as to avoid the possibility that env mutants that enhanced usage of macaque receptors would drop out of the pool. Importantly, replication kinetics of wild-type BF520 viruses in this mixed cell context were indistinguishable from replication in 293Thu cells alone (Figure S1), indicating that a reduction in susceptible target cells during this step was unlikely to impose a bottleneck on the diversity of the mutant virus pool. Following the low MOI passage, genotype–phenotype linked viruses were used to infect 293Trhm cells in order to select for mutants that enhance entry via macaque receptors. In order to investigate whether mutations enriched by this infection represented macaque-specific adaptations or instead improved infection by more general mechanisms, we simultaneously performed parallel infections of 293Thu cells with the same mutant virus pool. Importantly, expression of CCR5 and CD4 receptors was similar between 293Trhm and 293Thu cell types (Figure S2).
Following infection of 293Trhm and 293Thu cells, we deep sequenced env DNA isolated from infected cells, as well as from the initial mutant plasmid library. Similar to previous DMS experiments, the frequency of stop codons in post-selected samples was reduced to 19–24% of their frequency in the mutant plasmid library (Figure 1b), evidence that our DMS strategy effectively selected for functional Env variants [17,27]. Frequencies of nonsynonymous mutations were also depleted in post-selected samples compared to pre-selected plasmid libraries, consistent with previous reports that amino acid substitutions in Env are generally deleterious (Figure 1b) [14,15].
To quantify the effect of each Env mutation on infection of either 293Trhm or 293Thu cells, we first inferred the relative preference for all 20 amino acids at each Env site, for each infection, as previously described [26]. We next calculated a statistic, Smut/wt, which describes selection on each Env mutant as the log2 of the ratio of the mutant residue preference to the wild-type residue preference at a given site (Figure 1a). Smut/wt values from replicate 293Thu infections were well correlated (Figure 1c) (Pearson R = 74), whereas those from replicate 293Trhm infections were more moderately correlated (Figure 1d) (Pearson R = 53), although still comparable to previous Env DMS experiments [14]. The weaker correlation between 293Trhm replicates is suggestive of a highly selective bottleneck in which strongly beneficial mutations were enriched against a backdrop of stochastic infection by neutral and slightly deleterious variants (Figure 1d, compare upper right quadrant to other quadrants). As expected, the majority of Env mutations were selected against during both 293Trhm and 293Thu infections, evidence that these were generally detrimental to Env function (Figure 1c–e).
Infection of 293Trhm and 293Thu cells exerted similar overall selective pressures on Env mutants, as evidenced by the strong correlation between the average Smut/wt values obtained from these infections (Figure 1e). However, the magnitude of enrichment for strongly beneficial mutations was far greater following 293Trhm infection (Figure 1e, points above grey line in upper right quadrant). For example, the most highly enriched mutation following 293Trhm infection (N656Y) was enriched 34-fold relative to the wild-type residue, compared to a 12-fold enrichment for the most strongly selected mutation during 293Thu infection (E482K). This finding likely reflects the fact that wild-type BF520 Env is poorly adapted to infect macaque cells, and therefore mutants in this background have the potential to greatly enhance usage of macaque receptors.
Env mutations that had positive Smut/wt values following infection of 293Trhm cells are shown in Figure 2. Surprisingly, while previous studies have predominately identified mutations in Env’s gp120 subunit that improve infection of macaque cells [4,11,12,28], the majority of mutations highly enriched in our experiments were located in the gp41 subunit. In fact, all amino acid mutations enriched greater than 10-fold over wild-type residues in our DMS were located at just seven sites, all within gp41. Six of these sites were located within either the N terminal heptad repeat (NHR) or C-terminal heptad repeat (CHR), while the seventh was in the intervening loop region.
While the most strongly enriched mutations were in N/CHR regions, we observed modest enrichment of mutations in several other Env regions. At Env sites 58–63, 20 mutations were selected that improved 293Trhm infection 3- to 8-fold compared to wild-type residues (Figure 2). These sites lie within a short helix that was recently found to make contact with a CD4 molecule bound to the inner domain of a neighboring gp120 protomer [29]. Our DMS also revealed that 293Trhm infection is modestly improved by Env mutations that abolish the putative glycosylation motif at sites 611–613 (Figure 2).
Given the surprising result that previously uncharacterized mutations in the N/CHR showed the greatest enrichment in the DMS, we sought to validate the effect of these mutations on macaque receptor usage through an independent assay. We generated pseudoviruses bearing a subset of these gp41 mutations in the BF520 Env background, then tested the ability of the pseudovirus mutants to infect 293Trhm cells. Compared to wild-type BF520, all seven gp41 mutants assayed in these experiments improved entry into 293Trhm cells (Figure 3a). Similar to their relative enrichment in the DMS, mutations in the CHR were most beneficial to 293Trhm infection. For example, CHR mutations Q653L and Q652F improved infection of 293Trhm cells 9- and 7-fold, respectively, whereas the most advantageous NHR mutation, I573L, improved infection 4-fold (Figure 3a). These data thus validate that N/CHR mutations enriched in our DMS following infection of 293Trhm cells do indeed enhance entry of BF520 Env into cells bearing macaque receptors.

3.2. Similar Mutations Improve Usage of Both Macaque and Human Receptors

Infection of 293Thu cells with the mutant virus pools selected for beneficial mutations at many of the same sites where mutations were enriched following 293Trhm infection (Figure 4). Mutations at sites 58–61, as well as those in the N611 glycosylation motif, were enriched up to 8-fold over wild-type residues following 293Thu infection, similar to their enrichment following 293Trhm infection.
Within gp41, 293Thu infection selected for mutations at many of the same sites as 293Trhm infection, but with a lower magnitude of enrichment. Within the NHR and CHR, no mutations were enriched greater than 10-fold over the wild-type-residue following 293Thu infection, compared to 24 such mutations following 293Trhm infection (Table S1). This difference in enrichment was especially pronounced at CHR sites 653 and 656, where the most beneficial mutations were enriched 27-and 34-fold, respectively, following 293Trhm infection, but only 5- and 9-fold following 293Thu infection.
Few mutations outside of the CD4-contacting helix (sites 54–74), the N611 glycan motif, and N/CHR regions were selected for during either 293Trhm or 293Thu infection. However, two notable exceptions were sites 588 and 589, located in the loop region between the NHR and CHR. Following 293Trhm infection, many mutants were enriched up to 21-fold over WT residues at site 588, with a few mutants weakly enriched at site 589. The inverse was true for 293Thu selection, which enriched for mutants up to 6-fold at site 589, but weakly at site 588.
Pseudovirus infection assays qualitatively validated the effects of mutations enriched in our DMS following 293Thu selection (Figure 3b). Consistent with the DMS data, mutations R557L and K588L, which improved infection of 293Trhm cells, were mildly detrimental in the context of 293Thu infection. Five other N/CHR mutations that were enriched following both 293Thu and 293Trhm DMS selections improved entry of pseudovirus into 293Thu cells in this assay. The magnitude of the benefit to 293Thu infection in this assay was less than in the context of 293Trhm infection (Figure 3a,b), again mirroring the increased effect size of these mutations observed in the DMS. Taken together, these data validate that 293Trhm and 293Thu infection selected for beneficial mutations within similar regions of Env, often at identical sites.

3.3. Mutations Identified by DMS have Smaller Effects than Previously Identified Macaque Adaptive Mutations

Surprisingly, our DMS did not select for Env mutations at sites 204, 281, 312 or 375 that have previously been identified to enhance infection of macaque cells [4,11,12]. Mutations at these sites were adequately represented in the initial mutant plasmid library thus ruling out the possibility that they were not present at the start of our DMS experiment. Several possible hypotheses could explain the apparent lack of selection for these mutants in the DMS. First, these mutations may not have been tolerated in the BF520 Env strain background. If this had been the case, they would have dropped out of the DMS following passage of the mutant plasmid library. A second possibility is that these mutations may not have improved usage of macaque receptors in the BF520 strain background.
To determine whether either of these possibilities explained the lack of enrichment of these previously identified mutants, we first tested whether we could generate infectious pseudovirus bearing two of these mutations, A204E and S375Y, in the BF520 background. The titers of pseudoviruses bearing the S375Y Env mutation, as determined by infecting TZM-bl cells, were similar to those of pseudoviruses bearing wild-type BF520 (2.6 × 106 and 2.8 × 106 infectious units/mL, respectively). In contrast, replicate transfection supernatants containing A204E mutations were not infectious. We thus conclude that the A204E mutation is incompatible with the BF520 Env, thereby explaining its lack of enrichment in our DMS.
We next assessed whether the S375Y mutation in the BF520 background would enhance infection of 293Trhm and 293Thu cells in pseudovirus reporter assays. As shown in Figure 3a, the S375Y mutations increased infection of 293Trhm cells 34-fold compared to wild-type BF520. This was a higher fold increase in infection than any of the mutations that were most enriched in the DMS. In the context of 293Thu infection, S375Y enhanced wild-type infection by 1.7-fold (Figure 3b). These data indicate that the lack of enrichment of S375Y in our DMS cannot be attributed to this mutation being incompatible with the BF520 Env, nor to it impeding 293Trhm infection in this strain background. Additionally, the fact that S375Y was not detrimental to 293Thu infection indicates that the mixture of 293Thu and 293Trhm cells used to establish a mutant genotype–phenotype link in our DMS was likely not responsible for selection against this mutant.
Why then, was 375Y not enriched in our DMS? An examination of our data revealed that the frequency of S375Y, as well as previously identified macaque adaptive mutations G312V and A281T, was lower following 293Trhm and 293Thu than in the mutant plasmid libraries, evidence of selection against these mutations during the experiment. To further explore this finding, we revisited data from a previously conducted DMS experiment in which the same BF520 mutant Env library was passaged and selected in a human cell line expressing human CD4 and CCR5 (SupT1 T cell line; [17]). In these data, we observed a similar drop in frequencies of S375Y, G312V, and A281T that mirrored the depletion of these mutations in the current DMS experiment. Selection against these mutations in both DMS experiments suggests that properties inherent to the mutant plasmid library or the DMS protocol are likely responsible for their lack of enrichment following 293Trhm infection. One potential explanation is that while these mutations enhance entry of pseudovirus into 293Trhm and 293Thu cells, they are detrimental in the context of multi-cycle replication and purged from the mutant pool during the low MOI infection necessary to generate genotype–phenotype linked virions.

3.4. NHR/CHR Mutants Enhance Usage of Macaque Receptors by Env Strain BG505

We next sought to determine whether the mutations enriched in the DMS were advantageous solely in the context BF520 Env, or instead, reflected a more widely accessible mutational strategy through which Env could evolve to utilize macaque receptors. To investigate this, we introduced seven N/CHR mutants enriched in our DMS into the BG505.W6M.C2.T332N env (hereafter referred to as BG505) [30], then tested the ability of these mutants to enhance pseudovirus infection of 293Trhm cells. BG505 is a subtype A, T/F Env [30], 15.8% diverged from BF520 at the amino acid level, and used extensively as a SOSIP trimer for Env structural studies [20,31,32,33,34], and as a vaccine immunogen [20,35]. Recent DMS efforts have also revealed that the site-specific amino acid preferences and mutational tolerance of BG505 Env have diverged from BF520 [15].
Infection of 293Trhm cells with wild-type BG505 also revealed this strain to be ~7-fold less effective at entering cells bearing macaque receptors compared to wild-type BF520. At an MOI of 5, only 0.12% of 293Trhm cells were infected following wild-type BG505 infection (Figure 3c) compared to 0.87% following wild-type BF520 infection (Figure 3a).
Four of the seven N/CHR mutants introduced into BG505 improved the ability of this Env to infect 293Trhm cells (Figure 3c). The magnitude of this effect was smaller than in the BF520 strain, however, with the most advantageous mutant (Q653L) conferring a 2.5-fold improvement over wild-type BG505. Generally, mutations that conferred the greatest improvement to macaque receptor usage in the BF520 strain background also had the largest effect in BG505. For example, mutations Q652F and Q653L, which conferred the greatest benefit to infection in the BF520 background were also the most advantageous in the BG505 strain. Only one mutation, R557L, benefited BF520 infection of 293Trhm cells, but was detrimental to infection of these cells in BG505. Taken together, these data indicate that while some N/CHR mutations have strain-specific effects, the benefit conferred by larger-effect mutations is not restricted to the BF520 background.

3.5. N/CHR Mutations Increase Susceptibility to MPER Antibodies but Minimally Affect Other Epitopes

We next sought to investigate the mechanism through which the N/CHR mutations enhanced usage of macaque receptors. First, we tested whether these mutants conferred a CD4-independent phenotype by assaying whether they allowed entry of BF520 pseudovirus into 293T cells expressing only the macaque CCR5 receptor (293TrhmCCR5). No infection was observed by either wild-type or mutants Envs in this assay, thereby ruling out CD4-independence as the mechanism through which N/CHR mutants enhance 293Trhm infection (Figure 3d).
We next used a panel of monoclonal antibodies (mAbs) and inhibitory mimetic drugs to determine whether N/CHR mutations caused alterations in Env trimer structure that would inform on a mechanism through which they improved entry into 293Trhm cells. Drugs and mAbs in the panel were chosen so as to detect alterations to trimer structure in regions both proximal and distal to the mutant N/CHR residues, with particular interest in the CD4 binding site. In these experiments, we assessed neutralization susceptibility of three mutants that were highly enriched in our DMS: I573L (NHR), Q653L (CHR) and N656Y (CHR). Neutralization of each mutant was assessed in both the BF520 and BG505 Env strains so as to inform on the background dependence of any structural effects. Compared to wild-type Envs, introduction of N/CHR mutants had little effect on neutralization sensitivity to antibodies targeting the CD4bs, trimer apex, V3 loop, or six-helical-bundle epitopes (Figure 5). All tested mutants did, however, increase sensitivity to the MPER antibodies 4E10 and 10E8 in both BG505 and BF520 backgrounds (Figure 5). This effect was strongest for CHR mutants Q653L and N656Y, which increased sensitivity to MPER antibodies by 15- and 8-fold, in BF520 and BG505 backgrounds, respectively (Figure 5a,b). In the BF520 strain, NHR/CHR mutants also modestly increased sensitivity to T-20, a fusion-inhibitor drug that binds a hydrophobic pocket at the N-terminus of the NHR helix (Figure 5a,c). Again, the effect was largest for CHR mutants, which increased T-20 sensitivity 2- to 4-fold.

3.6. N/CHR Mutations Increase Sensitivity to Cold Inactivation

Prior to engagement of host receptors, it is thought that Env predominately maintains a closed conformation in which MPER and NHR residues are poorly exposed and not accessible to interactions with mAbs or small molecule inhibitors [36]. Receptor binding initiates formation of the pre-hairpin intermediate conformation, in which both the MPER epitopes, as well as the T-20 binding site in the NHR are exposed [37,38,39]. Several recent publications have described Env mutations that promote unliganded trimer sampling of an intermediate conformation, termed “state 2”, between closed (state 1) and receptor bound (state 3) conformations [40,41,42,43]. State 2 is characterized by a heightened intrinsic reactivity to ligand binding, increased capacity to infect cells expressing low levels of CD4 or CCR5 receptors, and increased exposure of MPER and NHR residues. Given that the N/CHR mutants we identified increase exposure of MPER epitopes, and in strain BF520, the NHR helix, we hypothesized that they may increase reactivity to macaque receptors by promoting sampling of a state 2 conformation.
Envs frequently sampling state 2 conformations display increased sensitivity to inactivation during cold incubation [40,44]. To investigate whether N/CHR mutants identified in our DMS also increased sensitivity to cold-inactivation, we assessed infectivity of wild-type and mutant BF520 pseudoviruses following incubations on ice. Introduction of Q653L, N656Y and I573L into BF520 decreased Env infectivity following ice incubation (Figure 6). This effect was more pronounced for Q653L and N656Y than for I573L, paralleling the relative effect sizes of these mutants observed during 293Trhm infection and in neutralization sensitivity assays. It thus appears that the degree to which N/CHR mutants enhance usage of macaque receptors may be correlated to increased cold-sensitivity and exposure of MPER epitopes on the unliganded Env trimer. We note that this is not uniformly the case, however, as the N656Y mutation in BG505 increased sensitivity to MPER antibodies but did not enhance infection of 293Trhm cells (Figure 4), evidencing strain-specific effects.

4. Discussion

We identified mutations in several regions of an Env from a transmitted HIV-1 virus that led to increased entry using macaque CD4 and CCR5 receptors. Mutations within gp120 modestly improved usage of macaque receptors, but mutations with the largest beneficial effect were located in gp41, most notably in the N/CHR regions. None of the tested mutations conferred a CD4-independent phenotype. Several N/CHR mutations also enhanced entry into cells bearing macaque receptors by a second transmitted Env strain, indicating that their effects were not strain specific. While N/CHR mutations have previously been identified that heighten fusogenic activity [45,46] and increase infectivity of cells expressing low levels of CD4 and CCR5 receptors [40,47], none have been shown to enhance entry of cells bearing macaque receptors. Many of the N/CHR mutations we identified also enhanced entry into cells bearing human CD4 and CCR5 receptors, but the magnitude of this effect was smaller than during infection of cells bearing macaque receptors.
Our observation that N/CHR mutations increase susceptibility to neutralization by MPER antibodies, and, in the BF520 background, to inhibition by the fusion inhibitor T-20, provide some insight into their effects on Env trimer conformation. The MPER epitope (sites 671–683) and T-20 binding pocket (sites 547–556) are located immediately N, and C terminal, respectively, to the NHR and CHR mutants enriched in our DMS [18,48]. Our neutralization data thus indicate that these mutations cause Env conformation changes that are confined to regions of the trimer directly adjacent to the mutated residues. The mechanism by which these mutants enhance usage of macaque receptors thus appears to be distinct from that of mutations A204E and G312V, which neutralization data suggest disrupt quaternary Env contacts and promote a more open trimer conformation [10]. The mechanism also appears to be distinct from mutations at site 375 and 281, which are located within the CD4 binding site region, and minimally impact Env susceptibility to neutralizing antibodies [11,12].
Viruses with mutations in NHR and CHR also demonstrated increased sensitivity to inactivation by cold. This phenotype has been linked to Envs frequently sampling state 2 conformations, in which the thermodynamic barrier between unliganded and CD4-bound states is lowered [40,41]. Envs displaying a state 2 phenotype are also known to display heightened sensitivity to MPER antibodies [40,41,42,43], and in some cases, increased exposure of the NHR helix [40]. Our observations that N/CHR mutations increase sensitivity to MPER antibodies, and in the BF520 background modestly increase T-20 sensitivity, are consistent with the hypothesis that these mutations promote a state 2 conformation.
Mapping of the N/CHR mutations identified here onto the structure of the unliganded BG505 SOSIP trimer reveals that NHR mutations are located within the central core of the trimer, while CHR mutations are positioned at the trimer base. Mutations in both regions are located at sites that make inter- and intra-protomer contacts with adjacent gp120 and gp41 subunits (Figure 7). Given that previous reports have described N/CHR mutations that modulate trimer reactivity through alterations of inter-protomer contacts [43,49,50], it is possible that the mutations identified here promote a state 2 phenotype by disrupting the stabilizing interactions between protomers within the unliganded trimer. If N/CHR identified here do indeed promote a state 2-like conformation, this would lower the thermodynamic barrier between unliganded and CD4 bound Env states, which could allow relatively weak interactions between Env and macaque CD4 to initiate the conformational changes leading to exposure of the co-receptor binding site. Subsequent binding of the co-receptor would lead to the cascade of Env rearrangements that results in membrane fusion and cell entry. While speculative, this model of macaque receptor usage is in line with a recent report of retroviral host-range expansion associated with envelope mutations that promote a receptor-bound conformation [51]. Alternatively, N/CHR mutations could enhance entry into macaque cells by a distinct mechanism. For example, mutations to Env’s gp41 and heptad repeat regions have previously been found to alter virus infectivity by altering Env expression [52], incorporation into virions [53,54], and trimer stability [55,56]. The N/CHR mutations identified in our DMS may enhance entry into 293Trhm though these or other mechanisms, either independent of, or in conjunction with their potential effects on trimer reactivity.
We note that although the NHR and CHR mutations identified in this work generally enhance usage of both macaque and human receptors, they occur at sites that are highly conserved among circulating HIV-1 strains and lentiviruses [50,57]. Selective pressures exerted by the human immune system could explain the discrepancy between the enrichment of these mutations in our experiments and their lack of representation in natural isolates. Specifically, our data suggest that N/CHR mutations that enhance entry also make the virus more susceptible to neutralization by MPER antibodies. As the MPER epitope is highly conserved [48] and functionally constrained [15], mutations that exposed this region are predicted to be strongly selected against. Thus, in an established HIV-1 infection, the advantage that N/CHR mutations have during entry may be balanced by the necessity of evading the immune response.
While the most strongly enriched mutations were in N/CHR regions, we observed modest enrichment of mutations in several other Env regions, including sites 58–63 and 611–613. Enriched mutations at sites 58–63 are located within the recently discovered CD4-contacting helix, in which several residues make contact with CD4 bound to the inner domain of an adjacent gp120 protomer [29]. Mutations in this region, when introduced into SOSIP trimers, have also been shown to affect reactivity to CD4 engagement and conformational dynamics [35]. The mutations enriched in our DMS may therefore improve entry by either strengthening contacts with CD4 or increasing Env reactivity to CD4 engagement. The role of mutations at sites 611–613, which abolish the putative N611 glycosylation motif, is less clear. Previous studies have found that removal of certain N-linked glycosylation sites in SIV gp120 enhance replication in macaque macrophages [58] and increase exposure of CD4-contact residues [59], but the relevance of these observations to the loss of the N611 glycan in the context of 293Trhm infection is unclear. Interestingly, the magnitude of mutant enrichment within both the N611 glycan region, and the CD4-contacting helix was similar between 293Trhm and 293Thu selections, suggesting that these mutations function by a mechanism that is similarly advantageous to both human and macaque infection.
Our DMS did not select for mutations that have previously been found to improve usage of macaque CD4 and CCR5. This observation can be attributed in part to incompatibilities between these mutants and the BF520 Env background, as was observed for A204E. For other mutations, such as S375Y that clearly improved entry of BF520 pseudovirus into 293Trhm cells, the lack of enrichment is harder to explain. One possibility is that these mutations are detrimental in the BF520 background during multi-cycle replication and were removed from the library during the initial low MOI passage. Regardless of the underlying cause, depletion of these mutations indicates that our DMS did not quantify the effects of every possible BF520 Env mutation on the usage of macaque receptors. Future studies that employ different experimental methods are therefore likely to discover additional Env mutants that enhance entry into macaque cells.
A notable limitation of the experiments we performed is that we characterized N/CHR mutants in engineered 293T cells lines that may express CD4 and CCR5 at different levels than human and macaque primary cells. The phenotypes of the Env mutants described here may therefore differ during infection of other cell types, and fitness tradeoffs associated with these mutations may arise during multi-cycle replication in a more clinically relevant context. Additionally, while our data suggest that N/CHR mutations enhance intrinsic Env reactivity, a mechanistic understanding of the effects of these mutations will benefit from further investigations into their effects on CD4 binding and fusion kinetics during infection of relevant cell types.
This work uncovered many novel determinants of entry into cells bearing both human and macaque receptors. Mutations within the N/CHRs had the largest benefit to usage of receptors from both species and may function by promoting a state 2 Env conformation with heightened reactivity to CD4 binding. Identification of these mutations provides insight into the mechanisms of viral entry and interactions with cellular receptors and may be of use in developing more robust SHIV models.

Supplementary Materials

The following are available online at https://www.mdpi.com/1999-4915/12/2/241/s1, Figure S1: Wild-type BF520 infection of 293T cells engineered to express either macaque or human CD4/CCR5 receptors, Figure S2: Expression of CD4 and CCR5 on 293Trhm and 293Thu cells. Table S1: Fold-change relative to the wild-type residue for all mutants that were enriched greater than 10-fold following either 293Trhm or 293Thu infection in the DMS.

Author Contributions

Conceptualization, J.I.R., J.D.B., A.S.D., and J.O.; software, J.I.R. and J.D.B.; formal analysis, J.I.R.; investigation, J.I.R., N.A.C., and A.S.D; writing—original draft preparation, J.I.R. and J.O.; writing—review and editing, J.I.R., J.O., A.S.D, and J.D.B. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by NIH DP1 DA039543, awarded to J.O., and Damon Runyon Cancer Research Institute Postdoctoral Fellowship 2306-17, awarded to J.I.R.

Acknowledgments

We thank Daryl Humes, Amit Sharma, and Hugh Haddox for helpful discussions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Marsh, M.; Helenius, A. Virus entry: Open sesame. Cell 2006, 124, 729–740. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Overbaugh, J.; Miller, A.D.; Eiden, M.V. Receptors and entry cofactors for retroviruses include single and multiple transmembrane-spanning proteins as well as newly described glycophosphatidylinositol-anchored and secreted proteins. Microbiol. Mol. Biol. Rev. 2001, 65, 371–389. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Humes, D.; Emery, S.; Laws, E.; Overbaugh, J. A species-specific amino acid difference in the macaque CD4 receptor restricts replication by global circulating HIV-1 variants representing viruses from recent infection. J. Virol. 2012, 86, 12472–12483. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  4. Humes, D.; Overbaugh, J. Adaptation of subtype a human immunodeficiency virus type 1 envelope to pig-tailed macaque cells. J. Virol. 2011, 85, 4409–4420. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Warren, C.J.; Meyerson, N.R.; Dirasantha, O.; Feldman, E.R.; Wilkerson, G.K.; Sawyer, S.L. Selective use of primate CD4 receptors by HIV-1. PLoS Biol. 2019, 17, e3000304. [Google Scholar] [CrossRef] [PubMed]
  6. Meyerson, N.R.; Sharma, A.; Wilkerson, G.K.; Overbaugh, J.; Sawyer, S.L. Identification of Owl Monkey CD4 Receptors Broadly Compatible with Early-Stage HIV-1 Isolates. J. Virol. 2015, 89, 8611–8622. [Google Scholar] [CrossRef] [Green Version]
  7. Sharma, A.; Boyd, D.F.; Overbaugh, J. Development of SHIVs with circulating, transmitted HIV-1 variants. J. Med. Primatol. 2015, 44, 296–300. [Google Scholar] [CrossRef]
  8. Hatziioannou, T.; Evans, D.T. Animal models for HIV/AIDS research. Nat. Rev. Microbiol. 2012, 10, 852–867. [Google Scholar] [CrossRef]
  9. Cayabyab, M.; Karlsson, G.B.; Etemad-Moghadam, B.A.; Hofmann, W.; Steenbeke, T.; Halloran, M.; Fanton, J.W.; Axthelm, M.K.; Letvin, N.L.; Sodroski, J.G. Changes in human immunodeficiency virus type 1 envelope glycoproteins responsible for the pathogenicity of a multiply passaged simian-human immunodeficiency virus (SHIV-HXBc2). J. Virol. 1999, 73, 976–984. [Google Scholar] [CrossRef] [Green Version]
  10. Boyd, D.F.; Peterson, D.; Haggarty, B.S.; Jordan, A.P.O.; Hogan, M.J.; Goo, L.; Hoxie, J.A.; Overbaugh, J. Mutations in HIV-1 envelope that enhance entry with the macaque CD4 receptor alter antibody recognition by disrupting quaternary interactions within the trimer. J. Virol. 2015, 89, 894–907. [Google Scholar] [CrossRef] [Green Version]
  11. Li, H.; Wang, S.; Kong, R.; Ding, W.; Lee, F.-H.; Parker, Z.; Kim, E.; Learn, G.H.; Hahn, P.; Policicchio, B.; et al. Envelope residue 375 substitutions in simian-human immunodeficiency viruses enhance CD4 binding and replication in rhesus macaques. Proc. Natl. Acad. Sci. USA 2016, 113, E3413. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Del Prete, G.Q.; Keele, B.F.; Fode, J.; Thumar, K.; Swanstrom, A.E.; Rodriguez, A.; Raymond, A.; Estes, J.D.; LaBranche, C.C.; Montefiori, D.C.; et al. A single gp120 residue can affect HIV-1 tropism in macaques. PLoS Pathog. 2017, 13, e1006572. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  13. Fowler, D.M.; Fields, S. Deep mutational scanning: A new style of protein science. Nat. Methods 2014, 11, 801–807. [Google Scholar] [CrossRef] [PubMed]
  14. Haddox, H.K.; Dingens, A.S.; Bloom, J.D. Experimental Estimation of the Effects of All Amino-Acid Mutations to HIV’s Envelope Protein on Viral Replication in Cell Culture. PLoS Pathog. 2016, 12, e1006114. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Haddox, H.K.; Dingens, A.S.; Hilton, S.K.; Overbaugh, J.; Bloom, J.D. Mapping mutational effects along the evolutionary landscape of HIV envelope. eLife 2018, 7, e34420. [Google Scholar] [CrossRef] [PubMed]
  16. Dingens, A.S.; Acharya, P.; Haddox, H.K.; Rawi, R.; Xu, K.; Chuang, G.-Y.; Wei, H.; Zhang, B.; Mascola, J.R.; Carragher, B.; et al. Complete functional mapping of infection- and vaccine-elicited antibodies against the fusion peptide of HIV. PLoS Pathog 2018, 14, e1007159. [Google Scholar] [CrossRef] [Green Version]
  17. Dingens, A.S.; Haddox, H.K.; Overbaugh, J.; Bloom, J.D. Comprehensive Mapping of HIV-1 Escape from a Broadly Neutralizing Antibody. Cell Host Microbe 2017, 21, 777–787. [Google Scholar] [CrossRef] [Green Version]
  18. Dingens, A.S.; Arenz, D.; Overbaugh, J.; Bloom, J.D. Massively Parallel Profiling of HIV-1 Resistance to the Fusion Inhibitor Enfuvirtide. Viruses 2019, 11, 439. [Google Scholar] [CrossRef] [Green Version]
  19. Simonich, C.A.; Williams, K.L.; Verkerke, H.P.; Williams, J.A.; Nduati, R.; Lee, K.K.; Overbaugh, J. HIV-1 Neutralizing Antibodies with Limited Hypermutation from an Infant. Cell 2016, 166, 77–87. [Google Scholar] [CrossRef] [Green Version]
  20. Sanders, R.W.; Derking, R.; Cupo, A.; Julien, J.-P.; Yasmeen, A.; de Val, N.; Kim, H.J.; Blattner, C.; de la Peña, A.T.; Korzun, J.; et al. A next-generation cleaved, soluble HIV-1 Env trimer, BG505 SOSIP.664 gp140, expresses multiple epitopes for broadly neutralizing but not non-neutralizing antibodies. PLoS Pathog. 2013, 9, e1003618. [Google Scholar] [CrossRef] [Green Version]
  21. Korber, B.T.; Foley, B.T.; Kuiken, C.L.; Pillai, S.K.; Sodroski, J.G. Numbering Positions in HIV Relative to HXB2CG. Human Retrovir. AIDS Compend. 1999, 102–111. [Google Scholar]
  22. Goo, L.; Chohan, V.; Nduati, R.; Overbaugh, J. Early development of broadly neutralizing antibodies in HIV-1-infected infants. Nat. Med. 2014, 20, 655–658. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Sarzotti-Kelsoe, M.; Bailer, R.T.; Turk, E.; Lin, C.-L.; Bilska, M.; Greene, K.M.; Gao, H.; Todd, C.A.; Ozaki, D.A.; Seaman, M.S.; et al. Optimization and validation of the TZM-bl assay for standardized assessments of neutralizing antibodies against HIV-1. J. Immunol. Methods 2014, 409, 131–146. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Garces, F.; Lee, J.H.; de Val, N.; de la Peña, A.T.; Kong, L.; Puchades, C.; Hua, Y.; Stanfield, R.L.; Burton, D.R.; Moore, J.P.; et al. Affinity Maturation of a Potent Family of HIV Antibodies Is Primarily Focused on Accommodating or Avoiding Glycans. Immunity 2015, 43, 1053–1063. [Google Scholar] [CrossRef] [Green Version]
  25. Pettersen, E.F.; Goddard, T.D.; Huang, C.C.; Couch, G.S.; Greenblatt, D.M.; Meng, E.C.; Ferrin, T.E. UCSF Chimera—A visualization system for exploratory research and analysis. J. Comput. Chem. 2004, 25, 1605–1612. [Google Scholar] [CrossRef] [Green Version]
  26. Bloom, J.D. Software for the analysis and visualization of deep mutational scanning data. BMC Bioinform. 2015, 16, 168. [Google Scholar] [CrossRef] [Green Version]
  27. Lee, J.M.; Huddleston, J.; Doud, M.B.; Hooper, K.A.; Wu, N.C.; Bedford, T.; Bloom, J.D. Deep mutational scanning of hemagglutinin helps predict evolutionary fates of human H3N2 influenza variants. Proc. Natl. Acad. Sci. USA 2018, 115, E8276–E8285. [Google Scholar] [CrossRef] [Green Version]
  28. Asmal, M.; Lane, S.; Tian, M.; Nickel, G.; Venner, C.; Dirk, B.; Dikeakos, J.; Luedemann, C.; Mach, L.; Balachandran, H.; et al. Pathogenic infection of Rhesus macaques by an evolving SIV-HIV derived from CCR5-using envelope genes of acute HIV-1 infections. Virology 2016, 499, 298–312. [Google Scholar] [CrossRef]
  29. Liu, Q.; Acharya, P.; Dolan, M.A.; Zhang, P.; Guzzo, C.; Lu, J.; Kwon, A.; Gururani, D.; Miao, H.; Bylund, T.; et al. Quaternary contact in the initial interaction of CD4 with the HIV-1 envelope trimer. Nat. Struct. Mol. Biol. 2017, 24, 370–378. [Google Scholar] [CrossRef] [Green Version]
  30. Wu, X.; Parast, A.B.; Richardson, B.A.; Nduati, R.; John-Stewart, G.; Mbori-Ngacha, D.; Rainwater, S.M.J.; Overbaugh, J. Neutralization escape variants of human immunodeficiency virus type 1 are transmitted from mother to infant. J. Virol. 2006, 80, 835–844. [Google Scholar] [CrossRef] [Green Version]
  31. Lyumkis, D.; Julien, J.-P.; de Val, N.; Cupo, A.; Potter, C.S.; Klasse, P.-J.; Burton, D.R.; Sanders, R.W.; Moore, J.P.; Carragher, B.; et al. Cryo-EM structure of a fully glycosylated soluble cleaved HIV-1 envelope trimer. Science 2013, 342, 1484–1490. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Pancera, M.; Zhou, T.; Druz, A.; Georgiev, I.S.; Soto, C.; Gorman, J.; Huang, J.; Acharya, P.; Chuang, G.-Y.; Ofek, G.; et al. Structure and immune recognition of trimeric pre-fusion HIV-1 Env. Nature 2014, 514, 455–461. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Huang, J.; Kang, B.H.; Pancera, M.; Lee, J.H.; Tong, T.; Feng, Y.; Imamichi, H.; Georgiev, I.S.; Chuang, G.-Y.; Druz, A.; et al. Broad and potent HIV-1 neutralization by a human antibody that binds the gp41-gp120 interface. Nature 2014, 515, 138–142. [Google Scholar] [CrossRef] [PubMed]
  34. Stewart-Jones, G.B.E.; Soto, C.; Lemmin, T.; Chuang, G.-Y.; Druz, A.; Kong, R.; Thomas, P.V.; Wagh, K.; Zhou, T.; Behrens, A.-J.; et al. Trimeric HIV-1-Env Structures Define Glycan Shields from Clades A, B, and G. Cell 2016, 165, 813–826. [Google Scholar] [CrossRef] [Green Version]
  35. De Taeye, S.W.; Ozorowski, G.; Torrents de la Peña, A.; Guttman, M.; Julien, J.-P.; van den Kerkhof, T.L.G.M.; Burger, J.A.; Pritchard, L.K.; Pugach, P.; Yasmeen, A.; et al. Immunogenicity of Stabilized HIV-1 Envelope Trimers with Reduced Exposure of Non-neutralizing Epitopes. Cell 2015, 163, 1702–1715. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. Montefiori, D.C.; Roederer, M.; Morris, L.; Seaman, M.S. Neutralization tiers of HIV-1. Curr. Opin. HIV AIDS 2018, 13, 128–136. [Google Scholar] [CrossRef]
  37. Chakrabarti, B.K.; Walker, L.M.; Guenaga, J.F.; Ghobbeh, A.; Poignard, P.; Burton, D.R.; Wyatt, R.T. Direct Antibody Access to the HIV-1 Membrane-Proximal External Region Positively Correlates with Neutralization Sensitivity. J. Virol. 2011, 85, 8217–8226. [Google Scholar] [CrossRef] [Green Version]
  38. Alam, S.M.; Morelli, M.; Dennison, S.M.; Liao, H.X.; Zhang, R.; Xia, S.-M.; Rits-Volloch, S.; Sun, L.; Harrison, S.C.; Haynes, B.F.; et al. Role of HIV membrane in neutralization by two broadly neutralizing antibodies. Proc. Natl. Acad. Sci. USA 2009, 106, 20234–20239. [Google Scholar] [CrossRef] [Green Version]
  39. Frey, G.; Peng, H.; Rits-Volloch, S.; Morelli, M.; Cheng, Y.; Chen, B. A fusion-intermediate state of HIV-1 gp41 targeted by broadly neutralizing antibodies. Proc. Natl. Acad. Sci. USA 2008, 105, 3739–3744. [Google Scholar] [CrossRef] [Green Version]
  40. Haim, H.; Strack, B.; Kassa, A.; Madani, N.; Wang, L.; Courter, J.R.; Princiotto, A.; McGee, K.; Pacheco, B.; Seaman, M.S.; et al. Contribution of intrinsic reactivity of the HIV-1 envelope glycoproteins to CD4-independent infection and global inhibitor sensitivity. PLoS Pathog. 2011, 7, e1002101. [Google Scholar] [CrossRef] [Green Version]
  41. Herschhorn, A.; Ma, X.; Gu, C.; Ventura, J.D.; Castillo-Menendez, L.; Melillo, B.; Terry, D.S.; Smith, A.B.; Blanchard, S.C.; Munro, J.B.; et al. Release of gp120 Restraints Leads to an Entry-Competent Intermediate State of the HIV-1 Envelope Glycoproteins. Mbio 2016, 7, 1884. [Google Scholar] [CrossRef] [Green Version]
  42. Herschhorn, A.; Gu, C.; Moraca, F.; Ma, X.; Farrell, M.; Smith, A.B.; Pancera, M.; Kwong, P.D.; Schön, A.; Freire, E.; et al. The β20-β21 of gp120 is a regulatory switch for HIV-1 Env conformational transitions. Nat. Commun. 2017, 8, 1049. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Keller, P.W.; Morrison, O.; Vassell, R.; Weiss, C.D. HIV-1 gp41 Residues Modulate CD4-Induced Conformational Changes in the Envelope Glycoprotein and Evolution of a Relaxed Conformation of gp120. J. Virol. 2018, 92, e00583. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  44. Kassa, A.; Finzi, A.; Pancera, M.; Courter, J.R.; Smith, A.B.; Sodroski, J. Identification of a human immunodeficiency virus type 1 envelope glycoprotein variant resistant to cold inactivation. J. Virol. 2009, 83, 4476–4488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Sivaraman, V.; Zhang, L.; Meissner, E.G.; Jeffrey, J.L.; Su, L. The Heptad Repeat 2 Domain Is a Major Determinant for Enhanced Human Immunodeficiency Virus Type 1 (HIV-1) Fusion and Pathogenicity of a Highly Pathogenic HIV-1 Env. J. Virol. 2009, 83, 11715–11725. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  46. Yu, D.; Ding, X.; Liu, Z.; Wu, X.; Zhu, Y.; Wei, H.; Chong, H.; Cui, S.; He, Y. Molecular mechanism of HIV-1 resistance to sifuvirtide, a clinical trial-approved membrane fusion inhibitor. J. Biol. Chem. 2018. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  47. Taylor, B.M.; Foulke, J.S.; Flinko, R.; Heredia, A.; DeVico, A.; Reitz, M. An alteration of human immunodeficiency virus gp41 leads to reduced CCR5 dependence and CD4 independence. J. Virol. 2008, 82, 5460–5471. [Google Scholar] [CrossRef] [Green Version]
  48. Muster, T.; Steindl, F.; Purtscher, M.; Trkola, A.; Klima, A.; Himmler, G.; Rüker, F.; Katinger, H. A conserved neutralizing epitope on gp41 of human immunodeficiency virus type 1. J. Virol. 1993, 67, 6642–6647. [Google Scholar] [CrossRef] [Green Version]
  49. Espy, N.; Pacheco, B.; Sodroski, J. Adaptation of HIV-1 to cells with low expression of the CCR5 coreceptor. Virology 2017, 508, 90–107. [Google Scholar] [CrossRef]
  50. Pacheco, B.; Alsahafi, N.; Debbeche, O.; Prévost, J.; Ding, S.; Chapleau, J.-P.; Herschhorn, A.; Madani, N.; Princiotto, A.; Melillo, B.; et al. Residues in the gp41 Ectodomain Regulate HIV-1 Envelope Glycoprotein Conformational Transitions Induced by gp120-Directed Inhibitors. J. Virol. 2017, 91, e02219. [Google Scholar] [CrossRef] [Green Version]
  51. Lounková, A.; Kosla, J.; Přikryl, D.; Štafl, K.; Kučerová, D.; Svoboda, J. Retroviral host range extension is coupled with Env-activating mutations resulting in receptor-independent entry. Proc. Natl. Acad. Sci. USA 2017, 114, E5148–E5157. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. LaBranche, C.C.; Sauter, M.M.; Haggarty, B.S.; Vance, P.J.; Romano, J.; Hart, T.K.; Bugelski, P.J.; Marsh, M.; Hoxie, J.A. A single amino acid change in the cytoplasmic domain of the simian immunodeficiency virus transmembrane molecule increases envelope glycoprotein expression on infected cells. J. Virol. 1995, 69, 5217–5227. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. LaBranche, C.C.; Sauter, M.M.; Haggarty, B.S.; Vance, P.J.; Romano, J.; Hart, T.K.; Bugelski, P.J.; Hoxie, J.A. Biological, molecular, and structural analysis of a cytopathic variant from a molecularly cloned simian immunodeficiency virus. J. Virol. 1994, 68, 5509–5522. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Sauter, M.M.; Pelchen-Matthews, A.; Bron, R.; Marsh, M.; LaBranche, C.C.; Vance, P.J.; Romano, J.; Haggarty, B.S.; Hart, T.K.; Lee, W.M.; et al. An internalization signal in the simian immunodeficiency virus transmembrane protein cytoplasmic domain modulates expression of envelope glycoproteins on the cell surface. J. Cell Biol. 1996, 132, 795–811. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Leaman, D.P.; Zwick, M.B. Increased functional stability and homogeneity of viral envelope spikes through directed evolution. PLoS Pathog. 2013, 9, e1003184. [Google Scholar] [CrossRef] [Green Version]
  56. Agrawal, N.; Leaman, D.P.; Rowcliffe, E.; Kinkead, H.; Nohria, R.; Akagi, J.; Bauer, K.; Du, S.X.; Whalen, R.G.; Burton, D.R.; et al. Functional stability of unliganded envelope glycoprotein spikes among isolates of human immunodeficiency virus type 1 (HIV-1). PLoS ONE 2011, 6, e21339. [Google Scholar] [CrossRef] [Green Version]
  57. Araújo, L.A.L.; Almeida, S.E.M. HIV-1 diversity in the envelope glycoproteins: Implications for viral entry inhibition. Viruses 2013, 5, 595–604. [Google Scholar] [CrossRef] [Green Version]
  58. Yen, P.-J.; Herschhorn, A.; Haim, H.; Salas, I.; Gu, C.; Sodroski, J.; Gabuzda, D. Loss of a conserved N-linked glycosylation site in the simian immunodeficiency virus envelope glycoprotein V2 region enhances macrophage tropism by increasing CD4-independent cell-to-cell transmission. J. Virol. 2014, 88, 5014–5028. [Google Scholar] [CrossRef] [Green Version]
  59. Pikora, C.; Wittish, C.; Desrosiers, R.C. Identification of two N-linked glycosylation sites within the core of the simian immunodeficiency virus glycoprotein whose removal enhances sensitivity to soluble CD4. J. Virol. 2005, 79, 12575–12583. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Deep mutational scanning of BF520 Env. (a) Schematic of the deep mutational scanning method used. Mutant virus pools were generated by transfection of mutant proviral plasmids followed by low MOI passage. Virus pools were then used to infect either 293Trhm or 293Thu cells. Following infection, amino acid preferences at each Env site were inferred by deep sequencing Env variants from pre- and post-selected samples. Preferences were used to infer Smut/wt values for each mutant as log2 (mutant preference / wild-type preference), as depicted by the logoplot. (b) The bar graph reports the frequencies of nonsynonymous, stop, and synonymous mutations, as indicated in the key to the right. Labels below each bar indicate the stage of library selection (DNA-mut: mutant plasmid DNA; HU-mut and Rh-mut: mutant virus following selection in 293Thu or 293Trhm; WT-DNA: wild-type plasmid DNA control; WT-inf: wild-type virus control following selecting in 293Thu cells). (c,d) Pearson correlations between Smut/wt values from replicate 293Thu infections (c) and 293Trhm infections (d). (e) Pearson correlation between average Smut/wt values from 293Trhm and 293Thu infections. Points above the diagonal grey line indicate mutants that were more advantageous during 293Trhm infection than during 293Thu infection.
Figure 1. Deep mutational scanning of BF520 Env. (a) Schematic of the deep mutational scanning method used. Mutant virus pools were generated by transfection of mutant proviral plasmids followed by low MOI passage. Virus pools were then used to infect either 293Trhm or 293Thu cells. Following infection, amino acid preferences at each Env site were inferred by deep sequencing Env variants from pre- and post-selected samples. Preferences were used to infer Smut/wt values for each mutant as log2 (mutant preference / wild-type preference), as depicted by the logoplot. (b) The bar graph reports the frequencies of nonsynonymous, stop, and synonymous mutations, as indicated in the key to the right. Labels below each bar indicate the stage of library selection (DNA-mut: mutant plasmid DNA; HU-mut and Rh-mut: mutant virus following selection in 293Thu or 293Trhm; WT-DNA: wild-type plasmid DNA control; WT-inf: wild-type virus control following selecting in 293Thu cells). (c,d) Pearson correlations between Smut/wt values from replicate 293Thu infections (c) and 293Trhm infections (d). (e) Pearson correlation between average Smut/wt values from 293Trhm and 293Thu infections. Points above the diagonal grey line indicate mutants that were more advantageous during 293Trhm infection than during 293Thu infection.
Viruses 12 00241 g001
Figure 2. Smut/wt values of amino acids enriched during 293Trhm infection. BF520 Env wild-type amino acid identities are shown below the axes for residues 32–702. Logoplots above the axes indicate mutants that were enriched relative to wild-type residues (Smut/wt > 0) following infection of 293Trhm cells in our DMS. The height of each letter is proportional to the Smut/wt value for that amino acid mutant following 293Trhm infection, averaged over both replicate infections. Site numbers correspond to HXB2 numbering. The underlay bar indicates Env region.
Figure 2. Smut/wt values of amino acids enriched during 293Trhm infection. BF520 Env wild-type amino acid identities are shown below the axes for residues 32–702. Logoplots above the axes indicate mutants that were enriched relative to wild-type residues (Smut/wt > 0) following infection of 293Trhm cells in our DMS. The height of each letter is proportional to the Smut/wt value for that amino acid mutant following 293Trhm infection, averaged over both replicate infections. Site numbers correspond to HXB2 numbering. The underlay bar indicates Env region.
Viruses 12 00241 g002
Figure 3. Infection of 293Trhm, 293TCCR5, and 293Thu cells by point mutant pseudoviruses. Each panel reports the percentage of the indicated cells that were GFP positive 48 h post-infection (MOI of 5) with the indicated wild-type or mutant viruses. The Env strain in which mutants were engineered is shown at the bottom of each panel. Every point in each panel represents one experimental replicate, and the median, indicated as a horizontal black bar, is representative of at least two independent experiments. Asterisks indicate mutants whose infection was significantly different than the wild-type (p < 0.05, Welch’s t-test). (a) Infection of 293Trhm by point mutants in the BF520 Env background. (b) Infection of 293Thu by point mutants in the BF520 Env background. (c) Infection of 293Trhm by point mutants in the BG505 Env background. Note that infection of these cells by wild-type BG505 is nearly 10-fold less efficient than by wild-type BF520. (d) Infection of 293TrhmCCR5 by point mutants in the BF520 Env background. QA255 CD4iB is a CD4-independent Env described in [10] used here as a positive control.
Figure 3. Infection of 293Trhm, 293TCCR5, and 293Thu cells by point mutant pseudoviruses. Each panel reports the percentage of the indicated cells that were GFP positive 48 h post-infection (MOI of 5) with the indicated wild-type or mutant viruses. The Env strain in which mutants were engineered is shown at the bottom of each panel. Every point in each panel represents one experimental replicate, and the median, indicated as a horizontal black bar, is representative of at least two independent experiments. Asterisks indicate mutants whose infection was significantly different than the wild-type (p < 0.05, Welch’s t-test). (a) Infection of 293Trhm by point mutants in the BF520 Env background. (b) Infection of 293Thu by point mutants in the BF520 Env background. (c) Infection of 293Trhm by point mutants in the BG505 Env background. Note that infection of these cells by wild-type BG505 is nearly 10-fold less efficient than by wild-type BF520. (d) Infection of 293TrhmCCR5 by point mutants in the BF520 Env background. QA255 CD4iB is a CD4-independent Env described in [10] used here as a positive control.
Viruses 12 00241 g003
Figure 4. Smut/wt values of amino acids enriched during 293Thu infection. BF520 Env wild-type amino acid identities are shown below the axes for residues 32–702. Logoplots above the axes indicate mutants that were enriched relative to wild-type residues (Smut/wt > 0) following infection of 293Thu cells in our deep mutational scanning (DMS). The height of each letter is proportional to the Smut/wt value for that amino acid mutant following 293Thu infection, averaged over both replicate infections. Site numbers correspond to HXB2 numbering. The underlay bar indicates Env region.
Figure 4. Smut/wt values of amino acids enriched during 293Thu infection. BF520 Env wild-type amino acid identities are shown below the axes for residues 32–702. Logoplots above the axes indicate mutants that were enriched relative to wild-type residues (Smut/wt > 0) following infection of 293Thu cells in our deep mutational scanning (DMS). The height of each letter is proportional to the Smut/wt value for that amino acid mutant following 293Thu infection, averaged over both replicate infections. Site numbers correspond to HXB2 numbering. The underlay bar indicates Env region.
Viruses 12 00241 g004
Figure 5. Neutralization of BF520 and BG505 point mutants. (a) IC50 values of the indicated drugs and mAbs with BF520 and BG505 Env point mutants. Values report the mean of two independent replicates and are representative of at least two distinct experiments. (b,c) Neutralization curves illustrating the effects of Env point mutants on neutralization by 10E8 (b) and T-20 (c). Each data point indicates the mean of two independent replicates and error bars report standard deviation.
Figure 5. Neutralization of BF520 and BG505 point mutants. (a) IC50 values of the indicated drugs and mAbs with BF520 and BG505 Env point mutants. Values report the mean of two independent replicates and are representative of at least two distinct experiments. (b,c) Neutralization curves illustrating the effects of Env point mutants on neutralization by 10E8 (b) and T-20 (c). Each data point indicates the mean of two independent replicates and error bars report standard deviation.
Viruses 12 00241 g005
Figure 6. Cold inactivation of BF520 mutant viruses. Viruses were incubated on ice for the indicated period of time then used to infect TZM-bl reporter cells. Percent infectivity of wild-type and mutant strains are shown as a function of time incubated on ice. Each data point indicates the average of four independent replicates and error bars report the standard deviation. Infectivity of viruses not incubated on ice was set to 100%.
Figure 6. Cold inactivation of BF520 mutant viruses. Viruses were incubated on ice for the indicated period of time then used to infect TZM-bl reporter cells. Percent infectivity of wild-type and mutant strains are shown as a function of time incubated on ice. Each data point indicates the average of four independent replicates and error bars report the standard deviation. Infectivity of viruses not incubated on ice was set to 100%.
Viruses 12 00241 g006
Figure 7. Location and predicted interactions of N-terminal heptad repeat (NHR) and C-terminal heptad repeat (CHR) mutations in the unliganded trimer. The BG505 SOSIP trimer structure (PDB 5CEZ) is shown. Gp120 and gp41 residues of one protomer are shown in bright green and blue respectively. Residues from other protomers are shown in dark green. Upper and lower insets display the regions surrounding the NHR and CHR of one gp41 subunit, respectively. In each inset, residues strongly enriched during the DMS are shown in red. N/CHR residue side chains are shown for amino acids that form predicted H-bond contacts with neighboring residues, and H-bonds are depicted with dotted lines. Side chains are shown for all residues on adjacent subunits that are within 4 angstroms of highly enriched N/CHR residues.
Figure 7. Location and predicted interactions of N-terminal heptad repeat (NHR) and C-terminal heptad repeat (CHR) mutations in the unliganded trimer. The BG505 SOSIP trimer structure (PDB 5CEZ) is shown. Gp120 and gp41 residues of one protomer are shown in bright green and blue respectively. Residues from other protomers are shown in dark green. Upper and lower insets display the regions surrounding the NHR and CHR of one gp41 subunit, respectively. In each inset, residues strongly enriched during the DMS are shown in red. N/CHR residue side chains are shown for amino acids that form predicted H-bond contacts with neighboring residues, and H-bonds are depicted with dotted lines. Side chains are shown for all residues on adjacent subunits that are within 4 angstroms of highly enriched N/CHR residues.
Viruses 12 00241 g007

Share and Cite

MDPI and ACS Style

Roop, J.I.; Cassidy, N.A.; Dingens, A.S.; Bloom, J.D.; Overbaugh, J. Identification of HIV-1 Envelope Mutations that Enhance Entry Using Macaque CD4 and CCR5. Viruses 2020, 12, 241. https://doi.org/10.3390/v12020241

AMA Style

Roop JI, Cassidy NA, Dingens AS, Bloom JD, Overbaugh J. Identification of HIV-1 Envelope Mutations that Enhance Entry Using Macaque CD4 and CCR5. Viruses. 2020; 12(2):241. https://doi.org/10.3390/v12020241

Chicago/Turabian Style

Roop, Jeremy I., Noah A. Cassidy, Adam S. Dingens, Jesse D. Bloom, and Julie Overbaugh. 2020. "Identification of HIV-1 Envelope Mutations that Enhance Entry Using Macaque CD4 and CCR5" Viruses 12, no. 2: 241. https://doi.org/10.3390/v12020241

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop