Online Inquiry

  • *
  •   
  • *
  •   
  •   
  • *
  • Please input "lifetech" as verification code

Ordering Information

Evolution of the SARS-CoV-2 spike protein in the human host

Abstract

Recently emerged variants of SARS-CoV-2 contain in their surface spike glycoproteins multiple substitutions associated with increased transmission and resistance to neutralising antibodies. We have examined the structure and receptor binding properties of spike proteins from the B.1.1.7 (Alpha) and B.1.351 (Beta) variants to better understand the evolution of the virus in humans. Spikes of both variants have the same mutation, N501Y, in the receptor-binding domains. This substitution confers tighter ACE2 binding, dependent on the common earlier substitution, D614G. Each variant spike has acquired other key changes in structure that likely impact virus pathogenesis. The spike from the Alpha variant is more stable against disruption upon binding ACE2 receptor than all other spikes studied. This feature is linked to the acquisition of a more basic substitution at the S1-S2 furin site (also observed for the variants of concern Delta, Kappa, and Omicron) which allows for near-complete cleavage. In the Beta variant spike, the presence of a new substitution, K417N (also observed in the Omicron variant), in combination with the D614G, stabilises a more open spike trimer, a conformation required for receptor binding. Our observations suggest ways these viruses have evolved to achieve greater transmissibility in humans.

 

Introduction

The SARS-CoV-2 spike glycoprotein is the major surface antigen of the virus. Its function is to bind the host receptor ACE2 and mediate the subsequent membrane fusion required for cell entry1,2,3,4,5,6,7,8. The virus has evolved in the human host during the pandemic9,10,11 and we and others have demonstrated that the predominant D614G substitution, located in a monomer-monomer interface of the spike trimer, increases its propensity to adopt the open conformation that is competent to bind receptor12,13,14. The D614G substitution has been shown to decrease shedding of S1 from spike on virions, consistent with increased stability of the pre-fusion conformation13,15. Recently emerging variants of SARS-CoV-2 have acquired other substitutions in the spike including a number located at the monomer-monomer interfaces, at the receptor-binding site, and near the furin-cleavage site (Supplementary Fig. 1). Here we have examined the structures and receptor-binding properties of spikes of the B.1.1.7 variant first described in Kent, United Kingdom (now termed Alpha)10,16,17,18,19 and the B.1.351 variant first described in South Africa (now termed Beta)20,21. We used pre-fusion-stabilised spikes22 which have been established as a tool to study the receptor-binding properties in the absence of the conformational change associated with fusion. This enabled us to directly compare the pre-fusion spikes of the new variants with those of the original strain (first identified in Wuhan)5,8 and the D614G-only variant12 we described in previous studies, findings that agree with reports that used non-stabilised spikes6,7,13,23.

Two recent studies showed increased RBD erectability and enhanced receptor binding by the Alpha and Beta spikes24,25. Our findings confirm and extend these observations, providing insights into their structures with bound receptor; and revealing several characteristics that these variant spikes evolved to optimise their interaction with the host, which could explain their increased infectivity.

Results and discussion

First, examination of the 2D class averages (Supplementary Fig. 2) of ACE2-bound Alpha spike with its furin site intact reveals that all spikes in this dataset are present as trimers (Fig. 1a, Supplementary Table 1). This is the first cleaved SARS-CoV-2 spike protein in complex with ACE2 we have observed to remain fully trimeric upon receptor binding. A substantial proportion of all other spike/ACE2 complexes are dissociated into monomers5. For example, in our previous study of the furin-cleaved spike of the original Wuhan strain (Wuhan) complexed with ACE2 (Fig. 1a), we observed that more than 70% of the particles present were monomeric S1-S2/ACE2 complexes5. There are many factors that prohibit quantitative estimates of solution equilibria from electron micrographs. However, given that the present study was done under the same conditions as our earlier studies, it is reasonable to conclude that the trimeric Alpha spike is more stable in complex with the receptor than the Wuhan5 or Beta spikes (this study). We also observed that the Alpha variant spike is almost fully cleaved into S1 and S2 (Supplementary Fig. 3a), which was not observed for the Wuhan spike7,8. This observation is consistent with one of the changes in the Alpha spike being the substitution P681H, which generates an even more basic furin-cleavage site (HRRAR). A similar observation has been made in studies of spike material isolated directly from the Alpha variant virus26 and the same P681H substitution has also been observed in the recent B.1.1.529 (Omicron) variant spike.

