Binding and molecular basis of the bat coronavirus RaTG13 virus to ACE2 in humans and other species

A SARS-CoV and SARS-CoV-2 RBD Heterodimer Vaccine Candidate

The continuous evolution of SARS-CoV-2 through accumulating mutations, combined with the persistent risk of zoonotic sarbecovirus transmission events, highlights the critical demand for broadly protective vaccines. Building on our previous findings that a heterodimeric receptor-binding domain (RBD) design substantially improves cross-reactive immunogenicity in vaccine candidates, we propose this strategy as a foundation for developing pan-sarbecovirus vaccines with cross-neutralizing capacity against diverse and emerging variants. In this study, we developed a sarbecovirus immunogen, utilizing a heterodimeric strategy incorporating the RBDs from both SARS-CoV and SARS-CoV-2. Pseudovirus neutralization assays revealed that mice immunized with the SARS-CoV-2 prototype (PT)-SARS-CoV heterodimer (PT-SARS) developed 39.9- to 305.6-fold higher neutralizing antibody (NAb) titers against SARS-CoV-2 sub-variants compared to the SARS-CoV RBD homodimer (SARS-SARS). Furthermore, PT-SARS elicited 17.6- and 31.2-fold enhanced neutralization against WIV1 and SARS-CoV, respectively, relative to the SARS-CoV-2 PT homodimer (PT-PT). To address evolving Omicron sub-variants, we further updated BA.1-SARS and BA.2-SARS immunogens. Notably, BA.2-SARS exhibited a 6.2-fold increase in neutralizing potency against BA.2.86 compared to PT-SARS. Crucially, the heterodimeric immunogen induced balanced and broadly reactive NAbs against multiple sarbecoviruses, including RaTG13, Pangolin GD, SARS-CoV, and SARS-CoV-2 variants/sub-variants, demonstrating its potential as a sarbecovirus immunogen candidate.

Molecular mechanisms of RaTG13 and SARS-CoV-2 RBD bound to Rhinolophus affinis bat ACE2

The discovery of the RaTG13 coronavirus in Rhinolophus affinis bats in 2013, sharing 96.3% genome homology with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), suggested bats as the origin of SARS-CoV-2. Although both human angiotensin-converting enzyme 2 (hACE2) and Rhinolophus affinis bat ACE2 (bACE2-Ra, seven polymorphic variants named 01-07) are known to serve as entry receptors for these coronaviruses, the binding mechanism of RaTG13 receptor binding domain (RBD) bound to bACE2-Ra remains poorly understood. Here, we found that RaTG13 RBD bound to bACE2-Ra-07 with a weaker affinity (2.42 μM) compared with SARS-CoV-2 RBD (372 nM). Additional glycosylation at residue N370 of RaTG13 had little influence on bACE2-Ra-07 binding by RaTG13 RBD. Crystal structures of the SARS-CoV-2 and RaTG13 N370Q RBD bound to bACE2-Ra-07 were solved. Interface analysis and surface plasmon resonance (SPR) assay indicated that residue substitutions at 493, 498, 501, and 505 may play a more important role in the cross-species recognition of bACE2-Ra-07 by the SARS-CoV-2 RBD. Besides, the N370Q mutation enhanced the binding affinity between the RBD of pangolin coronavirus isolated from Guangxi (PCoV-GX) and the bACE2-Ra-07 receptor by over 10-fold. Furthermore, the recently prevalent SARS-CoV-2 variant RBDs extensively retained the interaction with the bACE2-Ra-07 receptor. Our findings give new lights on the cross-species evolution of SARS-CoV-2 and prompt the urgency to monitor the circulation of coronaviruses in bats to better prevent future spillover.

Multiple independent acquisitions of ACE2 usage in MERS-related coronaviruses

The angiotensin-converting enzyme 2 (ACE2) receptor is shared by various coronaviruses with distinct receptor-binding domain (RBD) architectures, yet our understanding of these convergent acquisition events remains elusive. Here, we report that two bat MERS-related coronaviruses (MERSr-CoVs) infecting Pipistrellus nathusii (P.nat)-MOW15-22 and PnNL2018B-use ACE2 as their receptor, with narrow ortholog specificity. Cryoelectron microscopy structures of the MOW15-22/PnNL2018B RBD-ACE2 complexes unveil an unexpected and entirely distinct binding mode, mapping >45 Å away from that of any other known ACE2-using coronaviruses. Functional profiling of ACE2 orthologs from 105 mammalian species led to the identification of host tropism determinants, including an ACE2 N432-glycosylation restricting viral recognition, and the design of a soluble P.nat ACE2 mutant with potent viral neutralizing activity. Our findings reveal convergent acquisition of ACE2 usage for merbecoviruses found in European bats, underscoring the extraordinary diversity of ACE2 recognition modes among coronaviruses and the promiscuity of this receptor.

SARS-related coronavirus S-protein structures reveal synergistic RBM interactions underpinning high-affinity human ACE2 binding

High-affinity and specific binding toward the human angiotensin-converting enzyme 2 (hACE2) receptor by severe acute respiratory syndrome coronavirus (SARS)-related coronaviruses (SARSr-CoVs) remains incompletely understood. We report cryo-electron microscopy structures of eight different S-proteins from SARSr-CoVs found across Asia, Europe, and Africa. These S-proteins all adopt tightly packed, locked, prefusion conformations. These structures enable the classification of SARSr-CoV S-proteins into three types, based on their receptor-binding motif (RBM) structures and ACE2 binding characteristics. Type-2 S-proteins often preferentially bind bat ACE2 (bACE2) over hACE2. We report a structure of a type-2 BtKY72-RBD in complex with bACE2 to understand ACE2 specificity. Structure-guided mutagenesis of BtKY72-RBD reveals that multiple synergistic mutations in four different regions of RBM are required to achieve high-affinity hACE2 binding. Similar RBM changes can also confer hACE2 binding to another type-2 BM48-31 S-protein, which is primarily non-ACE2 binding. These results provide an understanding of how high-affinity hACE2 binding may be acquired by SARSr-CoV S-proteins.

Rapid preparation and characterization of pan-sarbecovirus mRNA vaccine candidates based on the receptor binding domain

•We developed a rapid approach for small-scale mRNA vaccine preparation against pan-sarbecovirus infection. •An mRNA vaccine developed on this platform effectively elicits the production of neutralizing antibodies in mice. •Our results highlight the potential for rapid and universal vaccine development against emerging coronavirus.

Cross-species recognition of two porcine coronaviruses to their cellular receptor aminopeptidase N of dogs and seven other species

Porcine deltacoronavirus (PDCoV) and transmissible gastroenteritis coronavirus (TGEV), the two causative agents of porcine diarrhea, have been reported to be at risk of cross-species transmission, including to humans. However, the potential host range in which these two CoVs interact remains unclear. We screened 16 animal counterparts for porcine aminopeptidase N (APN), the receptor of PDCoV and TGEV, and found that APNs from eight of 17 animals could bind to the receptor-binding domains (RBDs) of PDCoV and TGEV. Furthermore, the animal APNs that could bind to the RBDs could mediate cellular infection by both viruses. Dog APN (dAPN) has been identified as the animal receptor with the highest capability to mediate the virus infection. We further resolved the complex structures of dAPN bound to the PDCoV RBD/TGEV RBD, respectively, establishing its divergent receptor-binding modes. We identified R325 of dAPN as an important residue in the PDCoV RBD-dAPN interaction, and found the central role of Q746 and T749 in dAPN in the interaction with the TGEV RBD. These findings provide the molecular basis of the potential cross-species transmission of these two porcine CoVs and shed light on future surveillance of these CoVs.

Interfacial subregions of SARS-CoV-2 spike RBD to hACE2 affect intermolecular affinity by their distinct roles played in association and dissociation kinetics

