SARS-CoV-2 spike N-terminal domain modulates TMPRSS2-dependent viral entry and fusogenicity

Dual-Locking the SARS-CoV-2 Spike Trimer: An Amphipathic Molecular "Bolt" Stabilizes Conserved Druggable Interfaces for Coronavirus Inhibition

The SARS-CoV-2 spike (S) protein, a trimeric structure comprising three receptor binding domains (RBDs) and three N-terminal domains (NTDs), undergoes substantial conformational changes to a fusion-prone open state for angiotensin-converting enzyme 2 (ACE2) binding and host cell infection. Stabilizing its closed state is a key antiviral strategy but remains challenging. Here, we introduce S416, a novel amphipathic molecule acting as a "molecular bolt". Cryo-EM study reveals that S416 binds concurrently to six sites across two distinct druggable interfaces: three molecules at the RBD-RBD interfaces and three at the NTD-RBD interfaces. This unique "dual-locking" mechanism, driven by S416's polar carboxyl head and nonpolar phenylthiazole tail, robustly stabilizes the spike trimer in a locked, closed conformation through strong inter-domain interactions, reducing structural flexibility and atomic fluctuations compared to the apo structure resolved synchronously. Crucially, these RBD-RBD and NTD-RBD interfaces are conserved across human-infecting coronaviruses, suggesting potential as broad-spectrum antiviral targets. Our findings demonstrate that the highly dynamic spike trimer can be effectively stabilized by an amphipathic molecular bolt targeting both the inter- and intra-monomer interfaces, offering a promising strategy against emerging coronaviruses.

Genomic Evolution of the SARS-CoV-2 Omicron Variant in Córdoba, Argentina (2021-2022): Analysis of Uncommon and Prevalent Spike Mutations

Understanding the evolutionary patterns and geographic spread of SARS-CoV-2 variants, particularly Omicron, is essential for effective public health responses. This study focused on the genomic analysis of the Omicron variant in Cordoba, Argentina from 2021 to 2022. Phylogenetic analysis revealed the dominant presence of BA.1 and BA.2 lineages, with BA.5 emerging earlier than BA.4, aligning with observations from other regions. Haplotype network analysis showed significant genetic divergence within Omicron samples, forming distinct clusters. In comparison to global datasets, we identified mutations in the Omicron genomes (A27S, Y145D, and L212I) situated within the NTD region of the Spike protein. These mutations, while not widespread globally, showed higher prevalence in our region. Of particular interest were the Y145D and L212I substitutions, previously unreported in Argentina. In silico analysis revealed that both mutations impact the binding affinity of T-cell epitopes to HLA type I and II alleles. Notably, these alleles are among the most common in the Argentinian population, with some associated with protection against and others with susceptibility to SARS-CoV-2 infection. These findings strongly suggest that these prevalent mutations likely influence the immunogenicity of the Spike protein and contribute to immune evasion mechanisms. This study provides valuable insights into the genomic dynamics of the Omicron variant in Cordoba, Argentina and highlights unique mutations with potential implications for COVID-19 vaccines.

Spatiotemporal evolution and transmission dynamics of Alpha and Delta SARS-CoV-2 variants contributing to sequential outbreaks in Cambodia during 2021

BACKGROUND

Tracking the emergence, introduction and spread of SARS-CoV-2 variants of concern are essential for informing public health strategies. In 2021, Cambodia faced two major epidemic waves of SARS-CoV-2 triggered by the successive rise of the Alpha and Delta variants.

METHODS

Phylodynamic analysis of 1,163 complete SARS-CoV-2 genomes from Cambodia, along with global sequences, were conducted between February and September 2021 to infer viral introductions, molecular epidemiology and population dynamics. The relationship between epidemic trends and control strategies were evaluated. Bayesian phylogeographic reconstruction was employed to estimate and contrast the spatiotemporal dynamics of the Alpha and Delta variants over time.

RESULTS

Here we reveal that the Alpha variant displays rapid lineage diversification, accompanied by the acquisition of a spike E484K mutation that coincides with the national implementation of mass COVID-19 vaccination. Despite nationwide control strategies and increased vaccination coverage, the Alpha variant was quickly displaced by Delta variants that exhibits a higher effective reproductive number. Phylogeographic inference indicates that the Alpha variant was introduced through south-central region of Cambodia, with strong diffusion rates from the capital of Phnom Penh to other provinces, while the Delta variant likely entered the country via the northern border provinces.

CONCLUSIONS

Continual genomic surveillance and sequencing efforts, in combination with public health strategies, play a vital role in effectively tracking and responding to the emergence, evolution and dissemination of future emerging variants.

The accomplices: Heparan sulfates and N-glycans foster SARS-CoV-2 spike:ACE2 receptor binding and virus priming