a Surface representations of predominant molecular species present when ACE2 is mixed with the Alpha variant spike expressed in the presence (left panel) or absence (middle panel) of furin inhibitor I, compared to the predominant species when ACE2 is mixed with the original (Wuhan) spike from our previous study (right panel) (PDB ID 7A915). ACE2 is coloured in green, with monomers of the spike coloured in blue, goldenrod, and rosy brown in the middle panel and rosy brown for S monomers in the right and left panels. b Comparison of the domain movements of S complexed to ACE2 in its trimeric (middle) compared to monomeric (left) form. The monomeric and trimeric ACE2 S1 subunits can be aligned very closely on the ACE2/RBD components but the NTD subdomain (NTD-s, pallid blue) and the NTD (blue) are rotated by 95° and their centre of mass translated by 25 Å in the monomer compared with the trimer, resulting in a conformation incompatible with maintenance of the trimeric state. The furin-cleavage site lies on a loop (shown as dotted line) between the C-terminus of the NTDG and S2 and thus cleavage appears to be able to release sufficient strain to allow the resulting ACE2 complex to remain stable as a trimer. The S1 subunits are shown as surface representation and the S2 core as cartoons. c The ring formed by S1 subunits of the trimer upon binding to three ACE2 molecules (surface representation) remains attached to the exposed S2 core through interactions between the S2′ subunit of one monomer (peach) and NTD-s (pallid blue) and RBD-s (plum) domains of the adjacent one. Domains are coloured: RBD in rosy brown, RBD subdomain (RBD-s) in plum, NTD subdomain (NTD-s) in cyan, NTD in navy, and S2 of the same chain as coloured S1 domains in red.

To test if the cleavage of the S1-S2 subunits is directly responsible for the greater stability of the trimeric ACE2-bound S1/S2 complex, we expressed the Alpha spike in the presence of a furin inhibitor, decanoyl-Arg-Val-Lys-Arg-chloromethylketone, which resulted in an almost completely uncleaved protein (Supplementary Fig. 3b). Incubation of this uncleaved spike with ACE2 resulted in more than 50% of particles being monomeric S1-S2/ACE2 (Fig. 1a, Supplementary Fig. 2). This conformation, like that of the Wuhan strain S1-S2/ACE2 complex, is incompatible with the trimeric structure (Fig. 1b). This might be the result of some local effect not immediately evident, such as the recently reported allosteric effects of ACE2 binding on the local dynamics of the S1-S2 linker27, or a consequence of the (disordered) S2 moiety carried along with the dissociated S1:ACE2 from the uncleaved protein, since preponderance of monomer S/ACE2 species in other SARS-CoV-2 variants, such as Wuhan (Fig. 1a) and Beta (Supplementary Fig. 4), correlates with the lower levels of cleavage of these spikes (Supplementary Fig. 3). In addition our structures of unbound, uncleaved Alpha spike (Supplementary Figs. 5 and 6) suggest the trimeric state of the receptor-bound form of the Alpha variant spike might be further stabilised by the substitutions D1118H and A570D on the inter-monomer interfaces (Supplementary Fig. 5), as also recently suggested by Yang et al.28.

Second, in line with a recent study by Cai et al.24, we found enhanced binding of ACE2 to Alpha and Beta variant spikes when compared to Wuhan by surface biolayer interferometry. The data (Fig. 2a, b; Supplementary Fig. 7) show a sixfold increase in binding strength for Alpha spike, and a twofold increase for Beta spike (compared to Wuhan) arising from the shared substitution N501Y in the RBD. The substitution of the asparagine at position 501 in Wuhan for a tyrosine residue in both Alpha & Beta variants (Fig. 2c) leads to an increase in hydrophobic interactions between the aromatic ring of Y501(RBD) and the aromatic ring of Tyr-41(ACE2) and the aliphatic moiety of Lys-353(ACE2), in addition to a charged hydrogen bond between the phenolic hydroxyl of Tyr-501(RBD) and Lys-353(ACE2) (Fig. 2d). The smaller increase in affinity for ACE2 of the Beta spike versus Alpha spike, is consistent with the finding that, whereas Alpha has retained the same salt bridge as Wuhan between Lys-417(RBD) and Asp-30(ACE2), the RBD of Beta has acquired the additional substitution of an asparagine at residue 417 which cannot make the equivalent salt bridge (Fig. 2e).

a Kd of variant spikes binding to ACE2 measured using biolayer interferometry and calculated from koff/kon analysis (see Supplementary Fig. 7 for details). b Plots of fractional saturation binding measurements with data for the Alpha variant shown in red, Alpha variant with G614D substitution in pink, and Wuhan (D614) spike in blue. Wuhan (*the data for which shown here are adapted from our previous work47) and Alpha D614 spike show almost identical affinity towards ACE2. Similar results were obtained for G614 vs D614 mink (Y453F) spike (Supplementary Figs. 4 and 7). c cryoEM density of the complex of Beta variant S1 with ACE2, with ACE2 coloured in green, RBD in rosy brown, RBD subdomain in plum, with the remaining S1 disseminated density in cream. de Detail of changes in the binding interfaces present in variants (left column) compared to the Wuhan strain (right column). d The N501Y substitution present in both the Beta and Alpha variants allows formation of a new hydrogen bond or a salt bridge. e The K417N substitution present in the Beta variant eliminates a salt bridge between the RBD and ACE2. Source data are provided as a Source Data file.

To understand the evolution of receptor binding, we also expressed Alpha spike with an aspartic acid, rather than glycine residue, at position 614. The substitution D614G (relative to Wuhan) occurred earlier in the evolution of SARS-CoV-2, became the predominant global form of the virus9 and continues to be present in the Alpha and Beta variant forms of the virus. The engineered G614D Alpha spike (Y501, D614) shows the same binding affinity as Wuhan (N501, D614) (Fig. 2a, b; Supplementary Fig. 7). Similarly, we also show that the Y453F substitution in the RBD of mink spike, a structure of which we also report here (Supplementary Figures 46), only increases affinity for human ACE2 if residue 614 is a glycine but not if it is an aspartic acid (Supplementary Fig. 7).

Third, the most striking feature of the structure of Beta spike revealed by cryoEM is that all the trimers adopt an open conformation. This is in contrast to our earlier study, under the same conditions, that showed 83% of Wuhan spike particles were in the closed form (Fig. 3a)8. Again, solution equilibria cannot be extracted from these data but we interpret this near-complete switch in states as indicative of the open form of the Beta variant being more stable than the open form of Wuhan. Two recent studies also report Beta spike showing a higher proportion of open conformation relative to earlier viruses24,25. Inspection of the sequence of the Beta spike and comparison of its structure with that of the closed form of Wuhan spike and the open and closed forms of G614 spikes, suggests that the opening of Beta spike could be driven by the substitution K417N on the background of G614 (both also observed in the Omicron variant spike). In Wuhan, Lysine at residue 417 on the RBD not only makes an aliphatic packing interaction with Tyr-369 of the neighbouring subunit, it also forms a salt bridge/hydrogen bond network with Glu-406 and Arg-403 of the RBD and with Ser-373 of the neighbouring RBD that stabilises the closed conformation. In contrast, the substitution of an asparagine at position 417 in Beta removes an intramolecular salt bridge and would generate a steric clash with Tyr-369 of the neighbouring RBD leading to destabilisation of the closed form and thus promotion of the open form (Fig. 3b). In the same way that G614 is a prerequisite for realising tighter receptor binding by the substitution N501Y described above, it may also enable spike protein from the Beta virus to achieve an open conformation as a result of the K417N substitution.

a (Left and centre) Surface representations of cryoEM structures of the Beta variant spike, with monomers coloured in blue, goldenrod and rosy brown, compared to (right) the closed form of the Wuhan strain determined in our previous study (PDB ID 6ZGE8) shown in paler hues of the same colours. The Beta variant adopts only open conformations with either one (64%) or two (36%) of the RBDs in the erect, receptor-binding-competent form. b The K417N substitution present in the Beta variant likely destabilises the closed form of the protein by introducing a steric clash and interfering with the network of electrostatic interactions present at the trimer interface in the closed form of the Wuhan protein. Strongly charged R417 in the fully-closed Pangolin-CoV spike seems to stabilise the same network of RBD-RBD interactions47.

We have shown that the almost fully cleaved Alpha spike is more stable as a receptor-bound trimer than the Wuhan or Beta spikes which are subject to less efficient cleavage. We cannot explain why near-complete cleavage produces more stable spike/ACE2 complexes but similar increases in stability have been observed following cleavage of the influenza haemagglutinin precursor29,30 and we suggest that such enhanced stability increases the likelihood of productive membrane fusion events. The recent variants of concern B.1.617 (Kappa) and B.1.617.2 (Delta)31,32,33 contain the substitution P681R, which also results in full cleavage (Supplementary Fig. 3c). Our data suggest that this substitution P681R in Kappa and Delta spikes, as well as P681H in Omicron spike, will also increase stability of the receptor-bound form of these spikes, accounting at least in part for their increased transmissibility. Our binding data on the variants and engineered constructs show that the D614G substitution is a prerequisite for the tighter receptor binding of changes in RBD, like N501Y, but do not explain how it facilitates the increase in affinity. We suggest that one of effects of the more open conformation promoted by the D614G substitution is to increase the proportion of accessible RBD-binding sites, thus enhancing the avidity of virus binding to host cells.

We have demonstrated that the emergence of Alpha spike, which is completely cleaved, and Beta spike, which we have not observed to adopt the closed conformation, represent two related recent steps in viral adaptation to the human host. They follow D614G substitution which was acquired early in the pandemic and similarly acted to increase the spike stability12,13,14,15. We suggest that modifications in the spike glycoprotein during evolution of the SARS-CoV-2 virus in humans may have made the virus more infectious by promotion of the open forms of spike, increasing the stability of the pre-fusion-conformation of the receptor-bound trimeric spike, and by tighter receptor binding.