SARS-CoV-2's rapid global transmission depends on spike RBD's strong affinity to hACE2. In the context of binding hot spots well defined, the work investigated how interfacial subregions of SARS-CoV-2 spike RBD to hACE2 affect intermolecular affinity and their potential distinct roles involved in association and dissociation kinetics due to their local structural characteristics. Three spatially consecutive subregions of SARS-CoV-2 RBD were structurally partitioned across RBD's receptor binding motif (RBM). Their impacts on binding affinity and kinetics were differentiated through a comprehensive SPR measurement of hACE2 binding by chimeric swap mutants of respective subdomains from SARS-CoV-2 VOCs & phylogenetically close sarbecoviruses, and further compared with those of included single mutations across RBM and around the RBD core. The data supports that the intermediate interfacial subregion of RBD involving key residue at 417 is the rate-limiting effector of association kinetics and the subregion encompassing residues at 501/498/449 is the key binding energy contributor dictating dissociation kinetics, both of which relate to SARS-CoV-2's adaptive mutational evolution and host tropism closely. The kinetic data and structural analysis of local mutations' impact on spike RBD's binding and thermal stability provide a new perspective in evaluating SARS-CoV-2 evolution and other sarbecoviruses' evolvable binding to hACE2. The inherent binding mode offers direct clues of valid epitope in designing new antibodies that the coronavirus can't elude.

Receptor binding and structural basis of raccoon dog ACE2 binding to SARS-CoV-2 prototype and its variants

Raccoon dog was proposed as a potential host of SARS-CoV-2, but no evidence support such a notion. In our study, we investigated the binding affinities of raccoon dog ACE2 (rdACE2) to the spike (S) protein receptor binding domain (RBD) of SARS-CoV-2 prototype (PT) and its variants. It revealed that the binding affinities of RBD from SARS-CoV-2 variants were generally lower than that of the PT RBD. Through structural and functional analyses, we found amino acids H34 and M82 play pivotal roles in maintaining the binding affinity of ACE2 to different SARS-CoV-2 sub-variants. These results suggest that raccoon dogs exhibit lower susceptibility to SARS-CoV-2 compared to those animal species with a high prevalence of SARS-CoV-2 transmission.

Establishment of a human organoid-based evaluation system for assessing interspecies infection risk of animal-borne coronaviruses

The COVID-19 pandemic presents a major threat to global public health. Several lines of evidence have shown that the Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2), along with two other highly pathogenic coronaviruses, SARS-CoV and Middle East Respiratory Syndrome (MERS-CoV) originated from bats. To prevent and control future coronavirus outbreaks, it is necessary to investigate the interspecies infection and pathogenicity risks of animal-related coronaviruses. Currently used infection models, including in vitro cell lines and in vivo animal models, fail to fully mimic the primary infection in human tissues. Here, we employed organoid technology as a promising new model for studying emerging pathogens and their pathogenic mechanisms. We investigated the key host-virus interaction patterns of five human coronaviruses (SARS-CoV-2 original strain, Omicron BA.1, MERS-CoV, HCoV-229E, and HCoV-OC43) in different human respiratory organoids. Five indicators, including cell tropism, invasion preference, replication activity, host response and virus-induced cell death, were developed to establish a comprehensive evaluation system to predict coronavirus interspecies infection and pathogenicity risks. Using this system, we further examined the pathogenicity and interspecies infection risks of three SARS-related coronaviruses (SARSr-CoV), including WIV1 and rRsSHC014S from bats, and MpCoV-GX from pangolins. Moreover, we found that cannabidiol, a non-psychoactive plant extract, exhibits significant inhibitory effects on various coronaviruses in human lung organoid. Cannabidiol significantly enhanced interferon-stimulated gene expression but reduced levels of inflammatory cytokines. In summary, our study established a reliable comprehensive evaluation system to analyse infection and pathogenicity patterns of zoonotic coronaviruses, which could aid in prevention and control of potentially emerging coronavirus diseases.

Cellular sialoglycans are differentially required for endosomal and cell-surface entry of SARS-CoV-2 in lung cell lines

Cell entry of severe acute respiratory coronavirus-2 (SARS-CoV-2) and other CoVs can occur via two distinct routes. Following receptor binding by the spike glycoprotein, membrane fusion can be triggered by spike cleavage either at the cell surface in a transmembrane serine protease 2 (TMPRSS2)-dependent manner or within endosomes in a cathepsin-dependent manner. Cellular sialoglycans have been proposed to aid in CoV attachment and entry, although their functional contributions to each entry pathway are unknown. In this study, we used genetic and enzymatic approaches to deplete sialic acid from cell surfaces and compared the requirement for sialoglycans during endosomal and cell-surface CoV entry using lentiviral particles pseudotyped with the spike proteins of different sarbecoviruses. We show that entry of SARS-CoV-1, WIV1-CoV and WIV16-CoV, like the SARS-CoV-2 omicron variant, depends on endosomal cathepsins and requires cellular sialoglycans for entry. Ancestral SARS-CoV-2 and the delta variant can use either pathway for entry, but only require sialic acid for endosomal entry in cells lacking TMPRSS2. Binding of SARS-CoV-2 spike protein to cells did not require sialic acid, nor was sialic acid required for SARS-CoV-2 entry in TMRPSS2-expressing cells. These findings suggest that cellular sialoglycans are not strictly required for SARS-CoV-2 attachment, receptor binding or fusion, but rather promote endocytic entry of SARS-CoV-2 and related sarbecoviruses. In contrast, the requirement for sialic acid during entry of MERS-CoV pseudoparticles and authentic HCoV-OC43 was not affected by TMPRSS2 expression, consistent with a described role for sialic acid in merbecovirus and embecovirus cell attachment. Overall, these findings clarify the role of sialoglycans in SARS-CoV-2 entry and suggest that cellular sialoglycans mediate endosomal, but not cell-surface, SARS-CoV-2 entry.

Distinct pathways for evolution of enhanced receptor binding and cell entry in SARS-like bat coronaviruses

Understanding the zoonotic risks posed by bat coronaviruses (CoVs) is critical for pandemic preparedness. Herein, we generated recombinant vesicular stomatitis viruses (rVSVs) bearing spikes from divergent bat CoVs to investigate their cell entry mechanisms. Unexpectedly, the successful recovery of rVSVs bearing the spike from SHC014-CoV, a SARS-like bat CoV, was associated with the acquisition of a novel substitution in the S2 fusion peptide-proximal region (FPPR). This substitution enhanced viral entry in both VSV and coronavirus contexts by increasing the availability of the spike receptor-binding domain to recognize its cellular receptor, ACE2. A second substitution in the S1 N-terminal domain, uncovered through the rescue and serial passage of a virus bearing the FPPR substitution, further enhanced spike:ACE2 interaction and viral entry. Our findings identify genetic pathways for adaptation by bat CoVs during spillover and host-to-host transmission, fitness trade-offs inherent to these pathways, and potential Achilles' heels that could be targeted with countermeasures.

Sarbecovirus RBD indels and specific residues dictating multi-species ACE2 adaptiveness

Our comprehensive understanding of the multi-species ACE2 adaptiveness of sarbecoviruses remains elusive, particularly for those with various receptor binding motif (RBM) insertions/deletions (indels). Here, we analyzed RBM sequences from 268 sarbecoviruses categorized into four RBM indel types. We examined the ability of 20 representative sarbecovirus Spike glycoproteins (S) and derivatives in utilizing ACE2 from various bats and several other mammalian species. We reveal that sarbecoviruses with long RBMs (type-I) can achieve broad ACE2 tropism, whereas viruses with single deletions in Region 1 (type-II) or Region 2 (type-III) exhibit narrower ACE2 tropism. Sarbecoviruses with double region deletions (type-IV) completely lost ACE2 usage, which is restricted by clade-specific residues within and outside RBM. Lastly, we propose the evolution of sarbecovirus RBM indels and illustrate how loop lengths, disulfide, and residue determinants shape multi-species ACE2 adaptiveness. This study provides profound insights into the mechanisms governing ACE2 usage and spillover risks of sarbecoviruses.

Cross-Species Susceptibility of Emerging Variants of SARS-CoV-2 Spike

BACKGROUND

The continuous evolution of SARS-CoV-2 and the emergence of novel variants with numerous mutations have heightened concerns surrounding the possibility of cross-species transmission and the establishment of natural animal reservoirs for the virus, but the host range of emerging SARS-CoV-2 variants has not been fully explored yet.

METHODS

We employed an in vitro model comprising VSV∆G* pseudotyped viruses bearing SARS-CoV-2 spike proteins to explore the plausible host range of SARS-CoV-2 emerging variants.

RESULTS

