Structural and Molecular Evidence Suggesting Coronavirus-driven Evolution of Mouse Receptor

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.

Lys417 acts as a molecular switch that regulates the conformation of SARS-CoV-2 spike protein

SARS-CoV-2 spike protein plays a key role in mediating viral entry and inducing host immune responses. It can adopt either an open or closed conformation based on the position of its receptor-binding domain (RBD). It is yet unclear what causes these conformational changes or how they influence the spike's functions. Here, we show that Lys417 in the RBD plays dual roles in the spike's structure: it stabilizes the closed conformation of the trimeric spike by mediating inter-spike-subunit interactions; it also directly interacts with ACE2 receptor. Hence, a K417V mutation has opposing effects on the spike's function: it opens up the spike for better ACE2 binding while weakening the RBD's direct binding to ACE2. The net outcomes of this mutation are to allow the spike to bind ACE2 with higher probability and mediate viral entry more efficiently, but become more exposed to neutralizing antibodies. Given that residue 417 has been a viral mutational hotspot, SARS-CoV-2 may have been evolving to strike a balance between infection potency and immune evasion, contributing to its pandemic spread.

Kathryn V. Holmes: A Career of Contributions to the Coronavirus Field

Over the past two years, scientific research has moved at an unprecedented rate in response to the COVID-19 pandemic. The rapid development of effective vaccines and therapeutics would not have been possible without extensive background knowledge on coronaviruses developed over decades by researchers, including Kathryn (Kay) Holmes. Kay's research team discovered the first coronavirus receptors for mouse hepatitis virus and human coronavirus 229E and contributed a wealth of information on coronaviral spike glycoproteins and receptor interactions that are critical determinants of host and tissue specificity. She collaborated with several research laboratories to contribute knowledge in additional areas, including coronaviral pathogenesis, epidemiology, and evolution. Throughout her career, Kay was an extremely dedicated and thoughtful mentor to numerous graduate students and post-doctoral fellows. This article provides a review of her contributions to the coronavirus field and her exemplary mentoring.

Human CEACAM1 N-domain dimerization is independent from glycan modifications