References

  1. Walls, A. C. et al. Cryo-electron microscopy structure of a coronavirus spike glycoprotein trimer. Nature 531, 114–117 (2016).

    ADS CAS Article Google Scholar

  2. Walls, A. C. et al. Tectonic conformational changes of a coronavirus spike glycoprotein promote membrane fusion. PNAS 114, 11157–11162 (2017).

    CAS Article Google Scholar

  3. Walls, A. C. et al. Structure, function, and antigenicity of the SARS-CoV-2 spike glycoprotein. Cell 181, 281–292 (2020).

    CAS Article Google Scholar

  4. Wrapp, D. et al. Cryo-EM structure of the 2019-nCoV spike in the prefusion conformation. Science 367, 1260–1263 (2020).

    ADS CAS Article Google Scholar

  5. Benton, D. J. et al. Receptor binding and priming of the spike protein of SARS-CoV-2 for membrane fusion. Nature 588, 327–330 (2020).

    ADS CAS Article Google Scholar

  6. Cai, Y. et al. Distinct conformational states of SARS-CoV-2 spike protein. Science 369, 1586–1592 (2020).

    ADS CAS Article Google Scholar

  7. Ke, Z. et al. Structures and distributions of SARS-CoV-2 spike proteins on intact virions. Nature 588, 498–502 (2020).

    ADS CAS Article Google Scholar

  8. Wrobel, A. G. et al. SARS-CoV-2 and bat RaTG13 spike glycoprotein structures inform on virus evolution and furin-cleavage effects. Nat. Struct. Mol. Biol. 27, 763–767 (2020).

    CAS Article Google Scholar

  9. Korber, B. et al. Tracking changes in SARS-CoV-2 spike: evidence that D614G increases infectivity of the COVID-19 virus. Cell 182, 812–827.e19 (2020).

    CAS Article Google Scholar

  10. Challen, R. et al. Risk of mortality in patients infected with SARS-CoV-2 variant of concern 202012/1: Matched cohort study. BMJ 372, n579 (2021).

    Article Google Scholar

  11. Davies, N. G. et al. Increased mortality in community-tested cases of SARS-CoV-2 lineage B.1.1.7. Nature 593, 270–274 (2021).

    ADS CAS Article Google Scholar

  12. Benton, D. J. et al. The effect of the D614G substitution on the structure of the spike glycoprotein of SARS-CoV-2. PNAS 118, e2022586118 (2021).

    CAS Article Google Scholar

  13. Zhang, J. et al. Structural impact on SARS-CoV-2 spike protein by D614G substitution. Science 372, 525–530 (2021).

    ADS CAS Article Google Scholar

  14. Gobeil, S. M. C. et al. D614G mutation alters SARS-CoV-2 spike conformation and enhances protease cleavage at the S1/S2 junction. Cell Rep. 34, 108630 (2021).

    CAS Article Google Scholar

  15. Zhang, L. et al. SARS-CoV-2 spike-protein D614G mutation increases virion spike density and infectivity. Nat. Commun. 11, 1–9 (2020).

    ADS Google Scholar

  16. Horby, P. et al. NERVTAG - Presented to SAGE on 21/1/21. https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/961037/NERVTAG_note_on_B.1.1.7_severity_for_SAGE_77__1_.pdf.

  17. Public Health England. Investigation of novel SARS-CoV-2 variant: variant of Concern 202012/01. https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/959360/Variant_of_Concern_VOC_202012_01_Technical_Briefing_3.pdf.

  18. Davies, N. G. et al. Estimated transmissibility and impact of SARS-CoV-2 lineage B.1.1.7 in England. Science. 372, eabg3055 (2021).

    CAS Article Google Scholar

  19. Rambaut, A. et al. Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations. https://virological.org/t/preliminary-genomic-characterisation-of-an-emergent-sars-cov-2-lineage-in-the-uk-defined-by-a-novel-set-of-spike-mutations/563.

  20. Tegally, H. et al. Emergence of a SARS-CoV-2 variant of concern with mutations in spike glycoprotein. Nature 592, 438–443 (2021).

    ADS CAS Article Google Scholar

  21. Tegally, H. et al. Sixteen novel lineages of SARS-CoV-2 in South Africa. Nat. Med. 27, 440–446 (2021).

    CAS Article Google Scholar

  22. Pallesen, J. et al. Immunogenicity and structures of a rationally designed prefusion MERS-CoV spike antigen. PNAS 114, E7348–E7357 (2017).

    CAS Article Google Scholar

  23. Turoňová, B. et al. In situ structural analysis of SARS-CoV-2 spike reveals flexibility mediated by three hinges. Science 370, 203–208 (2020).

    ADS Article Google Scholar

  24. Cai, Y. et al. Structural basis for enhanced infectivity and immune evasion of SARS-CoV-2 variants. Science. eabi9745. https://doi.org/10.1126/SCIENCE.ABI9745 (2021).

  25. Gobeil, S. M.-C. et al. Effect of natural mutations of SARS-CoV-2 on spike structure, conformation, and antigenicity. Science https://doi.org/10.1126/SCIENCE.ABI6226 (2021).

    Article PubMed PubMed Central Google Scholar

  26. Brown, J. C. et al. Increased transmission of SARS-CoV-2 lineage B.1.1.7 (VOC 2020212/01) is not accounted for by a replicative advantage in primary airway cells or antibody escape. Preprint at bioRxiv https://doi.org/10.1101/2021.02.24.432576 (2021).

  27. Raghuvamsi, P. V. et al. Sars-cov-2 s protein:Ace2 interaction reveals novel allosteric targets. Elife 10, 1–47 (2021).

    Article Google Scholar

  28. Yang, T.-J. et al. Effect of SARS-CoV-2 B.1.1.7 mutations on spike protein structure and function. Nat. Struct. Mol. Biol. 28, 731–739 (2021).

    CAS Article Google Scholar

  29. Boulay, F., Doms, R., Wilson, I. & Helenius, A. The influenza hemagglutinin precursor as an acid-sensitive probe of the biosynthetic pathway. EMBO J. 6, 2643–2650 (1987).

    CAS Article Google Scholar

  30. Skehel, J. J., Wharton, S. A., Calder, L. J. & Stevens, D. J. Novel and re-emerging respiratory viral diseases. 56–62 (Novartis Foundation and Wiley, 2008).

  31. Threat Assessment Brief: Emergence of SARS-CoV-2 B.1.617 variants in India and situation in the EU/EEA. https://www.ecdc.europa.eu/en/publications-data/threat-assessment-emergence-sars-cov-2-b1617-variants.

  32. Nervtag. SARS-CoV-2 variants of concern and variants under investigation in England. https://assets.publishing.service.gov.uk/government/uploads/system/uploads/attachment_data/file/984274/Variants_of_Concern_VOC_Technical_Briefing_10_England.pdf (2021).

  33. Wall, E. C. et al. Neutralising antibody activity against SARS-CoV-2 VOCs B.1.617.2 and B.1.351 by BNT162b2 vaccination. Lancet 397, 2331–2333 (2021).

    CAS Article Google Scholar

  34. Ecdc. Detection of new SARS-CoV-2 variants related to mink. https://www.ecdc.europa.eu/sites/default/files/documents/RRA-SARS-CoV-2-in-mink-12-nov-2020.pdf (2020).

  35. Esposito, D. et al. Optimizing high-yield production of SARS-CoV-2 soluble spike trimers for serology assays. Protein Expr. Purif. 174, 105686 (2020).

    CAS Article Google Scholar

  36. Zheng, S. Q. et al. MotionCor2: anisotropic correction of beam-induced motion for improved cryo-electron microscopy. Nat. Methods 14, 331–332 (2017).

    CAS Article Google Scholar

  37. Scheres, S. H. W. RELION: implementation of a Bayesian approach to cryo-EM structure determination. J. Struct. Biol. 180, 519–530 (2012).

    CAS Article Google Scholar

  38. Rohou, A. & Grigorieff, N. CTFFIND4: Fast and accurate defocus estimation from electron micrographs. J. Struct. Biol. 192, 216–221 (2015).

    Article Google Scholar

  39. Wagner, T. et al. SPHIRE-crYOLO is a fast and accurate fully automated particle picker for cryo-EM. Commun. Biol. 2, 218 (2019).

    Article Google Scholar

  40. Punjani, A., Rubinstein, J. L., Fleet, D. J. & Brubaker, M. A. cryoSPARC: algorithms for rapid unsupervised cryo-EM structure determination. Nat. Methods 14, 290–296 (2017).

    CAS Article Google Scholar

  41. Cardone, G., Heymann, J. B. & Steven, A. C. One number does not fit all: mapping local variations in resolution in cryo-EM reconstructions. J. Struct. Biol. 184, 226–236 (2013).

    Article Google Scholar

  42. Rosenthal, P. B. & Henderson, R. Optimal determination of particle orientation, absolute hand, and contrast loss in single-particle electron cryomicroscopy. J. Mol. Biol. 333, 721–745 (2003).

    CAS Article Google Scholar

  43. Emsley, P., Lohkamp, B., Scott, W. G. & Cowtan, K. & IUCr. Features Dev. Coot. Acta Crystallogr. Sect. D. Biol. Crystallogr. 66, 486–501 (2010).

    CAS Article Google Scholar

  44. Adams, P. D. et al. PHENIX: a comprehensive Python-based system for macromolecular structure solution. Acta Crystallogr. Sect. D. Biol. Crystallogr. 66, 213–221 (2010).

    CAS Article Google Scholar

  45. McNicholas, S., Potterton, E., Wilson, K. S. & Noble, M. E. M. Presenting your structures: The CCP4mg molecular-graphics software. Acta Crystallogr. Sect. D. Biol. Crystallogr. 67, 386–394 (2011).

    CAS Article Google Scholar

  46. Pettersen, E. F. et al. UCSF Chimera - a visualization system for exploratory research and analysis. J. Comput. Chem. 25, 1605–1612 (2004).

    CAS Article Google Scholar

  47. Wrobel, A. G. et al. Structure and binding properties of Pangolin-CoV spike glycoprotein inform the evolution of SARS-CoV-2. Nat. Comm. 12, 1–6 (2021).