Experimental Evidence for Enhanced Receptor Binding by Rapidly Spreading SARS-CoV-2 Variants

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.     

Disulfide Bonds Play a Critical Role in the Structure and Function of the Receptor-binding Domain of the SARS-CoV-2 Spike Antigen

The current coronavirus pandemic is exerting a tremendously detrimental impact on global health. The Spike proteins of coronaviruses, responsible for cell receptor binding and viral internalization, possess multiple and frequently conserved disulfide bonds raising the question about their role in these proteins. Here, we present a detailed structural and functional investigation of the disulfide bonds of the SARS-CoV-2 Spike receptor-binding domain (RBD). Molecular dynamics simulations of the RBD predict increased flexibility of the surface loops when the four disulfide bonds of the domain are reduced. This flexibility is particularly prominent for the disulfide bond-containing surface loop (residues 456–490) that participates in the formation of the interaction surface with the Spike cell receptor ACE2. In vitro, disulfide bond reducing agents affect the RBD secondary structure, lower its melting temperature from 52 °C to 36–39 °C and decrease its binding affinity to ACE2 by two orders of magnitude at 37 °C. Consistent with these in vitro findings, the reducing agents tris(2-carboxyethyl)phosphine (TCEP) and dithiothreitol (DTT) were able to inhibit viral replication at low millimolar levels in cell-based assays. Our research demonstrates the mechanism by which the disulfide bonds contribute to the molecular structure of the RBD of the Spike protein, allowing the RBD to execute its viral function.     

Structural Modeling of the SARS-CoV-2 Spike/Human ACE2 Complex Interface can Identify High-Affinity Variants Associated with Increased Transmissibility

The COVID-19 pandemic has triggered concerns about the emergence of more infectious and pathogenic viral strains. As a public health measure, efficient screening methods are needed to determine the functional effects of new sequence variants. Here we show that structural modeling of SARS-CoV-2 Spike protein binding to the human ACE2 receptor, the first step in host-cell entry, predicts many novel variant combinations with enhanced binding affinities. By focusing on natural variants at the Spike-hACE2 interface and assessing over 700 mutant complexes, our analysis reveals that high-affinity Spike mutations (including N440K, S443A, G476S, E484R, G502P) tend to cluster near known human ACE2 recognition sites (K31 and K353). These Spike regions are structurally flexible, allowing certain mutations to optimize interface interaction energies. Although most human ACE2 variants tend to weaken binding affinity, they can interact with Spike mutations to generate high-affinity double mutant complexes, suggesting variation in individual susceptibility to infection. Applying structural analysis to highly transmissible variants, we find that circulating point mutations S477N, E484K and N501Y form high-affinity complexes (~40% more than wild-type). By combining predicted affinities and available antibody escape data, we show that fast-spreading viral variants exploit combinatorial mutations possessing both enhanced affinity and antibody resistance, including S477N/E484K, E484K/N501Y and K417T/E484K/N501Y. Thus, three-dimensional modeling of the Spike/hACE2 complex predicts changes in structure and binding affinity that correlate with transmissibility and therefore can help inform future intervention strategies.     

In vitro evolution predicts emerging SARS-CoV-2 mutations with high affinity for ACE2 and cross-species binding

Emerging SARS-CoV-2 variants are creating major challenges in the ongoing COVID-19 pandemic. Being able to predict mutations that could arise in SARS-CoV-2 leading to increased transmissibility or immune evasion would be extremely valuable in development of broad-acting therapeutics and vaccines, and prioritising viral monitoring and containment. Here we use in vitro evolution to seek mutations in SARS-CoV-2 receptor binding domain (RBD) that would substantially increase binding to ACE2. We find a double mutation, S477N and Q498H, that increases affinity of RBD for ACE2 by 6.5-fold. This affinity gain is largely driven by the Q498H mutation. We determine the structure of the mutant-RBD:ACE2 complex by cryo-electron microscopy to reveal the mechanism for increased affinity. Addition of Q498H to SARS-CoV-2 RBD variants is found to boost binding affinity of the variants for human ACE2 and confer a new ability to bind rat ACE2 with high affinity. Surprisingly however, in the presence of the common N501Y mutation, Q498H inhibits binding, due to a clash between H498 and Y501 side chains. To achieve an intermolecular bonding network, affinity gain and cross-species binding similar to Q498H alone, RBD variants with the N501Y mutation must acquire instead the related Q498R mutation. Thus, SARS-CoV-2 RBD can access large affinity gains and cross-species binding via two alternative mutational routes involving Q498, with route selection determined by whether a variant already has the N501Y mutation. These mutations are now appearing in emerging SARS-CoV-2 variants where they have the potential to influence human-to-human and cross-species transmission.     

