Deep mutational scans for ACE2 binding, RBD expression, and antibody escape in the SARS-CoV-2 Omicron BA.1 and BA.2 receptor-binding domains

Epistasis meets pleiotropy in shaping biophysical protein subspaces associated with antimicrobial resistance

Protein space is a rich analogy for genotype-phenotype maps, where amino acid sequence is organized into a high-dimensional space that highlights the connectivity between protein variants. It is a useful abstraction for understanding the process of evolution, and for efforts to engineer proteins towards desirable phenotypes. Few framings of protein space consider how higher-level protein phenotypes can be described in terms of their biophysical dimensions, nor do they rigorously interrogate how forces like epistasis-describing the nonlinear interaction between mutations and their phenotypic consequences-manifest across these dimensions. In this study, we deconstruct a low-dimensional protein space of a bacterial enzyme (dihydrofolate reductase; DHFR) into "subspaces" corresponding to a set of kinetic and thermodynamic traits [(kcat, KM, Ki, and Tm (melting temperature)]. We then examine how three mutations (eight alleles in total) display pleiotropy in their interactions across these subspaces. We extend this approach to examine protein spaces across three orthologous DHFR enzymes (Escherichia coli, Listeria grayi, and Chlamydia muridarum), adding a genotypic context dimension through which epistasis occurs across subspaces. In doing so, we reveal that protein space is a deceptively complex notion, and that the process of protein evolution and engineering should consider how interactions between amino acid substitutions manifest across different phenotypic subspaces.

Persistent topological Laplacian analysis of SARS-CoV-2 variants

Topological data analysis (TDA) is an emerging field in mathematics and data science. Its central technique, persistent homology, has had tremendous success in many science and engineering disciplines. However, persistent homology has limitations, including its inability to handle heterogeneous information, such as multiple types of geometric objects; being qualitative rather than quantitative, e.g., counting a 5-member ring the same as a 6-member ring, and a failure to describe non-topological changes, such as homotopic changes in protein-protein binding. Persistent topological Laplacians (PTLs), such as persistent Laplacian and persistent sheaf Laplacian, were proposed to overcome the limitations of persistent homology. In this work, we examine the modeling and analysis power of PTLs in the study of the protein structures of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike receptor binding domain (RBD). First, we employ PTLs to study how the RBD mutation-induced structural changes of RBD-angiotensin-converting enzyme 2 (ACE2) binding complexes are captured in the changes of spectra of the PTLs among SARS-CoV-2 variants. Additionally, we use PTLs to analyze the binding of RBD and ACE2-induced structural changes of various SARS-CoV-2 variants. Finally, we explore the impacts of computationally generated RBD structures on a topological deep learning paradigm and predictions of deep mutational scanning datasets for the SARS-CoV-2 Omicron BA.2 variant. Our results indicate that PTLs have advantages over persistent homology in analyzing protein structural changes and provide a powerful new TDA tool for data science.

SARS-CoV-2 Omicron (B.1.1.529) shows minimal neurotropism in a double-humanized mouse model

Although severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) initially infects the respiratory tract, it also directly or indirectly affects other organs, including the brain. However, little is known about the relative neurotropism of SARS-CoV-2 variants of concern (VOCs), including Omicron (B.1.1.529), which emerged in November 2021 and has remained the dominant pathogenic lineage since then. To address this gap, we examined the relative ability of Omicron, Beta (B.1.351), and Delta (B.1.617.2) to infect the brain in the context of a functional human immune system by using human angiotensin-converting enzyme 2 (hACE2) knock-in triple-immunodeficient NGC mice with or without reconstitution with human CD34+ stem cells. Intranasal inoculation of huCD34+-hACE2-NCG mice with Beta and Delta resulted in productive infection of the nasal cavity, lungs, and brain on day 3 post-infection, but Omicron was surprisingly unique in its failure to infect either the nasal tissue or brain. Moreover, the same infection pattern was observed in hACE2-NCG mice, indicating that antiviral immunity was not responsible for the lack of Omicron neurotropism. In independent experiments, we demonstrate that nasal inoculation with Beta or with D614G, an ancestral SARS-CoV-2 with undetectable replication in huCD34+-hACE2-NCG mice, resulted in a robust response by human innate immune cells, T cells, and B cells, confirming that exposure to SARS-CoV-2, even without detectable infection, is sufficient to induce an antiviral immune response. Collectively, these results suggest that modeling of the neurologic and immunologic sequelae of SARS-CoV-2 infection requires careful selection of the appropriate SARS-CoV-2 strain in the context of a specific mouse model.

Therapeutic and vaccine-induced cross-reactive antibodies with effector function against emerging Omicron variants

