Identification of S glycoprotein mutations that constrain the motility of RBD.
The native SARS-CoV-2 S trimer possess some structural flexibility that affects its stability and immunogenicity10,14,15. In addition to pre- and post-fusion S conformations, each RBD displays a dynamic equilibrium between open (up) and closed (down) configurations. In this regard, we aimed to design S variants with a preference for adopting the closed state, and thus, showing limited opening motion and RBD exposure. To this end, we envisioned a computational pipeline involving the three-dimensional modeling of all possible single mutations for both open and closed states, followed by the estimation of changes in their Gibbs free energy (ΔΔG) (Fig. 1a). We focused on all single mutations showing a strong predicted preference for the closed-state (ΔΔG<-1 kcal/mol) and among them, only those that clearly generated well-defined interactions (hydrogen bonds, ionic interactions of filling hydrophobic pockets) between the RBDs of the trimer were screened. We selected a total of 11 single mutations (A372W, K386R, G416R, D420R, D420Y, D427I, L517Y, S982F, D985L, V987H, and V987W) (Fig. 1a). We also included one double mutation (A372W-D420R) referred as 2M, and a combined quintuple mutation (D198F-G232L-A372W-N394Q-D420R) named 5M. Locations of the selected mutations are represented in Fig. 1b. Recombinant mutant proteins were expressed by transient transfection in Expi293F cells, and their production was evaluated by ELISA (Fig. 1c). Most variants displayed a substantial decrease in production when compared to the S-2P trimer (Fig. 1c). Then, we analyzed the exposure of the RBD by ELISA using a Fc fusion protein containing the extracellular portion of the human ACE2 receptor fused to the human IgG1-Fc domain. The results confirmed that most of the variants were in fact promoting a closed trimer conformation (Fig. 1d) as it was predicted by our in silico pipeline. Moreover, variants associated with a reduced exposure of the RBD also resulted in a very low production. In contrast, the V987H mutation promoted the exposure of the RBD (Fig. 1d) and showed higher production than the S-2P protein (2.5-fold). Thus, our results suggest that the RBD exposure may be associated with S production levels.
S-V987H trimer vaccination protects K18-hACE2 mice from SARS-CoV-2 infection-associated disease.
It has been described that K986P and V987P mutations stabilize and increase the expression and immunogenicity of the Spike glycoprotein14,15. Since the V987H mutation improved Spike trimer production and RBD exposure, we evaluated whether it could impact the Spike antigenicity in vivo. Thus, we compared the immunogenicity and protective capacity of the recombinant S-2P, S-V987H and RBD (Fig. 2a) after SARS-CoV-2 D614G challenge in K18-hACE2 mice (Fig. 2b).
Forty-five K18-hACE2 mice were immunized using a prime-boost immunization strategy (Fig. 2b). In addition to S-V987H (n = 14), we determined the immunogenicity of S-2P (n = 16), and a recombinant monomeric RBD protein (n = 15). Mice were first immunized by DNA electroporation with 40 µg of plasmid. Two weeks later, animals were boosted with the corresponding purified recombinant protein (15µg) formulated with aluminum phosphate as adjuvant. Prior to challenge, four mice from S-2P and control groups, two mice from S-V987H and three mice from RBD groups were euthanized to collect tissue samples. Then, 12 vaccinated mice for each group, and 16 unvaccinated controls were intranasally challenged with SARS-CoV-2 D614G (Fig. 2b). Four mice (two male/two female) from each group were euthanized on days 2, 4 and 7 (end of the experiment) post-challenge (Fig. 2b) to analyze the humoral immune responses, viral replication in target organs, and tissue damage. Mice that developed severe SARS-CoV-2 induced disease and/or showed a weight reduction higher than 20% of the initial weight were euthanized before the end of the experiment (day 7) as a humane endpoint. Four additional unvaccinated mice were used as uninfected controls.
The humoral responses elicited against the RBD (Fig. 2c), and the S protein (Supplementary Fig. 1a) were evaluated before each immunization (days − 28 and − 14) and viral challenge (day − 2), and in the mice euthanized on days 2, 4 and 7 after infection. At all-time points, mice immunized with S-V987H and S-2P showed similar levels of anti-RBD and anti-Spike IgG antibodies, which were greater than those observed in RBD vaccinated animals (p < 0.001, Conover-Iman test) (Fig. 2c and Supplementary Fig. 1a, b and c). For simplification purposes, challenged animals were grouped as a “post-challenge” group (Fig. 2c). The levels of anti-RBD and anti-Spike IgG antibodies in the S-2P and S-V987H groups increased after boost and viral challenge (p < 0.05, Conover-Iman test) (Fig. 2c and Supplementary Fig. 1a), while mice immunized with RBD only presented increased levels of these antibodies after viral challenge (p < 0.05, Conover-Iman test), indicating that infection may further boost humoral responses in vaccinated mice (Fig. 2c and Supplementary Fig. 1a, b and c).
In addition, we evaluated level of NAbs against the Wuhan-Hu-1 (WH1) strain and Beta (B.1.351) variant after SARS-CoV-2 challenge. Both S-2P or S-V987H mice groups developed equivalent titers of NAbs against SARS-CoV-2 WH1, but significantly higher than those measured in the RBD immunized group (Fig. 2d p < 0.05, Conover-Iman test). However, only S-2P and S-V987H vaccinated mice had systemic neutralizing activity against the Beta VoC prior to challenge (Fig. 2e). Despite this, the titers of NAbs against WH1 strain were higher than those targeting the Beta VoC (Supplementary Fig. 1d).
The progressive increase in the levels of total IgG antibodies and/or NAbs observed in both RBD and control-challenged groups after infection (Fig. 2c-d) also supports the idea of a boosting of the humoral response after virus challenge.
To assess the ability of each immunogen to prevent SARS-CoV-2 infection-associated disease, we measured weight evolution in all groups as an indicator of disease progression in this model16. We identified weight reduction on day 5 post-infection in mice belonging to infected-control and RBD groups, which is opposed to S-2P and S-V987H vaccinated groups (p < 0.05, Kruskal-Wallis corrected by Dunn’s test) (Fig. 2f). On day 6 and 7, all animals from the inoculated-control group (n = 4), one out of four S-2P vaccinated mice, and two out of four RBD immunized mice had to be euthanized due to disease development (Fig. 2g; p < 0.05 compared to control-infected group, Log-rank test) Mice from the S-V987H group did not show clinical signs of disease during the entire experimental procedure (p < 0.01 compared to control-infected group, Log-rank test) (Fig. 2f and Fig. 2g).
In addition to the clinical course, we measured the levels of viral replication by RT-qPCR in four samples (Fig. 3a-d): oropharyngeal swab, nasal turbinate, lung, and brain. Although the levels of genomic RNA (gRNA) in oropharyngeal swabs decreased in all groups over time, only S-V987H vaccinated mice became undetectable on day 7. On the contrary, S-2P and RBD vaccinated mice remained positive (Fig. 3a, p < 0.05 Peto & Peto Left-censored k sample test). Similar results were observed in nasal turbinates (Fig. 3b), where gRNA was lower in S-2P and S-V987H groups compared to control-infected mice or the RBD group on days 6–7 after viral challenge, and the S-V987H group showed the lowest values by day 7. Interestingly, low levels of gRNA were detected in lung and brain from S-V987H immunized mice (lung p = 0.055, brain p < 0.05, Peto & Peto Left-censored k sample test) (Fig. 3c and d), whereas a progressive increase was observed in brains from the infected-control and RBD groups, and in the only S-2P vaccinated mouse that developed disease (Fig. 3c and 3d). Conversely, subgenomic RNA (sgRNA), which indicates active viral replication, was not detected at any time point in lung or brain of mice immunized with any of the S trimers, or in oropharyngeal swabs on day 4 and 7, except for the single S-2P mouse that developed illness (Supplementary Fig. 1e).
To confirm active viral replication, we analyzed Nucleoprotein (NP) levels in tissues by immunohistochemistry (IHC). NP was detected in lung and brain of both control and RBD groups; and in one animal from S-2P group that developed the disease, but not in S-V987H or disease-free S-2P-vaccinated mice after challenge (Fig. 3e, p < 0.05 asymptotic generalized Pearson Chi-Squared test corrected for multiple comparison using FDR). Low IHC scores were observed in nasal turbinates on days 2 and 4 with no major differences among study groups (Fig. 3e). Tissue damage was in line with the levels of viral antigens detected by IHC (Fig. 3f). No tissue damage was observed in lung or brain of mice vaccinated with S-2P and S-V987H, except for the S-2P mouse that became sick (Fig. 3f). A low lesion score was observed at early time points after challenge in nasal turbinate of all infected mice (Fig. 3f).
Overall, the immunogenicity of both S-2P and S-V987H trimers was equivalent in K18-hACE2 mice, and greater than the produced by the monomeric RBD immunogen. However, S-V987H vaccination improved mice protection against SARS-CoV-2 D614G variant over the S-2P immunogen, since all mice in S-V987H group were disease free and showed a faster viral clearance in tissues.
S-V987H trimer vaccination protects golden Syrian hamsters from SARS-CoV-2 infection-associated disease.
To confirm the results obtained in the transgenic mouse model, we performed a second experiment using golden Syrian hamster (GSH). Similar to K18-hACE2 mice, GSH were immunized using a prime-boost strategy, and intranasally challenged with SARS-CoV-2 D614G (Fig. 4a). Animals were monitored until day 7 post-inoculation, since it has been described that animals start spontaneously recovering around a week after viral infection17,18.
The magnitude of the humoral responses elicited against the S and the RBD by both S-2P and S-V987H trimers was similar and greater than those elicited by the RBD immunogen (Fig. 4b and Supplementary Fig. 2a). The levels of anti-RBD and anti-S IgG antibodies increased after each immunization and after viral challenge (p < 0.05, Friedman test) (Fig. 4b and Supplementary Fig. 2a), confirming the results obtained in K18-hACE2 mice. However, unlike mice, infected-control GSHs rapidly developed an anti-S humoral response after challenge, showing similar levels of anti-S and anti-RBD antibodies on day 7 to those observed in animals immunized with the RBD protein (Fig. 4b and Supplementary Fig. 2b). When the neutralizing activity of serum samples was analyzed, we observed that GSHs immunized with S-2P or S-V987H proteins neutralized the WH1 variant and, to a lesser extent, the Beta VoC (Fig. 4c and d). The neutralizing activity against WH1 increased overtime after challenge in all study groups (p < 0.05, Conover-Iman test). Neutralization of WH1 was also detected in sera from infected control animals by day 4 after challenge, and their titers rapidly increased, becoming similar to the ones observed in S-V987H and RBD groups, and higher than those observed in S-2P vaccinated animals by day 7 (p < 0.05, Conover-Iman test). Intriguingly, despite all groups showed similar titers of NAbs on day 7 after challenge, the levels of anti-RBD and anti-S binding antibodies (Fig. 4b and Supplementary Fig. 2. a and b) were higher in the S-2P and S-V987H immunized groups than in infected-controls GSH. These results support that SARS-CoV-2 infection induced a rapid humoral response against SARS-CoV-2 in GSH that may be qualitatively different to the one elicited by immunization.
We then evaluated the clinical course after challenge. Animals in both control and RBD groups showed a progressive weight reduction until day 7 (end of the experiment) indicative of disease progression (% of weight in infected controls = 87.3 ± 3.1; RBD group = 84.4±1.4) (Fig. 4e). Such weight loss was not observed in S-2P (98.9 ± 1.3) or S-V987H (98.76 ± 2.4) vaccinated GSH (p < 0.05 Kruskal-Wallis corrected by Dunn´s test). Thus, both S trimers generated equivalent protection from disease development in vaccinated GSH (Fig. 4e).
The presence of SARS-CoV-2 was determined by RT-qPCR in oropharyngeal swabs and respiratory tissue samples (nasal turbinate and lung). Brain was not evaluated in GSHs since SARS-CoV-2 does not affect the brain in this animal model 19. Despite the levels of gRNA decreased over time in all analyzed samples, we detected lower gRNA levels in nasal turbinate of both S-2P and S-V987H groups on day 4 compared to the RBD and infected-control groups (Fig. 5; p < 0.05, Peto & Peto Left-censored k sample test). Remarkably, contrary to the RBD or infected-control groups, gRNA was undetectable on day 7 post-infection in S-2P and S-V987H immunized GSHs (Fig. 5a). No major differences were detected in the levels of gRNA in oropharyngeal swabs among the study groups (Fig. 5a). Interestingly, lower levels of gRNA were detected in lungs of S-2P and S-V987H vaccinated animals than in the RBD and infected-control groups on days 2 and 4 (Fig. 5a). While we did not observe significant differences among groups when sgRNA was analyzed in nasal turbinate or oropharyngeal swabs, we detected lower sgRNA levels in the lungs of S-V987H and S-2P groups compared to RBD and infected-control groups on days 2 and 4 after challenge (Supplementary Fig. 2c).
Viral RNA detection results were in line with the levels of NP detection by IHC (Fig. 5b). No major differences in NP levels were observed among study groups at any time points in nasal turbinate, becoming undetectable by day 7 (Fig. 5b) (p < 0.05, Asymptotic Generalized Pearson Chi-Squared Test). However, lower NP levels were detected in lungs of both S-2P and S-V987H vaccinated groups when compared with RBD and infected controls on days 2, 4 and 7. Interestingly, NP was not detected in lungs on day 7 in S-2P and S-V987H groups (Fig. 5b). All study groups showed a similar lesion degree in nasal turbinate, which decreased by day 7 after challenge (p < 0.05). By contrast, a lower tissue damage was observed in lung from S-V987H (on days 4 and 7) and in S-2P (on day 7) groups compared to RBD and infected control groups (p < 0.05) (Fig. 5c).
Overall, our results showed that the immunogenicity and protective efficacy of both S-2P and S-V987H trimers are equivalent in GSHs, and higher than the one conferred by RBD vaccination.
S-V987H trimer vaccination protects K18-hACE2 mice from the SARS-CoV-2 Beta-variant challenge.
From the beginning of the COVID-19 pandemic, several SARS-CoV-2 VoC have emerged. These VoC have shown different transmissibility, pathogenic potential and resistance to antibodies previously elicited by vaccination or natural infection20. The results described above have shown that S-V987H-vaccinated animals were protected from COVID-19 development after SARS-CoV-2 D614G strain challenge. Additionally, vaccinated animals showed low sera neutralizing activity against the SARS-CoV-2 Beta variant. Since the Beta VoC is one of the most resistant to antibodies elicited by natural infection and the currently available vaccines21, and also induces severe disease in K18-hACE2 mice22, we evaluated whether the immune responses induced by S-V987H could protect against disease development after challenge with the SARS-CoV-2 Beta variant. Thus, we immunized twenty-one K18-hACE2 mice with S-V987H or S-2P, using AddVax as adjuvant in this homologous prime-boost experiment (Fig. 6a). Unvaccinated mice were used as negative (n = 10) and positive (n = 16) controls of infection. Two weeks after receiving the protein boost, mice were challenged with the SARS-CoV-2 Beta variant (Fig. 6a). Six mice from each challenged group were euthanized on days 3 (n = 6) and 6 (n = 6) after infection. The remaining animals were euthanized on day 14 after challenge, excepting those mice that developed severe disease after day 3 (10 in the infected-control group and one in the S-2P group) that were euthanized before day 14 following the humane endpoints of the protocol and analyzed separately. Of note, both S-2P and S-V987H recombinant proteins induced similar levels of IgG antibodies against the S and the RBD, which increased after each boost and after viral challenge (p < 0.05, Conover-Iman test) (Fig. 6b and Supplementary Fig. 3a and b). Interestingly, three days after challenge, S-V987H immunized mice showed higher sera neutralizing activity against the WH1 (n = 6; 15376±9203) (Fig. 6c), and the Delta VoC (n = 6; 7750±8403) (Fig. 6d) than mice immunized with the S-2P (n = 6; WH1: 2913±3524; Delta: 1505±4773) (WH1: p < 0.01; Delta: p = 0.055; Conover-Iman test). Neutralizing activity against the Beta VoC increased after challenge (p < 0.05, Conover-Iman test) (Fig. 6e). In addition, we identified an increasing trend in sera neutralizing activity against Omicron over time (p = 0.055, Conover-Iman test) (Fig. 6f). These differences suggest that the humoral responses elicited after S-2P or S-V987H immunization evolved after challenge with the SARS-CoV-2 Beta variant, increasing neutralizing activity against Beta and Omicron VoC, as well as against Delta in the case of the S-2P group. Interestingly, neutralizing activity against the Beta VoC was detected in control-infected mice at clinical endpoint (Fig. 6e) with little or no cross-neutralization activity with other variants (Fig. 6c, d and f). No statistical differences in neutralizing activity were observed on days 6 and 14 between S-V987H and S-2P groups for any of four SARS-CoV-2 variants evaluated (Fig. 6c-f).
A reduction of body weight associated with disease progression was observed in mice from the infected control group starting on day 2 after challenge compared to mice vaccinated with S-2P and S-V987H (Fig. 6g) (p < 0.05, Kruskal-Wallis corrected by Dunn’s test). Mice in both Spike-vaccinated groups maintained their weight until day 14 [percentage of weight: S-2P = 99±4 (n = 9); S-V987H = 98±5 (n = 9)]. Contrarily to the S-2P and infected control groups, no mice from the S-V987H group (n = 9) showed any clinical signs of disease (Fig. 6g and 6h) during the experiment (day 14) (p < 0.001, Long rank test).
The analysis of viral load in tissues by RT-qPCR showed that both S-2P and S-V987H vaccinated groups had a progressive decrease in gRNA levels in oropharyngeal swabs and lung over time (Fig. 7a) (p < 0.05; Peto & Peto Left-censored k sample test). Interestingly, the S-V987H group displayed lower viral loads in nasal turbinate than S-2P and infected control animals on day 3, and also in oropharyngeal swab compared to the infected controls (Fig. 7a) (p < 0.05). However, these differences were not maintained over time and both S-trimer immunized groups showed low but equivalent values of gRNA on day 14 in all analyzed tissues (Fig. 7a). In addition, these groups displayed lower viral load in lung and brain compared to the infected control group at day 3 after challenge (Fig. 7a, p < 0.05).
Remarkably, NP was hardly detected in lung and brain from S-2P and S-V987H groups by IHC (Fig. 7b), which was in line with the low levels of gRNA detected in these animals. Despite that, S-2P vaccinated mice showed a higher lesion score in lung at day 14 than the S-V987H group (p < 0.01; Asymptotic Generalized Pearson Chi-Squared test) (Fig. 7c), indicating that these mice presented a severe lung damage. Interestingly, both Spike-based immunogens protected from viral dissemination to the brain (Fig. 7).
To summarize, the immunogenicity of both S-2P and S-V987H trimers was similar in K18-hACE2 SARS-CoV-2 Beta-infected mice, although S-V987H promoted the development of higher serum neutralization, which might explain the increase in protection observed in S-V987H vaccinated animals, compared to the S-2P group.