Although it is well established that the SARS-CoV-2 spike glycoprotein binds to the host cell ACE2 receptor to initiate infection, far less is known about the tissue tropism and host cell susceptibility to the virus. Differential expression across different cell types of heparan sulfate (HS) proteoglycans, with variably sulfated glycosaminoglycans (GAGs), and their synergistic interactions with host and viral N-glycans may contribute to tissue tropism and host cell susceptibility. Nevertheless, their contribution remains unclear since HS and N-glycans evade experimental characterization. We, therefore, carried out microsecond-long all-atom molecular dynamics simulations, followed by random acceleration molecular dynamics simulations, of the fully glycosylated spike:ACE2 complex with and without highly sulfated GAG chains bound. By considering the model GAGs as surrogates for the highly sulfated HS expressed in lung cells, we identified key cell entry mechanisms of spike SARS-CoV-2. We find that HS promotes structural and energetic stabilization of the active conformation of the spike receptor-binding domain (RBD) and reorientation of ACE2 toward the N-terminal domain in the same spike subunit as the RBD. Spike and ACE2 N-glycans exert synergistic effects, promoting better packing, strengthening the protein:protein interaction, and prolonging the residence time of the complex. ACE2 and HS binding trigger rearrangement of the S2' functional protease cleavage site through allosteric interdomain communication. These results thus show that HS has a multifaceted role in facilitating SARS-CoV-2 infection, and they provide a mechanistic basis for the development of GAG derivatives with anti-SARS-CoV-2 potential.

SARS-CoV-2 Variants from Long-Term, Persistently Infected Immunocompromised Patients Have Altered Syncytia Formation, Temperature-Dependent Replication, and Serum Neutralizing Antibody Escape

SARS-CoV-2 infection of immunocompromised individuals often leads to prolonged detection of viral RNA and infectious virus in nasal specimens, presumably due to the lack of induction of an appropriate adaptive immune response. Mutations identified in virus sequences obtained from persistently infected patients bear signatures of immune evasion and have some overlap with sequences present in variants of concern. We characterized virus isolates obtained greater than 100 days after the initial COVID-19 diagnosis from two COVID-19 patients undergoing immunosuppressive cancer therapy, wand compared them to an isolate from the start of the infection. Isolates from an individual who never mounted an antibody response specific to SARS-CoV-2 despite the administration of convalescent plasma showed slight reductions in plaque size and some showed temperature-dependent replication attenuation on human nasal epithelial cell culture compared to the virus that initiated infection. An isolate from another patient-who did mount a SARS-CoV-2 IgM response-showed temperature-dependent changes in plaque size as well as increased syncytia formation and escape from serum-neutralizing antibodies. Our results indicate that not all virus isolates from immunocompromised COVID-19 patients display clear signs of phenotypic change, but increased attention should be paid to monitoring virus evolution in this patient population.

Glycan masking of NTD loops with a chimeric RBD of the spike protein as a vaccine design strategy against emerging SARS-CoV-2 Omicron variants

The N-terminal domain (NTD) of the SARS-CoV-2 S protein comprises five exposed protruding loops. Deletions, insertions, and substitutions within these NTD loops play a significant role in viral evolution and contribute to immune evasion. We reported previously that introducing the glycan masking mutation R158N/Y160T in the NTD loop led to increased titers of neutralizing antibodies against the SARS-CoV-2 Wuhan-Hu-01 strain, as well as the Alpha, Beta, and Delta variants. In this study, we conducted further investigations on 10 additional glycan-masking sites in the NTD loops. Our findings indicate that the introduction of glycan masking mutations, specifically N87/G89T, H146N/N148T, N185/K187T, and V213N/D215T significantly enhanced neutralizing antibody titers against the Delta variant. The combination of dual glycan-masking mutations R158N/Y160T+V213N/D215T and R158N/Y160T+G219N results in a shift toward the Omicron BA.1. Furthermore, the introduction of the Omicron receptor binding domain (RBD) alongside these two dual glycan masking mutations of Wuhan-Hu-1 and XBB.1 NTD sequences resulted in a noticeable shift in antigenic distances, aligning with the Omicron BA.4/5, BA.2.75.2, BQ.1.1, and XBB.1 subvariants on the antigenic map. This strategic combination, which involves the dual glycan masking mutations R158N/Y160T+V213N/D215T and R158N/Y160T+G219N in the NTD loops, along with the domain swap incorporating the Omicron RBD, emerges as a promising vaccine design strategy for the continuous development of next-generation SARS-CoV-2 vaccines.

Primary SARS-CoV-2 variant of concern infections elicit broad antibody Fc-mediated effector functions and memory B cell responses