Currently circulating SARS-CoV-2 variants acquired convergent mutations at receptor-binding domain (RBD) hot spots. Their impact on viral infection, transmission, and efficacy of vaccines and therapeutics remains poorly understood. Here, we demonstrate that recently emerged BQ.1.1. and XBB.1 variants bind ACE2 with high affinity and promote membrane fusion more efficiently than earlier Omicron variants. Structures of the BQ.1.1 and XBB.1 RBDs bound to human ACE2 and S309 Fab (sotrovimab parent) explain the altered ACE2 recognition and preserved antibody binding through conformational selection. We show that sotrovimab binds avidly to all Omicron variants, promotes Fc-dependent effector functions and protects mice challenged with BQ.1.1, the variant displaying the greatest loss of neutralization. Moreover, in several donors vaccine-elicited plasma antibodies cross-react with and trigger effector functions against Omicron variants despite reduced neutralizing activity. Cross-reactive RBD-directed human memory B cells remained dominant even after two exposures to Omicron spikes, underscoring persistent immune imprinting. Our findings suggest that this previously overlooked class of cross-reactive antibodies, exemplified by S309, may contribute to protection against disease caused by emerging variants through elicitation of effector functions.

Imprinted SARS-CoV-2 humoral immunity induces convergent Omicron RBD evolution

Continuous evolution of Omicron has led to a rapid and simultaneous emergence of numerous variants that display growth advantages over BA.5 (ref. 1). Despite their divergent evolutionary courses, mutations on their receptor-binding domain (RBD) converge on several hotspots. The driving force and destination of such sudden convergent evolution and its effect on humoral immunity remain unclear. Here we demonstrate that these convergent mutations can cause evasion of neutralizing antibody drugs and convalescent plasma, including those from BA.5 breakthrough infection, while maintaining sufficient ACE2-binding capability. BQ.1.1.10 (BQ.1.1 + Y144del), BA.4.6.3, XBB and CH.1.1 are the most antibody-evasive strains tested. To delineate the origin of the convergent evolution, we determined the escape mutation profiles and neutralization activity of monoclonal antibodies isolated from individuals who had BA.2 and BA.5 breakthrough infections2,3. Owing to humoral immune imprinting, BA.2 and especially BA.5 breakthrough infection reduced the diversity of the neutralizing antibody binding sites and increased proportions of non-neutralizing antibody clones, which, in turn, focused humoral immune pressure and promoted convergent evolution in the RBD. Moreover, we show that the convergent RBD mutations could be accurately inferred by deep mutational scanning profiles4,5, and the evolution trends of BA.2.75 and BA.5 subvariants could be well foreseen through constructed convergent pseudovirus mutants. These results suggest that current herd immunity and BA.5 vaccine boosters may not efficiently prevent the infection of Omicron convergent variants.

Whole-genome single molecule real-time sequencing of SARS-CoV-2 Omicron

New variants and genetic mutations of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) genome can only be identified using accurate sequencing methods. Single molecule real-time (SMRT) sequencing has been used to characterize Alpha and Delta variants, but not Omicron variants harboring numerous mutations in the SARS-CoV-2 genome. This study assesses the performance of a target capture SMRT sequencing protocol for whole genome sequencing (WGS) of SARS-CoV-2 Omicron variants and compared it to that of an amplicon SMRT sequencing protocol optimized for Omicron variants. The failure rate of the target capture protocol (6%) was lower than that of the amplicon protocol (34%, p < 0.001) on our data set, and the median genome coverage with the target capture protocol (98.6% [interquartile range (IQR): 86-99.4]) was greater than that with the amplicon protocol (76.6% [IQR: 66-89.6], [p < 0.001]). The percentages of samples with >95% whole genome coverage were 64% with the target capture protocol and 19% with the amplicon protocol (p < 0.05). The clades of 96 samples determined with both protocols were 93% concordant and the lineages of 59 samples were 100% concordant. Thus, target capture SMRT sequencing appears to be an efficient method for WGS, genotyping and detecting mutations of SARS-CoV-2 Omicron variants.

Emerging Dominant SARS-CoV-2 Variants

Accurate and reliable forecasting of emerging dominant severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants enables policymakers and vaccine makers to get prepared for future waves of infections. The last three waves of SARS-CoV-2 infections caused by dominant variants, Omicron (BA.1), BA.2, and BA.4/BA.5, were accurately foretold by our artificial intelligence (AI) models built with biophysics, genotyping of viral genomes, experimental data, algebraic topology, and deep learning. On the basis of newly available experimental data, we analyzed the impacts of all possible viral spike (S) protein receptor-binding domain (RBD) mutations on the SARS-CoV-2 infectivity. Our analysis sheds light on viral evolutionary mechanisms, i.e., natural selection through infectivity strengthening and antibody resistance. We forecast that BP.1, BL*, BA.2.75*, BQ.1*, and particularly BN.1* have a high potential to become the new dominant variants to drive the next surge. Our key projection about these variants dominance made on Oct. 18, 2022 (see arXiv:2210.09485) became reality in late November 2022.

[Neutralizing antibodies against SARS-CoV-2]