Mechanism and evolution of human ACE2 binding by SARS-CoV-2 spike

Spike glycoprotein of SARS-CoV-2 mediates viral entry into host cells by facilitating virus attachment and membrane fusion. ACE2 is the main receptor of SARS-CoV-2 and its interaction with spike has shaped the virus’ emergence from an animal reservoir and subsequent evolution in the human host. Many structural studies on the spike:ACE2 interaction have provided insights into mechanisms driving viral evolution during the on-going pandemic. This review describes the molecular basis of spike binding to ACE2, outlines mechanisms that have optimised this interaction during viral evolution, and suggests directions for future research.     

SARS-CoV-2 beta variant substitutions alter spike glycoprotein receptor binding domain structure and stability

The emergence of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) and the subsequent COVID-19 pandemic have visited a terrible cost on the world in the forms of disease, death, and economic turmoil. The rapid development and deployment of extremely effective vaccines against SARS-CoV-2 have seemingly brought within reach the end of the pandemic. However, the virus has acquired mutations. and emerging variants of concern are more infectious and reduce the efficacy of existing vaccines. Although promising efforts to combat these variants are underway, the evolutionary pressures leading to these variants are poorly understood. To that end, here we have studied the effects on the structure and function of the SARS-CoV-2 spike glycoprotein receptor-binding domain of three amino-acid substitutions found in several variants of concern, including alpha (B.1.1.7), beta (B.1.351), and gamma (P.1). We found that these substitutions alter the receptor-binding domain structure, stability, and ability to bind to angiotensin converting enzyme 2, in such a way as to possibly have opposing and compensatory effects. These findings provide new insights into how these variants of concern may have been selected for infectivity while maintaining the structure and stability of the receptor binding domain.

The emergence of severe acute respiratory syndrome-related coronavirus 2 (SARS-CoV-2) in late 2019 and its subsequent spread around the world have caused the deadliest airborne pandemic in the United States, recently surpassing the 1918 influenza pandemic nearly a century ago (1). The international scientific community has risen to the challenge of combating SARS-CoV-2 and COVID-19. The year 2020 ended with the fastest development of vaccine candidates, starting with the genetic sequence of the virus being reported (2) to human trials of novel mRNA-based vaccines within 3 months. Now, there are three SARS-CoV-2 vaccines approved for use within the United States and many more next-generation and pan-coronavirus vaccines currently in development. These advances have made substantial contributions to the control of the COVID-19 pandemic within the United States. Despite multiple manufacturers receiving emergency use authorization and an unprecedented vaccination campaign, significant challenges remain including uncertainty regarding durability, vaccination hesitancy, limited access to healthcare among disadvantaged persons, as well as the continued emergence of variants of concern (VOC). Our ultimate success in quelling this pandemic may lie in our ability, not only to characterize new variants, but also to be able to predict the emergence of new variants. Such advances will require an increased understanding of evolutionary pressures and constraints on viral variation.Three SARS-CoV-2 lineages, the alpha variant lineage B.1.1.7 (or 501Y.V1) first identified within the United Kingdom, the beta variant lineage B.1.351 (or 501Y.V2) identified in South Africa, and the gamma variant lineage P.1 (or 501Y.V3) identified in Brazil, have been demonstrated to possess increased infectivity (3) and in the case, beta and gamma exhibit reduced neutralization by antibodies reacting with the cognate regions of the spike protein within the original Wuhan strain of SARS-CoV-2 (4, 5, 6). The alpha variant possesses the N501Y substitution within the spike glycoprotein receptor-binding domain (RBD) which has been shown to enhance binding to angiotensin converting enzyme 2 (ACE2), the entry receptor for SARS-CoV-2 (7, 8, 9). The beta and gamma variants possess N501Y as well as substitutions at two other sites within the RBD, E484K, and K417N in beta and K417T in gamma (10). These RBD substitutions present in the spike protein of the B.1.351 and P.1 variants have been shown to reduce the binding and neutralization of mRNA vaccine-induced antibodies as well as potent human monoclonal antibodies (11).The consequences of the K417N, E484K, and N501Y substitutions on RBD-ACE2 interactions have also been examined, with the increased infectivity of the alpha variant resulting from the enhanced binding to ACE2 when the RBD N501Y substitution is present (9). The E484K substitution has been shown to enhance ACE2 binding (12) and reduce the efficacy of neutralizing antibodies (13). A recent study examined the effects of the K417N substitution on ACE2 binding and antibody interactions using molecular dynamics and found that K417N disrupts RBD-ACE2 interactions, as well as interactions with a monoclonal antibody (14). However, the effects of these substitutions on the structure of the RBD itself have not been examined. Based on the nature of these substitutions, including residue changes in charge or polar to nonpolar substitutions, we hypothesized that the K417N, E484K, and N501Y substitutions alter the RBD structure and stability as well as ACE2 binding interactions. We studied those changes in single-substitution RBD variants as well as in the RBD containing all three substitutions using molecular dynamics and biophysical approaches. Our data suggest that these VOC substitutions significantly alter RBD structure and stability, with consequences for ACE2 binding and proteolytic susceptibility, having potentially opposing consequences for the fitness of new variants. These findings have implications for viral evolution and the design of subunit vaccine candidates.     

