Epitope-focused immunogens against the CD4-binding site of HIV-1 envelope protein induce neutralizing antibodies against auto- and heterologous viruses

Hua Wang (汪桦), Xiangjun Chen (陈相军), Dianhong Wang (王点红), Chen Yao (姚辰), Qian Wang (王茜), Jiayu Xie (谢家瑀), Xuanling Shi (史宣玲), Ye Xiang (向烨), Wanli Liu (刘万里), and Linqi Zhang (张林琦) From the Comprehensive AIDS Research Center, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, School of Life Sciences, and School of Medicine, MOE Key Laboratory of Protein Sciences, Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, School of Life Sciences, Institute for Immunology, and Beijing Advanced Innovation Center for Structural Biology, and Collaborative Innovation Center for Diagnosis and Treatment of Infectious Diseases, School of Medicine, Tsinghua University, Beijing 100084, China


Recent discoveries of broadly neutralizing antibodies (bnAbs)
in HIV-1-infected individuals have led to the identification of several major "vulnerable sites" on the HIV-1 envelope (Env) glycoprotein. These sites have provided precise targets for HIV-1 vaccine development, but identifying and utilizing many of these targets remain technically challenging. Using a yeast surface display-based approach, we sought to identify epitopefocused antigenic domains (EADs) containing one of the "vulnerable sites," the CD4-binding site (CD4bs), through screening and selection of a combinatorial antigen library of the HIV-1 envelope glycoprotein with the CD4bs bnAb VRC01. We isolated multiple EADs and found that their trimeric forms have biochemical and structural features that preferentially bind and activate B cells that express VRC01 in vitro. More importantly, these EADs could induce detectable levels of neutralizing antibodies against genetically related autologous and heterologous subtype B viruses in guinea pigs. Our results demonstrate that an epitope-focused approach involving a screen of a combinatorial antigen library is feasible. The EADs identified here represent a promising collection of possible targets in the rational design of HIV-1 vaccines and lay the foundation for harnessing the specific antigenicity of CD4bs for protective immunogenicity in vivo.
The global impact of a protective HIV-1 vaccine can hardly be overstated. Although both cell-mediated and humoral immune responses are believed to play critical roles in combating infection, it has been widely theorized that a successful vaccine should elicit a potent, sustainable neutralizing antibody response against a broad range of circulating HIV-1 strains (1,2). Recent scientific advances have led to the identification of a growing number of broadly neutralizing antibodies (bnAbs) 3 against HIV-1, which in turn have provided critical insights into the protective mechanisms of the human immune system and precise targets for immunogen design (3,4). Broadly speaking, these bnAbs recognized five major "vulnerable sites" on the HIV-1 envelope glycoprotein gp160 as follows: 1) the CD4binding site (CD4bs); 2) the glycan-associated V1V2; 3) V3 subdomains of gp120; 4) the membrane proximal external region (MPER) of gp41; and 5) the interface between gp120 and gp41 (5)(6)(7). Identifying antigens that can stimulate production of these bnAbs have become a major focus for HIV-1 vaccine researchers.
Traditionally, approaches to vaccine design have begun with the pathogen or its associated antigens aiming to manipulate their capacity to induce a protective antibody response in vivo. Although it has been remarkably successful against a large array of pathogens with little variability, it has become increasingly clear that a much more sophisticated approach is required to develop broadly protective vaccines against HIV-1 as it demonstrates continued sequence and structural evolution (8 -10). Common mechanisms of HIV-1 immune evasion, such as sequence diversity, conformation multiplicity, and the glycan shield, enable viral escape from immune recognition and antiretroviral agents, and they pose major challenges for HIV-1 vaccine development (11)(12)(13). Recent progress in the HIV-1 vaccine field has resulted in the development of antibody ontogeny-based HIV-1 subunits or trimeric Env immunogens with emphasis on triggering the specific antibody germ line ancestors (14 -19). These immunogens are able to trigger appreciable levels of autologous neutralizing antibodies in rabbits or to stimulate affinity maturation in the transgenic mouse models but failed to induce considerable heterologous bnAbs (20 -30).
The idea of reverse vaccinology, in which vaccine design works backward from the host bnAbs and the antigenic domains they recognized, has recently become an exciting new avenue for HIV-1 research (31). But how exactly this idea will be translated into vaccine design has yet to be clearly defined and tested. Several methodologies are currently being investigated such as the epitope-focused scaffolding and engineering as well as stabilizing the envelope glycoprotein in the pre-fusion closed conformation using the structure-based design (17)(18)(19)(20)(32)(33)(34). These approaches aim to characterize the antigenic regions closest to the bnAb epitope, thereby focusing the immune response on the key vulnerable sites while excluding regions that may elicit an ineffective or deleterious immune response. However, as many of the epitopes are conformational in nature, carving them out from the entire envelope glycoprotein while maintaining and stabilizing their native structures remains a technical challenge.
Our laboratory has recently developed a novel combinatorial antigen library approach for identifying antigenic domains containing the linear or conformational epitopes recognized by bnAbs (35)(36)(37). Using this system, we demonstrated here the successful isolation of the epitope-focused antigenic domain recognized by CD4bs bnAb VRC01 (EAD-VRC01). Structure and function analysis showed that EAD-VRC01 contains the large majority of CD4bs and is able to bind and attenuate the neutralizing activity of VRC01 in a pseudovirus-based assay. EAD-VRC01 is also able to trigger B cell activation through the formation of microclusters of the membrane-bound VRC01-BCRs and downstream signalings, the critical steps for B cells to differentiate into antibody-producing plasma cells. More importantly, immunization of guinea pigs with EAD-VRC01 induced detectable levels of neutralizing antibodies against genetically related autologous and heterologous subtype B viruses. Profiling of the neutralizing sera indicated the antibody specificity was largely directed to the CD4bs as well as V3 peptide on the envelope glycoprotein. We believe that the study reported here represents a novel and promising strategy for identifying EAD containing the desired bnAb epitopes. The EAD-VRC01 identified here could serve as a promising prototype for the future HIV-1 vaccines aiming to translate the specific antigenicity of CD4bs into the protective immunogenicity in vivo.

