Antibody escape and cryptic cross-domain stabilization in the SARS-CoV-2 Omicron spike protein

Mammalian Antigen Display for Pandemic Countermeasures

Pandemic countermeasures require the rapid design of antigens for vaccines, profiling patient antibody responses, assessing antigen structure-function landscapes, and the surveillance of emerging viral lineages. Cell surface display of a viral antigen or its subdomains can facilitate these goals by coupling the phenotypes of protein variants to their DNA sequence. Screening surface-displayed proteins via flow cytometry also eliminates time-consuming protein purification steps. Prior approaches have primarily relied on yeast as a display chassis. However, yeast often cannot express large viral glycoproteins, requiring their truncation into subdomains. Here, we describe a method to design and express antigens on the surface of mammalian HEK293T cells. We discuss three use cases, including screening of stabilizing mutations, deep mutational scanning, and epitope mapping. The mammalian antigen display platform described herein will accelerate ongoing and future pandemic countermeasures.

Evolution of the SARS-CoV-2 Omicron spike

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron variant of concern, first identified in November 2021, rapidly spread worldwide and diversified into several subvariants. The Omicron spike (S) protein accumulated an unprecedented number of sequence changes relative to previous variants. In this review, we discuss how Omicron S protein structural features modulate host cell receptor binding, virus entry, and immune evasion and highlight how these structural features differentiate Omicron from previous variants. We also examine how key structural properties track across the still-evolving Omicron subvariants and the importance of continuing surveillance of the S protein sequence evolution over time.

SARS-COV-2 Omicron variants conformationally escape a rare quaternary antibody binding mode

The ongoing evolution of SARS-CoV-2 into more easily transmissible and infectious variants has provided unprecedented insight into mutations enabling immune escape. Understanding how these mutations affect the dynamics of antibody-antigen interactions is crucial to the development of broadly protective antibodies and vaccines. Here we report the characterization of a potent neutralizing antibody (N3-1) identified from a COVID-19 patient during the first disease wave. Cryogenic electron microscopy revealed a quaternary binding mode that enables direct interactions with all three receptor-binding domains of the spike protein trimer, resulting in extraordinary avidity and potent neutralization of all major variants of concern until the emergence of Omicron. Structure-based rational design of N3-1 mutants improved binding to all Omicron variants but only partially restored neutralization of the conformationally distinct Omicron BA.1. This study provides new insights into immune evasion through changes in spike protein dynamics and highlights considerations for future conformationally biased multivalent vaccine designs.

Impact of BA.1, BA.2, and BA.4/BA.5 Omicron mutations on therapeutic monoclonal antibodies

The emergence of Omicron SARS-CoV-2 subvariants (BA.1, BA.2, BA.4, and BA.5), with an unprecedented number of mutations in their receptor-binding domain (RBD) of the spike-protein, has fueled a resurgence of COVID-19 infections, posing a major challenge to the efficacy of existing vaccines and monoclonal antibody (mAb) therapeutics. We conducted a systematic molecular dynamics (MD) simulation to investigate how the RBD mutations of these subvariants affect the interactions with broad mAbs including AstraZeneca (COV2-2196 and COV2-2130), Brii Biosciences (BRII-196), Celltrion (CT-P59), Eli Lilly (LY-CoV555 and LY-CoV016), Regeneron (REGN10933 and REGN10987), Vir Biotechnology (S309), and S2X259. Our results show a complete loss of binding for COV2-2196, BRII-196, CT-P59, and LY-CoV555 with all Omicron RBDs. Additionally, REGN10987 totally loses its binding with BA.1, but retains a partial binding with BA.2 and BA.4/5. The binding reduction is significant for LY-CoV016 and REGN10933 but moderate for COV2-2130. S309 and S2X259 retain their binding with BA.1 but exhibit decreased binding with other subvariants. We introduce a mutational escape map for each mAb to identify the key RBD sites and the corresponding critical mutations. Overall, our findings suggest that the majority of therapeutic mAbs have diminished or missing activity against Omicron subvariants, indicating the urgent need for a new therapeutic mAb with a better design.

Graphene Field Effect Biosensor for Concurrent and Specific Detection of SARS-CoV-2 and Influenza

The SARS-CoV-2 pandemic has highlighted the need for devices capable of carrying out rapid differential detection of viruses that may manifest similar physiological symptoms yet demand tailored treatment plans. Seasonal influenza may be exacerbated by COVID-19 infections, increasing the burden on healthcare systems. In this work, we demonstrate a technology based on liquid-gated graphene field-effect transistors (GFETs), for rapid and ultraprecise sensing and differentiation of influenza and SARS-CoV-2 surface protein. Most distinctively, the device consists of 4 onboard GFETs arranged in a quadruple architecture, where each quarter is functionalized individually (with either antibodies or chemically passivated control) but measured jointly. The sensor platform was tested against a range of concentrations of viral surface proteins from both viruses with the lowest tested and detected concentration at ∼50 ag/mL, or 88 zM for COVID-19 and 227 zM for Flu, which is 5-fold lower than the values reported previously on a similar platform. Unlike the classic real-time polymerase chain reaction test, which has a turnaround time of a few hours, the graphene technology presents an ultrafast response time of ∼10 s even in complex and clinically relevant media such as saliva. Thus, we have developed a multianalyte, highly sensitive, and fault-tolerant technology for rapid diagnostic of contemporary, emerging, and future pandemics.

Deep mutational scans of XBB.1.5 and BQ.1.1 reveal ongoing epistatic drift during SARS-CoV-2 evolution

Substitutions that fix between SARS-CoV-2 variants can transform the mutational landscape of future evolution via epistasis. For example, large epistatic shifts in mutational effects caused by N501Y underlied the original emergence of Omicron variants, but whether such large epistatic saltations continue to define ongoing SARS-CoV-2 evolution remains unclear. We conducted deep mutational scans to measure the impacts of all single amino acid mutations and single-codon deletions in the spike receptor-binding domain (RBD) on ACE2-binding affinity and protein expression in the recent Omicron BQ.1.1 and XBB.1.5 variants, and we compared mutational patterns to earlier viral strains that we have previously profiled. As with previous RBD deep mutational scans, we find many mutations that are tolerated or even enhance binding to ACE2 receptor. The tolerance of sites to single-codon deletion largely conforms with tolerance to amino acid mutation. Though deletions in the RBD have not yet been seen in dominant lineages, we observe many tolerated deletions including at positions that exhibit indel variation across broader sarbecovirus evolution and in emerging SARS-CoV-2 variants of interest, most notably the well-tolerated Δ483 deletion in BA.2.86. The substitutions that distinguish recent viral variants have not induced as dramatic of epistatic perturbations as N501Y, but we identify ongoing epistatic drift in SARS-CoV-2 variants, including interaction between R493Q reversions and mutations at positions 453, 455, and 456, including mutations like F456L that define the newly emerging EG.5 lineage. Our results highlight ongoing drift in the effects of mutations due to epistasis, which may continue to direct SARS-CoV-2 evolution into new regions of sequence space.

SARS-CoV-2 spike variants differ in their allosteric responses to linoleic acid

The SARS-CoV-2 spike protein contains a functionally important fatty acid (FA) binding site, which is also found in some other coronaviruses, e.g. SARS-CoV and MERS-CoV. The occupancy of the FA site by linoleic acid (LA) reduces infectivity by 'locking' the spike in a less infectious conformation. Here, we use dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to compare the allosteric responses of spike variants to LA removal. D-NEMD simulations show that the FA site is coupled to other functional regions of the protein, e.g. the receptor-binding motif (RBM), N-terminal domain (NTD), furin cleavage site, and regions surrounding the fusion peptide. D-NEMD simulations also identify the allosteric networks connecting the FA site to these functional regions. The comparison between the wild-type spike and four variants (Alpha, Delta, Delta plus, and Omicron BA.1) shows that the variants differ significantly in their responses to LA removal. The allosteric connections to the FA site on Alpha are generally similar to those on the wild-type protein, with the exception of the RBM and the S71-R78 region, which show a weaker link to the FA site. In contrast, Omicron is the most different variant, exhibiting significant differences in the RBM, NTD, V622-L629, and furin cleavage site. These differences in the allosteric modulation may be of functional relevance, potentially affecting transmissibility and virulence. Experimental comparison of the effects of LA on SARS-CoV-2 variants, including emerging variants, is warranted.