Neutralization of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) by human sera is a strong correlate of protection against symptomatic and severe Coronavirus Disease 2019 (COVID-19). The emergence of antigenically distinct SARS-CoV-2 variants of concern (VOCs) and the relatively rapid waning of serum antibody titers, however, raises questions about the sustainability of serum protection. In addition to serum neutralization, other antibody functionalities and the memory B cell (MBC) response are suggested to help maintaining this protection. In this study, we investigate the breadth of spike (S) protein-specific serum antibodies that mediate effector functions by interacting with Fc-gamma receptor IIa (FcγRIIa) and FcγRIIIa, and of the receptor binding domain (RBD)-specific MBCs, following a primary SARS-CoV-2 infection with the D614G, Alpha, Beta, Gamma, Delta, Omicron BA.1 or BA.2 variant. Irrespectively of the variant causing the infection, the breadth of S protein-specific serum antibodies that interact with FcγRIIa and FcγRIIIa and the RBD-specific MBC responses exceeded the breadth of serum neutralization, although the Alpha-induced B cell response seemed more strain-specific. Between VOC groups, both quantitative and qualitative differences in the immune responses were observed, suggesting differences in immunogenicity. Overall, this study contributes to the understanding of protective humoral and B cell responses in the light of emerging antigenically distinct VOCs, and highlights the need to study the immune system beyond serum neutralization to gain a better understanding of the protection against emerging variants.

Single-molecule imaging reveals allosteric stimulation of SARS-CoV-2 spike receptor binding domain by host sialic acid

Conformational dynamics of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike glycoprotein (S) mediate exposure of the binding site for the cellular receptor, angiotensin-converting enzyme 2 (ACE2). The N-terminal domain (NTD) of S binds terminal sialic acid (SA) moieties on the cell surface, but the functional role of this interaction in virus entry is unknown. Here, we report that NTD-SA interaction enhances both S-mediated virus attachment and ACE2 binding. Through single-molecule Förster resonance energy transfer imaging of individual S trimers, we demonstrate that SA binding to the NTD allosterically shifts the S conformational equilibrium, favoring enhanced exposure of the ACE2-binding site. Antibodies that target the NTD block SA binding, which contributes to their mechanism of neutralization. These findings inform on mechanisms of S activation at the cell surface.

Spike N354 glycosylation augments SARS-CoV-2 fitness for human adaptation through structural plasticity

Selective pressures have given rise to a number of SARS-CoV-2 variants during the prolonged course of the COVID-19 pandemic. Recently evolved variants differ from ancestors in additional glycosylation within the spike protein receptor-binding domain (RBD). Details of how the acquisition of glycosylation impacts viral fitness and human adaptation are not clearly understood. Here, we dissected the role of N354-linked glycosylation, acquired by BA.2.86 sub-lineages, as a RBD conformational control element in attenuating viral infectivity. The reduced infectivity is recovered in the presence of heparin sulfate, which targets the 'N354 pocket' to ease restrictions of conformational transition resulting in a 'RBD-up' state, thereby conferring an adjustable infectivity. Furthermore, N354 glycosylation improved spike cleavage and cell-cell fusion, and in particular escaped one subset of ADCC antibodies. Together with reduced immunogenicity in hybrid immunity background, these indicate a single spike amino acid glycosylation event provides selective advantage in humans through multiple mechanisms.

Nanomolar anti-SARS-CoV-2 Omicron activity of the host-directed TMPRSS2 inhibitor N-0385 and synergistic action with direct-acting antivirals

SARS-CoV-2 Omicron subvariants with increased transmissibility and immune evasion are spreading globally with alarming persistence. Whether the mutations and evolution of spike (S) Omicron subvariants alter the viral hijacking of human TMPRSS2 for viral entry remains to be elucidated. This is particularly important to investigate because of the large number and diversity of mutations of S Omicron subvariants reported since the emergence of BA.1. Here we report that human TMPRSS2 is a molecular determinant of viral entry for all the Omicron clinical isolates tested in human lung cells, including ancestral Omicron subvariants (BA.1, BA.2, BA.5), contemporary Omicron subvariants (BQ.1.1, XBB.1.5, EG.5.1) and currently circulating Omicron BA.2.86. First, we used a co-transfection assay to demonstrate the endoproteolytic cleavage by TMPRSS2 of spike Omicron subvariants. Second, we found that N-0385, a highly potent TMPRSS2 inhibitor, is a robust entry inhibitor of virus-like particles harbouring the S protein of Omicron subvariants. Third, we show that N-0385 exhibits nanomolar broad-spectrum antiviral activity against live Omicron subvariants in human Calu-3 lung cells and primary patient-derived bronchial epithelial cells. Interestingly, we found that N-0385 is 10-20 times more potent than the repositioned TMPRSS2 inhibitor, camostat, against BA.5, EG.5.1, and BA.2.86. We further found that N-0385 shows broad synergistic activity with clinically approved direct-acting antivirals (DAAs), i.e., remdesivir and nirmatrelvir, against Omicron subvariants, demonstrating the potential therapeutic benefits of a multi-targeted treatment based on N-0385 and DAAs.

SARS-CoV-2 biology and host interactions

The zoonotic emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and the ensuing coronavirus disease 2019 (COVID-19) pandemic have profoundly affected our society. The rapid spread and continuous evolution of new SARS-CoV-2 variants continue to threaten global public health. Recent scientific advances have dissected many of the molecular and cellular mechanisms involved in coronavirus infections, and large-scale screens have uncovered novel host-cell factors that are vitally important for the virus life cycle. In this Review, we provide an updated summary of the SARS-CoV-2 life cycle, gene function and virus-host interactions, including recent landmark findings on general aspects of coronavirus biology and newly discovered host factors necessary for virus replication.

Therapeutic development targeting host heparan sulfate proteoglycan in SARS-CoV-2 infection

The global pandemic caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has led to an urgent need for effective therapeutic options. SARS-CoV-2 is a novel coronavirus responsible for the COVID-19 pandemic that has resulted in significant morbidity and mortality worldwide. The virus is known to enter host cells by binding to the angiotensin-converting enzyme 2 (ACE2) receptor, and emerging evidence suggests that heparan sulfate proteoglycans (HSPGs) play a crucial role in facilitating this process. HSPGs are abundant cell surface proteoglycan present in many tissues, including the lung, and have been shown to interact directly with the spike protein of SARS-CoV-2. This review aims to summarize the current understanding of the role of HSPGs in SARS-CoV-2 infection and the potential of developing new therapies targeting HSPGs.

Molecular characterization of SARS-CoV-2 Omicron clade and clinical presentation in children

Since its emergence, SARS-CoV-2 Omicron clade has shown a marked degree of variability and different clinical presentation compared with previous clades. Here we demonstrate that at least four Omicron lineages circulated in children since December 2021, and studied until November 2022: BA.1 (33.6%), BA.2 (40.6%), BA.5 (23.7%) and BQ.1 (2.1%). At least 70% of infections concerned children under 1 year, most of them being infected with BA.2 lineages (n = 201, 75.6%). Looking at SARS-CoV-2 genetic variability, 69 SNPs were found to be significantly associated in pairs, (phi <  - 0.3 or > 0.3 and p-value < 0.001). 16 SNPs were involved in 4 distinct clusters (bootstrap > 0.75). One of these clusters (A23040G, A27259C, T23617G, T23620G) was also positively associated with moderate/severe COVID-19 presentation (AOR [95% CI] 2.49 [1.26-4.89] p-value: 0.008) together with comorbidities (AOR [95% CI] 2.67 [1.36-5.24] p-value: 0.004). Overall, these results highlight the extensive SARS-CoV-2 Omicron circulation in children, mostly aged < 1 year, and provide insights on viral diversification even considering low-abundant SNPs, finally suggesting the potential contribution of viral diversification in affecting disease severity.

AlphaFold2 Reveals Structural Patterns of Seasonal Haplotype Diversification in SARS-CoV-2 Spike Protein Variants

The slow experimental acquisition of high-quality atomic structures of the rapidly changing proteins of the COVID-19 virus challenges vaccine and therapeutic drug development efforts. Fortunately, deep learning tools such as AlphaFold2 can quickly generate reliable models of atomic structure at experimental resolution. Current modeling studies have focused solely on definitions of mutant constellations of Variants of Concern (VOCs), leaving out the impact of haplotypes on protein structure. Here, we conduct a thorough comparative structural analysis of S-proteins belonging to major VOCs and corresponding latitude-delimited haplotypes that affect viral seasonal behavior. Our approach identified molecular regions of importance as well as patterns of structural recruitment. The S1 subunit hosted the majority of structural changes, especially those involving the N-terminal domain (NTD) and the receptor-binding domain (RBD). In particular, structural changes in the NTD were much greater than just translations in three-dimensional space, altering the sub-structures to greater extents. We also revealed a notable pattern of structural recruitment with the early VOCs Alpha and Delta behaving antagonistically by suppressing regions of structural change introduced by their corresponding haplotypes, and the current VOC Omicron behaving synergistically by amplifying or collecting structural change. Remarkably, haplotypes altering the galectin-like structure of the NTD were major contributors to seasonal behavior, supporting its putative environmental-sensing role. Our results provide an extensive view of the evolutionary landscape of the S-protein across the COVID-19 pandemic. This view will help predict important regions of structural change in future variants and haplotypes for more efficient vaccine and drug development.

SARS-CoV-2 and innate immunity: the good, the bad, and the "goldilocks"

An ancient conflict between hosts and pathogens has driven the innate and adaptive arms of immunity. Knowledge about this interplay can not only help us identify biological mechanisms but also reveal pathogen vulnerabilities that can be leveraged therapeutically. The humoral response to SARS-CoV-2 infection has been the focus of intense research, and the role of the innate immune system has received significantly less attention. Here, we review current knowledge of the innate immune response to SARS-CoV-2 infection and the various means SARS-CoV-2 employs to evade innate defense systems. We also consider the role of innate immunity in SARS-CoV-2 vaccines and in the phenomenon of long COVID.

Cooperative and structural relationships of the trimeric Spike with infectivity and antibody escape of the strains Delta (B.1.617.2) and Omicron (BA.2, BA.5, and BQ.1)

Herein, we conducted simulations of trimeric Spike from several SARS-CoV-2 variants of concern (Delta and Omicron sub-variants BA.2, BA.5, and BQ.1) and investigated the mechanisms by which specific mutations confer resistance to neutralizing antibodies. We observed that the mutations primarily affect the cooperation between protein domains within and between protomers. The substitutions K417N and L452R expand hydrogen bonding interactions, reducing their interaction with neutralizing antibodies. By interacting with nearby residues, the K444T and N460K mutations in the SpikeBQ.1 variant potentially reduces solvent exposure, thereby promoting resistance to antibodies. We also examined the impact of D614G, P681R, and P681H substitutions on Spike protein structure that may be related to infectivity. The D614G substitution influences communication between a glycine residue and neighboring domains, affecting the transition between up- and -down RBD states. The P681R mutation, found in the Delta variant, enhances correlations between protein subunits, while the P681H mutation in Omicron sub-variants weakens long-range interactions that may be associated with reduced fusogenicity. Using a multiple linear regression model, we established a connection between inter-protomer communication and loss of sensitivity to neutralizing antibodies. Our findings underscore the importance of structural communication between protein domains and provide insights into potential mechanisms of immune evasion by SARS-CoV-2. Overall, this study deepens our understanding of how specific mutations impact SARS-CoV-2 infectivity and shed light on how the virus evades the immune system.

Molecular dynamics simulations of the spike trimeric ectodomain of the SARS-CoV-2 Omicron variant: structural relationships with infectivity, evasion to immune system and transmissibility

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) Omicron is currently the most prevalent SARS-CoV-2 variant worldwide. Herein, we calculated molecular dynamics simulations of the trimeric spikeWT and SpikeBA.1 for 300 ns. Our results show that SpikeBA.1 has more conformational flexibility than SpikeWT. Our principal component analysis (PCA) allowed us to observe a broader spectrum of different conformations for SpikeBA.1, mainly at N-terminal domain (NTD) and receptor-binding domain (RBD). Such increased flexibility could contribute to decreased neutralizing antibody recognition of this variant. Our molecular dynamics data show that the RBDBA.1 easily visits an up-conformational state and the prevalent D614G mutation is pivotal to explain molecular dynamics results for this variant because to lost hydrogen bonding interactions between the residue pairs K854SC/D614SC, Y837MC/D614MC, K835SC/D614SC, T859SC/D614SC. In addition, SpikeBA.1 residues near the furin cleavage site are more flexible than in SpikeWT, probably due to P681H and D614G substitutions. Finally, dynamical cross-correlation matrix (DCCM) analysis reveals that D614G and P681H may allosterically affect the cleavage site S1/S2. Conversely, S2' site may be influenced by residues located between NTD and RBD of a neighboring protomer of the SpikeWT. Such communication may be lost in SpikeBA.1, explaining the changes of the cell tropism in the viral infection. In addition, the movements of the NTDWT and NTDBA.1 may modulate the RBD conformation through allosteric effects. Taken together, our results explain how the structural aspects may explain the observed gains in infectivity, immune system evasion and transmissibility of the Omicron variant.Communicated by Ramaswamy H. Sarma.

Performance of SARS COV-2 IgG Anti-N as an Independent Marker of Exposure to SARS COV-2 in an Unvaccinated West African Population

Determination of previous SARS-COV-2 infection is hampered by the absence of a standardized test. The marker used to assess previous exposure is IgG antibody to the nucleocapsid (IgG anti-N), although it is known to wane quickly from peripheral blood. The accuracies of seven antibody tests (virus neutralization test, IgG anti-N, IgG anti-spike [anti-S], IgG anti-receptor binding domain [anti-RBD], IgG anti-N + anti-RBD, IgG anti-N + anti-S, and IgG anti-S + anti-RBD), either singly or in combination, were evaluated on 502 cryopreserved serum samples collected before the COVID-19 vaccination rollout in Kumasi, Ghana. The accuracy of each index test was measured using a composite reference standard based on a combination of neutralization test and IgG anti-N antibody tests. According to the composite reference, 262 participants were previously exposed; the most sensitive test was the virus neutralization test, with 95.4% sensitivity (95% CI: 93.6-97.3), followed by 79.0% for IgG anti-N + anti-S (95% CI: 76.3-83.3). The most specific tests were virus neutralization and IgG anti-N, both with 100% specificity. Viral neutralization and IgG anti-N + anti-S were the overall most accurate tests, with specificity/sensitivity of 100/95.2% and 79.0/92.1%, respectively. Our findings indicate that IgG anti-N alone is an inadequate marker of prior exposure to SARS COV-2 in this population. Virus neutralization assay appears to be the most accurate assay in discerning prior infection. A combination of IgG anti-N and IgG anti-S is also accurate and suited for assessment of SARS COV-2 exposure in low-resource settings.

Drug repurposing and sequence analysis in S-glycoprotein variants reveals critical signature patterns and destabilization of receptor-binding domain in omicron variant