Selection of the epitope-focused antigenic domain by VRC01
We have previously shown that the combinatorial antigen library displayed on the yeast surface was able to express both linear and conformational epitopes and was successfully applied to the profiling studies of antibody response during HIV-1 and influenza infection (35)(36)(37). Taking the advantage of the scale and diversity of displayed antigens, we sought to identify the epitope-focused antigenic domain (EAD) for HIV-1 immunogen design and development. To this end, we first constructed a combinatorial HIV-1 envelope antigen library based on the full-length envelope gene of CNE11 derived from an HIV-1 subtype B chronically infected patient in China (38). The pseudoviruses bearing the CNE11 envelope were able to infect the CCR5 ϩ CD4 ϩ TZM-bl cell line and demonstrated sensitivity to the majority of well-characterized bnAbs targeting the key vulnerable sites on the envelope (38). The actual protocol for construction, validation, and selection of the combinatorial antigen library has been illustrated in greater detail elsewhere (35)(36)(37). Specifically, through multiple rounds of fluorescenceactivated cell sorting (FACS) and cloning, we were able to obtain 33 yeast clones that were highly reactive to VRC01, and a substantial number of clones to the control monoclonal antibodies (mAbs) such as PG9 (50 clones), 2F5 (33 clones), 4E10 (23 clones), and 10E8 (29 clones) ( Fig. 1 and Fig. S1). Sequence analysis of the coding sequences revealed that they were variable in length but invariably contained the epitope specificity defined by each of the corresponding antibodies ( Fig. 1B and Fig. S2) (5)(6)(7). This suggests that the selected EADs were correctly folded and displayed on the surface of the yeast. In particular, those selected by VRC01 (EAD-VRC01) largely spanned the V3, V4, and V5 regions of gp120, and those by PG9 (EAD-PG9) encompassed the V1V2 region of gp120, whereas those by 2F5 (EAD-2F5), 4E10 (EAD-4E10), and 10E8 (EAD-10E8) covered the MPER region of gp41 (Fig. 1, A-C). Of note, as PG9 is a glycan-dependent antibody, the EAD-PG9 must also contain an appropriate glycan moiety for specific recognition. The average length of the selected EAD-VRC01 is 282 aa, which is relatively longer than the EAD-PG9 (155 aa), EAD-2F5 (32 aa), EAD-4E10 (27 aa), and EAD-10E8 (29 aa). Furthermore, all of the 33 EAD-VRC01 sequences (Fig. S2), exemplified by the shortest (V61, 227 aa), the longest (V93, 372 aa), and the two intermediates (V51, 254 aa and V17, 260 aa), contained the critical residues as well as the conformational structure for VRC01 binding (Fig. 1, C and D). In fact, when superimposed onto the monomeric gp120, EAD-VRC01 significantly overlapped with the outer domain of gp120 and covered over 94% of the VRC01 footprint (Fig. 1D), in agreement with our initial design and expectation for EAD-VRC01.

Binding kinetics of monomeric and trimeric EAD-VRC01 to VRC01 class mAbs
To further characterize the binding properties of EAD-VRC01 to a spectrum of mAbs with well-defined specificities, we expressed the monomeric and trimeric forms of three representatives EAD-VRC01 with different lengths (V61, 227 aa; V51, 254 aa; and V17, 260 aa) in Sf9 cells using the Bac-to-Bac baculovirus expression system. A foldon trimerization domain of T4 fibritin was fused to the C terminus of EAD-VRC01 for the production of trimers. The monomeric and trimeric V61, V51, and V17 in the cell supernatant were concentrated and purified by nickel-affinity chromatography and size-exclusion chromatography before being subjected to the binding kinetic studies by surface plasmon resonance (SPR) on a Biacore T200 (GE Healthcare). As shown in Fig. 2, the monomeric and trimeric V61, V51, and V17 were able to bind to the antigenic fragment (Fab) of all members of the VRC01 class mAbs, suggesting that the VRC01 epitope on the EAD-VRC01 was well-maintained and accessible. Furthermore, despite their differences in length, V61, V51, and V17 demonstrated comparable binding affinities (K D ) to the Fab of VRC01 class mAbs, although the trimeric form demonstrated an overall average close to 10-fold stronger binding than the monomeric ones. Among all the VRC01 class mAbs tested, 12A12 bound the strongest whereas VRC-PG20 was the weakest (Fig. 2A). Furthermore, both the monomeric and trimeric forms of V61, V51, and V17 demonstrated similar levels of binding to the Fab of non-VRC01 class The consensus of amino acid sequences was numbered on the top of each graph according to HXB2 and highlighted by black background. The length and sequences of each consensus are shown on the far right wherever possible. C, sequence alignment of four EAD-VRC01 together with their consensus (con.) and corresponding sequences from CNE11 and HXB2. The EAD-VRC01 sequences represented the shortest (V61, 227 aa), the longest (V93, 372 aa), and the two intermediate ones (V51, 254 aa and V17, 260 aa). The contact residues with VRC01 were indicated with the symbol E for main-chain only contacts; ¤ for side-chain only contacts, and • for both main-and side-chain contacts. Those colored in yellow were located outside the selected EAD-VRC01 fragments, in purple among the consensus EAD-VRC01 fragment, and in dark green in the shortest EAD-VRC01 V61 fragment. Potential glycosylation sites with the signature sequence NX(T/S) were highlighted in cyan. Dashes represented the gaps introduced to the preserve alignment. The hypervariable regions V1-V5, loop D, CD4-binding loop, and ␤20/␤21 in gp120 were indicated. D, superimposing EAD-VRC01 sequences onto the monomeric gp120 (PDB code 3NGB) in ribbon and surface representations. The longest EAD-VRC01 was highlighted in light blue, the shortest in light green, and the consensus in red. The contact residues with VRC01 were colored as in C, and the contact surface of CD4 on the outer domain of gp120 was outlined by the white line. A, binding affinity between EAD-VRC01 and HIV-1 mAbs was measured by SPR and presented by K on (on-rate constant), K off (off-rate constant), and K D (equilibrium dissociation constant). The tested mAbs were CD4bs-targeted mAbs from multiple donors, including VRC01 class mAbs (VRC01, 12A12, VRC-CH31, 3BNC60, 3BNC117, NIH45-46, VRC-PG04, and VRC-PG20) and non-VRC01 class mAbs (CH103 and b12), V3 glycan-targeted mAbs (PGT121, PGT128, and PGT135), a V1V2-targeted mAb (PG9), CD4-IgG, and an irrelevant influenza hemagglutinin-targeted mAb (3C11) as the negative controls. The sensorgrams were fitted globally in 1:1 Langmuir model, and several K D values (denoted by asterisk symbols) were determined by the steady-state affinity model as the interactions were too fast to provide kinetic information. B, comparison of K D values within and between each mAb or combined group. The median value in each group was represented by the horizontal red line.

Epitope-focused HIV-1 vaccine design
mAbs CH103 and b12, although the overall average was about 10-fold less compared with that of VRC01 class mAbs (Fig. 2, A  and B). In addition, strong binding to the Fab of V3 glycan mAbs PGT121, PGT128, and PGT135 was also observed. In particular, the trimers demonstrated more than 2 log stronger binding than the monomers and also exceeded those of the VRC01 class mAbs (Fig. 2, A and B). This could be explained by the inclusion of the V3 region in the EAD-VRC01, which may become more exposed and accessible in the trimeric form. It is worthwhile to note that among the anti-V3 narrowly neutralizing mAbs tested (447-52D, 19b, and 39F), only 39F was able to bind to EAD-VRC01 (V61, V51, and V17) with significant preference for the trimers over the monomers (Fig. S3), suggesting that at least the N-terminal side of the V3 loop is exposed and accessible. These results clearly indicate that the three representatives EAD-VRC01 V61, V51, and V17 contain CD4bs and V3 glycan vulnerable sites identified on the full-length envelope glycoprotein of HIV-1, and the trimeric forms appear to adopt superior conformational structure than their monomeric counterparts for mAbs binding. Finally, none of the EAD-VRC01 V61, V51, and V17 had detectable bindings to the control mAbs PG9 and 3C11 or to soluble CD4-IgG ( Fig. 2A). This is consistent with earlier reports where the outer domain of gp120 also failed to engage with the CD4 molecule despite its strong interaction with VRC01 (33).

Attenuation of VRC01 Fab neutralization by monomeric and trimeric EAD-VRC01
Although the CD4bs on EAD-VRC01 V61, V51, and V17 are detectable on the solid surface by SPR, it is uncertain whether they are well-preserved in solution. This question is critical for their engagement with the membrane-bound VRC01 in cell cultures and their capacity in inducing the preferred antibody specificities in experimental animals. To this end, we mixed the Fab of VRC01 with the monomeric and trimeric V61, V51, and V17 in solution before applying the pseudovirus-based neutralization assay. We hypothesized that if the CD4bs was indeed well-preserved in solution, such mixing would result in a spontaneous absorption and therefore a dramatic depletion of the neutralizing activity of VRC01 Fab. In fact, this is what we have found. As shown in Fig. 3, both monomeric and trimeric V61, V51, and V17 were able to attenuate neutralizing activity of VRC01 Fab in a dose-dependent manner against the three representative HIV-1 subtypes or circulating recombinant forms (CRFs) from China (38). In general, the trimer was more effective than the monomer in depleting the neutralizing activity but such effect appears to be EAD-VRC01 as well as viral strain-dependent. For instance, rather small differences were identified between EAD-VRC01 V51 and V51tri compared with that between EAD-VRC01 V61 and V61tri as well as between EAD-VRC01 V17 and V17tri (Fig. 3B). In contrast, the viral strain CNE11 was relatively insensitive to VRC01 Fab depletion than CNE55 and CNE54 (Fig. 3). Nevertheless, these results indicate that both the monomeric and trimeric forms of EAD-VRC01 V61, V51, and V17 preserved the CD4bs in solution, although the latter appeared to adopt superior conformational structure for VRC01 binding than their monomeric counterparts, consis-tent with the binding kinetic analysis conducted on the solid surface using SPR (see above).

Structural analysis of EAD-VRC01 trimer and its complex with VRC01 Fab by electron microscopy
To further characterize the structural features of EAD-VRC01 and its binding with VRC01 Fab, we went further to analyze V51tri, one of the EAD-VRC01s, and its complex with VRC01 Fab using negative staining electron microscopy (EM). As shown in Fig. 4, trimeric features of the V51tri and V51tri-VRC01 Fab complex particles were clearly observed in the raw EM images (Fig. 4A). Reference-free 2D classification showed more details of the overall trimeric shapes of V51tri and V51tri-VRC01 Fab complex in both top and side views (Fig.  4B). In particular, the 2D classification of the V51tri-VRC01 Fab complex suggested that the majority of V51tri proteins were complexed with VRC01 Fab. Specifically, the reconstructed 3D EM maps of the V51tri and V51tri-VRC01 Fab complex showed the three monomeric subunits surrounding a cylindrical central core (Fig. 4C). Compared with V51tri, the V51tri-VRC01 Fab complex has additional densities protruding radially out from the lateral surface of each inner V51tri subunit. To further interpret the EM results, the V51tri structure model generated from the BG505 SOSIP.664 gp140 trimer (PDB code 4NCO) (39) and the VRC01 Fab crystal structure (PDB code 3NGB) (40) were fitted into the reconstructed 3D EM map (Fig. 4C). The fitted V51tri-VRC01 Fab complex model shows that the CD4bs was exposed to the outside and was spatially accessible to VRC01 Fab. Furthermore, the comparison between V51tri 3D EM map and its equivalent map generated from BG505 SOSIP.664 gp140 trimer (PDB code 4NCO) showed that the orientation and the overall arrangement of the monomers closely resemble those in the native gp140 trimers (Fig. 4D). Taken together, these results reinforced the notion that the trimeric form of V51, and most likely that of V61 and V17, contained the properly positioned and oriented epitope of VRC01 and are suitable for downstream biological and immunological studies.

EAD-VRC01 triggers surface BCR clustering and downstream signalings in B cells
Although EAD-VRC01 V61, V51, and V17 are able to bind to the soluble VRC01, it is uncertain whether they could activate B cells by triggering surface VRC01-BCR clustering and downstream signalings, the critical steps for B cells to differentiate into antibody-producing plasma cells. To address this question, we constructed the human Ramos B cells to express human IgM-BCRs with the variable regions of the heavy and light chain from VRC01 (Ramos-VRC01 B cells) and loaded onto a planar lipid bilayer containing the monomeric or trimeric V61, V51, and V17 for 10 min at 37°C. The formation and accumulation of BCR microclusters on the antigen contact area were examined using total internal reflection fluorescence microscopy (TIRFM) as reported previously (41)(42)(43). As shown in Fig. 5, A-C, V61, V51, and V17 were all able to trigger the formation of VRC01-BCR microclusters within B cell immunological synapses, but the trimeric forms were clearly more potent than the monomeric ones. Among the monomeric forms, V51 was the Epitope-focused HIV-1 vaccine design weakest, whereas V61 and V17 were comparable in terms of the contact area on the antigen-containing surface (Fig. 5B) and the total fluorescence intensity of the VRC01-BCR microclusters (Fig. 5C). The differences among the three trimeric forms were quite minimum but significantly more potent than their respective monomeric counterparts. In fact, they were equivalent if not more potent than the positive controls such as the monomeric full-length gp120 glycoprotein (HXB2, CNE54, and CNE55) and resurfaced stabilized core 3 of gp120 (RSC3) (Fig.  5, B and C). The negative control P01 (a EAD-PG9 from the same library) and RSC3⌬371I had either negligible or barely detectable activity in inducing the formation of VRC01-BCR microclusters.
The ability of V61, V51, and V17 to activate BCR downstream signalings in the Ramos-VRC01 B cells was further analyzed using intracellular immunofluorescence staining of the signaling molecules, phosphorylated spleen tyrosine kinase (pSyk), and downstream signaling molecules of phosphorylated tyrosine (pTyr) as reported previously (41)(42)(43). The phosphorylation of Syk and total tyrosine of BCR downstream signaling molecules is the direct outcome of BCR aggregation and therefore has been widely used as the downstream marker for B cell activation. Consistent with the effect on VRC01-BCR microclusters, the monomeric and trimeric V61, V51, and V17 were able to trigger the recruitment of pSyk and downstream signaling molecules of pTyr into the immunological synapses dem-onstrated by the accumulation and significant colocalization with VRC01-BCR in the TIRFM images ( Fig. 5D and Fig. S4A). Quantitative measurement of the contact area and total fluorescence intensity of pSyk and pTyr in the immunological synapse showed significantly more potency for the trimeric than the monomeric form, except for V17 and V17tri where the comparable effect on pSyk was found (Fig. 5, E and F, and Fig.  S4, B and C). Furthermore, the trimeric forms were at least equivalent to if not more potent than the positive controls such as the monomeric full-length gp120 glycoprotein (HXB2, CNE54, and CNE55). The negative control P01 showed insufficient ability to induce the phosphorylation of Syk and the total tyrosine of BCR downstream signaling molecules (Fig. 5, E and F, and Fig. S4, B and C).
We also examined the ability of V61, V51, and V17 to induce Ca 2ϩ mobilization in the Ramos-VRC01 B cells, an event following the formation of VRC01-BCR microclusters and the recruitment of pSyk and downstream signaling molecules of pTyr into the immunological synapses. Specifically, Ramos-VRC01-Gcamp5 B cells were generated to express the intracellular calcium indicator Gcamp5, stimulated with the monomeric or trimeric V61, V51, and V17, and then monitored calcium dynamics by flow cytometry. As shown in Fig. 5G, whereas the monomeric V61, V51, and V17 failed to induce Ca 2ϩ flux in the Ramos-VRC01-Gcamp5 B cells, the corresponding trimeric forms were quite capable, although V61 and

Epitope-focused HIV-1 vaccine design
V17 were more robust than V51. Taken together, these results demonstrate that EAD-VRC01 V61, V51, and V17 can activate B cells through the VRC01-BCR and are therefore appropriate for downstream immunological characterization.

EAD-VRC01 induces neutralizing antibody against autologous and heterologous HIV-1
As the trimeric forms of V61, V51, and V17 are superior to the monomeric ones in binding to the soluble VRC01 and acti-vating B cell signaling through the membrane-bound VRC01, we have chosen the trimeric forms to immunize guinea pigs. Two different immunization strategies were explored with the intention to test whether the immunization method itself has significant impact on immunogenicity in inducing neutralizing antibodies. One strategy involved a sequential immunization from the longest to the shortest, starting with the V17tri, followed by V51tri and then V61tri, and with a final boost of V17tri. The other involved immunization with the equal mix- The scale bars, 100 nm in the raw micrographs and represented 10 nm in the boxed particle images. B, reference-free class averages of the V51tri (left) and V51tri-VRC01 Fab complex (right). C, three-dimensional (3D) reconstructions of the V51tri (left) and V51tri-VRC01 Fab complex (right) were shown in the top (top) and side (bottom) views with the fitted models. The V51tri and V51tri-VRC01 Fab models were generated from the crystal structures of BG505 SOSIP.664 gp140 trimer (PDB code 4NCO) and gp120 -VRC01 Fab complex (PDB code 3NGB). The fitted V51 monomers were indicated by the cyan ribbon and the fitted VRC01 Fabs were indicated by the orange ribbon. The N and C termini of the fitted V51 monomers were marked with the blue and red balls, respectively. Positions of the foldon tags were indicated by the red triangles and the scale bars, 5 nm. D, comparison of the V51tri EM map (in gray) and its equivalent map (in magenta) generated from the crystal structure of BG505 SOSIP.664 gp140 (left, top view; right, side view). The position of the 3-fold axis was indicated by a black dotted line.

Epitope-focused HIV-1 vaccine design
ture of V61tri, V51tri, and V17tri throughout the course. For each strategy, six animals received four rounds of subcutaneous immunization at weeks 0, 4, 12, and 20, and blood samples were collected at weeks Ϫ2 for prebleed and at weeks 14 and 22 after the first immunization (Fig. S5, A and B). Serum-binding activity to a broad array of the envelope recombinant proteins was evaluated using standard ELISA (Fig. 6A). It is clear that both immunization strategies induced high levels of binding antibodies to the autologous trimeric V61, V51, and V17. Binding activity to the heterologous monomeric gp120 (HXB2, CM235, CNE54, and CNE55) and trimeric NL4-3 gp140 as well as RSC3 was also significant, although the actual binding titer was about 1 log less compared with the autologous trimers. The serum sample of Mix 05 in the mixture of the immunized group had equivalent high levels of binding to the autologous trimers but unusually low levels to the heterologous recombinant envelopes for reasons that are currently unknown (Fig. 6A). Serumneutralizing activity was tested against a panel of autologous and heterologous CNE pseudoviruses originally from China (38) and the global panel and other standard controls from the AIDS Reagent Program, National Institutes of Health (Fig. 6B) (44,45). It needs to be noted that all of the CNE pseudoviruses as well as those in the global panel were properly tier-phenotyped by Dr. David Montefiori and co-workers at Duke University Medical Center using widely accepted tier phenotype definition (44). Results from multiple independent experiments showed appreciable levels of neutralizing antibodies against the genetically related autologous and heterologous subtype B CNE pseudoviruses, although no detectable neutralizing activity was found against distantly related subtype B or other subtypes or CRFs (Fig. 6B). In general, neutralizing activity in the mixture group was slightly higher than that in the sequentially immunized group, but no significant difference in potency or breadth was noticed. Among all of the animals tested, Mix 02 demonstrated the strongest activity against the autologous (CNE11) as well as the heterologous (CNE4, CNE10, CNE14, and CNE57) strains. However, the prebleed sera did not have any detectable neutralization activity against the same virus panel (Fig. S5C). It is important to note phylogenetic analysis based on the amino acid sequences corresponding to the EAD-VRC01 region revealed that the sensitive viruses were clustered in a single branch, whereas those insensitive ones were rather distant (Fig.  S6). These results clearly indicate that the trimeric forms of V61, V51, and V17 were able to induce appreciable levels of neutralizing antibodies in vivo against the autologous and heterologous subtype B CNE strains that are genetically more related than other divergent subtypes or CRFs. The immunization strategy per se does not appear to have significant impact on the potency and breadth of the neutralizing antibodies induced in the experimental animals.