The overall host tropism of emerging SARS-CoV-2 variants are consistent with that of the SARS-CoV-2 wuhan-hu-1 strain with minor difference. Pseudotyped viruses bearing spike protein from RaTG13 and RmYN02 can enter cell cultures from a broad range of mammalian species, revealing that mink and hamsters may act as potential intermediate hosts. We further investigated 95 potential site-specific mutations in the SARS-CoV-2 spike protein that could impact viral infectivity across different species. The results showed that 13 of these mutations notably increased the transduction rates by more than two-fold when compared to the wild-type spike protein. Further examination of these 13 mutations within cell cultures from 31 different species revealed heightened sensitivity in cells derived from palm civets, minks, and Chinese horseshoe bats to the VSV∆G*-SARS2-S mutants. Specific mutations, such as L24F, R158G, and L212I, were seen to significantly enhance the capacity for SARS-CoV-2 of cross-species transmission.

CONCLUSIONS

This study offers critical insights for the ongoing surveillance and monitoring efforts of SARS-CoV-2 evolution, emphasizing the need for the vigilant monitoring of specific mutations in both human and animal populations.

Structural characteristics of BtKY72 RBD bound to bat ACE2 reveal multiple key residues affecting ACE2 usage of sarbecoviruses

Two different sarbecoviruses, severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2, have caused serious challenges to public health. Certain sarbecoviruses utilize angiotensin-converting enzyme 2 (ACE2) as their cellular receptor, whereas some do not, speculatively due to the two deletions in their receptor-binding domain (RBD). However, it remains unclear whether sarbecoviruses with one deletion in the RBD can still bind to ACE2. Here, we showed that two phylogenetically related sarbecoviruses with one deletion, BtKY72 and BM48-31, displayed a different ACE2-usage range. The cryo-electron microscopy structure of BtKY72 RBD bound to bat ACE2 identified a key residue important for the interaction between RBD and ACE2. In addition, we demonstrated that the mutations involving four types of core residues enabled the sarbecoviruses with deletion(s) to bind to human ACE2 (hACE2) and broadened the ACE2 usage of SARS-CoV-2. Our findings help predict the potential hACE2-binding ability to emerge sarbecoviruses and develop pan-sarbecovirus therapeutic agents.

IMPORTANCE

Many sarbecoviruses, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), possess the ability to bind to receptor angiotensin-converting enzyme 2 (ACE2) through their receptor-binding domain (RBD). However, certain sarbecoviruses with deletion(s) in the RBD lack this capability. In this study, we investigated two closely related short-deletion sarbecoviruses, BtKY72 and BM48-31, and revealed that BtKY72 exhibited a broader ACE2-binding spectrum compared to BM48-31. Structural analysis of the BtKY72 RBD-bat ACE2 complex identifies a critical residue at position 493 contributing to these differences. Furthermore, we demonstrated that the mutations involving four core residues in the RBD enabled the sarbecoviruses with deletion(s) to bind to human ACE2 and expanded the ACE2 usage spectra of SARS-CoV-2. These findings offer crucial insights for accurately predicting the potential threat of newly emerging sarbecoviruses to human health.

Resistance to SARS-CoV-2 infection in camelid nasal organoids is associated with lack of ACE2 expression

The severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infects a variety of animal species. Susceptibility to SARS-CoV-2 is primarily determined by the utilization of the viral receptor, ACE2. SARS-CoV-2 can utilize a broad range of animal ACE2 isoforms in vitro, including the ACE2 from various camelid species. However, experimental infection of these animals does not lead to productive infection or seroconversion. In this study, we investigate the susceptibility of camelids to SARS-CoV-2 using novel well-differentiated camelid nasal organoids. We show that camelid nasal organoids are highly susceptible to Middle East respiratory syndrome coronavirus (MERS-CoV) infection, but not to infection with different SARS-CoV-2 variants (614G, BA.1 or EG.5.1.1). All viruses efficiently infected human airway organoids. Immunohistochemistry analysis revealed the absence of ACE2 on camelid nasal organoids and dromedary camel upper respiratory tract. In contrast, DPP4 was expressed in both camelid nasal organoids and the camel upper respiratory tract, which correlates with MERS-CoV infection. This study indicates that the camelid upper respiratory tract lacks expression of ACE2, which is associated with resistance to SARS-CoV-2 infection.

Spike structures, receptor binding, and immune escape of recently circulating SARS-CoV-2 Omicron BA.2.86, JN.1, EG.5, EG.5.1, and HV.1 sub-variants

The recently emerged BA.2.86, JN.1, EG.5, EG.5.1, and HV.1 variants have a growth advantage. In this study, we explore the structural bases of receptor binding and immune evasion for the Omicron BA.2.86, JN.1, EG.5, EG.5.1, and HV.1 sub-variants. Our findings reveal that BA.2.86 exhibits strong receptor binding, whereas its JN.1 sub-lineage displays a decreased binding affinity to human ACE2 (hACE2). Through complex structure analyses, we observed that the reversion of R493Q in BA.2.86 receptor binding domain (RBD) plays a facilitating role in receptor binding, while the L455S substitution in JN.1 RBD restores optimal affinity. Furthermore, the structure of monoclonal antibody (mAb) S309 complexed with BA.2.86 RBD highlights the importance of the K356T mutation, which brings a new N-glycosylation motif, altering the binding pattern of mAbs belonging to RBD-5 represented by S309. These findings emphasize the importance of closely monitoring BA.2.86 and its sub-lineages to prevent another wave of SARS-CoV-2 infections.

Structural basis for mouse receptor recognition by bat SARS2-like coronaviruses

The animal origin of SARS-CoV-2 remains elusive, lacking a plausible evolutionary narrative that may account for its emergence. Its spike protein resembles certain segments of BANAL-236 and RaTG13, two bat coronaviruses considered possible progenitors of SARS-CoV-2. Additionally, its spike contains a furin motif, a common feature of rodent coronaviruses. To explore the possible involvement of rodents in the emergence of SARS-CoV-2 spike, we examined the crystal structures of the spike receptor-binding domains (RBDs) of BANAL-236 and RaTG13 each complexed with mouse receptor ACE2. Both RBDs have residues at positions 493 and 498 that align well with two virus-binding hotspots on mouse ACE2. Our biochemical evidence supports that both BANAL-236 and RaTG13 spikes can use mouse ACE2 as their entry receptor. These findings point to a scenario in which these bat coronaviruses may have coinfected rodents, leading to a recombination of their spike genes and a subsequent acquisition of a furin motif in rodents, culminating in the emergence of SARS-CoV-2.

The binding and structural basis of fox ACE2 to RBDs from different sarbecoviruses

Foxes are susceptible to SARS-CoV-2 in laboratory settings, and there have also been reports of natural infections of both SARS-CoV and SARS-CoV-2 in foxes. In this study, we assessed the binding capacities of fox ACE2 to important sarbecoviruses, including SARS-CoV, SARS-CoV-2, and animal-origin SARS-CoV-2 related viruses. Our findings demonstrated that fox ACE2 exhibits broad binding capabilities to receptor-binding domains (RBDs) of sarbecoviruses. We further determined the cryo-EM structures of fox ACE2 complexed with RBDs of SARS-CoV, SARS-CoV-2 prototype (PT), and Omicron BF.7. Through structural analysis, we identified that the K417 mutation can weaken the ability of SARS-CoV-2 sub-variants to bind to fox ACE2, thereby reducing the susceptibility of foxes to SARS-CoV-2 sub-variants. In addition, the Y498 residue in the SARS-CoV RBD plays a crucial role in forming a vital cation-π interaction with K353 in the fox ACE2 receptor. This interaction is the primary determinant for the higher affinity of the SARS-CoV RBD compared to that of the SARS-CoV-2 PT RBD. These results indicate that foxes serve as potential hosts for numerous sarbecoviruses, highlighting the critical importance of surveillance efforts.

Distinct pathway for evolution of enhanced receptor binding and cell entry in SARS-like bat coronaviruses

Understanding the zoonotic risks posed by bat coronaviruses (CoVs) is critical for pandemic preparedness. Herein, we generated recombinant vesicular stomatitis viruses (rVSVs) bearing spikes from divergent bat CoVs to investigate their cell entry mechanisms. Unexpectedly, the successful recovery of rVSVs bearing the spike from SHC014, a SARS-like bat CoV, was associated with the acquisition of a novel substitution in the S2 fusion peptide-proximal region (FPPR). This substitution enhanced viral entry in both VSV and coronavirus contexts by increasing the availability of the spike receptor-binding domain to recognize its cellular receptor, ACE2. A second substitution in the spike N-terminal domain, uncovered through forward-genetic selection, interacted epistatically with the FPPR substitution to synergistically enhance spike:ACE2 interaction and viral entry. Our findings identify genetic pathways for adaptation by bat CoVs during spillover and host-to-host transmission, fitness trade-offs inherent to these pathways, and potential Achilles' heels that could be targeted with countermeasures.

Variation in structural motifs within SARS-related coronavirus spike proteins

SARS-CoV-2 is the third known coronavirus (CoV) that has crossed the animal-human barrier in the last two decades. However, little structural information exists related to the close genetic species within the SARS-related coronaviruses. Here, we present three novel SARS-related CoV spike protein structures solved by single particle cryo-electron microscopy analysis derived from bat (bat SL-CoV WIV1) and civet (cCoV-SZ3, cCoV-007) hosts. We report complex glycan trees that decorate the glycoproteins and density for water molecules which facilitated modeling of the water molecule coordination networks within structurally important regions. We note structural conservation of the fatty acid binding pocket and presence of a linoleic acid molecule which are associated with stabilization of the receptor binding domains in the "down" conformation. Additionally, the N-terminal biliverdin binding pocket is occupied by a density in all the structures. Finally, we analyzed structural differences in a loop of the receptor binding motif between coronaviruses known to infect humans and the animal coronaviruses described in this study, which regulate binding to the human angiotensin converting enzyme 2 receptor. This study offers a structural framework to evaluate the close relatives of SARS-CoV-2, the ability to inform pandemic prevention, and aid in the development of pan-neutralizing treatments.

Pseudotyped virus infection of multiplexed ACE2 libraries reveals SARS-CoV-2 variant shifts in receptor usage

Pairwise compatibility between virus and host proteins can dictate the outcome of infection. During transmission, both inter- and intraspecies variabilities in receptor protein sequences can impact cell susceptibility. Many viruses possess mutable viral entry proteins and the patterns of host compatibility can shift as the viral protein sequence changes. This combinatorial sequence space between virus and host is poorly understood, as traditional experimental approaches lack the throughput to simultaneously test all possible combinations of protein sequences. Here, we created a pseudotyped virus infection assay where a multiplexed target-cell library of host receptor variants can be assayed simultaneously using a DNA barcode sequencing readout. We applied this assay to test a panel of 30 ACE2 orthologs or human sequence mutants for infectability by the original SARS-CoV-2 spike protein or the Alpha, Beta, Gamma, Delta, and Omicron BA1 variant spikes. We compared these results to an analysis of the structural shifts that occurred for each variant spike's interface with human ACE2. Mutated residues were directly involved in the largest shifts, although there were also widespread indirect effects altering interface structure. The N501Y substitution in spike conferred a large structural shift for interaction with ACE2, which was partially recreated by indirect distal substitutions in Delta, which does not harbor N501Y. The structural shifts from N501Y greatly influenced the set of animal orthologs the variant spike was capable of interacting with. Out of the thirteen non-human orthologs, ten exhibited unique patterns of variant-specific compatibility, demonstrating that spike sequence changes during human transmission can toggle ACE2 compatibility and potential susceptibility of other animal species, and cumulatively increase overall compatibilities as new variants emerge. These experiments provide a blueprint for similar large-scale assessments of protein compatibility during entry by diverse viruses. This dataset demonstrates the complex compatibility relationships that occur between variable interacting host and virus proteins.

Unveiling hidden structural patterns in the SARS-CoV-2 genome: Computational insights and comparative analysis

SARS-CoV-2, the causative agent of COVID-19, is known to exhibit secondary structures in its 5' and 3' untranslated regions, along with the frameshifting stimulatory element situated between ORF1a and 1b. To identify additional regions containing conserved structures, we utilized a multiple sequence alignment with related coronaviruses as a starting point. We applied a computational pipeline developed for identifying non-coding RNA elements. Our pipeline employed three different RNA structural prediction approaches. We identified forty genomic regions likely to harbor structures, with ten of them showing three-way consensus substructure predictions among our predictive utilities. We conducted intracomparisons of the predictive utilities within the pipeline and intercomparisons with four previously published SARS-CoV-2 structural datasets. While there was limited agreement on the precise structure, different approaches seemed to converge on regions likely to contain structures in the viral genome. By comparing and combining various computational approaches, we can predict regions most likely to form structures, as well as a probable structure or ensemble of structures. These predictions can be used to guide surveillance, prophylactic measures, or therapeutic efforts. Data and scripts employed in this study may be found at https://doi.org/10.5281/zenodo.8298680.

Key mechanistic features of the trade-off between antibody escape and host cell binding in the SARS-CoV-2 Omicron variant spike proteins

Since SARS-CoV-2 Omicron variant emerged, it is constantly evolving into multiple sub-variants, including BF.7, BQ.1, BQ.1.1, XBB, XBB.1.5 and the recently emerged BA.2.86 and JN.1. Receptor binding and immune evasion are recognized as two major drivers for evolution of the receptor binding domain (RBD) of the SARS-CoV-2 spike (S) protein. However, the underlying mechanism of interplay between two factors remains incompletely understood. Herein, we determined the structures of human ACE2 complexed with BF.7, BQ.1, BQ.1.1, XBB and XBB.1.5 RBDs. Based on the ACE2/RBD structures of these sub-variants and a comparison with the known complex structures, we found that R346T substitution in the RBD enhanced ACE2 binding upon an interaction with the residue R493, but not Q493, via a mechanism involving long-range conformation changes. Furthermore, we found that R493Q and F486V exert a balanced impact, through which immune evasion capability was somewhat compromised to achieve an optimal receptor binding. We propose a "two-steps-forward and one-step-backward" model to describe such a compromise between receptor binding affinity and immune evasion during RBD evolution of Omicron sub-variants.

Rapid Degradation of the Human ACE2 Receptor Upon Binding and Internalization of SARS-Cov-2-Spike-RBD Protein

It is widely accepted that the SARS-CoV-2 betacoronavirus infects humans through binding the human Angiotensin Receptor 2 (ACE2) that lines the nasal cavity and lungs, followed by import into a cell utilizing the Transmembrane Protease, Serine 2 (TMPRSS2) cofactor. ACE2 binding is mediated by an approximately 200-residue portion of the SARS-CoV-2 extracellular spike protein, the receptor binding domain (RBD). Robust interactions are shown using a novel cell-based assay between an RBD membrane tethered-GFP fusion protein and the membrane bound ACE2-Cherry fusion protein. Several observations were not predicted including, quick and sustained interactions leading to internalization of RBD fusion protein into the ACE2 cells and rapid downregulation of the ACE2-Cherry fluorescence. Targeted mutation in the RBD disulfide Loop 4 led to a loss of internalization for several variants tested. However, a secreted RBD did not cause ACE2 downregulation of ACE2-Cherry fluorescence. Thus, the membrane associated form of RBD found on the viral coat may have long-term system wide consequences on ACE2 expressing cells.

Pseudotyped virus infection of multiplexed ACE2 libraries reveals SARS-CoV-2 variant shifts in receptor usage

Pairwise compatibility between virus and host proteins can dictate the outcome of infection. During transmission, both inter- and intraspecies variabilities in receptor protein sequences can impact cell susceptibility. Many viruses possess mutable viral entry proteins and the patterns of host compatibility can shift as the viral protein sequence changes. This combinatorial sequence space between virus and host is poorly understood, as traditional experimental approaches lack the throughput to simultaneously test all possible combinations of protein sequences. Here, we created a pseudotyped virus infection assay where a multiplexed target-cell library of host receptor variants can be assayed simultaneously using a DNA barcode sequencing readout. We applied this assay to test a panel of 30 ACE2 orthologs or human sequence mutants for infectability by the original SARS-CoV-2 spike protein or the Alpha, Beta, Gamma, Delta, and Omicron BA1 variant spikes. We compared these results to an analysis of the structural shifts that occurred for each variant spike's interface with human ACE2. Mutated residues were directly involved in the largest shifts, although there were also widespread indirect effects altering interface structure. The N501Y substitution in spike conferred a large structural shift for interaction with ACE2, which was partially recreated by indirect distal substitutions in Delta, which does not harbor N501Y. The structural shifts from N501Y greatly influenced the set of animal orthologs the variant spike was capable of interacting with. Out of the thirteen non-human orthologs, ten exhibited unique patterns of variant-specific compatibility, demonstrating that spike sequence changes during human transmission can toggle ACE2 compatibility and potential susceptibility of other animal species, and cumulatively increase overall compatibilities as new variants emerge. These experiments provide a blueprint for similar large-scale assessments of protein compatibility during entry by diverse viruses. This dataset demonstrates the complex compatibility relationships that occur between variable interacting host and virus proteins.

Emergence of SARS and COVID-19 and preparedness for the next emerging disease X

Severe acute respiratory syndrome (SARS) and Coronavirus disease 2019 (COVID-19) are two human Coronavirus diseases emerging in this century, posing tremendous threats to public health and causing great loss to lives and economy. In this review, we retrospect the studies tracing the molecular evolution of SARS-CoV, and we sort out current research findings about the potential ancestor of SARS-CoV-2. Updated knowledge about SARS-CoV-2-like viruses found in wildlife, the animal susceptibility to SARS-CoV-2, as well as the interspecies transmission risk of SARS-related coronaviruses (SARSr-CoVs) are gathered here. Finally, we discuss the strategies of how to be prepared against future outbreaks of emerging or re-emerging coronaviruses.

Recombination-aware phylogenetic analysis sheds light on the evolutionary origin of SARS-CoV-2

SARS-CoV-2 can infect human cells through the recognition of the human angiotensin-converting enzyme 2 receptor. This affinity is given by six amino acid residues located in the variable loop of the receptor binding domain (RBD) within the Spike protein. Genetic recombination involving bat and pangolin Sarbecoviruses, and natural selection have been proposed as possible explanations for the acquisition of the variable loop and these amino acid residues. In this study we employed Bayesian phylogenetics to jointly reconstruct the phylogeny of the RBD among human, bat and pangolin Sarbecoviruses and detect recombination events affecting this region of the genome. A recombination event involving RaTG13, the closest relative of SARS-CoV-2 that lacks five of the six residues, and an unsampled Sarbecovirus lineage was detected. This result suggests that the variable loop of the RBD didn't have a recombinant origin and the key amino acid residues were likely present in the common ancestor of SARS-CoV-2 and RaTG13, with the latter losing five of them probably as the result of recombination.

ACE2-Fc and DPP4-Fc decoy receptors against SARS-CoV-2 and MERS-CoV variants: a quick therapeutic option for current and future coronaviruses outbreaks

The severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) and the Middle East respiratory syndrome coronavirus (MERS-CoV) are highly pathogenic human coronaviruses (CoVs). Anti-CoVs mAbs and vaccines may be effective, but the emergence of neutralization escape variants is inevitable. Angiotensin-converting enzyme 2 and dipeptidyl peptidase 4 enzyme are the getaway receptors for SARS-CoV-2 and MERS-CoV, respectively. Thus, we reformatted these receptors as Fc-fusion decoy receptors. Then, we tested them in parallel with anti-SARS-CoV (ab1-IgG) and anti-MERS-CoV (M336-IgG) mAbs against several variants using pseudovirus neutralization assay. The generated Fc-based decoy receptors exhibited a strong inhibitory effect against all pseudotyped CoVs. Results showed that although mAbs can be effective antiviral drugs, they might rapidly lose their efficacy against highly mutated viruses. We suggest that receptor traps can be engineered as Fc-fusion proteins for highly mutating viruses with known entry receptors, for a faster and effective therapeutic response even against virus harboring antibodies escape mutations.

Cross-species recognition of bat coronavirus RsYN04 and cross-reaction of SARS-CoV-2 antibodies against the virus

Following the outbreak of coronavirus disease 2019 (COVID-19), several severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2)-related coronaviruses have been discovered. Previous research has identified a novel lineage of SARS-CoV-2-related CoVs in bats, including RsYN04, which recognizes human angiotensin-converting enzyme 2 (ACE2) and thus poses a potential threat to humans. Here, we screened the binding of the RsYN04 receptor-binding domain (RBD) to ACE2 orthologs from 52 animal species and found that the virus showed a narrower ACE2-binding spectrum than SARS-CoV-2. However, the presence of the T484W mutation in the RsYN04 RBD broadened its range. We also evaluated 44 SARS-CoV-2 antibodies targeting seven epitope communities in the SARS-CoV-2 RBD, together with serum obtained from COVID-19 convalescents and vaccinees, to determine their cross-reaction against RsYN04. Results showed that no antibodies, except for the RBD-6 and RBD-7 classes, bound to the RsYN04 RBD, indicating substantial immune differences from SARS-CoV-2. Furthermore, the structure of the RsYN04 RBD in complex with cross-reactive antibody S43 in RBD-7 revealed a potently broad epitope for the development of therapeutics and vaccines. Our findings suggest RsYN04 and other viruses belonging to the same clade have the potential to infect several species, including humans, highlighting the necessity for viral surveillance and development of broad anti-coronavirus countermeasures.

Determination of the factors responsible for the tropism of SARS-CoV-2-related bat coronaviruses to Rhinolophus bat ACE2

IMPORTANCE

The efficiency of infection receptor use is the first step in determining the species tropism of viruses. After the coronavirus disease 2019 pandemic, a number of SARS-CoV-2-related coronaviruses (SC2r-CoVs) were identified in Rhinolophus bats, and some of them can use human angiotensin converting enzyme 2 (ACE2) for the infection receptor without acquiring additional mutations. This means that the potential of certain SC2r-CoVs to cause spillover from bats to humans is "off-the-shelf." However, both SC2r-CoVs and Rhinolophus bat species are highly diversified, and the host tropism of SC2r-CoVs remains unclear. Here, we focus on two Laotian SC2r-CoVs, BANAL-20-236 and BANAL-20-52, and determine how the tropism of SC2r-CoVs to Rhinolophus bat ACE2 is determined at the amino acid resolution level.

Determinants of species-specific utilization of ACE2 by human and animal coronaviruses

Utilization of human ACE2 allowed several bat coronaviruses (CoVs), including the causative agent of COVID-19, to infect humans directly or via intermediate hosts. However, the determinants of species-specific differences in ACE2 usage and the frequency of the ability of animal CoVs to use human ACE2 are poorly understood. Here we applied VSV pseudoviruses to analyze the ability of Spike proteins from 26 human or animal CoVs to use ACE2 receptors across nine reservoir, potential intermediate and human hosts. We show that SARS-CoV-2 Omicron variants evolved towards more efficient ACE2 usage but mutation of R493Q in BA.4/5 and XBB Spike proteins disrupts utilization of ACE2 from Greater horseshoe bats. Variations in ACE2 residues 31, 41 and 354 govern species-specific differences in usage by coronaviral Spike proteins. Mutation of T403R allows the RaTG13 bat CoV Spike to efficiently use all ACE2 orthologs for viral entry. Sera from COVID-19 vaccinated individuals neutralize the Spike proteins of various bat Sarbecoviruses. Our results define determinants of ACE2 receptor usage of diverse CoVs and suggest that COVID-19 vaccination may protect against future zoonoses of bat coronaviruses.

A pan-coronavirus peptide inhibitor prevents SARS-CoV-2 infection in mice by intranasal delivery

Coronaviruses (CoVs) have brought serious threats to humans, particularly severe acute respiratory syndrome Coronavirus 2 (SARS-CoV-2), which continually evolves into multiple variants. These variants, especially Omicron, reportedly escape therapeutic antibodies and vaccines, indicating an urgent need for new antivirals with pan-SARS-CoV-2 inhibitory activity. We previously reported that a peptide fusion inhibitor, P3, targeting heptad repeated-1 (HR1) of SARS-CoV-2 spike (S) protein, could inhibit viral infections. Here, we further designed multiple derivatives of the P3 based on structural analysis and found that one derivative, the P315V3, showed the most efficient antiviral activity against SARS-CoV-2 variants and several other sarbecoviruses, as well as other human-CoVs (HCoVs). P315V3 also exhibited effective prophylactic efficacy against the SARS-CoV-2 Delta and Omicron variants in mice via intranasal administration. These results suggest that P315V3, which is in Phase II clinical trial, is promising for further development as a nasal pan-SARS-CoV-2 or pan-CoVs inhibitor to prevent or treat CoV diseases.

Host range, transmissibility and antigenicity of a pangolin coronavirus

The pathogenic and cross-species transmission potential of SARS-CoV-2-related coronaviruses (CoVs) remain poorly characterized. Here we recovered a wild-type pangolin (Pg) CoV GD strain including derivatives encoding reporter genes using reverse genetics. In primary human cells, PgCoV replicated efficiently but with reduced fitness and showed less efficient transmission via airborne route compared with SARS-CoV-2 in hamsters. PgCoV was potently inhibited by US Food and Drug Administration approved drugs, and neutralized by COVID-19 patient sera and SARS-CoV-2 therapeutic antibodies in vitro. A pan-Sarbecovirus antibody and SARS-CoV-2 S2P recombinant protein vaccine protected BALB/c mice from PgCoV infection. In K18-hACE2 mice, PgCoV infection caused severe clinical disease, but mice were protected by a SARS-CoV-2 human antibody. Efficient PgCoV replication in primary human cells and hACE2 mice, coupled with a capacity for airborne spread, highlights an emergence potential. However, low competitive fitness, pre-immune humans and the benefit of COVID-19 countermeasures should impede its ability to spread globally in human populations.

Structural basis of white-tailed deer, Odocoileus virginianus, ACE2 recognizing all the SARS-CoV-2 variants of concern with high affinity

SARS-CoV-2 has been expanding its host range, among which the white-tailed deer (WTD), Odocoileus virginianus, became the first wildlife species infected on a large scale and might serve as a host reservoir for variants of concern (VOCs) in case no longer circulating in humans. In this study, we comprehensively assessed the binding of the WTD angiotensin-converting enzyme 2 (ACE2) receptor to the spike (S) receptor-binding domains (RBDs) from the SARS-CoV-2 prototype (PT) strain and multiple variants. We found that WTD ACE2 could be broadly recognized by all of the tested RBDs. We further determined the complex structures of WTD ACE2 with PT, Omicron BA.1, and BA.4/5 S trimer. Detailed structural comparison revealed the important roles of RBD residues on 486, 498, and 501 sites for WTD ACE2 binding. This study deepens our understanding of the interspecies transmission mechanisms of SARS-CoV-2 and further addresses the importance of constant monitoring on SARS-CoV-2 infections in wild animals. IMPORTANCE Even if we manage to eliminate the virus among humans, it will still circulate among wildlife and continuously be transmitted back to humans. A recent study indicated that WTD may serve as reservoir for nearly extinct SARS-CoV-2 strains. Therefore, it is critical to evaluate the binding abilities of SARS-CoV-2 variants to the WTD ACE2 receptor and elucidate the molecular mechanisms of binding of the RBDs to assess the risk of spillback events.

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.

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.

Structural basis for receptor binding and broader interspecies receptor recognition of currently circulating Omicron sub-variants

Multiple SARS-CoV-2 Omicron sub-variants, such as BA.2, BA.2.12.1, BA.4, and BA.5, emerge one after another. BA.5 has become the dominant strain worldwide. Additionally, BA.2.75 is significantly increasing in some countries. Exploring their receptor binding and interspecies transmission risk is urgently needed. Herein, we examine the binding capacities of human and other 28 animal ACE2 orthologs covering nine orders towards S proteins of these sub-variants. The binding affinities between hACE2 and these sub-variants remain in the range as that of previous variants of concerns (VOCs) or interests (VOIs). Notably, R493Q reverse mutation enhances the bindings towards ACE2s from humans and many animals closely related to human life, suggesting an increased risk of cross-species transmission. Structures of S/hACE2 or RBD/hACE2 complexes for these sub-variants and BA.2 S binding to ACE2 of mouse, rat or golden hamster are determined to reveal the molecular basis for receptor binding and broader interspecies recognition.

Broad Sarbecovirus Neutralizing Antibodies Obtained by Computational Design and Synthetic Library Screening

Members of the Sarbecovirus subgenus of Coronaviridae have twice caused deadly threats to humans. There is increasing concern about the rapid mutation of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which has evolved into multiple generations of epidemic variants in 3 years. Broad neutralizing antibodies are of great importance for pandemic preparedness against SARS-CoV-2 variants and divergent zoonotic sarbecoviruses. Here, we analyzed the structural conservation of the receptor-binding domain (RBD) from representative sarbecoviruses and chose S2H97, a previously reported RBD antibody with ideal breadth and resistance to escape, as a template for computational design to enhance the neutralization activity and spectrum. A total of 35 designs were purified for evaluation. The neutralizing activity of a large proportion of these designs against multiple variants was increased from several to hundreds of times. Molecular dynamics simulation suggested that extra interface contacts and enhanced intermolecular interactions between the RBD and the designed antibodies are established. After light and heavy chain reconstitution, AI-1028, with five complementarity determining regions optimized, showed the best neutralizing activity across all tested sarbecoviruses, including SARS-CoV, multiple SARS-CoV-2 variants, and bat-derived viruses. AI-1028 recognized the same cryptic RBD epitope as the parental prototype antibody. In addition to computational design, chemically synthesized nanobody libraries are also a precious resource for rapid antibody development. By applying distinct RBDs as baits for reciprocal screening, we identified two novel nanobodies with broad activities. These findings provide potential pan-sarbecovirus neutralizing drugs and highlight new pathways to rapidly optimize therapeutic candidates when novel SARS-CoV-2 escape variants or new zoonotic coronaviruses emerge. IMPORTANCE The subgenus Sarbecovirus includes human SARS-CoV, SARS-CoV-2, and hundreds of genetically related bat viruses. The continuous evolution of SARS-CoV-2 has led to the striking evasion of neutralizing antibody (NAb) drugs and convalescent plasma. Antibodies with broad activity across sarbecoviruses would be helpful to combat current SARS-CoV-2 mutations and longer term animal virus spillovers. The study of pan-sarbecovirus NAbs described here is significant for the following reasons. First, we established a structure-based computational pipeline to design and optimize NAbs to obtain more potent and broader neutralizing activity across multiple sarbecoviruses. Second, we screened and identified nanobodies from a highly diversified synthetic library with a broad neutralizing spectrum using an elaborate screening strategy. These methodologies provide guidance for the rapid development of antibody therapeutics against emerging pathogens with highly variable characteristics.

Genomic screening of 16 UK native bat species through conservationist networks uncovers coronaviruses with zoonotic potential

There has been limited characterisation of bat-borne coronaviruses in Europe. Here, we screened for coronaviruses in 48 faecal samples from 16 of the 17 bat species breeding in the UK, collected through a bat rehabilitation and conservationist network. We recovered nine complete genomes, including two novel coronavirus species, across six bat species: four alphacoronaviruses, a MERS-related betacoronavirus, and four closely related sarbecoviruses. We demonstrate that at least one of these sarbecoviruses can bind and use the human ACE2 receptor for infecting human cells, albeit suboptimally. Additionally, the spike proteins of these sarbecoviruses possess an R-A-K-Q motif, which lies only one nucleotide mutation away from a furin cleavage site (FCS) that enhances infectivity in other coronaviruses, including SARS-CoV-2. However, mutating this motif to an FCS does not enable spike cleavage. Overall, while UK sarbecoviruses would require further molecular adaptations to infect humans, their zoonotic risk warrants closer surveillance.

Broad host tropism of ACE2-using MERS-related coronaviruses and determinants restricting viral recognition

Recently, two Middle East respiratory syndrome coronavirus (MERS-CoV) closely related to bat merbecoviruses, NeoCoV and PDF-2180, were discovered to use angiotensin-converting enzyme 2 (ACE2) for entry. The two viruses cannot use human ACE2 efficiently, and their host range and cross-species transmissibility across a wide range of mammalian species remain unclear. Herein, we characterized the species-specific receptor preference of these viruses by testing ACE2 orthologues from 49 bats and 53 non-bat mammals through receptor-binding domain (RBD)-binding and pseudovirus entry assays. Results based on bat ACE2 orthologues revealed that the two viruses were unable to use most, but not all, ACE2 from Yinpterochiropteran bats (Yin-bats), which is distinct from NL63 and SARS-CoV-2. Besides, both viruses exhibited broad receptor recognition spectra across non-bat mammals. Genetic and structural analyses of bat ACE2 orthologues highlighted four crucial host range determinants, all confirmed by subsequent functional assays in human and bat cells. Notably, residue 305, participating in a critical viral receptor interaction, plays a crucial role in host tropism determination, particularly in non-bat mammals. Furthermore, NeoCoV and PDF-2180 mutants with enhanced human ACE2 recognition expanded the potential host range, especially by enhancing their interaction with an evolutionarily conserved hydrophobic pocket. Our results elucidate the molecular basis for the species-specific ACE2 usage of MERS-related viruses and shed light on their zoonotic risks.

Genome characterization based on the Spike-614 and NS8-84 loci of SARS-CoV-2 reveals two major possible onsets of the COVID-19 pandemic

The global COVID-19 pandemic has lasted for 3 years since its outbreak, however its origin is still unknown. Here, we analyzed the genotypes of 3.14 million SARS-CoV-2 genomes based on the amino acid 614 of the Spike (S) and the amino acid 84 of NS8 (nonstructural protein 8), and identified 16 linkage haplotypes. The GL haplotype (S_614G and NS8_84L) was the major haplotype driving the global pandemic and accounted for 99.2% of the sequenced genomes, while the DL haplotype (S_614D and NS8_84L) caused the pandemic in China in the spring of 2020 and accounted for approximately 60% of the genomes in China and 0.45% of the global genomes. The GS (S_614G and NS8_84S), DS (S_614D and NS8_84S), and NS (S_614N and NS8_84S) haplotypes accounted for 0.26%, 0.06%, and 0.0067% of the genomes, respectively. The main evolutionary trajectory of SARS-CoV-2 is DS→DL→GL, whereas the other haplotypes are minor byproducts in the evolution. Surprisingly, the newest haplotype GL had the oldest time of most recent common ancestor (tMRCA), which was May 1 2019 by mean, while the oldest haplotype DS had the newest tMRCA with a mean of October 17, indicating that the ancestral strains that gave birth to GL had been extinct and replaced by the more adapted newcomer at the place of its origin, just like the sequential rise and fall of the delta and omicron variants. However, the haplotype DL arrived and evolved into toxic strains and ignited a pandemic in China where the GL strains had not arrived in by the end of 2019. The GL strains had spread all over the world before they were discovered, and ignited the global pandemic, which had not been noticed until the virus was declared in China. However, the GL haplotype had little influence in China during the early phase of the pandemic due to its late arrival as well as the strict transmission controls in China. Therefore, we propose two major onsets of the COVID-19 pandemic, one was mainly driven by the haplotype DL in China, the other was driven by the haplotype GL globally.

The Nucleocapsid Proteins of SARS-CoV-2 and Its Close Relative Bat Coronavirus RaTG13 Are Capable of Inhibiting PKR- and RNase L-Mediated Antiviral Pathways

Coronaviruses (CoVs), including severe acute respiratory syndrome CoV (SARS-CoV), Middle East respiratory syndrome CoV (MERS-CoV), and SARS-CoV-2, produce double-stranded RNA (dsRNA) that activates antiviral pathways such as PKR and OAS/RNase L. To successfully replicate in hosts, viruses must evade such antiviral pathways. Currently, the mechanism of how SARS-CoV-2 antagonizes dsRNA-activated antiviral pathways is unknown. In this study, we demonstrate that the SARS-CoV-2 nucleocapsid (N) protein, the most abundant viral structural protein, is capable of binding to dsRNA and phosphorylated PKR, inhibiting both the PKR and OAS/RNase L pathways. The N protein of the bat coronavirus (bat-CoV) RaTG13, the closest relative of SARS-CoV-2, has a similar ability to inhibit the human PKR and RNase L antiviral pathways. Via mutagenic analysis, we found that the C-terminal domain (CTD) of the N protein is sufficient for binding dsRNA and inhibiting RNase L activity. Interestingly, while the CTD is also sufficient for binding phosphorylated PKR, the inhibition of PKR antiviral activity requires not only the CTD but also the central linker region (LKR). Thus, our findings demonstrate that the SARS-CoV-2 N protein is capable of antagonizing the two critical antiviral pathways activated by viral dsRNA and that its inhibition of PKR activities requires more than dsRNA binding mediated by the CTD. IMPORTANCE The high transmissibility of SARS-CoV-2 is an important viral factor defining the coronavirus disease 2019 (COVID-19) pandemic. To transmit efficiently, SARS-CoV-2 must be capable of disarming the innate immune response of its host efficiently. Here, we describe that the nucleocapsid protein of SARS-CoV-2 is capable of inhibiting two critical innate antiviral pathways, PKR and OAS/RNase L. Moreover, the counterpart of the closest animal coronavirus relative of SARS-CoV-2, bat-CoV RaTG13, can also inhibit human PKR and OAS/RNase L antiviral activities. Thus, the importance of our discovery for understanding the COVID-19 pandemic is 2-fold. First, the ability of SARS-CoV-2 N to inhibit innate antiviral activity is likely a factor contributing to the transmissibility and pathogenicity of the virus. Second, the bat relative of SARS-CoV-2 has the capacity to inhibit human innate immunity, which thus likely contributed to the establishment of infection in humans. The findings described in this study are valuable for developing novel antivirals and vaccines.

Spike-Independent Infection of Human Coronavirus 229E in Bat Cells

Bats are the reservoir for numerous human pathogens, including coronaviruses. Despite many coronaviruses having descended from bat ancestors, little is known about virus-host interactions and broader evolutionary history involving bats. Studies have largely focused on the zoonotic potential of coronaviruses with few infection experiments conducted in bat cells. To determine genetic changes derived from replication in bat cells and possibly identify potential novel evolutionary pathways for zoonotic virus emergence, we serially passaged six human 229E isolates in a newly established Rhinolophus lepidus (horseshoe bat) kidney cell line. Here, we observed extensive deletions within the spike and open reading frame 4 (ORF4) genes of five 229E viruses after passaging in bat cells. As a result, spike protein expression and infectivity of human cells was lost in 5 of 6 viruses, but the capability to infect bat cells was maintained. Only viruses that expressed the spike protein could be neutralized by 229E spike-specific antibodies in human cells, whereas there was no neutralizing effect on viruses that did not express the spike protein inoculated on bat cells. However, one isolate acquired an early stop codon, abrogating spike expression but maintaining infection in bat cells. After passaging this isolate in human cells, spike expression was restored due to acquisition of nucleotide insertions among virus subpopulations. Spike-independent infection of human coronavirus 229E may provide an alternative mechanism for viral maintenance in bats that does not rely on the compatibility of viral surface proteins and known cellular entry receptors. IMPORTANCE Many viruses, including coronaviruses, originated from bats. Yet, we know little about how these viruses switch between hosts and enter human populations. Coronaviruses have succeeded in establishing in humans at least five times, including endemic coronaviruses and the recent severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In an approach to identify requirements for host switches, we established a bat cell line and adapted human coronavirus 229E viruses by serial passage. The resulting viruses lost their spike protein but maintained the ability to infect bat cells, but not human cells. Maintenance of 229E viruses in bat cells appears to be independent of a canonical spike receptor match, which in turn might facilitate cross-species transmission in bats.

Using machine learning to detect coronaviruses potentially infectious to humans

Establishing the host range for novel viruses remains a challenge. Here, we address the challenge of identifying non-human animal coronaviruses that may infect humans by creating an artificial neural network model that learns from spike protein sequences of alpha and beta coronaviruses and their binding annotation to their host receptor. The proposed method produces a human-Binding Potential (h-BiP) score that distinguishes, with high accuracy, the binding potential among coronaviruses. Three viruses, previously unknown to bind human receptors, were identified: Bat coronavirus BtCoV/133/2005 and Pipistrellus abramus bat coronavirus HKU5-related (both MERS related viruses), and Rhinolophus affinis coronavirus isolate LYRa3 (a SARS related virus). We further analyze the binding properties of BtCoV/133/2005 and LYRa3 using molecular dynamics. To test whether this model can be used for surveillance of novel coronaviruses, we re-trained the model on a set that excludes SARS-CoV-2 and all viral sequences released after the SARS-CoV-2 was published. The results predict the binding of SARS-CoV-2 with a human receptor, indicating that machine learning methods are an excellent tool for the prediction of host expansion events.

Exposure to diverse sarbecoviruses indicates frequent zoonotic spillover in human communities interacting with wildlife

BACKGROUND

Sarbecoviruses are a subgenus of Coronaviridae that mostly infect bats with known potential to infect humans (SARS-CoV and SARS-CoV-2). Populations in Southeast Asia, where these viruses are most likely to emerge, have been undersurveyed to date.

METHODS

We surveyed communities engaged in extractive industries and bat guano harvesting from rural areas in Myanmar. Participants were screened for exposure to sarbecoviruses, and their interactions with wildlife were evaluated to determine the factors associated with exposure to sarbecoviruses.

RESULTS

Of 693 people screened between July 2017 and February 2020, 12.1% were seropositive for sarbecoviruses. Individuals were significantly more likely to have been exposed to sarbecoviruses if their main livelihood involved working in extractive industries (logging, hunting, or harvesting of forest products; odds ratio [OR] = 2.71, P = 0.019) or had been hunting/slaughtering bats (OR = 6.09, P = 0.020). Exposure to a range of bat and pangolin sarbecoviruses was identified.

CONCLUSION

Exposure to diverse sarbecoviruses among high-risk human communities provides epidemiologic and immunologic evidence that zoonotic spillover is occurring. These findings inform risk mitigation efforts needed to decrease disease transmission at the bat-human interface, as well as future surveillance efforts warranted to monitor isolated populations for viruses with pandemic potential.

A comprehensive survey of bat sarbecoviruses across China in relation to the origins of SARS-CoV and SARS-CoV-2

SARS-CoV and SARS-CoV-2 have been thought to originate from bats. In this study, we screened pharyngeal and anal swabs from 13 064 bats collected between 2016 and 2021 at 703 locations across China for sarbecoviruses, covering almost all known southern hotspots, and found 146 new bat sarbecoviruses. Phylogenetic analyses of all available sarbecoviruses show that there are three different lineages-L1 as SARS-CoV-related CoVs (SARSr-CoVs), L2 as SARS-CoV-2-related CoVs (SC2r-CoVs) and novel L-R (recombinants of L1 and L2)-present in Rhinolophus pusillus bats, in the mainland of China. Among the 146 sequences, only four are L-Rs. Importantly, none belong in the L2 lineage, indicating that circulation of SC2r-CoVs in China might be very limited. All remaining 142 sequences belong in the L1 lineage, of which YN2020B-G shares the highest overall sequence identity with SARS-CoV (95.8%). The observation suggests endemic circulations of SARSr-CoVs, but not SC2r-CoVs, in bats in China. Geographic analysis of the collection sites in this study, together with all published reports, indicates that SC2r-CoVs may be mainly present in bats of Southeast Asia, including the southern border of Yunnan province, but absent in all other regions within China. In contrast, SARSr-CoVs appear to have broader geographic distribution, with the highest genetic diversity and sequence identity to human sarbecoviruses along the southwest border of China. Our data provide the rationale for further extensive surveys in broader geographical regions within, and beyond, Southeast Asia in order to find the most recent ancestors of human sarbecoviruses.

Targetable elements in SARS-CoV-2 S2 subunit for the design of pan-coronavirus fusion inhibitors and vaccines

The ongoing global pandemic of coronavirus disease 2019 (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), has caused devastating impacts on the public health and the global economy. Rapid viral antigenic evolution has led to the continual generation of new variants. Of special note is the recently expanding Omicron subvariants that are capable of immune evasion from most of the existing neutralizing antibodies (nAbs). This has posed new challenges for the prevention and treatment of COVID-19. Therefore, exploring broad-spectrum antiviral agents to combat the emerging variants is imperative. In sharp contrast to the massive accumulation of mutations within the SARS-CoV-2 receptor-binding domain (RBD), the S2 fusion subunit has remained highly conserved among variants. Hence, S2-based therapeutics may provide effective cross-protection against new SARS-CoV-2 variants. Here, we summarize the most recently developed broad-spectrum fusion inhibitors (e.g., nAbs, peptides, proteins, and small-molecule compounds) and candidate vaccines targeting the conserved elements in SARS-CoV-2 S2 subunit. The main focus includes all the targetable S2 elements, namely, the fusion peptide, stem helix, and heptad repeats 1 and 2 (HR1-HR2) bundle. Moreover, we provide a detailed summary of the characteristics and action-mechanisms for each class of cross-reactive fusion inhibitors, which should guide and promote future design of S2-based inhibitors and vaccines against new coronaviruses.

Perspective Chapter: Real-Time Genomic Surveillance for SARS-CoV-2 on Center Stage

The course of the COVID-19 pandemic depends not only on how the SARS-CoV-2 virus mutates but on the actions taken to respond to it. Important public health decisions can only be taken if we know viral dynamics, viral variants distribution, and whether new variants are emerging that may be more transmissible or/and more virulent, displaying evasion to vaccines or antiviral treatments. This situation has put the use of different approaches, such as molecular techniques and real-time genomic sequencing, to support public health decision-making on center stage. To achieve this, robust programs based on: (i) diagnostic capacity; (ii) high-throughput sequencing technologies; and (iii) high-performance bioinformatic resources, need to be established. This chapter focuses on how SARS-CoV-2 evolved since its discovery and it summarizes the scientific efforts to obtain genomic data as the virus spread throughout the globe.

Mutation-driven parallel evolution in emergence of ACE2-utilizing sarbecoviruses

Mutation and recombination are two major genetic mechanisms that drive the evolution of viruses. They both exert an interplay during virus evolution, in which mutations provide a first ancestral source of genetic diversity for subsequent recombination. Sarbecoviruses are a group of evolutionarily related β-coronaviruses including human severe acute respiratory syndrome coronavirus (SARS-CoV) and SARS-CoV-2 and a trove of related animal viruses called SARS-like CoVs (SL-CoVs). This group of members either use or not use angiotensin-converting enzyme 2 (ACE2) as their entry receptor, which has been linked to the properties of their spike protein receptor binding domains (RBDs). This raises an outstanding question regarding how ACE2 binding originated within sarbecoviruses. Using a combination of analyses of phylogenies, ancestral sequences, structures, functions and molecular dynamics, we provide evidence in favor of an evolutionary scenario, in which three distinct ancestral RBDs independently developed the ACE2 binding trait via parallel amino acid mutations. In this process, evolutionary intermediate RBDs might be firstly formed through loop extensions to offer key functional residues accompanying point mutations to remove energetically unfavorable interactions and to change the dynamics of the functional loops, all required for ACE2 binding. Subsequent optimization in the context of evolutionary intermediates led to the independent emergence of ACE2-binding RBDs in the SARS-CoV and SARS-CoV-2 clades of Asian origin and the clade comprising SL-CoVs of European and African descent. These findings will help enhance our understanding of mutation-driven evolution of sarbecoviruses in their early history.

Two pan-SARS-CoV-2 nanobodies and their multivalent derivatives effectively prevent Omicron infections in mice

With the widespread vaccinations against coronavirus disease 2019 (COVID-19), we are witnessing gradually waning neutralizing antibodies and increasing cases of breakthrough infections, necessitating the development of drugs aside from vaccines, particularly ones that can be administered outside of hospitals. Here, we present two cross-reactive nanobodies (R14 and S43) and their multivalent derivatives, including decameric ones (fused to the immunoglobulin M [IgM] Fc) that maintain potent neutralizing activity against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) after aerosolization and display not only pan-SARS-CoV-2 but also varied pan-sarbecovirus activities. Through respiratory administration to mice, monovalent and decameric R14 significantly reduce the lung viral RNAs at low dose and display potent pre- and post-exposure protection. Furthermore, structural studies reveal the neutralizing mechanisms of R14 and S43 and the multiple inhibition effects that the multivalent derivatives exert. Our work demonstrates promising convenient drug candidates via respiratory administration against SARS-CoV-2 infection, which can contribute to containing the COVID-19 pandemic.