The evolution of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) virus since its emergence in 2019 has yielded several new viral variants with varied infectivity, disease severity, and antigenicity. Although most mutations are expected to be relatively neutral, mutations at the Spike region of the genome have shown to have a major impact on the viral transmission and infection in humans. Therefore, it is crucial to survey the structures of spike protein across the global virus population to contextualize the rate of therapeutic success against these variants. In this study, high-frequency mutational variants from different geographic regions were pooled in order to study the structural evolution of the spike protein through drug docking and MD simulations. We investigated the mutational burden in the spike subregions and have observed that the different variants harbour unique signature patterns in the spike subregions, with certain domains being highly prone to mutations. Further, the MD simulations and docking study revealed that different variants show differential stability when docked for the same set of drug targets. This work sheds light on the mutational burden and the stability landscape of the spike protein across the variants from different geographical regions.Communicated by Ramaswamy H. Sarma.

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.

Host susceptibility and structural and immunological insight of S proteins of two SARS-CoV-2 closely related bat coronaviruses

The bat coronaviruses (CoV) BANAL-20-52 and BANAL-20-236 are two newly identified severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) closely related coronaviruses (SC2r-CoV) and the genome of BANAL-20-52 shares the highest homology with SARS-CoV-2. However, the risk of their potential zoonotic transmission has not been fully evaluated. Here, we determined their potential host susceptibility among 13 different bat species and 26 different animal species, and found that both might have extensive host ranges, indicating high zoonotic transmission potential. We also determined the cryo-EM structures of BANAL-20-52 and BANAL-20-236 S proteins at pH 5.5 and the complex of BANAL-20-236 S1 and Rhinolophus affinis ACE2, and found that both trimeric S proteins adopt all three receptor binding domains (RBDs) in "closed" conformation and are more compact than SARS-CoV-2. Strikingly, the unique sugar moiety at N370 of bat SC2r-CoVs acts like a "bolt" and crosses over two neighboring subunits, facilitating the S proteins in the locked conformation and underpinning the architecture stability. Removal of the glycosylation at N370 by a T372A substitution substantially enhances virus infectivity but becomes highly sensitive to trypsin digestion at pH 5.5, a condition roughly mimicking the insectivorous bat's stomach digestion. In contrast, WT S proteins of SC2r-CoVs showed considerable resistance to trypsin digestion at pH 5.5, indicating that the highly conserved T372 in bat CoVs might result from the selective advantages in stability during the fecal-oral transmission over A372. Moreover, the results of cross-immunogenicity among S proteins of SARS-CoV-2, BANAL-20-52, and BANAL-20-236 showed that A372 pseudoviruses are more sensitive to anti-S sera than T372, indicating that immune evasion might also play a role in the natural selection of T372 over A372 during evolution. Finally, residues 493 and 498 of the S protein affect host susceptibility, and residue 498 also influences the immunogenicity of the S protein. Together, our findings aid a better understanding of the molecular basis of CoV entry, selective evolution, and immunogenicity and highlight the importance of surveillance of susceptible hosts of these viruses to prevent potential outbreaks.

High-throughput identification of prefusion-stabilizing mutations in SARS-CoV-2 spike

Designing prefusion-stabilized SARS-CoV-2 spike is critical for the effectiveness of COVID-19 vaccines. All COVID-19 vaccines in the US encode spike with K986P/V987P mutations to stabilize its prefusion conformation. However, contemporary methods on engineering prefusion-stabilized spike immunogens involve tedious experimental work and heavily rely on structural information. Here, we establish a systematic and unbiased method of identifying mutations that concomitantly improve expression and stabilize the prefusion conformation of the SARS-CoV-2 spike. Our method integrates a fluorescence-based fusion assay, mammalian cell display technology, and deep mutational scanning. As a proof-of-concept, we apply this method to a region in the S2 domain that includes the first heptad repeat and central helix. Our results reveal that besides K986P and V987P, several mutations simultaneously improve expression and significantly lower the fusogenicity of the spike. As prefusion stabilization is a common challenge for viral immunogen design, this work will help accelerate vaccine development against different viruses.

Dissection of Antibody Responses of Gam-COVID-Vac-Vaccinated Subjects Suggests Involvement of Epitopes Outside RBD in SARS-CoV-2 Neutralization

Millions of people have been vaccinated with Gam-COVID-Vac but fine specificities of induced antibodies have not been fully studied. Plasma from 12 naïve and 10 coronavirus disease 2019 (COVID-19) convalescent subjects was obtained before and after two immunizations with Gam-COVID-Vac. Antibody reactivity in the plasma samples (n = 44) was studied on a panel of micro-arrayed recombinant folded and unfolded severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) proteins and 46 peptides spanning the spike protein (S) and by immunoglobulin G (IgG) subclass enzyme-linked immunosorbent assay (ELISA). The ability of Gam-COVID-Vac-induced antibodies to inhibit binding of the receptor-binding domain (RBD) to its receptor angiotensin converting enzyme 2 (ACE2) was investigated in a molecular interaction assay (MIA). The virus-neutralizing capacity of antibodies was studied by the pseudo-typed virus neutralization test (pVNT) for Wuhan-Hu-1 and Omicron. We found that Gam-COVID-Vac vaccination induced significant increases of IgG1 but not of other IgG subclasses against folded S, spike protein subunit 1 (S1), spike protein subunit 2 (S2), and RBD in a comparable manner in naïve and convalescent subjects. Virus neutralization was highly correlated with vaccination-induced antibodies specific for folded RBD and a novel peptide (i.e., peptide 12). Peptide 12 was located close to RBD in the N-terminal part of S1 and may potentially be involved in the transition of the pre- to post-fusion conformation of the spike protein. In summary, Gam-COVID-Vac vaccination induced S-specific IgG1 antibodies in naive and convalescent subjects in a comparable manner. Besides the antibodies specific for RBD, the antibodies induced against a peptide close to the N-terminus of RBD were also associated with virus-neutralization.

SARS-CoV-2 variant biology: immune escape, transmission and fitness

In late 2020, after circulating for almost a year in the human population, severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) exhibited a major step change in its adaptation to humans. These highly mutated forms of SARS-CoV-2 had enhanced rates of transmission relative to previous variants and were termed 'variants of concern' (VOCs). Designated Alpha, Beta, Gamma, Delta and Omicron, the VOCs emerged independently from one another, and in turn each rapidly became dominant, regionally or globally, outcompeting previous variants. The success of each VOC relative to the previously dominant variant was enabled by altered intrinsic functional properties of the virus and, to various degrees, changes to virus antigenicity conferring the ability to evade a primed immune response. The increased virus fitness associated with VOCs is the result of a complex interplay of virus biology in the context of changing human immunity due to both vaccination and prior infection. In this Review, we summarize the literature on the relative transmissibility and antigenicity of SARS-CoV-2 variants, the role of mutations at the furin spike cleavage site and of non-spike proteins, the potential importance of recombination to virus success, and SARS-CoV-2 evolution in the context of T cells, innate immunity and population immunity. SARS-CoV-2 shows a complicated relationship among virus antigenicity, transmission and virulence, which has unpredictable implications for the future trajectory and disease burden of COVID-19.

COVID-19 Biogenesis and Intracellular Transport

SARS-CoV-2 is responsible for the COVID-19 pandemic. The structure of SARS-CoV-2 and most of its proteins of have been deciphered. SARS-CoV-2 enters cells through the endocytic pathway and perforates the endosomes' membranes, and its (+) RNA appears in the cytosol. Then, SARS-CoV-2 starts to use the protein machines of host cells and their membranes for its biogenesis. SARS-CoV-2 generates a replication organelle in the reticulo-vesicular network of the zippered endoplasmic reticulum and double membrane vesicles. Then, viral proteins start to oligomerize and are subjected to budding within the ER exit sites, and its virions are passed through the Golgi complex, where the proteins are subjected to glycosylation and appear in post-Golgi carriers. After their fusion with the plasma membrane, glycosylated virions are secreted into the lumen of airways or (seemingly rarely) into the space between epithelial cells. This review focuses on the biology of SARS-CoV-2's interactions with cells and its transport within cells. Our analysis revealed a significant number of unclear points related to intracellular transport in SARS-CoV-2-infected cells.

Towards Quantum-Chemical Level Calculations of SARS-CoV-2 Spike Protein Variants of Concern by First Principles Density Functional Theory

The spike protein (S-protein) is a crucial part of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), with its many domains responsible for binding, fusion, and host cell entry. In this review we use the density functional theory (DFT) calculations to analyze the atomic-scale interactions and investigate the consequences of mutations in S-protein domains. We specifically describe the key amino acids and functions of each domain, which are essential for structural stability as well as recognition and fusion processes with the host cell; in addition, we speculate on how mutations affect these properties. Such unprecedented large-scale ab initio calculations, with up to 5000 atoms in the system, are based on the novel concept of amino acid-amino acid-bond pair unit (AABPU) that allows for an alternative description of proteins, providing valuable information on partial charge, interatomic bonding and hydrogen bond (HB) formation. In general, our results show that the S-protein mutations for different variants foster an increased positive partial charge, alter the interatomic interactions, and disrupt the HB networks. We conclude by outlining a roadmap for future computational research of biomolecular virus-related systems.

SARS-CoV-2 Spike Protein Mutation at Cysteine-488 Impairs Its Golgi Localization and Intracellular S1/S2 Processing

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) spike protein binds to the cellular receptor-angiotensin-converting enzyme-2 (ACE2) as the first step in viral cell entry. SARS-CoV-2 spike protein expression in the ACE2-expressing cell surface induces cell-cell membrane fusion, thus forming syncytia. To exert its fusogenic activity, the spike protein is typically processed at a specific site (the S1/S2 site) by cellular proteases such as furin. The C488 residue, located at the spike-ACE2 interacting surface, is critical for the fusogenic and infectious roles of the SARS-CoV-2 spike protein. We have demonstrated that the C488 residue of the spike protein is involved in subcellular targeting and S1/S2 processing. C488 mutant spike localization to the Golgi apparatus and cell surface were impaired. Consequently, the S1/S2 processing of the spike protein, probed by anti-Ser-686-cleaved spike antibody, markedly decreased in C488 mutant spike proteins. Moreover, brefeldin-A-mediated endoplasmic-reticulum-to-Golgi traffic suppression also suppressed spike protein S1/S2 processing. As brefeldin A treatment and C488 mutation inhibited S1/S2 processing and syncytia formation, the C488 residue of spike protein is required for functional spike protein processing.