Profiling neutralizing sera by competitive ELISA and neutralization depletion
To dissect the target specificity of the neutralizing sera, we first assessed their ability to compete with various mAbs for binding to the immunogens through a previously published ELISA (46). As the immunogens were recognized by CD4bs as well as V3 glycan mAbs, representative mAbs from both classes were biotinylated and tested. Biotinylated PG9 and serum samples from the control group injected with PBS were used as negative controls. Competitions were considered as significant when the binding titer of biotinylated mAbs in the immune sera was reduced to less than half of that in the control sera. Specifically, all serum samples were diluted 1:100 and mixed with a serial dilution of the biotinylated mAbs before being applied to the ELISA plate coated with the immunogen V61tri. The competitive ability of the immune sera was defined as the residual binding percentage of the biotinylated mAbs compared with the controls. As shown in Fig. 7A, all serum samples showed the significant competition with CD4bs mAbs (VRC01, NIH45-46, 12A12, 3BNC60, and b12). In fact, except for the Mix 05, all samples reduced the binding activity of CD4bs mAbs for more than 97%, suggesting antibodies in the immune sera are largely directed to the epitope initially focused in the immunogen design. In addition, all serum samples also demonstrated the significant competition with V3 glycan mAbs, especially a more profound effect on PGT128 than PGT135 (Fig. 7A). The Mix 02 serum with the strongest neutralizing activity (see above) did not seem to have any special competitive edge over other samples. The Mix 05 serum, however, demonstrated the lowest competitive activity consistent with its genuinely low levels of binding as well as neutralizing activity (see above). As expected, no competition effect was detected against the control PG9 mAb. These findings are consistent with the antigenicity of the immunogens analyzed by SPR (Fig. 2) and suggested that some of their antigenic features were translated into the immunogenicity in the experimental animals. It needs to be noted, however, that the competition activity observed here could represent the direct overlaps with the epitope and/or steric interference with the tested mAbs. We need to be cautious when interpreting these results.

Epitope-focused HIV-1 vaccine design
To further profile the neutralizing sera, we also performed ELISA analysis against a 15-mer peptide library derived from the autologous CNE11 EAD-VRC01 envelope sequence. All serum samples were 1:1000 diluted before applied to the ELISA plate and their binding activity was measured as the absorbance (A) at 450 nm. Fig. S7 shows the binding profiles of all serum samples across the entire HIV-1 CNE11 EAD-VRC01 sequence. Each individual profile is presented together with the average values from the sequential and mixed immunization group. For comparative purposes, regions that have important structural and functional relevance are also indicated along the top of each profile. It is clear that binding profiles of all serum samples were confined to the region starting from the upstream of loop D all the way to the downstream of the V5 loop, matching the sequences of EAD-VRC01 immunogens ( Fig. 1C and  Fig. S7). However, the binding activity was not evenly distributed but was largely focused on six stretches of amino acid residues in the upstream of loop D, the V3 loop, the CD4-bind- Figure 6. Immunogenicity of EAD-VRC01 trimers in guinea pigs. A, ELISA binding titers of the immunized guinea pig sera against the autologous EAD-VRC01 trimer immunogens (V61tri, V51tri, V17tri, and clade B), HXB2 gp120 (clade B), NL4-3 SOSIP.664 gp140 trimer (clade B), RSC3 (clade B), CM235 gp120 (clade AE), CNE55 gp120 (clade CRF01_AE), and CNE54 gp120 (clade CRF08_BC). The sera tested here were collected 2 weeks (week ϭ 22) after the fourth inoculation. The ELISAs were performed on two independent experiments. B, neutralizing activity of the immune sera against a panel of pseudoviruses as well as 10 strains from the global panel in the TZM-bl assay. Moloney murine leukemia virus (MMLV) was used as a negative control. The sera tested here were collected 2 weeks (week ϭ 22) after the fourth inoculation. The neutralization assays were conducted on at least two independent occasions, and each was in duplicate. The color code was used to define neutralizing ID 50 titers, and those Ͻ20 were in black, 20 -100 in orange, 100 -300 in red, and Ͼ300 in purple.

Epitope-focused HIV-1 vaccine design
ing loop, the V4, the ␤20/␤21, and the V5 regions. In fact, these residues overlap somewhat with those in the CD4bs and V3 region recognized by VRC01 class and V3 glycan mAbs previously identified through the crystal structure analysis (40,47,48). Furthermore, no significant difference was found in the average profiles between the sequential and mixed immunization groups demonstrated by the almost overlapping orange dashed and green dotted lines in Fig. S7. Taken together, these results suggest that the targeted specificity of the neutralizing sera is confined to the intended regions of our initial design and overlaps with the CD4bs and V3 epitope recognized by VRC01 class and V3 glycan mAbs, respectively.
To further map the neutralizing sera, we conducted neutralization depletion experiments in the presence of recombinant RSC3, RSC3⌬371I, and the autologous V3 peptide. We hypothesized that if the sera contained the CD4bs-and V3-directed neutralizing antibodies, such mixing would result in a spontaneous absorption and therefore a reduction of the neutralizing activity of immune sera. To this end, the three relatively potent neutralizing sera (Mix 01, Mix 02, and Mix 04) were studied for their ability to neutralize the three sensitive pseudoviruses (CNE14, CNE4, and CNE57) in the presence or absence of RSC3, RSC3⌬371I, and the autologous V3 peptide in the TZM-bl assay. The fold of reduction in sera ID 50 was then compared and recorded. Throughout the entire analysis, RSC3 was on average more than twice as efficient as RSC3⌬371I in reducing the serum-neutralizing activity, although apparent differences existed for different samples and different viruses (Fig.  7B). For instance, the serum sample Mix 04 was the most sensitive and Mix 02 was the least sensitive to RSC3 depletion reflected by the actual fold of reduction in serum ID 50 . Among the tested viruses, however, CNE14 (tier 1B) was the most responsive to RSC3 depletion (RSC3/RSC3⌬371I ϭ 2.60/ 0.96 ϭ 2.71) followed by CNE57 (tier 2) (4.80/2.00 ϭ 2.40) and then by CNE4 (tier 2) (2.29/1.64 ϭ 1.39). Similarly, the V3 peptide was also able to reduce neutralizing activity of all serum samples, although Mix 04 was the most responsive and Mix 02 was the least responsive. These results clearly indicate that the immune sera from Mix 01, Mix 02, and Mix 04 contained CD4bs-and V3-directed neutralizing antibodies that can be absorbed and reduced by RSC3 and the V3 peptide, respectively. However, as sequence and structure differences exist between the EAD-VRC01 immunogens and RSC3 and the V3 peptide, the observed reduction in neutralizing activity and its extrapolated target specificity could only be interpreted with caution. More definitive answer for target specificity would have to rely on isolation and characterization of a sizable number of monoclonal antibodies from the immune animals.

