key: cord-0994760-zhywodm4
authors: Laffeber, Charlie; de Koning, Kelly; Kanaar, Roland; HG Lebbink, Joyce
title: Experimental evidence for enhanced receptor binding by rapidly spreading SARS-CoV-2 variants
date: 2021-05-21
journal: J Mol Biol
DOI: 10.1016/j.jmb.2021.167058
sha: adf4d5ff6c59952e2e0a0fe1d6b97964821d27d8
doc_id: 994760
cord_uid: zhywodm4

Rapidly spreading new variants of SARS-CoV-2 carry multiple mutations in the viral spike protein which attaches to the angiotensin converting enzyme 2 (ACE2) receptor on host cells. Among these mutations are amino acid changes N501Y (lineage B.1.1.7, first identified in the UK), and the combination N501Y, E484K, K417N (B.1.351, first identified in South Africa), all located at the interface on the receptor binding domain (RBD). We experimentally establish that RBD containing the N501Y mutation results in 7-fold stronger binding to the hACE2 receptor than wild type RBD. The E484K mutation only slightly enhances the affinity for the receptor, while K417N attenuates affinity. As a result, RBD from B.1.351 containing all three mutations binds 3-fold stronger to hACE2 than wild type RBD but 2-fold weaker than N501Y. However, the recently emerging double mutant E484K/N501Y binds even stronger than N501Y. The independent evolution of lineages containing mutations with different effects on receptor binding affinity, viral transmission and immune evasion underscores the importance of global viral genome surveillance and functional characterization.

Since its emergence in late 2019, SARS-CoV-2 has rapidly spread across the globe, resulting in more than 160 million confirmed COVID-19 cases and more than 3 million confirmed casualties as of May 15 th , 2021 (coronavirus.jhu.edu). As revealed by the GISAID initiative (gisaid.org), SARS-CoV-2 slowly but continuously mutates, resulting in some instances in variants that become dominant in the population due to increased transmission and/or immune evasion. A number of these variants carry mutations in the spike protein, which is located on the viral surface and interacts with the angiotensin-converting enzyme (ACE2) on host cells, resulting in membrane fusion and viral entry. One such mutation is D614G, which confers increased infectivity and transmissibility and has rapidly become the dominant global variant [1, 2] . Another common mutation is N439K, which is located in the receptor binding domain and enhances affinity for the hACE2 receptor by creating a new salt-bridge across the binding interface [3] . SARS-CoV-2 N439K retains fitness and causes infections with similar clinical outcome, but also shows immune evasion.

In December 2020, new variants of concern have been identified in the UK (B.1.1.7; [4] ), South Africa (B.1.351; [5] ) and Brazil (P.1 and P.2, both descendants from B.1.1.28; [6, 7] ).

These strains carry multiple mutations in the spike protein and form the dominant variants in multiple countries. B.1.1.7 contains a change from asparagine to tyrosine at position 501 (N501Y) in the receptor binding motif of the receptor binding domain (RBD) (Fig. 1C) . A high-throughput deep mutational scan using yeast, as well as recent protein interaction analysis, shows this change increases binding affinity for the hACE2 receptor [8] [9] [10] .

Furthermore, N501Y was identified as adaptive mutation during serial passaging of a clinical SARS-CoV-2 isolate in mice [11] . The N501Y mutation is accompanied by two additional changes at the receptor binding interface in two of the other strains (Fig. 1C) ; glutamate to lysine at position 484 (E484K) and lysine to asparagine or threonine at position 417 (K417N in B.1.351; K417T in P.1). P.2 contains only the E484K change in its receptor binding domain. While for B.1.1.7 increased transmission has been established [12] , increased prevalence of the other lineages may (also) be due to immune escape [5, 13] .

Here we experimentally analyze the effects of different combinations of these RBD mutations on hACE2 receptor binding affinities . Residues N501, E484 and K417 are located relatively far apart on the receptor binding motif (Fig. 1A ) and mutations will directly impact molecular interactions across the interface (Fig. 1B) . Asparagine 501 in the spike RBD forms a single hydrogen bond across the interface with hACE2 tyrosine 41. Mutation of the asparagine into tyrosine will result in loss of this hydrogen bond. In silico structure analysis suggests that the aromatic tyrosine 501 sidechain might be able to stack onto hACE2 tyrosine 41 and form favorable van der Waals interactions using its pi-electron orbitals ( Fig. S1A-B) . Furthermore, the Y501 hydroxyl group may form a hydrogen bond with hACE2 lysine 353. Verification of these predictions awaits the determination of a high-resolution structure of the mutant complex.

