Identification of VHHs binding to different epitopes of SARS-CoV-2 S1 protein RBD
Recently, we reported the identification of llama VHHs that bind to the SARS-CoV-2 S1 protein RBD 19. Briefly, we used two llama VHH libraries (one naïve library and another humanized synthetic library derived from the naïve library) to screen for VHHs that bind to the SARS-CoV-2 S1 protein in-vitro 19. We identified a total of 89 S1 protein binders, 64 from the naïve and 25 from the humanized libraries, respectively. Out of the S1 protein binders, 19 VHHs blocked the interaction between SARS-CoV-2 S1 RBD and ACE2 receptor, with 12 S/ACE2 blockers identified from the naïve library and 7 identified from the humanized library (data not shown). Furthermore, we observed that the pairwise addition of some of the VHHs caused synergistic effects on SARS-CoV-2 S/ACE2 blocking 19. We hypothesized that this synergistic effect is caused by binding of the VHHs to different epitopes within the S1 RBD. To test this, we performed epitope binning assays by biolayer interferometry (Figure 1a-c) or ELISA (Figure 1d) on a selected number of candidates.
In the initial epitope binning assay (Figure 1a-c), we used an S1 RBD sensor to capture 2A-Fc, 1B-Fc, or 3F-Fc separately, followed by the incubation with our lead candidates 1B-Fc, 3F-Fc or 2A-Fc. The VHHs were fused to human IgG1 Fc domains to render the Fc effector functions against SARS-CoV-2 19. This analysis showed that with the 2A-Fc-loaded probe, the addition of 3F-Fc further increased the signal compared to the 2A-Fc control, while the addition of 1B-Fc decreased the signal compared to the control (Figure 1a). This indicates that 3F-Fc does not compete with the 2A-Fc site and it is likely that they bind to different S1 RBD epitopes. In contrast, 1B-Fc competed with 2A-Fc, indicating that they compete for binding to the same S1 RBD epitope (Figure 1a). Similarly, with the 1B-Fc-loaded probe (Figure 1b), 3F-Fc increased the signal compared to the 1B-Fc control. This shows that 3F-Fc does not compete with 1B-Fc. Interestingly, 2A-Fc also increased its signal compared to 1B-Fc control. This suggests that although having a common binding region, 2A binds to a wider epitope than 1B (Figure 1b). With 3F-Fc-loaded probe, both 2A-Fc and 1B-Fc showed an increase of the signal compared to the 3F-Fc control. This further shows that 3F-Fc does not compete with either 1B-Fc or 2A-Fc, and likely bind to a different epitope (Figure 1c). These results confirm our hypothesis and show that S/ACE2 blocking VHHs bind to at least two separate unique epitopes within the S1 RBD.
Next, we performed an ELISA-based epitope binning assay to assess five additional VHHs (1C, 1F, 3A, 4F, and G4) unfused to Fc, but previously assessed to block the SARS-CoV- 2 S/ACE2 interaction 19. The assessment of more VHHs would allow us to categorize several of our other VHHs into binding groups, which will aid in multi-specific antibody design and construction. In this ELISA, wells were coated with SARS-CoV-2 S1 and incubated with bi- specific VHH-Fc 1B-2A (based on previous data, 1B and 2A likely bind the same epitope) or monoclonal VHH-Fc 3F-Fc (based on previous data, this binds a different epitope than 1B or 2A) premixed with the VHH candidates. The resulting relative fluorescence signals obtained for each sample were calculated to reflect the percent difference from 1C, 1F, 3A, 4F, G4, and controls (3F-Fc and 1B-2A-Fc) signals, when the VHHs are combined with 1B-2A-Fc or 3F-Fc (Figure 1d). The results show that VHH-Fcs 1C, 1F, 4F, as well as the 1B-2A-Fc control have almost 100% difference from 1B-2A-Fc, which highly suggest that they compete for the same epitope (Highlighted in Red). However, G4 (Highlighted in Light Red) may partially compete with 1B-2A-Fc, whereas 3A does not likely compete for the same epitope (Highlighted in Green). Additionally, these results show that 3A and the 1B-2A-Fc control may compete with 3F (Figure 1d), while other VHHs, including the 1B-2A-Fc control resulted in a lower percent differences We also performed additional epitope binning assays using biolayer interferometry to assess the competition of the VHH-Fcs 1C, G4, and 3A to bind to S1 RBD. The VHH-Fcs 1F and 4F poorly bound to the biolayer interferometry probes used for this assay and were excluded from analysis. This approach confirmed the results that we obtained by ELISA and showed that 1C and G4 likely belong to Group 1, and 3A belongs to Group 2 in terms of the binding competition (Supplementary Fig. 1). Interestingly, G4-Fc shows competition with either 1B-Fc and 2A-Fc when it is loaded onto the probe first (Data not shown). In contrast, reversal the of loading further increased its signal compared to both 1B-Fc control and 2A-Fc control, suggesting that the epitope for G4 is wider than that of both 1B and 2A (Supplementary Fig. 1). Taken together, we could categorize 8 VHH blockers of S/ACE2 interaction into 2 major groups based on their epitope binding; Group 1 consist of 6 VHHs, whereas Group 2 consist of 2 VHHs (Figure 1c).
Elucidation of epitopes on S1 RBD that bind to VHH-Fcs
In an effort to elucidate the structural basis of the newly discovered epitope binding groups, we computationally generated structural models for 1B, 3F, and 2A VHHs and docked them with SARS-CoV-2 S1 RBD structure exported from PDB 6M0J using Schrodinger BioLuminate software. For context, Figure 2a shows the SARS-CoV-2 S1 protein with the ACE2 binding residues red font. This approach generated an array of poses of each S1 RBD/VHH complex structure, which allowed us to further analyze the interfaces of those poses with a good PIPER cluster size and led us to identify five regions in the RBD which may interact with VHH 1B, 2A, and 3F, respectively (Figure 2a and b). Next, we generated 5 different S1 RBD deletion mutants to validate the computationally mapped epitopes in-vitro to select the best docking model for molecular analysis. Interestingly, these S1 RBD deletion regions have been shown to block the S1 RBD/ACE2 interaction in recently published literature 10-13,20 (Table 1). We tested wild-type and all the S1 deletion mutants for their ability to bind to a tri-specific VHH-Fc to check whether the proteins are folded and expressed on the cells. The results show that they are indeed expressed and folded as they all bind to the tri-specific VHH-Fc, although the level of expression and/or folding might be different across the mutants based on the strength of the binding signals. The wild-type S1 RBD and the deletion 2 (del2) shows stronger binding, whereas the deletions 1 (del1), 3 (del3), 4 (del4) and 5 (del5) show weaker binding to the tri- specific VHH-Fc (Supplementary Fig. 2). Then we assessed the binding profiles of the S1 RBD wild-type and the deletion mutants with selected VHH-Fcs from Group 1 and Group 2, as well as ACE2 (Figure 2c and d). The binding of VHH-Fc candidates from both Group 1 and Group 2, as well as ACE2 to S1 RBD are affected following the removal of del1. It is possible that this results is due to a conformational change or decrease of S1 protein expression following its deletion because based on crystal structure of RBD/ACE2 complex (PDB 6M0J), the deleted domain is not part of the S1 RBD/ACE2 interface. The del2 which is adjacent to a computationally- predicted epitope domain in region 1 does not prevent the binding of both Group 1 and Group 2 VHH-Fcs to S1 RBD.. In addition, it does not prevent the binding of ACE2 to S1 RBD.. The removal of del3, 4, and 5 all decrease the binding of both Group 1 and Group 2 VHH-Fcs to S1 RBD. However, these regions are more critical for Group 1 than the Group 2 for their binding. In addition, these regions are critical for ACE2 to bind to S1 RBD. Taken together, the binding epitopes for Group 1 is more associated with del3, 4 and 5 regions which are located at the interface of S1 RBD/ACE2, while at least part of the epitopes for Group 2 are shifted farther away from the S1 RBD/ACE2 interface relative to the epitopes for Group 1 VHHs (Figure 2c and d). Based on the binding and epitope binning data, we constructed 3D docking models that predicted the interactions between SARS-CoV-2 S1 RBD, ACE2 and lead VHH-Fcs (Figure 2e). These models show that predicted binding epitopes for Group 1 VHHs 1B and 2A are located at the S1 RBD/ACE2 interface. In contrast, the epitope for Group 2 VHH 3F is located away from the S1 RBD/ACE2 interface (Figure 2e). . Interestingly, there are binding variations seen within Group 1. The binding of 2A to del1, del3, del4 and del5 have decreased more than that of 1B. This shows that epitopes for 2A and 1B are not the same even though they compete with each other and were initially characterized to be within the same binding Group 1 (Figure 2c and d). Taken together, our analysis confirms that there are two major binding groups (Group 1 and Group 2) and we show the likely binding regions on the SARS-CoV-2 S1 protein for each VHH.
Tri-specific VHH-Fcs show potent S1 RBD binding and S/ACE2 blocking activity
Next, we tested whether the combination of individual VHHs binding to different S1 RBD epitopes into bi-specific antibody molecules would yield synergistic effects in SARS-CoV-2 binding and S/ACE2 blocking. As expected, the resulting bi-specific VHH-Fc 1B-3F showed superior binding to S1 RBD and S/ACE2 blocking compared to individual component VHH-Fcs 19. Since SARS-CoV-2 S proteins formed trimers, we started to study whether tri-specific antibodies with two binding units from Group 1 and another binding unit from Group 2 or vice versa would have better binding and blocking function than the bi-specific antibody 21-23. Here, we only focused on tri-specific, as any larger multi-specific molecules will likely affect developability with Fc fusion proteins. We selected the VHHs from both Group 1 and 2 with the most favorable binding, functional and developability features, and constructed tri-specific VHH- Fcs with the computer-aided antibody design using the software BioLuminate (Schrodinger) that enabled their effective construction and optimization. Then, we tested the tri-specific, bi-specific and mono-specific VHH-Fcs for their ability in-vitro for SARS-CoV-2 S1 protein binding and S/ACE2 blocking (Figure 3a and d). As expected, the multi-specific antibodies showed higher binding affinities to SARS-CoV-2 S1 protein RBD in-vitro, with the tri-specific VHH-Fcs 3F-1B- 2A (KD ~ 0.047 nM) and 1B-3F-2A (KD ~ 0.095 nM) showing more potent binding than bi- specific VHH-Fc 1B-3F (Figure 3a, b, c and e). The binding affinities for tri-specific VHH-Fcs were higher than that of individual component VHH-Fcs 1B, 3F and 2A used in combination, and the binding affinity for 1B-3F-Fc was higher than that of individual component VHH-Fcs 1B, and 3F used in combination (Figure 3a). In addition, 3F-1B-2A and 1B-3F-2A showed potent blocking of the SARS-CoV-2 S/ACE2 interaction, with IC50 values of 0.71 nM and 0.74 nM, and full inhibition around 10 nM for both, respectively, that were far superior to using individual component VHH-Fcs as combinations (IC50 of 2.21 nM and full inhibition around 100 nM). In addition, 3F-1B-2A and 1B-3F-2A were more potent than bi-specific VHH-Fc 1B-3F in blocking SARS-CoV-2 S/ACE2 interaction (Figure 3d). Interestingly, the tri-specific VHH-Fc 2A-1B-3F had lower S/ACE2 blocking ability showing the physical arrangement and/or binding orientation of the VHHs in a multi-specific antibody is important for its binding and blocking (Figure 3d). Taken together, this data indicates that the tri-specific VHH-Fcs have a higher synergistic potency in both binding and blocking the S1 or S1/ACE2 interaction than bi-specific or mono- specific antibodies.
Tri-specific VHH-Fcs have favorable developability features
During the computer-aided design process, we incorporated several development- enhancing features in the structures of our VHH-Fcs. Therefore, we analyzed the physico- chemical properties, using DLS and DSF/SLS methods, of our lead bi- and tri-specific antibodies to determine whether they possess favorable characteristics for large-scale manufacturing that is essential for the commercial development of the antibodies (Figure 3e). Our data revealed that the lead tri-specific VHH-Fc 3F-1B-2A has lower aggregation potential based on the DLS method and is thermostable based on the DSF/SLS method (Figure 3e).
Tri-specific VHH-Fc 3F-1B-2A neutralizes SARS-CoV-2 infection in cells
We tested the multi-specific VHH-Fcs for their ability to target SARS-CoV-2 in cell biological functional assays. First, we analyzed the virus neutralizing ability of our antibodies using a pseudovirus that expresses the SARS-CoV-2 S1 protein. The tri-specific VHH-Fcs 3F- 1B-2A, 1B-3F-2A, and the mono-specific combinations of VHHs (1B-Fc + 3F-Fc + 2A-Fc) prevented the infection of human cells by the pseudoviruses (Figure 4a). In accordance with the SARS-CoV-2 S/ACE2 blocking data, the tri-specific VHH-Fcs were more effective in neutralizing the pseudovirus infection than the combination treatment of VHH-Fcs 1B, 3F and 2A, with IC50 values of 3.00 nM for 3F-1B-2A, 6.44 nM for 1B-3F-2A, and 29.19 nM for the combination treatment. (Figure 4a). This pseudovirus data presented here confirm the synergistic effect of the tri-specific antibodies and most importantly, it suggests that it is likely effective in preventing the SARS-CoV-2 infection.
As our VHH-Fcs contain the Fc domain of human IgG1, we expected it would be able to trigger the Fc-dependent functions to eliminate the viruses from the body. To test this, we used a cell line that transiently expresses the SARS-CoV-2 S1 protein. Then, we assessed the ability of our multi-specific VHH-Fcs to promote antibody‐dependent cellular cytotoxicity (ADCC) that is an Fc-dependent function of the antibodies. In addition to our lead tri-specific VHH-Fc antibody 3F-1B-2A, we also tested 3A-3F-2A-Fc another tri-specific antibody we constructed with similar S1 binding and S/ACE2 blocking potency (Supplementary Fig. 3). As expected, the VHH-Fcs were able to induce ADCC in the cells (Figure 4b). This suggests that these VHH-Fcs could bind to immune cells through their Fc domain and elicit Fc-dependent functions, thereby allowing multiple mechanisms of actions against SARS-CoV-2, including binding SARS-CoV-S1 and blocking S1/ACE2 interactions.
Generation of a structure docking model showing the interaction of 3F-1B-2A-Fc with SARS-CoV-2 S1 RBD
Using a computational approach, we generated a 3-dimensional docking model depicting the interaction of our lead tri-specific antibody 3F-1B-2A-Fc with the RBDs of SARS- CoV-2 S1 with BioLuminate software.This model predicted how the individual VHH-Fcs, belonging from two binding groups, can interact with a single RBD using different epitopes. In addition, this model suggests that individual VHH-Fcs of the tri-specific VHH-Fc interact with multiple RBDs in the SARS-CoV-2 S trimer (Figure 5a, b, and c)These modes of interaction are in line with our experimental findings and explain why the tri-specific VHH-Fcs are more potent than the bi-specific or monoclonal VHH-Fcs.