Discussion
Identification of vulnerable sites on the HIV-1 envelope glycoprotein by bnAbs has provided precise targets for HIV-1 vac-cine design and development. Among many other approaches being explored, the epitope-focused strategy aims to target the vulnerable sites while eliminating the surrounding antigenic regions to minimize unnecessary or deleterious immune responses in vivo. However, as many of the epitopes are conformational in nature, carving them out from the entire envelope glycoprotein while maintaining and stabilizing their native structures remains a technical challenge. In this study, we demonstrated a successful isolation of the epitope-focused antigenic domain recognized by VRC01 (EAD-VRC01) through FACS sorting of the combinatorial antigen library displayed on the surface of yeast. Sequences and structural analysis showed that EAD-VRC01 had a significant structural overlap with the outer domain of gp120, covering the upstream of loop D all the way to V3, V4, and the downstream of V5 in sequence (33,40). Biochemical characterizations of the three representatives EAD-VRC01 (V61, V51, and V17) revealed their strong binding to VRC01 class as well as to V3 glycan mAbs and therefore contained both CD4bs and V3 glycan-vulnerable sites. V61, V51, and V17 were also able to attenuate the neutralizing activity of VRC01 Fab in a pseudovirus-based assay and to trigger B cell activation through the formation of microclusters of the membrane-bound VRC01-BCRs and downstream signalings, the critical steps for B cells to differentiate into antibody-producing plasma cells. In all these settings, the trimeric forms of V61, V51, and V17 were invariably more potent than the monomeric ones. In particular, the sequential or mixed immunization of guinea pigs with the trimers was able to induce appreciable levels of neutralizing antibodies against genetically related autologous and heterologous subtype B CNE strains from China, although no detectable neutralizing activity was found against distantly related subtype B or other subtypes or CRFs. Systematic profiling of the immune sera by competitive ELISA, autologous peptide library, and neutralization depletion identified appreciable levels of CD4bs-and V3-directed neutralizing antibodies that can be absorbed by RSC3 and the V3 peptide, respectively. Taken together, these results demonstrated that the rational design of epitope-focused vaccine is not only theoretically viable but also technically feasible. The EADs identified here represent a promising collection of possible targets in the rational design of HIV-1 vaccines and lay the foundation for harnessing the specific antigenicity of CD4bs for protective immunogenicity in vivo.
Designing immunogens targeting the various sites of vulnerability and dissecting the antibody response in animal models will provide a critical guidance for the selection and optimization of vaccine candidates (20,(32)(33)(34). Studies on EAD-VRC01 here represent one of the first examples that a portion instead of a complete envelope immunogen can induce appreciable levels of neutralizing antibody in guinea pigs. We speculate that such immunogenicity could result from the structural and biochem- Figure 7. Profiling of the guinea pig immune sera by competitive ELISA and neutralization depletion. A, residual binding percentage (%) of biotinylated HIV-1 mAbs to V61tri measured by competitive ELISA. The binding titers of biotinylated mAbs with known epitopes in the presence of the immune guinea pig sera were evaluated and numerically presented. The sera from the NC group and PG9-biotin mAb were used as negative controls. B, reduction of serumneutralizing activity in the presence of RSC3, RSC3⌬371I, and the autologous V3 peptide. The three relatively potent neutralizing sera (Mix 01, Mix 02, and Mix 04) were studied for their ability to neutralize the three relatively sensitive pseudoviruses (CNE14, CNE4, and CNE57) in the presence of RSC3, RSC3⌬371I, and the autologous V3 peptide in the TZM-bl assay. The actual neutralization curves (top) and fold reduction in neutralization ID 50 (bottom) were showed. The Ͼ symbol indicates the complete depletion of neutralizing activity.

Epitope-focused HIV-1 vaccine design
ical properties of EAD-VRC01 where the CD4bs might have become more stabilized and exposed through the trimeric formation thereby to improve the accessibility and affinity for VRC01 class and the related antibodies. In other words, our selection criteria for EAD-VRC01, such as the structural mimicry and stability, the antigenic similarity and multiplicity, as well as the ability to trigger B cell activation in vitro, have served as a meaningful gatekeeper for the immunogen selection and animal immunization. In fact, these criteria are somewhat similar to those for evaluating the eOD-GT8 immunogen except it is precisely engineered to trigger the germ line version of VRC01 in the transgenic animals (14). It is certain that future optimization would be required for EAD-VRC01 to obtain additional desired features such as improved binding to the intermediate or germ line ancestors of the mature VRC01 to trigger and enhance the development process of VRC01-like antibodies in vivo. EAD-VRC01 without the V3 region would also be a viable alternative to see whether deletion of the V3 region could further drive the immune response toward the CD4bs. In addition, detailed analyses of the neutralizing antibody responses in the immunized guinea pigs, including isolating monoclonal antibodies and next-generation sequencing, will certainly provide a more concrete reference for antibody specificity and the subsequent rounds of optimization. Furthermore, as the different animal models contain the distinct antibody repertoire, it is imperative to test the immunogenicity of EAD-VRC01 in those closest to humans to provide the most relevant and precise reference for downstream human use.
It needs to be stressed that our combinatorial antigen library approach offers several advantages pertaining to the epitopefocused vaccine strategies. The sheer size and diversity of the library undoubtedly improve the likelihood of success in identifying those antigenic domains containing the vulnerable sites that may otherwise be impossible by manual design on a one by one basis. More importantly, the vulnerable sites in the antigenic domains maintain their native conformation and closely resemble those in the full-length envelope expressed in mammalian cells with the post-translational modifications such as glycosylation and disulfide isomerization (49). This is largely due to the fact that yeast is a eukaryotic species, and the selection of antigenic domains is conducted in solution by FACS, which likely reflects the intrinsic interaction between antigen and antibody within the host. Our ability to identify EAD-VRC01 and EAD-PG09 that contain the complex conformational epitopes as well as the sugar moiety is a clear testimony for these advantages. Furthermore, random mutagenesis could also be applied to the selected EAD to generate the novel library for further selection of EADs with the additional desired attributes such as improved binding to the mature mAb or to its intermediate and germ line ancestors. As all of these procedures are conducted through the yeast-surface display system, and there is no need to produce each individual EAD in either prokaryotic or eukaryotic system and hence to greatly improve the chance and efficiency in identifying EADs with the desired properties. The versatility with high throughput capacity has therefore made the combinatorial library approach potentially the method of choice and serve as the critical starting point toward the successful epitope-focused vaccine strategies.
In summary, our results demonstrated that the rational design of epitope-focused vaccine is not only theoretically viable but also technically feasible. The EADs identified here represent a promising collection of possible targets in the rational design of HIV-1 vaccines and lay the foundation for harnessing the specific antigenicity of CD4bs for protective immunogenicity in vivo. The yeast surface display system will offer a powerful tool for the epitope-focused vaccine strategies targeting other major vulnerable sites on the envelope glycoprotein of HIV-1.

Ethics statement
The guinea pig experimental protocol for immunization and blood collection was approved by the Institutional Animal Care and Use Committee (IACUC) at Tsinghua University. The actual experiments were designed and performed in strict accordance with the government rules and regulations for animal welfare and the principles of the 3Rs (Replacement, Refinement, and Reduction).

Construction of the HIV-1 gp160 combinatorial antigen library
The combinatorial antigen fragment library of HIV-1 gp160 was constructed on the surface of yeast Saccharomyces cerevisiae as described previously (35,36). In brief, the full-length envelope gene of CNE11 derived from a subtype B HIV-1 chronically infected patient in China (38) was amplified by PCR, purified by electrophoresis (QIAquick DNA purification kit, Qiagen), and then digested into about 50-bp fragments by DNase I. Next, the digested fragments were reassembled into larger fragments about 100 -800 bp in length through controlled rounds of PCR amplification, which were then A-tailed (DNA A-tailing kit, Takara) and ligated to the modified pCTCON2-T vector. The ligation products were transformed into the Escherichia coli competent cells, amplified, extracted, and then further transformed into the competent S. cerevisiae EBY100 by electroporation. Transformed yeast cells were partially spread on SDCAA Amp plates and incubated overnight at 30°C to estimate the number and sequences of recombinant colonies for quality control purpose. Conditions and procedures for yeast growing and induction of the surface antigen expression in solution have been previously described (50).

Identification of the epitope-focused antigenic domain by FACS and sequence analysis
Immunofluorescence staining and sorting of the yeast clones expressing the appropriate antigenic domains were performed as described previously (35,36). The EAD was defined as the antigenic regions closest to the mAb epitope previously identified by crystallography while excluding the surrounding irrelevant regions. Specifically, the induced yeast cells (10 6 -10 7 ) of the combinatorial antigen fragment library were stained with VRC01 or other HIV-1 mAbs (10 g/ml), and the positive clones were sorted by FACS using Aria II (BD Biosciences). The sorted yeast clone mixture was spread onto the SDCAA Amp plates before single colonies were picked, grown, and further confirmed for binding with mAbs. To ensure accurate determination and representation of the EAD, multiple yeast clones

Epitope-focused HIV-1 vaccine design
were analyzed in parallel for their reactivity with the selected HIV-1 mAb. Those with detectable and desired reactivity to the mAb were further grown in SDCAA before plasmids were extracted (yeast plasmid kit, Omega Bio-tek) for sequencing and sequence analysis (Sequencher 5.0, Gene Codes Corp.). The EAD selected by VRC01 (EAD-VRC01) was also modeled onto the HIV-1 gp120 monomer to determine its spatial location using the program PyMOL (DeLano Scientific).

Production of recombinant EAD-VRC01, Fab, and their complex
The EAD-VRC01 was further expressed by Bac-to-Bac Baculovirus Expression System (Invitrogen) in the monomeric and trimeric forms, according to the manufacturer's instructions. All EAD-VRC01 sequences were fused with a gp67 signal peptide at the N terminus and a His 6 tag at the C terminus. For expression of the trimeric EAD-VRC01, a foldon domain (GYIPEAPRDGQAYVRKDGEWVLLSTFL) of bacteriophage T4 fibritin was added upstream of the His 6 tag at the C terminus. Recombinant EAD-VRC01 expressed in the supernatants of Sf9 cells were concentrated and purified by the nickel-nitrilotriacetic acid metal-affinity chromatography (Qiagen) and size-exclusion chromatography using a Superdex 200 column (GE Healthcare).
Antigen-binding fragments (Fab) of mAbs were produced by the protease digestion of full-length antibodies using human rhinovirus 3C or Lys-C (Sigma) protease. The complex of the trimeric EAD-VRC01 and VRC01 Fab was generated by the on-column cleavage strategy. Additional details regarding the generation of Fab and the complex can be found in the supporting Methods.

Surface plasmon resonance
The binding kinetics of the EAD-VRC01 with various mAbs were determined by SPR on a Biacore T200 (GE Healthcare). All experiments were performed at 25°C with HBS-EP buffer (10 mM Hepes, 150 mM NaCl, 3 mM EDTA, and 0.005% Surfactant P20, pH 7.4). The monomeric or trimeric EAD-VRC01 was immobilized onto CM5 sensor chips at 200 -300 response units by an amine coupling kit (GE Healthcare). The Fab of HIV-1 mAbs or a control mAb 3C11 was injected through flow cells at the concentration ranging from 1 to 500 nM followed by a series of 2-fold increase in concentration. The flow rate was 30 l/min for 3 min for association and for 5-30 min for dissociation. Regenerations were conducted with 10 mM glycine HCl, pH 2.0, at a flow rate of 50 l/min for 1 min. Biacore T200 Evaluation Software (GE Healthcare) was performed to subtract blank references. The sensorgrams were fitted globally in a 1:1 Langmuir model, and several K D values (denoted by asterisk symbols in Fig. 2A) were determined by the steady-state affinity model as the interactions were too fast to provide kinetic information.

Neutralization and neutralization depletion experiments
Neutralizing activities of guinea pig immune sera and mAbs were evaluated using a single round infection of TZM-bl target cells by pseudoviruses as previously described (51). The pseudoviruses were selected in such a way that included the autologous and heterologous CNE viruses initially isolated from the Chinese patients (38) and the subtype B and C and the global panel from the AIDS Reagent Program, National Institutes of Health (44,45). It needs to be noted that all of the CNE viruses were properly tier-phenotyped by Dr. David Montefiori and co-workers at Duke University Medical Center using the same protocol for widely accepted standard strains, including those in the global panel (44). Moloney murine leukemia virus envelope was used as a negative control. Neutralization depletion assay was conducted for VRC01 Fab in the presence of the EAD-VRC01 to verify whether the EAD-VRC01 was able to bind and attenuate the neutralizing activity. Neutralizing depletion assay was also performed for the immune sera of guinea pigs in the presence of RSC3, RSC3⌬371I, and autologous V3 peptide to study whether neutralizing activity could be competed out by these specific domains of gp120. Specifically, the competitor proteins or peptides (30 g/ml) were incubated with a serial dilution of mAbs or sera at 37°C for 1 h before adding to 100 TCID 50 of pseudoviruses. V3 peptide (KSIHL-GQGRAWYTTG) derived from CNE11 was designed according to the format (catalog no. 1830) listed in the AIDS Reagent Program, National Institutes of Health, and synthesized by ChinaPeptides (Shanghai, China). The depletion effect of the competitor proteins or peptides was assessed by the changes in halfmaximal inhibitory concentrations (IC 50 ) or half-maximal inhibitory dilutions (ID 50 ) compared with the mock controls.

Electron microscopy and image reconstruction
Three microliters of the purified V51tri (0.01 mg/ml) or V51tri-VRC01 Fab complex (0.01 mg/ml) were applied to a pre-glow discharged carbon coated EM grid and then stained with 1.5% uranyl acetate. Micrographs of the V51tri (101 micrographs) and the V51tri-VRC01 Fab complex (115 micrographs) were collected on a 120-kV FEI Tecnai Spirit microscope using a charge-coupled device camera at a nominal magnification of 49,000 (the calibrated pixel size is 2.27 Å). The structure model of the V51tri was generated using the crystal structures of the BG505 SOSIP.664 gp140 trimer (PDB code 4NCO). Because the V3 loop was truncated in the gp120-VRC01 Fab crystal structure (PDB code 3NGB), the V51-VRC01 Fab complex model was generated as a chimera with the V51 fragment from the crystal structure of the BG505 SOSIP.664 gp140 trimer (PDB code 4NCO) and the VRC01 Fab from the crystal structure of the gp120 -VRC01 Fab complex (PDB code 3NGB). Model Fitting was performed by manual adjustment and the "fit in map" function in the UCSF Chimera software (52). Map comparison was accomplished with EMAN2 e2pdb2mrc.py package (53) and the fit in map function in Chimera. Additional details of data processing and model fitting can be found in the supporting Methods.

Triggering VRC01-BCR signaling by EAD-VRC01
The Ramos B cell line was first constructed to express VRC01-BCR on the cell membrane by lentiviral infection, according to the supporting Methods. The EAD-VRC01 triggering VRC01-BCR signaling was recorded by an Olympus IX-81 TIRFM, following our well-established protocols (41,42,54,55). Additional details related to the entire procedure for Epitope-focused HIV-1 vaccine design molecule staining and TIRFM imaging could be found in the supporting Methods.

Calcium mobilization analysis
Ca 2ϩ mobilization in Ramos-VRC01 B cells was measured using intracellular calcium indicator Gcamp5 after EAD-VRC01 stimulation. Briefly, Ramos-VRC01-Gcamp5 cells were washed with Hanks' balanced salt solution (HBSS) once and then incubated at 37°C for 3 min in HBSS buffer (with 1.26 mM Ca 2ϩ ). After the basal level of Ca 2ϩ concentration was monitored for 30 s, cells were directly stimulated with 100 nM EAD-VRC01 monomer, trimer, or control antigens. The fluorescence intensity of Gcamp5 was measured on BD Accuri C6 and analyzed in FlowJo software (Tree Star). The Gcamp5 mFI was normalized against the average basal level.

Animal immunization and sample collection
Approximately 8 -10-week-old Dunkin Hartley female guinea pigs (Beijing Vital River Laboratory) were immunized subcutaneously with the trimeric EAD-VRC01 either sequentially or together as illustrated in Fig. S5. Specifically, the sequential immunization group (n ϭ 6) was immunized sequentially with 300 g (in 300 l of PBS) each of EAD-VRC01 V17tri, V51tri, and then V61tri. The mix immunization group (n ϭ 6) was immunized with a 300-g (in 300 l PBS) mixture of the three EAD trimers and each contributed 100 g. The first immunization was conducted at week 0 with Freund's complete adjuvant (300 l), and the following boosts were at weeks 4, 12, and 20 with Freund's incomplete adjuvant (300 l) (Sigma). The negative control group (n ϭ 3) was injected with PBS and then adjuvant. Blood samples were collected from the anesthetized animals at week Ϫ2 for prebleed and at weeks 14 and 22 after the first immunization, and no sign of distress was noticed during the experiment.

ELISA, competitive ELISA, and sera profiling
ELISA was used to define the binding activities of mAbs or sera against various antigens as described previously (35,56). The end-point titers of the sera were determined at the dilution where A 450 was 3-fold higher than the negative control animal sera. To study the antibody specificity of immunized animals, competitive ELISA was conducted against various biotinylated mAbs with known epitopes, following the previous protocol (46). The sera profiling of the immunized guinea pigs was performed by ELISA against the autologous 15-mer peptides library, overlapping by 11 residues, derived from the CNE11 envelope sequence (clade B). These peptides (5 g/ml) were coated onto 96-well plates followed by incubation with 1:1000 dilution of the sera as described previously (57). Additional details for the profiling experiment could be found in the supporting Methods.

Sequence and statistical analysis
The sequences corresponding to EAD-VRC01 were downloaded from the GenBank TM and aligned using the ClustalW program. Phylogenetic analysis was performed by the neighbor-joining method, and the reliability of the branching orders was detected by bootstrap analysis of 1000 replicates. In the pseudovirus-based neutralization and neutralization depletion experiments, IC 50 and ID 50 values were calculated by the doseresponse inhibition model in GraphPad Prism 5.0 (GraphPad Software). In the BCR clustering and downstream signaling molecule imaging analysis, the contact area and total fluorescent intensity among the study groups were analyzed using a two-tailed unpaired t test (*** means p Ͻ 0.001). In the ELISA and competitive ELISA, data were also analyzed in GraphPad Prism 5.0.
Author contributions-H. W., Y. X., W. L., and L. Z. designed the research. H. W., X. C., D. W., C. Y., Q. W., J. X., and X. S. performed the experimental work. Specifically, H. W., C. Y., Q. W., J. X., and X. S. performed the EAD selection and expression, SPR, neutralization assay, ELISA, and the animal experiments. X. C. performed the B cell trigger assay. D. W. conducted the EM analysis. H. W., X. C., D. W., Y. X., W. L., and L. Z. analyzed the data and wrote the manuscript.