E484K is a charge reversal mutation, resulting in the loss of an ion-pair across the interface with hACE2 lysine 31 (Fig. 1B) . In silico predictions suggest that a new ion-pair might be formed with neighboring glutamate 35 in the same hACE2 alpha-helix ( Fig. S1C-D) , probably accompanied by local rearrangements of flexible side chains to avoid energetically unfavorable electrostatics with lysine 31. According to deep mutational scanning analysis, the polar rearrangements upon mutation of glutamate to lysine at position 484 slightly increase affinity for the receptor [8] . The K417N mutation is expected to reduce affinity as replacement of the lysine with a shorter asparagine (in B.1.351, or threonine in P.1) will disrupt the salt-bridge across the interface (Fig. S1E-F) . A similar mutation (K417V) was experimentally shown to reduce affinity 2-fold [3] .

To experimentally determine the effect of different combinations of these mutations on binding affinity for the hACE2 receptor, we purified the hACE2 ectodomain and different SARS-CoV-2 RBD variants from human cells (Fig S2A-D) , and determined rate and affinity constants for complex formation using surface plasmon resonance ( Fig. 1D -F). Complex formation is accurately described by a 1:1 binding model (Fig. S2E ). Obtained rate and affinity constants for wild type RBD (Table 1 ) are similar to previously published values [3, 9, 10] . The increase in affinity for the N501Y variant is 7.1-fold (K D = 2.4 nM instead of 17 nM; Fig. 1F , table 1). This is relatively large for a single amino acid change and typical for mutations that improve the hydrophobic effect, in this case possibly due to the ring stacking of the two tyrosine side chains across the interface. The change in affinity is predominantly caused by a reduction in the dissociation rate constant, indicating the N501Y spike protein remains bound to the receptor for a longer time period than wild type RBD, increasing the chance to undergo the proper conformational change and induce membrane fusion and cell entry.

Compared to the effect of N501Y, the effect of the single E484K mutation on binding affinity is minor (1.4-fold; Fig. 1F , table 1). Obviously, plasticity at the receptor binding interface allows rearrangements of the polar interactions into a new local conformation that is energetically at least as favorable as for the original RBD. Interestingly, while the increase in affinity for N501Y is due to decreased dissociation, increased affinity for E484K is accomplished through faster association. The single K417N mutation destabilizes the interaction with hACE2 4-fold through a combination of slower binding and faster dissociation (Fig. 1F, table 1 ).

The combination of all three mutations, as present in strain B.1.351 (first identified in South Africa), is predominantly additive and results in a 2.4-fold less stable complex than for N501Y alone due to the effect of K417N, but still 3-fold more stable than with wild type RBD Taken together, our results show that receptor binding domains from rapidly spreading variants of SARS-CoV-2 bind with increased affinity to the hACE2 receptor and that this is predominantly caused by the N501Y mutation. This is in agreement with recent experimental observations of a similar set of RBD mutations [9, 14, 15] . The large increase in receptor binding strength of variants carrying N501Y, and a reported overall positive correlation between the stabilizing effect of mutations on receptor binding and their incidence in the population [8, 16] , indicate that affinity for the receptor may be one of the factors that determine viral transmission.

The use of the minimal, highly purified system we describe here allows accurate determination of the intrinsic kinetic and thermodynamic parameters for 1:1 complex formation between a single RBD and single receptor subunit under controlled circumstances.

During viral infection, actual binding of virus particles to host cells will be further influenced by the biochemical environment in the host, as well as avidity effects due to the receptor being dimeric when embedded in the cell membrane, the spike protein being a trimer, and the presence of multiple copies of the receptor and the spike oligomers on the cell and viral surface. Synergy between the N501Y and E484K mutations in a cellular multivalent system seems to be more prominent [16] than in our 1:1 system, and it would be worthwhile to unravel if this is a result of differences in experimental conditions or protein homogeneity or if this is caused by avidity. In addition, within the spike trimer, the receptor binding domains are in a dynamic equilibrium between up and down conformations, with only the up conformation being capable of receptor binding [17] . If the mutations influence this equilibrium, this will affect the association rates for viral attachment.