Carcinoembryonic cellular adhesion molecules (CEACAMs) serve diverse roles in cell signaling, proliferation, and survival and are made up of one or several immunoglobulin (Ig)-like ectodomains glycosylated in vivo. The physiological oligomeric state and how it contributes to protein function are central to understanding CEACAMs. Two putative dimer conformations involving different CEACAM1 N-terminal Ig-like domain (CCM1) protein faces (ABED and GFCC'C″) were identified from crystal structures. GFCC'C″ was identified as the dominant CCM1 solution dimer, but ambiguity regarding the effect of glycosylation on dimer formation calls its physiological relevance into question. We present the first crystal structure of minimally glycosylated CCM1 in the GFCC'C″ dimer conformation and characterization in solution by continuous-wave and double electron-electron resonance electron paramagnetic resonance spectroscopy. Our results suggest the GFCC'C″ dimer is dominant in solution with different levels of glycosylation, and structural conservation and co-evolved residues support that the GFCC'C″ dimer is conserved across CEACAMs.

Molecular switches regulating the potency and immune evasiveness of SARS-CoV-2 spike protein

SARS-CoV-2 spike protein plays a key role in viral entry and host immune responses. The conformation of the spike protein can be either open or closed, yet it is unclear how the conformations affect the protein's functions or what regulate the conformational changes. Using SARS-CoV-1 and bat RaTG13-CoV as comparisons, we identified two molecular switches that regulate the conformations of SARS-CoV-2 spike protein: (i) a furin motif loop turns SARS-CoV-2 spike from a closed conformation to a mixture of open and closed conformations, and (ii) a K417V mutation turns SARS-CoV-2 spike from mixed conformations to an open conformation. We showed that the open conformation favors viral potency by exposing the RBD for receptor binding and viral entry, whereas the closed conformation supports viral immune evasion by hiding the RBD from neutralizing antibodies. Hence SARS-CoV-2 spike has evolved to reach a balance between potency and immune evasiveness, which may contribute to the pandemic spread of SARS-CoV-2. The dynamics between viral potency and invasiveness is likely to further evolve, providing insights into future evolution of SARS-CoV-2.

Novel virus-like nanoparticle vaccine effectively protects animal model from SARS-CoV-2 infection

The key to battling the COVID-19 pandemic and its potential aftermath is to develop a variety of vaccines that are efficacious and safe, elicit lasting immunity, and cover a range of SARS-CoV-2 variants. Recombinant viral receptor-binding domains (RBDs) are safe vaccine candidates but often have limited efficacy due to the lack of virus-like immunogen display pattern. Here we have developed a novel virus-like nanoparticle (VLP) vaccine that displays 120 copies of SARS-CoV-2 RBD on its surface. This VLP-RBD vaccine mimics virus-based vaccines in immunogen display, which boosts its efficacy, while maintaining the safety of protein-based subunit vaccines. Compared to the RBD vaccine, the VLP-RBD vaccine induced five times more neutralizing antibodies in mice that efficiently blocked SARS-CoV-2 from attaching to its host receptor and potently neutralized the cell entry of variant SARS-CoV-2 strains, SARS-CoV-1, and SARS-CoV-1-related bat coronavirus. These neutralizing immune responses induced by the VLP-RBD vaccine did not wane during the two-month study period. Furthermore, the VLP-RBD vaccine effectively protected mice from SARS-CoV-2 challenge, dramatically reducing the development of clinical signs and pathological changes in immunized mice. The VLP-RBD vaccine provides one potentially effective solution to controlling the spread of SARS-CoV-2.

Evolutionary analysis of the anti-viral STAT2 gene of primates and rodents: Signature of different stages of an arms race

STAT2 plays a strategic role in defending viral infection through the signaling cascade involving the immune system initiated after type I interferon release. Many flaviviruses target the inactivation or degradation of STAT2 as a strategy to impair this host's line of defense. Primates are natural reservoirs for a range of disease-causing flaviviruses (e.g., Zika, Dengue, and Yellow Fever virus), while rodents appear less susceptible. We analyzed the STAT2 coding sequence of 28 Rodentia species and 49 Primates species. Original data from 19 Platyrrhini species were sequenced for the SH2 domain of STAT2 and included in the analysis. STAT2 has many sites whose variation can be explained by positive selection, measurement by two methods (PALM indicated 12, MEME 61). Both evolutionary tests significantly marked sites 127, 731, 739, 766, and 780. SH2 is under evolutionary constraint but presents episodic positive selection events within Rodentia: in one of them, a moderately radical change (serine > arginine) at position 638 is found in Peromyscus species, and can be implicated in the difference in susceptibility to flaviviruses within Rodentia. Some other positively selected sites are functional such as 5, 95, 203, 251, 782, and 829. Sites 251 and 287 regulate the signaling mediated by the JAK-STAT2 pathway, while 782 and 829 create a stable tertiary structure of STAT2, facilitating its connection with transcriptional co-activators. Only three positively selected sites, 5, 95, and 203, are recognized members who act on the interface between STAT2 and flaviviruses NS5 protein. We suggested that due to the higher evolutionary rate, rodents are, at this moment, taking some advantage in the battle against infections for some well-known Flaviviridae, in particular when compared to primates. Our results point to dynamics that fit with a molecular evolutionary scenario shaped by a thought-provoking virus-host arms race.

Anti-SARS-CoV-2 Natural Products as Potentially Therapeutic Agents

Severe acute respiratory syndrome-related coronavirus-2 (SARS-CoV-2), a β-coronavirus, is the cause of the recently emerged pandemic and worldwide outbreak of respiratory disease. Researchers exchange information on COVID-19 to enable collaborative searches. Although there is as yet no effective antiviral agent, like tamiflu against influenza, to block SARS-CoV-2 infection to its host cells, various candidates to mitigate or treat the disease are currently being investigated. Several drugs are being screened for the ability to block virus entry on cell surfaces and/or block intracellular replication in host cells. Vaccine development is being pursued, invoking a better elucidation of the life cycle of the virus. SARS-CoV-2 recognizes O-acetylated neuraminic acids and also several membrane proteins, such as ACE2, as the result of evolutionary switches of O-Ac SA recognition specificities. To provide information related to the current development of possible anti-SARS-COV-2 viral agents, the current review deals with the known inhibitory compounds with low molecular weight. The molecules are mainly derived from natural products of plant sources by screening or chemical synthesis via molecular simulations. Artificial intelligence-based computational simulation for drug designation and large-scale inhibitor screening have recently been performed. Structure-activity relationship of the anti-SARS-CoV-2 natural compounds is discussed.

Antiviral and Anti‐Inflammatory Treatment with Multifunctional Alveolar Macrophage‐Like Nanoparticles in a Surrogate Mouse Model of COVID‐19

The pandemic of coronavirus disease 2019 (COVID‐19) is continually worsening. Clinical treatment for COVID‐19 remains primarily supportive with no specific medicines or regimens. Here, the development of multifunctional alveolar macrophage (AM)‐like nanoparticles (NPs) with photothermal inactivation capability for COVID‐19 treatment is reported. The NPs, made by wrapping polymeric cores with AM membranes, display the same surface receptors as AMs, including the coronavirus receptor and multiple cytokine receptors. By acting as AM decoys, the NPs block coronavirus from host cell entry and absorb various proinflammatory cytokines, thus achieving combined antiviral and anti‐inflammatory treatment. To enhance the antiviral efficiency, an efficient photothermal material based on aggregation‐induced emission luminogens is doped into the NPs for virus photothermal disruption under near‐infrared (NIR) irradiation. In a surrogate mouse model of COVID‐19 caused by murine coronavirus, treatment with multifunctional AM‐like NPs with NIR irradiation decreases virus burden and cytokine levels, reduces lung damage and inflammation, and confers a significant survival advantage to the infected mice. Crucially, this therapeutic strategy may be clinically applied for the treatment of COVID‐19 at early stage through atomization inhalation of the NPs followed by NIR irradiation of the respiratory tract, thus alleviating infection progression and reducing transmission risk.

Research progress on coronavirus S proteins and their receptors

Coronaviruses are a large family of important pathogens that cause human and animal diseases. At the end of 2019, a pneumonia epidemic caused by a novel coronavirus brought attention to coronaviruses. Exploring the interaction between the virus and its receptor will be helpful in developing preventive vaccines and therapeutic drugs. The coronavirus spike protein (S) plays an important role in both binding to receptors on host cells and fusion of the viral membrane with the host cell membrane. This review introduces the structure and function of the S protein and its receptor, focusing on the binding mode and binding region of both.

Evidence for the loss and recovery of SLAMF9 during human evolution: implications on Dollo's law

Signaling lymphocyte activation molecule family member 9 (SLAMF9) is a cell surface protein of the CD2/SLAM family of leukocyte surface receptors. It is conserved throughout mammals and has roles in the initiation of inflammatory responses and regulation of plasmacytoid dendritic cell function. Through comparison of reference sequences encoding SLAMF9 in human, mouse, and primate sequences, we have determined that the SLAMF9 gene underwent successive mutation events, resulting in the loss of the protein and subsequent recovery of a less stable version. The mutations included a single base pair deletion in the second exon and a change in the splice acceptor site of that same exon. These changes would have had the effect of creating and later repairing a frameshift in the coding sequence. These events took place since the divergence of the human lineage from the chimpanzee-human last common ancestor and represent the first known case of the functional loss and recovery of a gene within the human lineage.

SARS-CoV-2 Evolutionary Adaptation toward Host Entry and Recognition of Receptor O-Acetyl Sialylation in Virus-Host Interaction

The recently emerged SARS-CoV-2 is the cause of the global health crisis of the coronavirus disease 2019 (COVID-19) pandemic. No evidence is yet available for CoV infection into hosts upon zoonotic disease outbreak, although the CoV epidemy resembles influenza viruses, which use sialic acid (SA). Currently, information on SARS-CoV-2 and its receptors is limited. O-acetylated SAs interact with the lectin-like spike glycoprotein of SARS CoV-2 for the initial attachment of viruses to enter into the host cells. SARS-CoV-2 hemagglutinin-esterase (HE) acts as the classical glycan-binding lectin and receptor-degrading enzyme. Most β-CoVs recognize 9-O-acetyl-SAs but switched to recognizing the 4-O-acetyl-SA form during evolution of CoVs. Type I HE is specific for the 9-O-Ac-SAs and type II HE is specific for 4-O-Ac-SAs. The SA-binding shift proceeds through quasi-synchronous adaptations of the SA-recognition sites of the lectin and esterase domains. The molecular switching of HE acquisition of 4-O-acetyl binding from 9-O-acetyl SA binding is caused by protein-carbohydrate interaction (PCI) or lectin-carbohydrate interaction (LCI). The HE gene was transmitted to a β-CoV lineage A progenitor by horizontal gene transfer from a 9-O-Ac-SA-specific HEF, as in influenza virus C/D. HE acquisition, and expansion takes place by cross-species transmission over HE evolution. This reflects viral evolutionary adaptation to host SA-containing glycans. Therefore, CoV HE receptor switching precedes virus evolution driven by the SA-glycan diversity of the hosts. The PCI or LCI stereochemistry potentiates the SA-ligand switch by a simple conformational shift of the lectin and esterase domains. Therefore, examination of new emerging viruses can lead to better understanding of virus evolution toward transitional host tropism. A clear example of HE gene transfer is found in the BCoV HE, which prefers 7,9-di-O-Ac-SAs, which is also known to be a target of the bovine torovirus HE. A more exciting case of such a switching event occurs in the murine CoVs, with the example of the β-CoV lineage A type binding with two different subtypes of the typical 9-O-Ac-SA (type I) and the exclusive 4-O-Ac-SA (type II) attachment factors. The protein structure data for type II HE also imply the virus switching to binding 4-O acetyl SA from 9-O acetyl SA. Principles of the protein-glycan interaction and PCI stereochemistry potentiate the SA-ligand switch via simple conformational shifts of the lectin and esterase domains. Thus, our understanding of natural adaptation can be specified to how carbohydrate/glycan-recognizing proteins/molecules contribute to virus evolution toward host tropism. Under the current circumstances where reliable antiviral therapeutics or vaccination tools are lacking, several trials are underway to examine viral agents. As expected, structural and non-structural proteins of SARS-CoV-2 are currently being targeted for viral therapeutic designation and development. However, the modern global society needs SARS-CoV-2 preventive and therapeutic drugs for infected patients. In this review, the structure and sialobiology of SARS-CoV-2 are discussed in order to encourage and activate public research on glycan-specific interaction-based drug creation in the near future.

ACE2 diversity in placental mammals reveals the evolutionary strategy of SARS-CoV-2

The recent emergence of SARS-CoV-2 is responsible for the current pandemic of COVID-19, which uses the human membrane protein ACE2 as a gateway to host-cell infection. We performed a comparative genomic analysis of 70 ACE2 placental mammal orthologues to identify variations and contribute to the understanding of evolutionary dynamics behind this successful adaptation to infect humans. Our results reveal that 4% of the ACE2 sites are under positive selection, all located in the catalytic domain, suggesting possibly taxon-specific adaptations related to the ACE2 function, such as cardiovascular physiology. Considering all variable sites, we selected 30 of them located at the critical ACE2 binding sites to the SARS-CoV-like viruses for analysis in more detail. Our results reveal a relatively high diversity of ACE2 between placental mammal species, while showing no polymorphism within human populations, at least considering the 30 inter-species variable sites. A perfect scenario for natural selection favored this opportunistic new coronavirus in its trajectory of infecting humans. We suggest that SARS-CoV-2 became a specialist coronavirus for human hosts. Differences in the rate of infection and mortality could be related to the innate immune responses, other unknown genetic factors, as well as non-biological factors.

An Evolutionary Insertion in the Mxra8 Receptor-Binding Site Confers Resistance to Alphavirus Infection and Pathogenesis

Alphaviruses are emerging, mosquito-transmitted RNA viruses with poorly understood cellular tropism and species selectivity. Mxra8 is a receptor for multiple alphaviruses including chikungunya virus (CHIKV). We discovered that while expression of mouse, rat, chimpanzee, dog, horse, goat, sheep, and human Mxra8 enables alphavirus infection in cell culture, cattle Mxra8 does not. Cattle Mxra8 encodes a 15-amino acid insertion in its ectodomain that prevents Mxra8 binding to CHIKV. Identical insertions are present in zebu, yak, and the extinct auroch. As other Bovinae lineages contain related Mxra8 sequences, this insertion likely occurred at least 5 million years ago. Removing the Mxra8 insertion in Bovinae enhances alphavirus binding and infection, while introducing the insertion into mouse Mxra8 blocks CHIKV binding, prevents infection by multiple alphaviruses in cells, and mitigates CHIKV-induced pathogenesis in mice. Our studies on how this insertion provides resistance to CHIKV infection could facilitate countermeasures that disrupt Mxra8 interactions with alphaviruses.

Structure of mouse coronavirus spike protein complexed with receptor reveals mechanism for viral entry

Coronaviruses recognize a variety of receptors using different domains of their envelope-anchored spike protein. How these diverse receptor recognition patterns affect viral entry is unknown. Mouse hepatitis coronavirus (MHV) is the only known coronavirus that uses the N-terminal domain (NTD) of its spike to recognize a protein receptor, CEACAM1a. Here we determined the cryo-EM structure of MHV spike complexed with mouse CEACAM1a. The trimeric spike contains three receptor-binding S1 heads sitting on top of a trimeric membrane-fusion S2 stalk. Three receptor molecules bind to the sides of the spike trimer, where three NTDs are located. Receptor binding induces structural changes in the spike, weakening the interactions between S1 and S2. Using protease sensitivity and negative-stain EM analyses, we further showed that after protease treatment of the spike, receptor binding facilitated the dissociation of S1 from S2, allowing S2 to transition from pre-fusion to post-fusion conformation. Together these results reveal a new role of receptor binding in MHV entry: in addition to its well-characterized role in viral attachment to host cells, receptor binding also induces the conformational change of the spike and hence the fusion of viral and host membranes. Our study provides new mechanistic insight into coronavirus entry and highlights the diverse entry mechanisms used by different viruses.

Coronaviruses recognize many receptors using their envelope-anchored spike protein. The role of receptor binding in coronavirus entry into host cells is a fundamental question in virology. Mouse hepatitis coronavirus (MHV) is unique among all coronaviruses in that it uses the N-terminal domain (NTD) of its spike protein to bind a protein receptor CEACAM1a. While extensive research has been performed on the cell entry mechanisms of coronaviruses that use a different domain of their spike protein for receptor binding, the cell entry mechanism for MHV is still elusive. Here we determined the cryo-EM structure of MHV spike protein complexed with CEACAM1a. The structure reveals unique features of receptor binding by MHV spike that facilitate the structural changes of MHV spike and promote cell entry of MHV. We further confirmed the structural results with biochemical and negative-stain EM analyses. These results suggest that receptor binding plays dual roles in MHV entry: it promotes both viral attachment to host cells and the fusion of host and viral membranes. Our study provides insight into the molecular mechanism of MHV entry, demonstrating how cell entry of MHV has been adapted to its unique way of receptor binding.

Evolution: Decoy Receptors as Unique Weapons to Fight Pathogens

Many pathogens must bind to entry receptors on the surfaces of host cells yet avoid any closely-related phagocytic decoy receptors on granulocytes that evolved as a host defense mechanism. The discovery of decoy-receptor polymorphisms in human populations now points to an evolutionary process that allows the host to catch up with pathogens.

Glycine 29 Is Critical for Conformational Changes of the Spike Glycoprotein of Mouse Hepatitis Virus A59 Triggered by either Receptor Binding or High pH

Mouse hepatitis virus (MHV) uses its N-terminal domain (NTD) of the viral spike (S) protein to bind the host receptor mouse carcinoembryonic antigen-related cell adhesion molecule 1a (mCEACAM1a) and mediate virus entry. Our previous crystal structure study of the MHV NTD/mCEACAM1a complex (G. Peng, D. Sun, K. R. Rajashankar, Z. Qian, et al., Proc Natl Acad Sci U S A 108:10696-10701, 2011, https://doi.org/10.1073/pnas.1104306108) reveals that there are 14 residues in the NTD interacting with the receptor. However, their contribution to receptor binding and virus entry has not been fully investigated. Here we analyzed 13 out of 14 contact residues by mutagenesis and identified I22 as being essential for receptor binding and virus entry. Unexpectedly, we found that G29 was critical for the conformational changes of the S protein triggered by either receptor binding or high pH. Replacement of G29 with A, D, F, K, M, and T, to different extents, caused spontaneous dissociation of S1 from the S protein, resulting in an enhancement of high-pH-triggered receptor-independent syncytium (RIS) formation in HEK293T cells, compared to the wild type (WT). In contrast, replacement of G29 with P, a turn-prone residue with a strict conformation, hindered virus entry and conformational changes of the S protein triggered by either receptor binding or pH 8.0, suggesting that the structural turn around G29 and its flexibility are critical. Finally, stabilization of the NTD by G29P had almost no effect on pH-independent RIS induced by the Y320A mutation in the C-terminal domain (CTD) of the S1 subunit, indicating that there might be an absence of cross talk between the NTD and CTD during conformational changes of the S protein. Our study will aid in better understanding the mechanism of how conformational changes of the S protein are triggered.IMPORTANCE Binding of the MHV S protein to the receptor mCEACAM1a triggers conformational changes of S proteins, leading to the formation of a six-helix bundle and viral and cellular membrane fusion. However, the mechanism by which the conformational change of the S protein is initiated after receptor binding has not been determined. In this study, we showed that while replacement of G29, a residue at the edge of the receptor binding interface and the center of the structural turn after the β1-sheet of the S protein, with D or T triggered spontaneous conformational changes of the S protein and pH-independent RIS, the G29P mutation significantly impeded the conformational changes of S proteins triggered by either receptor binding or pH 8.0. We reason that this structural turn might be critical for conformational changes of the S protein and that altering this structural turn could initiate conformational changes of the S protein, leading to membrane fusion.

Cryo-EM structure of infectious bronchitis coronavirus spike protein reveals structural and functional evolution of coronavirus spike proteins

As cell-invading molecular machinery, coronavirus spike proteins pose an evolutionary conundrum due to their high divergence. In this study, we determined the cryo-EM structure of avian infectious bronchitis coronavirus (IBV) spike protein from the γ-genus. The trimeric IBV spike ectodomain contains three receptor-binding S1 heads and a trimeric membrane-fusion S2 stalk. While IBV S2 is structurally similar to those from the other genera, IBV S1 possesses structural features that are unique to different other genera, thereby bridging these diverse spikes into an evolutionary spectrum. Specifically, among different genera, the two domains of S1, the N-terminal domain (S1-NTD) and C-terminal domain (S1-CTD), diverge from simpler tertiary structures and quaternary packing to more complex ones, leading to different functions of the spikes in receptor usage and membrane fusion. Based on the above structural and functional comparisons, we propose that the evolutionary spectrum of coronavirus spikes follows the order of α-, δ-, γ-, and β-genus. This study has provided insight into the evolutionary relationships among coronavirus spikes and deepened our understanding of their structural and functional diversity.