Cryo-electron microscopy structures of the N501Y SARS-CoV-2 spike protein in complex with ACE2 and 2 potent neutralizing antibodies

The recently reported “UK variant” (B.1.1.7) of SARS-CoV-2 is thought to be more infectious than previously circulating strains as a result of several changes, including the N501Y mutation. We present a 2.9-Å resolution cryo-electron microscopy (cryo-EM) structure of the complex between the ACE2 receptor and N501Y spike protein ectodomains that shows Y501 inserted into a cavity at the binding interface near Y41 of ACE2. This additional interaction provides a structural explanation for the increased ACE2 affinity of the N501Y mutant, and likely contributes to its increased infectivity. However, this mutation does not result in large structural changes, enabling important neutralization epitopes to be retained in the spike receptor binding domain. We confirmed this through biophysical assays and by determining cryo-EM structures of spike protein ectodomains bound to 2 representative potent neutralizing antibody fragments.     

Purification and characterization of the receptor‐binding domain of SARS‐CoV‐2 spike protein from Escherichia coli

SARS‐CoV‐2 is responsible for a disruptive worldwide viral pandemic, and renders a severe respiratory disease known as COVID‐19. Spike protein of SARS‐CoV‐2 mediates viral entry into host cells by binding ACE2 through the receptor‐binding domain (RBD). RBD is an important target for development of virus inhibitors, neutralizing antibodies, and vaccines. RBD expressed in mammalian cells suffers from low expression yield and high cost. E. coli is a popular host for protein expression, which has the advantage of easy scalability with low cost. However, RBD expressed by E. coli (RBD‐1) lacks the glycosylation, and its antigenic epitopes may not be sufficiently exposed. In the present study, RBD‐1 was expressed by E. coli and purified by a Ni Sepharose Fast Flow column. RBD‐1 was structurally characterized and compared with RBD expressed by the HEK293 cells (RBD‐2). The secondary structure and tertiary structure of RBD‐1 were largely maintained without glycosylation. In particular, the major β‐sheet content of RBD‐1 was almost unaltered. RBD‐1 could strongly bind ACE2 with a dissociation constant (K_D) of 2.98 × 10^–8 M. Thus, RBD‐1 was expected to apply in the vaccine development, screening drugs and virus test kit.     

Neutralization of SARS‐CoV‐2 by highly potent,hyperthermostable, and mutation‐tolerant nanobodies

Monoclonal anti‐SARS‐CoV‐2 immunoglobulins represent a treatment option for COVID‐19. However, their production in mammalian cells is not scalable to meet the global demand. Single‐domain (VHH) antibodies (also called nanobodies) provide an alternative suitable for microbial production. Using alpaca immune libraries against the receptor‐binding domain (RBD) of the SARS‐CoV‐2 Spike protein, we isolated 45 infection‐blocking VHH antibodies. These include nanobodies that can withstand 95°C. The most effective VHH antibody neutralizes SARS‐CoV‐2 at 17–50 pM concentration (0.2–0.7 µg per liter), binds the open and closed states of the Spike, and shows a tight RBD interaction in the X‐ray and cryo‐EM structures. The best VHH trimers neutralize even at 40 ng per liter. We constructed nanobody tandems and identified nanobody monomers that tolerate the K417N/T, E484K, N501Y, and L452R immune‐escape mutations found in the Alpha, Beta, Gamma, Epsilon, Iota, and Delta/Kappa lineages. We also demonstrate neutralization of the Beta strain at low‐picomolar VHH concentrations. We further discovered VHH antibodies that enforce native folding of the RBD in the E. coli cytosol, where its folding normally fails. Such “fold‐promoting” nanobodies may allow for simplified production of vaccines and their adaptation to viral escape‐mutations.     

Structural basis of SARS-CoV-2 and its variants binding to intermediate horseshoe bat ACE2

Coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has caused a global pandemic. Intermediate horseshoe bats (Rhinolophus affinis) are hosts of RaTG13, the second most phylogenetically related viruses to SARS-CoV-2. We report the binding between intermediate horseshoe bat ACE2 (bACE2-Ra) and SARS-CoV-2 receptor-binding domain (RBD), supporting the pseudotyped SARS-CoV-2 viral infection. A 3.3 Å resolution crystal structure of the bACE2-Ra/SARS-CoV-2 RBD complex was determined. The interaction networks of Patch 1 showed differences in R34 and E35 of bACE2-Ra compared to hACE2 and big-eared horseshoe bat ACE2 (bACE2-Rm). The E35K substitution, existing in other species, significantly enhanced the binding affinity owing to its electrostatic attraction with E484 of SARS-CoV-2 RBD. Furthermore, bACE2-Ra showed extensive support for the SARS-CoV-2 variants. These results broaden our knowledge of the ACE2/RBD interaction mechanism and emphasize the importance of continued surveillance of intermediate horseshoe bats to prevent spillover risk.     

A Comparative Study between Spanish and British SARS-CoV-2 Variants

The study of the interaction between the SARS-CoV-2 spike protein and the angiotensin-converting enzyme 2 (ACE2) receptor is key to understanding binding affinity and stability. In the present report, we sought to investigate the differences between two already sequenced genome variants (Spanish and British) of SARS-CoV-2. Methods: In silico model evaluating the homology, identity and similarity in the genome sequence and the structure and alignment of the predictive spike by computational docking methods. Results: The identity results between the Spanish and British variants of the Spike protein were 28.67%. This close correspondence in the results between the Spanish and British SARS-CoV-2 variants shows that they are very similar (99.99%). The alignment obtained results in four deletions. There were 23 nucleotide substitutions also predicted which could affect the functionality of the proteins produced from this sequence. The interaction between the binding receptor domain from the spike protein and the ACE2 receptor produces some of the mutations found and, therefore, the energy of this ligand varies. However, the estimated antigenicity of the British variant is higher than its Spanish counterpart. Conclusions: Our results indicate that minimal mutations could interfere in the infectivity of the virus due to changes in the fitness between host cell recognition and interaction proteins. In particular, the N501Y substitution, situated in the RBD of the spike of the British variant, might be the reason for its extraordinary infective potential.     

SARS-CoV-2 receptor binding domain fusion protein efficiently neutralizes virus infection

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is responsible for the COVID-19 pandemic. Currently, as dangerous mutations emerge, there is an increased demand for specific treatments for SARS-CoV-2 infected patients. The spike glycoprotein on the virus envelope binds to the angiotensin converting enzyme 2 (ACE2) on host cells through its receptor binding domain (RBD) to mediate virus entry. Thus, blocking this interaction may inhibit viral entry and consequently stop infection. Here, we generated fusion proteins composed of the extracellular portions of ACE2 and RBD fused to the Fc portion of human IgG1 (ACE2-Ig and RBD-Ig, respectively). We demonstrate that ACE2-Ig is enzymatically active and that it can be recognized by the SARS-CoV-2 RBD, independently of its enzymatic activity. We further show that RBD-Ig efficiently inhibits in-vivo SARS-CoV-2 infection better than ACE2-Ig. Mechanistically, we show that anti-spike antibody generation, ACE2 enzymatic activity, and ACE2 surface expression were not affected by RBD-Ig. Finally, we show that RBD-Ig is more efficient than ACE2-Ig at neutralizing high virus titers. We thus propose that RBD-Ig physically blocks virus infection by binding to ACE2 and that RBD-Ig should be used for the treatment of SARS-CoV-2-infected patients.     

Differential binding of SARS-CoV-2 Spike protein variants to its cognate receptor hACE2 using molecular modeling based binding analysis

The current emergence of novel coronavirus, SARS-CoV-2 and its ceaseless expansion worldwide has posed a global health emergency that has adversely affected the humans. With the entire world striving to understand the newly emerged virus, differences in morbidity and infection rate of SARS-CoV-2 have been observed across varied geographic areas, which have been ascribed to viral mutation and evolution over time. The homotrimeric Spike (S) glycoprotein on the viral envelope surface is responsible for binding, priming, and initiating infection in the host. Our phylogeny analysis of 1947 sequences of S proteins indicated there is a change in amino acid (aa) from aspartate (Group-A) to glycine (Group-B) at position 614, near the receptor- binding domain (RBD; aa positions 331-524). The two variants are reported to be in circulation, disproportionately across the world, with Group-A dominant in Asia and Group-B in North America. The trimeric, monomeric, and RBD of S protein of both the variant groups (A & B) were modeled using the Swiss-Model server and were docked with the human receptor angiotensin-converting enzyme 2 (hACE2) employing the PatchDock server and visualized in PyMol. Group-A S protein''s RBD bound imperceptibly to the two binding clefts of the hACE2 protein, on the other hand, Group-B S protein''s RBD perfectly interacted inside the binding clefts of hACE2, with higher number of hydrogen and hydrophobic interactions. This implies that the S protein''s amino acid at position 614 near the core RBD influences its interaction with the cognate hACE2 receptor, which may induce its infectivity that should be explored further with molecular and biochemical studies.     

Screening of inhibitors against SARS-CoV-2 spike protein and their capability to block the viral entry mechanism: A viroinformatics study

SARS-CoV-2, previously named 2019 novel coronavirus (2019-nCoV), has been associated with the global pandemic of acute respiratory distress syndrome. First reported in December 2019 in the Wuhan province of China, this new RNA virus has several folds higher transmission among humans than its other family member (SARS-CoV and MERS-CoV). The SARS-CoV-2 spike receptor-binding domain (RBD) is the region mediating the binding of the virus to host cells via Angiotensin-converting enzyme 2 (ACE2), a critical step of viral. Here in this study, we have utilized in silico approach for the virtual screening of antiviral library extracted from the Asinex database against the Receptor binding domain (RBD) of the S1 subunit of the SARS-CoV-2 spike glycoprotein. Further, the molecules were ranked based on their binding affinity against RBD, and the top 15 molecules were selected. The affinity of these selected molecules to interrupt the ACE2-Spike interaction was also studied. It was found that the chosen molecules were demonstrating excellent binding affinity against spike protein, and these molecules were also very effectively interrupting the ACE2-RBD interaction.Furthermore, molecular dynamics (MD) simulation studies were utilized to investigate the top 3 selected molecules' stability in the ACE2-RBD complexes. To the best of our knowledge, this is the first study where molecules' inhibitory potential against the Receptor binding domain (RBD) of the S1 subunit of the SARS-CoV-2 spike glycoprotein and their inhibitory potential against the ACE2-Spike has been studied. We believe that these compounds can be further tested as a potential therapeutic option against COVID-19.     

Identification of lectin receptors for conserved SARS‐CoV‐2 glycosylation sites

New SARS‐CoV‐2 variants are continuously emerging with critical implications for therapies or vaccinations. The 22 N‐glycan sites of Spike remain highly conserved among SARS‐CoV‐2 variants, opening an avenue for robust therapeutic intervention. Here we used a comprehensive library of mammalian carbohydrate‐binding proteins (lectins) to probe critical sugar residues on the full‐length trimeric Spike and the receptor binding domain (RBD) of SARS‐CoV‐2. Two lectins, Clec4g and CD209c, were identified to strongly bind to Spike. Clec4g and CD209c binding to Spike was dissected and visualized in real time and at single‐molecule resolution using atomic force microscopy. 3D modelling showed that both lectins can bind to a glycan within the RBD‐ACE2 interface and thus interferes with Spike binding to cell surfaces. Importantly, Clec4g and CD209c significantly reduced SARS‐CoV‐2 infections. These data report the first extensive map and 3D structural modelling of lectin‐Spike interactions and uncovers candidate receptors involved in Spike binding and SARS‐CoV‐2 infections. The capacity of CLEC4G and mCD209c lectins to block SARS‐CoV‐2 viral entry holds promise for pan‐variant therapeutic interventions.     

NeutrobodyPlex—monitoring SARS‐CoV‐2 neutralizing immune responses using nanobodies

In light of the COVID‐19 pandemic, there is an ongoing need for diagnostic tools to monitor the immune status of large patient cohorts and the effectiveness of vaccination campaigns. Here, we present 11 unique nanobodies (Nbs) specific for the SARS‐CoV‐2 spike receptor‐binding domain (RBD), of which 8 Nbs potently inhibit the interaction of RBD with angiotensin‐converting enzyme 2 (ACE2) as the major viral docking site. Following detailed epitope mapping and structural analysis, we select two inhibitory Nbs, one of which binds an epitope inside and one of which binds an epitope outside the RBD:ACE2 interface. Based on these, we generate a biparatopic nanobody (bipNb) with viral neutralization efficacy in the picomolar range. Using bipNb as a surrogate, we establish a competitive multiplex binding assay (“NeutrobodyPlex”) for detailed analysis of the presence and performance of neutralizing RBD‐binding antibodies in serum of convalescent or vaccinated patients. We demonstrate that NeutrobodyPlex enables high‐throughput screening and detailed analysis of neutralizing immune responses in infected or vaccinated individuals, to monitor immune status or to guide vaccine design.     

Modelling conformational state dynamics and its role on infection for SARS-CoV-2 Spike protein variants

The SARS-CoV-2 Spike protein needs to be in an open-state conformation to interact with ACE2 to initiate viral entry. We utilise coarse-grained normal mode analysis to model the dynamics of Spike and calculate transition probabilities between states for 17081 variants including experimentally observed variants. Our results correctly model an increase in open-state occupancy for the more infectious D614G via an increase in flexibility of the closed-state and decrease of flexibility of the open-state. We predict the same effect for several mutations on glycine residues (404, 416, 504, 252) as well as residues K417, D467 and N501, including the N501Y mutation recently observed within the B.1.1.7, 501.V2 and P1 strains. This is, to our knowledge, the first use of normal mode analysis to model conformational state transitions and the effect of mutations on such transitions. The specific mutations of Spike identified here may guide future studies to increase our understanding of SARS-CoV-2 infection mechanisms and guide public health in their surveillance efforts.     

Biomechanical characterization of SARS-CoV-2 spike RBD and human ACE2 protein-protein interaction

The current COVID-19 pandemic has led to a devastating impact across the world. Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) (the virus causing COVID-19) is known to use the receptor-binding domain (RBD) at viral surface spike (S) protein to interact with the angiotensin-converting enzyme 2 (ACE2) receptor expressed on many human cell types. The RBD-ACE2 interaction is a crucial step to mediate the host cell entry of SARS-CoV-2. Recent studies indicate that the ACE2 interaction with the SARS-CoV-2 S protein has a higher affinity than its binding with the structurally identical S protein of SARS-CoV-1, the virus causing the 2002–2004 SARS outbreak. However, the biophysical mechanism behind such binding affinity difference is unclear. This study utilizes combined single-molecule force spectroscopy and steered molecular dynamics (SMD) simulation approaches to quantify the specific interactions between SARS-CoV-2 or SARS-CoV-1 RBD and ACE2. Depending on the loading rates, the unbinding forces between SARS-CoV-2 RBD and ACE2 range from 70 to 105 pN and are 30–40% higher than those of SARS-CoV-1 RBD and ACE2 under similar loading rates. SMD results indicate that SARS-CoV-2 RBD interacts with the N-linked glycan on Asn90 of ACE2. This interaction is mostly absent in the SARS-CoV-1 RBD-ACE2 complex. During the SMD simulations, the extra RBD-N-glycan interaction contributes to a greater force and prolonged interaction lifetime. The observation is confirmed by our experimental force spectroscopy study. After removing N-linked glycans on ACE2, its mechanical binding strength with SARS-CoV-2 RBD decreases to a similar level of the SARS-CoV-1 RBD-ACE2 interaction. Together, the study uncovers the mechanism behind the difference in ACE2 binding between SARS-CoV-2 and SARS-CoV-1 and could help develop new strategies to block SARS-CoV-2 entry.     

Free Energy Perturbation Calculations of Mutation Effects on SARS-CoV-2 RBD::ACE2 Binding Affinity

The strength of binding between human angiotensin converting enzyme 2 (ACE2) and the receptor binding domain (RBD) of viral spike protein plays a role in the transmissibility of the SARS-CoV-2 virus. In this study we focus on a subset of RBD mutations that have been frequently observed in infected individuals and probe binding affinity changes to ACE2 using surface plasmon resonance (SPR) measurements and free energy perturbation (FEP) calculations. Our SPR results are largely in accord with previous studies but discrepancies do arise due to differences in experimental methods and to protocol differences even when a single method is used. Overall, we find that FEP performance is superior to that of other computational approaches examined as determined by agreement with experiment and, in particular, by its ability to identify stabilizing mutations. Moreover, the calculations successfully predict the observed cooperative stabilization of binding by the Q498R N501Y double mutant present in Omicron variants and offer a physical explanation for the underlying mechanism. Overall, our results suggest that despite the significant computational cost, FEP calculations may offer an effective strategy to understand the effects of interfacial mutations on protein–protein binding affinities and, hence, in a variety of practical applications such as the optimization of neutralizing antibodies.     

Mutation Y453F in the spike protein of SARS-CoV-2 enhances interaction with the mink ACE2 receptor for host adaption

COVID-19 patients transmitted SARS-CoV-2 to minks in the Netherlands in April 2020. Subsequently, the mink-associated virus (miSARS-CoV-2) spilled back over into humans. Genetic sequences of the miSARS-CoV-2 identified a new genetic variant known as “Cluster 5” that contained mutations in the spike protein. However, the functional properties of these “Cluster 5” mutations have not been well established. In this study, we found that the Y453F mutation located in the RBD domain of miSARS-CoV-2 is an adaptive mutation that enhances binding to mink ACE2 and other orthologs of Mustela species without compromising, and even enhancing, its ability to utilize human ACE2 as a receptor for entry. Structural analysis suggested that despite the similarity in the overall binding mode of SARS-CoV-2 RBD to human and mink ACE2, Y34 of mink ACE2 was better suited to interact with a Phe rather than a Tyr at position 453 of the viral RBD due to less steric clash and tighter hydrophobic-driven interaction. Additionally, the Y453F spike exhibited resistance to convalescent serum, posing a risk for vaccine development. Thus, our study suggests that since the initial transmission from humans, SARS-CoV-2 evolved to adapt to the mink host, leading to widespread circulation among minks while still retaining its ability to efficiently utilize human ACE2 for entry, thus allowing for transmission of the miSARS-CoV-2 back into humans. These findings underscore the importance of active surveillance of SARS-CoV-2 evolution in Mustela species and other susceptible hosts in order to prevent future outbreaks.