Deep mutagenesis scanning using whole trimeric SARS-CoV-2 spike highlights the importance of NTD-RBD interactions in determining spike phenotype

New variants of SARS-CoV-2 are continually emerging with mutations in spike associated with increased transmissibility and immune escape. Phenotypic maps can inform the prediction of concerning mutations from genomic surveillance, however most of these maps currently derive from studies using monomeric RBD, while spike is trimeric, and contains additional domains. These maps may fail to reflect interdomain interactions in the prediction of phenotypes. To try to improve on this, we developed a platform for deep mutational scanning using whole trimeric spike. We confirmed a previously reported epistatic effect within the RBD affecting ACE2 binding, that highlights the importance of updating the base spike sequence for future mutational scanning studies. Using post vaccine sera, we found that the immune response of vaccinated individuals was highly focused on one or two epitopes in the RBD and that single point mutations at these positions can account for most of the immune escape mediated by the Omicron BA.1 RBD. However, unexpectedly we found that the BA.1 RBD alone does not account for the high level of antigenic escape by BA.1 spike. We show that the BA.1 NTD amplifies the immune evasion of its associated RBD. BA.1 NTD reduces neutralistion by RBD directed monoclonal antibodies, and impacts ACE2 interaction. NTD variation is thus an important mechanism of immune evasion by SARS-CoV-2. Such effects are not seen when pre-stabilized spike proteins are used, suggesting the interdomain effects require protein mobility to express their phenotype.

Spike-protein proteolytic antibodies in COVID-19 convalescent plasma contribute to SARS-CoV-2 neutralization

Understanding the mechanisms of antibody-mediated neutralization of SARS-CoV-2 is critical in combating the COVID-19 pandemic. Based on previous reports of antibody catalysis, we investigated the proteolysis of spike (S) by antibodies in COVID-19 convalescent plasma (CCP) and its contribution to viral neutralization. Quenched fluorescent peptides were designed based on S epitopes to sensitively detect antibody-mediated proteolysis. We observed epitope cleavage by CCP from different donors which persisted when plasma was heat-treated or when IgG was isolated from plasma. Further, purified CCP antibodies proteolyzed recombinant S domains, as well as authentic viral S. Cleavage of S variants suggests CCP antibody-mediated proteolysis is a durable phenomenon despite antigenic drift. We differentiated viral neutralization occurring via direct interference with receptor binding from that occurring by antibody-mediated proteolysis, demonstrating that antibody catalysis enhanced neutralization. These results suggest that antibody-catalyzed damage of S is an immunologically relevant function of neutralizing antibodies against SARS-CoV-2.

Genotype-phenotype landscapes for immune-pathogen coevolution

Our immune systems constantly coevolve with the pathogens that challenge them, as pathogens adapt to evade our defense responses, with our immune repertoires shifting in turn. These coevolutionary dynamics take place across a vast and high-dimensional landscape of potential pathogen and immune receptor sequence variants. Mapping the relationship between these genotypes and the phenotypes that determine immune-pathogen interactions is crucial for understanding, predicting, and controlling disease. Here, we review recent developments applying high-throughput methods to create large libraries of immune receptor and pathogen protein sequence variants and measure relevant phenotypes. We describe several approaches that probe different regions of the high-dimensional sequence space and comment on how combinations of these methods may offer novel insight into immune-pathogen coevolution.

Understanding Mutations in Human SARS-CoV-2 Spike Glycoprotein: A Systematic Review & Meta-Analysis

Genetic variant(s) of concern (VoC) of SARS-CoV-2 have been emerging worldwide due to mutations in the gene encoding spike glycoprotein. We performed comprehensive analyses of spike protein mutations in the significant variant clade of SARS-CoV-2, using the data available on the Nextstrain server. We selected various mutations, namely, A222V, N439K, N501Y, L452R, Y453F, E484K, K417N, T478K, L981F, L212I, N856K, T547K, G496S, and Y369C for this study. These mutations were chosen based on their global entropic score, emergence, spread, transmission, and their location in the spike receptor binding domain (RBD). The relative abundance of these mutations was mapped with global mutation D614G as a reference. Our analyses suggest the rapid emergence of newer global mutations alongside D614G, as reported during the recent waves of COVID-19 in various parts of the world. These mutations could be instrumentally imperative for the transmission, infectivity, virulence, and host immune system's evasion of SARS-CoV-2. The probable impact of these mutations on vaccine effectiveness, antigenic diversity, antibody interactions, protein stability, RBD flexibility, and accessibility to human cell receptor ACE2 was studied in silico. Overall, the present study can help researchers to design the next generation of vaccines and biotherapeutics to combat COVID-19 infection.