Evolutionary remodelling of N-terminal domain loops fine-tunes SARS-CoV-2 spike

The emergence of SARS-CoV-2 variants has exacerbated the COVID-19 global health crisis. Thus far, all variants carry mutations in the spike glycoprotein, which is a critical determinant of viral transmission being responsible for attachment, receptor engagement and membrane fusion, and an important target of immunity. Variants frequently bear truncations of flexible loops in the N-terminal domain (NTD) of spike; the functional importance of these modifications has remained poorly characterised. We demonstrate that NTD deletions are important for efficient entry by the Alpha and Omicron variants and that this correlates with spike stability. Phylogenetic analysis reveals extensive NTD loop length polymorphisms across the sarbecoviruses, setting an evolutionary precedent for loop remodelling. Guided by these analyses, we demonstrate that variations in NTD loop length, alone, are sufficient to modulate virus entry. We propose that variations in NTD loop length act to fine-tune spike; this may provide a mechanism for SARS-CoV-2 to navigate a complex selection landscape encompassing optimisation of essential functionality, immune-driven antigenic variation and ongoing adaptation to a new host.

High-throughput identification of prefusion-stabilizing mutations in SARS-CoV-2 spike

Designing prefusion-stabilized SARS-CoV-2 spike is critical for the effectiveness of COVID-19 vaccines. All COVID-19 vaccines in the US encode spike with K986P/V987P mutations to stabilize its prefusion conformation. However, contemporary methods on engineering prefusion-stabilized spike immunogens involve tedious experimental work and heavily rely on structural information. Here, we established a systematic and unbiased method of identifying mutations that concomitantly improve expression and stabilize the prefusion conformation of the SARS-CoV-2 spike. Our method integrated a fluorescence-based fusion assay, mammalian cell display technology, and deep mutational scanning. As a proof-of-concept, this method was applied to a region in the S2 domain that includes the first heptad repeat and central helix. Our results revealed that besides K986P and V987P, several mutations simultaneously improved expression and significantly lowered the fusogenicity of the spike. As prefusion stabilization is a common challenge for viral immunogen design, this work will help accelerate vaccine development against different viruses.

Receptor-Binding Domain (RBD) Antibodies Contribute More to SARS-CoV-2 Neutralization When Target Cells Express High Levels of ACE2

Neutralization assays are experimental surrogates for the effectiveness of infection- or vaccine-elicited polyclonal antibodies and therapeutic monoclonal antibodies targeting SARS-CoV-2. However, the measured neutralization can depend on the details of the experimental assay. Here, we systematically assess how ACE2 expression in target cells affects neutralization by antibodies to different spike epitopes in lentivirus pseudovirus neutralization assays. For high ACE2-expressing target cells, receptor-binding domain (RBD) antibodies account for nearly all neutralizing activity in polyclonal human sera. However, for lower ACE2-expressing target cells, antibodies targeting regions outside the RBD make a larger (although still modest) contribution to serum neutralization. These serum-level results are mirrored for monoclonal antibodies: N-terminal domain (NTD) antibodies and RBD antibodies that do not compete for ACE2 binding incompletely neutralize on high ACE2-expressing target cells, but completely neutralize on cells with lower ACE2 expression. Our results show that the ACE2 expression level in the target cells is an important experimental variable, and that high ACE2 expression emphasizes the role of a subset of RBD-directed antibodies.

Receptor binding domain (RBD) antibodies contribute more to SARS-CoV-2 neutralization when target cells express high levels of ACE2.. [PMID: 36093349 PMCID: PMC9460967 DOI: 10.1101/2022.08.29.505713] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

Neutralization assays are experimental surrogates for the effectiveness of infection- or vaccine-elicited polyclonal antibodies and therapeutic monoclonal antibodies targeting SARS-CoV-2. However, the measured neutralization can depend on details of the experimental assay. Here we systematically assess how ACE2 expression in target cells affects neutralization by antibodies to different spike epitopes in lentivirus pseudovirus neutralization assays. For high ACE2-expressing target cells, receptor binding domain (RBD) antibodies account for nearly all neutralizing activity in polyclonal human sera. But for lower ACE2-expressing target cells, antibodies targeting regions outside the RBD make a larger (although still modest) contribution to serum neutralization. These serum-level results are mirrored for monoclonal antibodies: N-terminal domain (NTD) antibodies and RBD antibodies that do not compete for ACE2 binding incompletely neutralize on high ACE2-expressing target cells, but completely neutralize on cells with lower ACE2 expression. Our results show that ACE2 expression level in the target cells is an important experimental variable, and that high ACE2 expression emphasizes the role of a subset of RBD-directed antibodies.