Severe acute respiratory syndrome (SARS) corona virus 2 (SARS-CoV-2) is a novel coronavirus that infects humans and causes respiratory symptoms, resulting in a global pandemic since its appearance in 2019. Neutralizing antibody production is an important immune response following SARS-CoV-2 infection, and a great deal of research has been performed regarding the immune response against SARS-CoV-2 infection. However, SARS-CoV-2 is constantly changing and multiple amino acid reconstitutions accumulated in the spike protein enabled viruses to escape from immune responses, especially from neutralizing antibodies. In this review, the antibody responses to SARS-CoV-2 and the emergence of escape variants, and the development of broadly neutralizing antibodies will be introduced.

Rational identification of potent and broad sarbecovirus-neutralizing antibody cocktails from SARS convalescents

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron sublineages have escaped most receptor-binding domain (RBD)-targeting therapeutic neutralizing antibodies (NAbs), which proves that previous NAb drug screening strategies are deficient against the fast-evolving SARS-CoV-2. Better broad NAb drug candidate selection methods are needed. Here, we describe a rational approach for identifying RBD-targeting broad SARS-CoV-2 NAb cocktails. Based on high-throughput epitope determination, we propose that broad NAb drugs should target non-immunodominant RBD epitopes to avoid herd-immunity-directed escape mutations. Also, their interacting antigen residues should focus on sarbecovirus conserved sites and associate with critical viral functions, making the antibody-escaping mutations less likely to appear. Following these criteria, a featured non-competing antibody cocktail, SA55+SA58, is identified from a large collection of broad sarbecovirus NAbs isolated from SARS-CoV-2-vaccinated SARS convalescents. SA55+SA58 potently neutralizes ACE2-utilizing sarbecoviruses, including circulating Omicron variants, and could serve as broad SARS-CoV-2 prophylactics to offer long-term protection, especially for individuals who are immunocompromised or with high-risk comorbidities.

Identification of broad, potent antibodies to functionally constrained regions of SARS-CoV-2 spike following a breakthrough infection

The antiviral benefit of antibodies can be compromised by viral escape especially for rapidly evolving viruses. Therefore, durable, effective antibodies must be both broad and potent to counter newly emerging, diverse strains. Discovery of such antibodies is critically important for SARS-CoV-2 as the global emergence of new variants of concern (VOC) has compromised the efficacy of therapeutic antibodies and vaccines. We describe a collection of broad and potent neutralizing monoclonal antibodies (mAbs) isolated from an individual who experienced a breakthrough infection with the Delta VOC. Four mAbs potently neutralize the Wuhan-Hu-1 vaccine strain, the Delta VOC, and also retain potency against the Omicron VOCs, including recently circulating BA.4/BA.5, in both pseudovirus-based and live virus assays, and one also potently neutralizes SARS-CoV-1. The potency of these mAbs was greater against Omicron VOCs than all but one of the mAbs that had been approved for therapeutic applications. The mAbs target distinct epitopes on the spike glycoprotein, three in the receptor binding domain (RBD) and one in an invariant region downstream of the RBD in subdomain 1 (SD1). The escape pathways we defined at single amino acid resolution with deep mutational scanning show they target conserved, functionally constrained regions of the glycoprotein, suggesting escape could incur a fitness cost. Overall, these mAbs are novel in their breadth across VOCs, their epitope specificity, and include a highly potent mAb targeting a rare epitope outside of the RBD in SD1.

A bias of Asparagine to Lysine mutations in SARS-CoV-2 outside the receptor binding domain affects protein flexibility

Introduction

COVID-19 pandemic has been threatening public health and economic development worldwide for over two years. Compared with the original SARS-CoV-2 strain reported in 2019, the Omicron variant (B.1.1.529.1) is more transmissible. This variant has 34 mutations in its Spike protein, 15 of which are present in the Receptor Binding Domain (RBD), facilitating viral internalization via binding to the angiotensin-converting enzyme 2 (ACE2) receptor on endothelial cells as well as promoting increased immune evasion capacity.

Methods

Herein we compared SARS-CoV-2 proteins (including ORF3a, ORF7, ORF8, Nucleoprotein (N), membrane protein (M) and Spike (S) proteins) from multiple ancestral strains. We included the currently designated original Variant of Concern (VOC) Omicron, its subsequent emerged variants BA.1, BA2, BA3, BA.4, BA.5, the two currently emerging variants BQ.1 and BBX.1, and compared these with the previously circulating VOCs Alpha, Beta, Gamma, and Delta, to better understand the nature and potential impact of Omicron specific mutations.

Results

Only in Omicron and its subvariants, a bias toward an Asparagine to Lysine (N to K) mutation was evident within the Spike protein, including regions outside the RBD domain, while none of the regions outside the Spike protein domain were characterized by this mutational bias. Computational structural analysis revealed that three of these specific mutations located in the central core region, contribute to a preference for the alteration of conformations of the Spike protein. Several mutations in the RBD which have circulated across most Omicron subvariants were also analysed, and these showed more potential for immune escape.

Conclusion

This study emphasizes the importance of understanding how specific N to K mutations outside of the RBD region affect SARS-CoV-2 conformational changes and the need for neutralizing antibodies for Omicron to target a subset of conformationally dependent B cell epitopes.