Host immunological responses facilitate development of SARS-CoV-2 mutations in patients receiving monoclonal antibody treatments

BackgroundThe role of host immunity in emergence of evasive SARS-CoV-2 Spike mutations under therapeutic monoclonal antibody (mAb) pressure remains to be explored.MethodsIn a prospective, observational, monocentric ORCHESTRA cohort study, conducted between March 2021 and November 2022, mild-to-moderately ill COVID-19 patients (n = 204) receiving bamlanivimab, bamlanivimab/etesevimab, casirivimab/imdevimab, or sotrovimab were longitudinally studied over 28 days for viral loads, de novo Spike mutations, mAb kinetics, seroneutralization against infecting variants of concern, and T cell immunity. Additionally, a machine learning-based circulating immune-related biomarker (CIB) profile predictive of evasive Spike mutations was constructed and confirmed in an independent data set (n = 19) that included patients receiving sotrovimab or tixagevimab/cilgavimab.ResultsPatients treated with various mAbs developed evasive Spike mutations with remarkable speed and high specificity to the targeted mAb-binding sites. Immunocompromised patients receiving mAb therapy not only continued to display significantly higher viral loads, but also showed higher likelihood of developing de novo Spike mutations. Development of escape mutants also strongly correlated with neutralizing capacity of the therapeutic mAbs and T cell immunity, suggesting immune pressure as an important driver of escape mutations. Lastly, we showed that an antiinflammatory and healing-promoting host milieu facilitates Spike mutations, where 4 CIBs identified patients at high risk of developing escape mutations against therapeutic mAbs with high accuracy.ConclusionsOur data demonstrate that host-driven immune and nonimmune responses are essential for development of mutant SARS-CoV-2. These data also support point-of-care decision making in reducing the risk of mAb treatment failure and improving mitigation strategies for possible dissemination of escape SARS-CoV-2 mutants.FundingThe ORCHESTRA project/European Union's Horizon 2020 research and innovation program.

The landscape of antibody binding affinity in SARS-CoV-2 Omicron BA.1 evolution

The Omicron BA.1 variant of SARS-CoV-2 escapes convalescent sera and monoclonal antibodies that are effective against earlier strains of the virus. This immune evasion is largely a consequence of mutations in the BA.1 receptor binding domain (RBD), the major antigenic target of SARS-CoV-2. Previous studies have identified several key RBD mutations leading to escape from most antibodies. However, little is known about how these escape mutations interact with each other and with other mutations in the RBD. Here, we systematically map these interactions by measuring the binding affinity of all possible combinations of these 15 RBD mutations (215=32,768 genotypes) to 4 monoclonal antibodies (LY-CoV016, LY-CoV555, REGN10987, and S309) with distinct epitopes. We find that BA.1 can lose affinity to diverse antibodies by acquiring a few large-effect mutations and can reduce affinity to others through several small-effect mutations. However, our results also reveal alternative pathways to antibody escape that does not include every large-effect mutation. Moreover, epistatic interactions are shown to constrain affinity decline in S309 but only modestly shape the affinity landscapes of other antibodies. Together with previous work on the ACE2 affinity landscape, our results suggest that the escape of each antibody is mediated by distinct groups of mutations, whose deleterious effects on ACE2 affinity are compensated by another distinct group of mutations (most notably Q498R and N501Y).

Deep mutational scanning: A versatile tool in systematically mapping genotypes to phenotypes

Unveiling how genetic variations lead to phenotypic variations is one of the key questions in evolutionary biology, genetics, and biomedical research. Deep mutational scanning (DMS) technology has allowed the mapping of tens of thousands of genetic variations to phenotypic variations efficiently and economically. Since its first systematic introduction about a decade ago, we have witnessed the use of deep mutational scanning in many research areas leading to scientific breakthroughs. Also, the methods in each step of deep mutational scanning have become much more versatile thanks to the oligo-synthesizing technology, high-throughput phenotyping methods and deep sequencing technology. However, each specific possible step of deep mutational scanning has its pros and cons, and some limitations still await further technological development. Here, we discuss recent scientific accomplishments achieved through the deep mutational scanning and describe widely used methods in each step of deep mutational scanning. We also compare these different methods and analyze their advantages and disadvantages, providing insight into how to design a deep mutational scanning study that best suits the aims of the readers' projects.

Inferring selection effects in SARS-CoV-2 with Bayesian Viral Allele Selection

The global effort to sequence millions of SARS-CoV-2 genomes has provided an unprecedented view of viral evolution. Characterizing how selection acts on SARS-CoV-2 is critical to developing effective, long-lasting vaccines and other treatments, but the scale and complexity of genomic surveillance data make rigorous analysis challenging. To meet this challenge, we develop Bayesian Viral Allele Selection (BVAS), a principled and scalable probabilistic method for inferring the genetic determinants of differential viral fitness and the relative growth rates of viral lineages, including newly emergent lineages. After demonstrating the accuracy and efficacy of our method through simulation, we apply BVAS to 6.9 million SARS-CoV-2 genomes. We identify numerous mutations that increase fitness, including previously identified mutations in the SARS-CoV-2 Spike and Nucleocapsid proteins, as well as mutations in non-structural proteins whose contribution to fitness is less well characterized. In addition, we extend our baseline model to identify mutations whose fitness exhibits strong dependence on vaccination status as well as pairwise interaction effects, i.e. epistasis. Strikingly, both these analyses point to the pivotal role played by the N501 residue in the Spike protein. Our method, which couples Bayesian variable selection with a diffusion approximation in allele frequency space, lays a foundation for identifying fitness-associated mutations under the assumption that most alleles are neutral.

Probing the biophysical constraints of SARS-CoV-2 spike N-terminal domain using deep mutational scanning

Increasing the expression level of the SARS-CoV-2 spike (S) protein has been critical for COVID-19 vaccine development. While previous efforts largely focused on engineering the receptor-binding domain (RBD) and the S2 subunit, the amino-terminal domain (NTD) has been long overlooked because of the limited understanding of its biophysical constraints. In this study, the effects of thousands of NTD single mutations on S protein expression were quantified by deep mutational scanning. Our results revealed that in terms of S protein expression, the mutational tolerability of NTD residues was inversely correlated with their proximity to the RBD and S2. We also identified NTD mutations at the interdomain interface that increased S protein expression without altering its antigenicity. Overall, this study not only advances the understanding of the biophysical constraints of the NTD but also provides invaluable insights into S-based immunogen design.

The Comparison of Mutational Progression in SARS-CoV-2: A Short Updated Overview

The COVID-19 pandemic has impacted the world population adversely, posing a threat to human health. In the past few years, various strains of SARS-CoV-2, each with different mutations in its structure, have impacted human health in negative ways. The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) mutations influence the virulence, antibody evasion, and Angiotensin-converting enzyme 2 (ACE2) affinity of the virus. These mutations are essential to understanding how a new strain of SARS-CoV-2 has changed and its possible effects on the human body. This review provides an insight into the spike mutations of SARS-CoV-2 variants. As the current scientific data offer a scattered outlook on the various type of mutations, we aimed to categorize the mutations of Beta (B.1.351), Gamma (P.1), Delta (B.1.612.2), and Omicron (B.1.1.529) systematically according to their location in the subunit 1 (S1) and subunit 2 (S2) domains and summarized their consequences as a result. We also compared the miscellany of mutations that have emerged in all four variants to date. The comparison shows that mutations such as D614G and N501Y have emerged in all four variants of concern and that all four variants have multiple mutations within the N-terminal domain (NTD), as in the case of the Delta variant. Other mutations are scattered in the receptor binding domain (RBD) and subdomain 2 (SD2) of the S1 domain. Mutations in RBD or NTD are often associated with antibody evasion. Few mutations lie in the S2 domain in the Beta, Gamma, and Delta variants. However, in the Omicron variant many mutations occupy the S2 domain, hinting towards a much more evasive virus.