While the E484K mutation does not have a large effect on affinity for the receptor, it seems to have a significant effect on immune response. Mutation of glutamate to lysine at this position abrogates binding to certain antibodies and results in immune evasion, reinfection and reduced efficacy of vaccines [9, 13, [18] [19] [20] [21] . It is therefore possible that the increased prevalence of lineage P.2, carrying only E484K in its RBD, is due to immune escape rather than increased transmissibility. The recently reported observation of the independent emergence of the E484K mutation into the more transmissible B.1.1.7 strain in for example the UK [9] is of particular concern, as it combines the immune evasion properties of the E484K mutation with N501Y's high affinity, as shown here for the E484K/N501Y double mutant. Additional adaptation seems to have occurred in certain global regions with relatively high previous exposure; mutation of K417 in P.1 and B.1.351 confers additional immune evasion [18] at the expense of RBD-receptor complex stability due to loss of the salt bridge across the interface. The K417N mutation in B.1.351 has occurred as an independent event onto the E484K/N501Y combination before rapidly spreading into the population [5] . As herd immunity will build up due to increased exposure and vaccination, continuous genomic and functional characterization will be crucial to tailor restrictive guidelines, vaccine composition and vaccination strategies towards optimal control of variants with increased receptor affinity, transmissibility and/or immune evasion.

The hACE2 ectodomain (residues 18-615) was fused with the HA signal sequence, a linker containing a 8*HIS-tag, a Twin-Strep-tag and TEV protease cleavage site into vector pCEP4 using Gibson assembly [22] . The same strategy was used to create the expression plasmid for wild type SARS-CoV-2 spike receptor binding domain (reference genome Wuhan-HU-1; residues 333-529). Templates used to create the PCR fragments were pTwist-EF1alpha- Protein purity and concentration were verified using 15% SDS-PAGE ( Fig. S2A-B) .

Glycosylation status of the purified hACE2 and RBD variants was investigated using PNGase Reactions were analyzed using 15% SDS-PAGE (Fig. S2B) .

Immediately prior to surface plasmon resonance (SPR) analysis, aggregates and oligomers were removed from stored protein fractions using size exclusion chromatography (Fig. S2C-D) . ACE2 was injected onto a Superdex 200 Increase 3.2/300, while RBD variants were injected onto a Superdex 75 Increase 3.2/300. Columns were equilibrated with assay buffer (25 mM Hepes-KOH pH 7.5, 150 mM NaCl, 10 µM ZnCl 2 to stabilize hACE2, 0.05% Tween 20) and operated using an Äkta Mikro (Cytiva, Marlborough, USA). The central peak fraction (50 µl) was used for SPR analysis. 

An optimized model for the complex between hACE2 receptor and the SARS-CoV-2 receptor binding domain was created from 6m0j.pdb [23] using PDB-redo [24] . This resulted, amongst others, in a side chain flip for asparagine 501 in the RBD such that it forms a hydrogen bond across the complex interface with Tyrosine 41 in ACE2. Models for mutant RBDs N501Y and E484K were created using Misssense3D [25] . Models were analyzed and figures were created using PyMOL (The PyMOL Molecular Graphics System, Version 2.0 Schrödinger, LLC). 

Tracking Changes in SARS-CoV-2 Spike: Evidence that D614G Increases Infectivity of the COVID-19 Virus

SARS-CoV-2 D614G variant exhibits efficient replication ex vivo and transmission in vivo

Circulating SARS-CoV-2 spike N439K variants maintain fitness while evading antibody-mediated immunity

Preliminary genomic characterisation of an emergent SARS-CoV-2 lineage in the UK defined by a novel set of spike mutations

Emergence and rapid spread of a new severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) lineage with multiple spike mutations in South Africa medRxiv

Genomic characterisation of an emergent SARS-CoV-2 lineage in Manuas: preliminary findings. virologicalorg

Genomic characterization of a novel SARS-CoV-2 lineage from Rio de Janeiro, Brazil

Deep Mutational Scanning of SARS-CoV-2 Receptor Binding Domain Reveals Constraints on Folding and ACE2 Binding

Sensitivity of SARS-CoV-2 B.1.1.7 to mRNA vaccine-elicited antibodies

The basis of a more contagious 501Y.V1 variant of SARS-CoV-2

Adaptation of SARS-CoV-2 in BALB/c mice for testing vaccine efficacy

Transmission of SARS-CoV-2 Lineage B.1.1.7 in England: Insights from linking epidemiological and genetic data

Genomic Evidence of a Sars-Cov-2 Reinfection Case With E484K Spike Mutation in Brazil

Mutation N501Y in RBD of Spike Protein Strengthens the Interaction between COVID-19 and its Receptor ACE2

V2 and 501Y.V3 variants of SARS-CoV-2 lose binding to Bamlanivimab in vitro

SARS-CoV-2 RBD in vitro evolution follows contagious mutations spread, yet generates an able infection inhibitor

D614G Mutation Alters SARS-CoV-2 Spike Conformation and Enhances Protease Cleavage at the S1/S2 Junction

Structural and functional ramifications of antigenic drift in recent SARS-CoV-2 variants

SARS-CoV-2 escape in vitro from a highly neutralizing COVID-19 convalescent plasma

Landscape analysis of escape variants identifies SARS-CoV-2 spike mutations that attenuate monoclonal and serum antibody neutralization

Comprehensive mapping of mutations in the SARS-CoV-2 receptor-binding domain that affect recognition by polyclonal human plasma antibodies

Enzymatic assembly of DNA molecules up to several hundred kilobases

Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor

The PDB_REDO server for macromolecular structure model optimization

Can Predicted Protein 3D Structures Provide Reliable Insights into whether Missense Variants Are Disease Associated?

created expression clones. K.K and C.L purified proteins. C.L. performed SPR analysis. C.L and J.L. analyzed SPR data. J.L performed structure analysis and visualized results. C.L and J.L wrote draft

☒ The authors declare that they have no known competing financial interests or personal relationships that could have appeared to influence the work reported in this paper.☐The authors declare the following financial interests/personal relationships which may be considered as potential competing interests: