Antigenic heterogeneity of Env spikes on PV
To dissect antigenic heterogeneity in the PV population, we conjugated bNAbs to Sepharose beads to deplete inocula of virions most reactive with the respective bNAb (Fig. 1C). The depletion is designed to be partial; the number and density of Env spikes with at least some antigenicity will determine the avidity of virion capture (Fig. 1C). Meaningful depletions were possible only for B41 PV, which showed the substantial PF; depletion of BG505 PV left negligible and insufficient infectivity for a mechanistic dissection.
Henceforth, we focus on B41 neutralization potency and efficacy, in particular seeking explanations in antigenic heterogeneity of its large PF with PGT151. Depletion of the B41 PV by PGT145- and PGT151-conjugated beads gave differentially neutralized fractions (Fig. 2). Neutralization by PGT145 was consistently but only moderately more potent for PGT151- than PGT145-depleted PV, mock-depleted PV falling in between; the neutralization efficacies were close to 100% after the three depletions (Fig. 2). Neither neutralization by 2G12 nor VRC01, directed to a CD4-binding-site (CD4bs) epitope, was affected by the depletions. In contrast, PGT151 neutralization was 2 orders of magnitude more potent against PGT145- than PGT151-depleted PV, the potency for mock-depleted PV falling between, somewhat closer to the former than the latter. The corresponding PF was reduced for PGT145- compared with mock-depleted PV, whereas the PGT151 neutralization of PGT151-depleted PV was so diminished that no PF plateau could accurately be extrapolated. ASC202, directed to an interface epitope overlapping that for PGT151, neutralized with potencies ranking as for PGT151 against the three depleted PV fractions but with smaller shifts and without yielding detectable PFs for any depletion. Finally, N123-VRC34.01 (hereafter VRC34.01), also specific for an interface epitope [48], weakly and partially neutralized the PGT145-depleted PV (40% at 50 µg/ml) but had no effect against the other two fractions.
Eight sera from B41 SOSIP.664-immunized rabbits [17, 49] yielded a range of PFs against mock-depleted PV: 5–60% (Fig. 3). Neutralization potency and efficacy were consistently lower against the PGT145-depleted than the PGT151-depleted PV, the curves for the mock-depleted PV falling in between. Just as the sizes of the PFs varied greatly among the sera, however, so did the differences in PFs with each serum for the three depletions.
Two monoclonal NAbs (mNAbs), 13A and 16D, were isolated from two of these B41-immunized rabbits (5713 and 5716 [50, 51]). Like the sera, they neutralize the autologous Tier-2 PV well, and the neutralization of a mutant with a knock-in N289-glycan is markedly reduced [50, 51]. Congruently with the results for the sera, neutralization potency and efficacy were lower against the PGT145-depleted than the PGT151-depleted PV, whereas curves for the mock-depleted PV fell in between; the corresponding PF differences were clear with both mNAbs (Fig. 4). Although it is not feasible to elute infectious PV from the bNAbs on the beads, the depletion results suggest that the more PGT151-reactive Env on the PV virions, which preferentially tethers them to the beads, exposes this off-target glycan-hole epitope less well than does the more PGT145-reactive Env.
Differential purification and antigenicity of soluble native-like Env SOSIP.664 trimers
The 2G12, PGT145, and PGT151 epitopes, as well as the N289 residue, which is part of a glycan-hole epitope on wild-type B41 Env, are shown in their oligomeric-structural contexts in Fig. 5A. The model is based on the crystal structure of the B41 SOSIP.664 trimer, with the addition of Man9 oligomannose glycans to all potential N-linked glycosylation sites (PNGS) [50, 52]. Figure 5A illustrates how the glycan-knock-in mutation at N289 fills a defect in the glycan shield, thereby blocking a prominent epitope for autologous NAbs [49]. That soluble trimers with structures similar to the native one shown in Fig. 5A can be obtained by 2G12- and PGT145-affinity purification has been shown multiple times for BG505 and B41 SOSIP.664 trimers by negative-stain electron microscopy (NS-EM) [17, 22–25, 30, 53]. As found here, the PGT151-purified B41 SOSIP.664 trimer also appeared trimeric by BN-PAGE (Fig. 5B). NS-EM 2D-class averages revealed intact trimer molecules with overall movement of the 3 protomeric lobes relative to the center, suggesting substantial conformational “breathing” within the basic native structure, as is typical for B41 trimers (Fig. 5C) [25]. Particles lacking a central triangular mass, a defect that characterizes non-native structures, such as those of uncleaved Env [54], were not observed. We have thus demonstrated that PGT151-affinity purification also yields native-like B41 SOSIP.664 trimer molecules, which enabled us to perform bNAb-binding analyses with differentially purified trimers and explore the antigenic heterogeneity further.
In pursuit of explanations of the large PF specifically with PGT151 against B41 PV, we first used ELISA to explore the binding of bNAbs to B41 SOSIP.664 trimer purified by 2G12-, PGT145-, or PGT151-affinity chromatography and thereafter by SEC (SI Fig. 1). 2G12 bound similarly to the three trimer preparations. PGT145 bound strongly to PGT145- and 2G12-purified trimers but considerably less well to PGT151-purified trimer. Conversely, PGT151 bound strongly only to the trimer purified with PGT151 itself and weakly to 2G12- and PGT145-purified trimer. The differential binding of VRC34.01, which like PGT151 is directed to an interface epitope [48], was similar to that of PGT151. The two autologous rabbit mNAbs resembled each other in binding profile: strong binding to PGT145- and 2G12-purified, and substantially weaker to PGT151-purified trimer. These results suggest that PGT151 has a high affinity for a subpopulation of the antigenically heterogeneous trimer molecules, which exposes the 289-glycan-hole epitope less well than the most PGT145-reactive subpopulation.
The bNAb binding to differentially purified B41 and BG505 SOSIP.664 trimers was then compared by SPR (Fig. 6). Binding to the BG505 SOSIP.664 trimer was indistinguishable after 2G12, PGT145, and PGT151 purification by the bNAbs VRC01, VRC34.01, 35O22 (gp120-gp41 interface [55]), and 3BC315 (gp41, inter-protomeric [56]), as well as by the two purification bNAbs, PGT145 and PGT151, whereas 2G12 and PGT121 binding subtly favored the 2G12-purifed form.
In contrast, the B41 SOSIP.664 trimer showed a wide range of distinct antigenicities resulting from the different purifications (Fig. 6A). VRC01 bound better to the PGT151-purified trimer than to the other two. 2G12 bound marginally better to trimers purified with itself than with the other two bNAbs. Although the PGT121 epitope includes the N332 glycan, which is central to the 2G12 epitope, PGT121 bound most strongly to PGT151-purifed, and least strongly to PGT145-purified trimer; binding to 2G12-purified trimer was intermediate. PGT145 binding to the trimer purified with itself was stronger than to the other two forms. The greatest differences occurred with PGT151 and VRC34.01 binding, the binding to the PGT151-purified trimer strongly dominating, whereas binding to neither of the other two forms was detectable with VRC34.01. Another interface bNAb, 35O22, bound distinctly better to the PGT151-purified form than to the other two, although the difference was smaller than for PGT151 itself and VRC34.01. 3BC315, binding inter-protomerically and closer to the base than the interface bNAbs, gave another distinct ranking: 2G12-purified trimer highest, PGT145-purified somewhat lower, and PGT151-purified markedly lower. Finally, the autologous mNAb 16D bound best to PGT145-purified, more weakly to 2G12-purified, and only to a low level to the PGT151-purified form, indicating, again, that most PGT151-reactive forms in the trimer population expose the 289-glycan-hole epitope less well or present it in less antigenic shape than does the PGT145-purified trimer. The exposure of the 289 epitope may require a flexibility that the most PGT151-reactive forms lack [38].
The affinity-purified fractions of trimer correspond to the eluate in the chromatography. Because of the distinct antigenic effects of the differential affinity-purification specifically on the B41 trimer SOSIP.664 trimer, we also depleted it with bNAb-affinity columns, collecting the effluent. Thus, 2G12 and SEC purification of the B41 SOSIP.664 trimer followed by PGT145, PGT151, or mock depletion gave further evidence of antigenic heterogeneity. Depletion with PGT145 reduced PGT145 binding to the trimer and depletion with PGT151 increased it. Conversely and more markedly, depletion with PGT151 reduced PGT151 binding to the trimer and depletion with PGT145 increased it. 2G12 binding was also affected by the depletions: PGT151 depletion enhanced it and PGT145 depletion reduced it (Fig. 6B); these effects suggest allosteric connections between non-overlapping epitopes, as has been described [36, 57]. They are intriguing since no marked corresponding effects were recorded for differential purifications, but the latter may involve not only selection but also persistent induction of conformations.
The sensorgrams for Fab titrations against BG505 and B41 SOSIP.664 trimers are shown in Fig. 7. Kinetic constants kon and koff, the dissociation constant, KD, and the stoichiometry, Sm, were first obtained by Langmuir modeling (Table 1). Langmuir modeling gave passable fits, χ2 ranging from 0.13 to 0.65; the values of the kinetic constants were significant: T values (= mean/(s.e.m.)) were > 10, except for koff of the highly stable PGT145-Fab binding to PGT145-purified B41 SOSIP.664, which fell below the level of detectability, 10− 5 (s− 1). For the other combinations the T values were in the range 51–545 (SI Table 1).
Table 1
SOSIP.664 trimer | Fab | Affinity purification | konb (1/Ms) | koffb (1/s) | KDb (nM) | Smb |
BG505 | PGT145 | 2G12 (n = 3) | 1.8. 104 ± 3.7. 102 | 5.1. 10− 4 ± 2.0. 10− 6 | 29 ± 0.69 | 0.74 ± 6.8. 10− 3 |
PGT151 | 2G12 (n = 2) | 8.7. 104 ± 3.5. 102 | 1.8. 10− 4 ± 5.0. 10− 7 | 2.0 ± 1.3. 10− 2 | 1.9 ± 7.5. 10− 3 |
B41 | PGT145 | 2G12 (n = 3) | 2.0. 104 ± 8.8. 102 | 1.2. 10− 4 ± 1.5. 10− 5 | 6.2 ± 0.28 | 0.63 ± 1.6. 10− 2 |
PGT145 (n = 3) | 3.2. 104 ± 2.3. 103 | < 10− 5 | < 0.5 | 0.87 ± 4.8. 10− 2 |
PGT151 | 2G12 (n = 3) | 4.0. 104 ± 4.1. 103 | 3.1. 10− 3 ± 2.4. 10− 4 | 77 ± 5.0 | 0.45 ± 1.5. 10− 2 |
PGT145 (n = 2) | Minimal binding: < 10 RU |
a Tabulated values are means ± S.E.M of n independent replicates |
bkon and koff are the on- and off- rate constants, respectively; KD is equilibrium dissociation constant = koff/ kon. Sm is stoichiometric number. Sm denotes the stoichiometry, i.e., number of Fab molecules per trimer. |
The kon value for PGT151 binding to BG505 SOSIP.664 was higher and the koff value lower than for PGT145 binding: the net effect was a 14-fold higher intrinsic affinity of PGT151 than PGT145 Fab, in agreement with the higher neutralization potency of PGT151; the partial bivalency of IgG binding to virions is not expected to change IC50 ratios but to enhance potency weakly or moderately [3, 30, 58, 59].
The kinetic analysis of binding to B41 SOSIP.664 showed that the PGT145 purification gave a somewhat higher (60%) kon, a lower koff, falling below detectability, and lower KD, whose upper limit only could therefore be determined, than the corresponding values for 2G12-purified trimer. Those differences suggest an antigenic heterogeneity in the B41 SOSIP.664 population, which is marked after 2G12 purification and diminished by PGT145 purification.
The most central explanatory findings were the differences in stoichiometry between PGT145 and PGT151 binding to BG505 and B41 SOSIP.664 trimers. The stoichiometric Sm value of PGT145 Fab binding to BG505 trimer, 0.74, approached the ideal maximum of 1.0 for this epitope [30, 36, 60]; that of PGT151 Fab binding to BG505 trimer was close to its previously described maximum of 2.0 [30, 38, 61]. PGT145 bound with distinct stoichiometries to 2G12-purified, and PGT145-purified B41 trimer, Sm = 0.63 and 0.87, respectively, showing that the 2G12-purified trimer population contains species lacking detectable affinity for this antibody. Most striking was the low PGT151 stoichiometry for 2G12-purified B41 trimer: 4.2-fold lower than for BG505. PGT151 Fab bound so poorly to PGT145-purified B41 trimer that the data could not be modeled. Taken together, these findings strongly suggest subpopulations among the B41 SOSIP.664-trimer molecules with distinct antigenicities. The heterogeneity manifesting itself as reduced Sm values comprises binding and non-binding forms.
We investigated potential heterogeneities within the binding populations further by applying a heterogeneous-ligand model. Four criteria can be applied to evaluate the meaningfulness of the more complex model: first, a marked reduction in χ2; secondly, T values > 10, except if the modeling suggests the existence of a site from which the Fab dissociates below the level of detection; thirdly, that the modeled kinetic parameters for the two sites are distinct; and fourthly, that the component Sm values are not highly distinct, i.e., that a minority site is not negligible in population size. The outcome was that meaningful heterogeneity was discernible within the population of BG505- but not the B41-trimer molecules that showed any detectable binding: for BG505, the χ2 values were reduced 2.9- to 6.2-fold; T values were high except for a component of extremely slow dissociation of each Fab; kon1/kon2 was 5.8 (PGT145) or 5.0 (PGT151) and koff1/koff2 was > 200 (PGT145) or > 76 (PGT151); Sm1/Sm2 was 1.0 (PGT145) or 1.7 (PGT151). The complex modeling also increased the stoichiometry for PGT145 from Sm = 0.74 for Langmuir to a cumulative Sm1+Sm2 = 0.92 for heterogeneous-ligand modeling, a value close to the ideal Sm = 1.0 (Table 2 and SI Table 2).
Table 2
Heterogeneous-ligand modelinga
SOSIP.664 trimer | Fab | Affinity purification | kon1b (1/Ms) | koff1b (1/s) | kon2b (1/Ms) | koff2b (1/s) | KD1b (nM) | KD2b (nM) | Sm1b | Sm2b | Sm1 +Sm2b |
BG505 | PGT145 | 2G12 (n = 3) | 2.9. 104 ± 7.5. 102 | < 10− 5 | 5.0. 103 ± 7.0. 102 | 2.0. 10− 3 ± 9.1. 10− 5 | < 0.5 | 4.1. 102 ± 64 | 0.46 ± 1.6. 10− 2 | 0.46 ± 2.4. 10− 2 | 0.92 ± 3.8. 10− 2 |
PGT151 | 2G12 (n = 2) | 1.4. 105 ± 1.3. 103 | < 10− 5 | 2.8. 104 ± 1.0. 103 | 7.6. 10− 4 ± 2.5. 10− 5 | < 0.1 | 27 ± 1.9 | 1.2 ± 1.9. 10− 2 | 0.70 ± 1.1. 10− 2 | 1.9 ± 7.7. 10− 3 |
B41 | PGT145 | 2G12 (n = 3) | 5.2. 104 ± 1.4. 104 | 1.9. 10− 4 ± 3.6. 10− 5 | 1.6. 104 ± 1.5. 103 | 1.1. 10− 4 ± 1.7. 10− 5 | 3.7 ± 0.23 | 6.6 ± 0.80 | 4.4. 10− 2 ± 4.4. 10− 2 | 0.59 ± 3.6. 10− 2 | 0.63 ± 1.8. 10− 2 |
PGT145 (n = 3) | 3.5. 104 ± 1.4. 104 | < 10− 5 | 4.5. 104 ± 1.1. 104 | < 10− 5 | < 0.5 | < 0.5 | 0.13 ± 0.13 | 0.74 ± 0.16 | 0.87 ± 5.1. 10− 2 |
PGT151 | 2G12 (n = 2) | 3.1. 104 ± 1.9. 103 | 8.5. 10− 4 ± 2.5. 10− 4 | 2.5. 104 ± 5.5. 103 | 7.8. 10− 4 ± 9.8. 10− 5 | 28 ± 9.8 | 34 ± 12 | 0.45 ± 2.9. 10− 2 | 4.1. 10− 3 ± 4.0. 10− 3 | 0.45 ± 2.5. 10− 2 |
a Tabulated values are means ± S.E.M of n independent replicates |
b The constants and Sm are explained in Table 1 but here heterogeneous model allocates distinct values for two different sites, each indicated by the subscript 1 or 2. Sm1 +Sm2= the cumulative stoichiometry. |
In contrast, only the T-value criterion was met for the B41 trimer: the suggested kinetic constants and affinities were close or indistinguishable, and the minority populations small or negligible: specifically, for PGT151 kon1/kon2 was 1.2, koff1/koff2 1.1, KD1/KD2 0.82, and Sm1/Sm2 110 (SI Table 2): the sites did not differ tangibly in kinetics and the minority site was negligible.
In conclusion, two kinds of antigenic heterogeneity were detected: an inclusive kind, prevailing within the binding population of epitopes, was kinetically discernable for BG505; an exclusive kind, dividing the binding from the non-binding population, manifested itself as reduced stoichiometry for B41, very moderately for PGT145 but prominently for PGT151 binding to 2G12-purified trimer.
Structural modeling of PGT151 binding to the B41 SOSIP.664 trimer
Using available structural data, we examined whether the differential PGT151 binding could be caused by pronounced conformational flexibility or “breathing” in B41 SOSIP.664 that would limit access to the PGT151 epitope (Fig. 8). As noted, we previously observed substantial conformational heterogeneity in B41 SOSIP.664 by NS-EM (Fig. 5, [23–26, 52]), as well as the ability of a CD4bs bNAb, b12, to bind to its conformationally dependent epitope on that trimer [26]. Earlier hydrogen-deuterium-exchange mass spectrometry (HDX-MS) experiments demonstrated that although b12 can bind BG505 SOSIP.664 trimers, it does so very slowly, and fails to neutralize the BG505 PV, suggesting that the required more open conformation is rarely sampled by BG505 SOSIP.664 and not triggered by the antibody in a timeframe conducive to neutralization [62]. Hence, the B41 trimer is intrinsically more flexible than the BG505 trimer (Figs. 5 and 8). Double-electron-electron resonance (DEER) spectroscopy has revealed multiple conformations in the trimer base and inner domain of both BG505 and B41 SOSIP trimers, a flexibility that is uncoupled from that of the conformationally more fixed trimer apex [57]. Less conformational heterogeneity in the apex could explain the absence of marked PFs in PGT145 neutralization of either virus (Fig. 1). The same DEER experiments suggested a degree of conformational homogeneity in B41 SOSIP.v4.1, which includes stabilizing mutations to limit exposure of the V3 region and to shield non-NAb epitopes. But we did not study B41 SOSIP.v4.1 or other hyper-stabilized variants here, because we sought to mimic the neutralization-relevant conformational flexibility of Env on virions. We note, however, that the rate of breathing, accordingly, is quite plausibly limited by stabilizing mutations [53, 57]. The conformations of epitopes from apex to base are interconnected in a complex network of long- and short-range effects. For example, pre-binding of PGT145 to the BG505 SOSIP.664 trimer markedly reduces subsequent PGT151 binding, whereas no converse effect is detectable [63]. This non-reciprocal allosteric effect may not be identical for B41, but it suggests an intricate relationship between the two epitopes pertinent to the non-overlapping antigenicity maxima in the trimer population (Figs. 6 and 7; Table 1; SI Fig. 1).
The PGT151 epitope is unusually complex: PGT151 binding strictly depends on the native quaternary structure of proteolytically cleaved gp120-gp41 protomers. The paratope closely interacts with the interface between gp120 and gp41 of one protomer and glycans on both subunits of another protomer, inserting itself inter-protomerically [38] (Fig. 8A). Most residues implicated in the interaction are, however, identical for BG505 and B41, including the PNGSs at gp41 positions 611 and 637 (Fig. 8B). Both the fusion peptide (FP) and FP-proximate region (FPPR), reached by the inter-protomeric insertion of the long complementarity-determining region 3 of the heavy chain of the antibody, CDR H3, are highly similar in the two Envs (Fig. 8B). Modeling of PGT151 (based on a complex of its Fab with JR-FL SOSIP.664 [38]) onto a structure of closed-conformation B41 SOSIP.664 does not suggest any clashes or unfavorable interactions that would explain the different stoichiometries (Fig. 8C). Partially open, intermediate conformations of B41 SOSIP.664 have been solved by cryo-EM, including the aforementioned b12-bound state, one in complex with sCD4 and the interface-specific bNAb 8ANC195, and a more rearranged sCD4- and 17b-bound one [26, 28]. Because conformational changes occur both in the gp120 and gp41 subunits in these states relative to the closed conformation, we aligned the PGT151-bound structure in three different ways for three parts of the epitope: to the part of the primary gp120 (gp1201), to part of the primary gp41 (gp411), which includes fusion-peptide residues, and to the component of the adjacent gp41 (gp412) at its interface with the paratope (Fig. 8C). In the first two cases, the alignment to the primary subunit component of the epitope results in a clash of PGT151 with the adjacent gp120, which has rotated in response to b12 or sCD4 and 8ANC195 binding (Fig. 8C). This clash alone would prevent binding of PGT151 at this binding angle. Alignment to gp412 relieves the clash but allows paratope contacts with the trimer exclusively by the CDR H3, which is unlikely to yield tangible binding. Lastly, it was shown that, in the b12- or sCD4- and 17b-bound state, the fusion peptide of B41 SOSIP.664 is not solvent-accessible and is instead sequestered in a newly formed pocket in gp41 in all three protomers [26]. Sequestration of the fusion peptide may be a typical response to trimer opening, and this would remove a key component of the PGT151 epitope during such intervals of breathing.
In conclusion, the binding data suggest heterogeneity in how PGT145 and PGT151 recognize the total population of BG505 trimer, but that their binding is sufficient in strength and extent for the antibody to neutralize potently and effectively. Some subtler heterogeneity of PGT145 binding specifically to 2G12-purified B41 trimer was also apparent. Notably, however, the stark difference in stoichiometry of PGT151 binding to the BG505 and B41 trimers explains the larger PF of the B41 than of the BG505 neutralization. The low stoichiometry of PGT151 binding to B41 SOSIP.664 is explained by reduced paratope access through steric hindrance and the sequestration of a key component of the epitope by the partial opening of the trimer, which the B41 SOSIP.664 trimer is prone to.