scholarly journals Conserved Sequence Features in the Spike Protein Provide Evidence Suggesting the Origin of SARS-CoV-2 (COVID-19) Virus by Recombination between SARS virus and Another Sarbecovirus

Author(s):  
Radhey S. Gupta ◽  
Bijendra Khadka

Both SARS-CoV-2 (COVID-19) and SARS coronaviruses (CoVs) are members of the subgenus Sarbecovirus. To understand the origin of SARS-CoV-2, protein sequences from sarbecoviruses were analyzed to identify highly-specific molecular markers consisting of conserved inserts or deletions (termed CSIs) in the spike (S) and nucleocapsid (N) proteins that are specific for either particular clusters/lineages of these viruses or are commonly shared by specific lineages. Three novel CSIs in the N-terminal domain of the spike protein S1-subunit (S1-NTD) are uniquely shared by the SARS-CoV-2, BatCoV-RaTG13 and most pangolin CoVs, distinguishing this cluster of viruses (SARS-CoV-2r) from all others. In the same positions, where these CSIs are found, related CSIs are also present in two other sarbecoviruses (viz. CoVZXC21 and CoVZC45 forming CoVZC cluster), which form an out group of the SARS-CoV-2r cluster. These three CSIs are not found in the SARS-CoVs. However, both SARS and SARS-CoV-2r CoVs contain two large CSIs in the C-terminal domain of S1 (S1-CTD), which binds the human ACE-2 receptor, that are absent in the CoVZC cluster of CoVs. These results indicate that while the S1-NTD of the SARS-CoV-2r viruses possesses the sequence characteristics of the CoVZC cluster of CoVs, their S1-CTD resembles the SARS viruses. Thus, the spike protein of SARS-CoV-2r viruses has likely originated from a recombination event between the S1-NTD of the CoVZC viruses and the S1-CTD of SARS viruses. This inference is also supported by the amino acid sequence similarity of the S1-NTD and S1-CTD from SARS-CoV-2 compared to the CoVZC and SARS CoVs. We also present evidence that one of the pangolin-CoV_MP789, whose receptor-binding domain is most similar to the SARS-CoV-2, is also derived by a recent recombination between the S1-NTD of the CoVZC CoVs and the S1-CTD of a SARS-CoV-2 related virus. Several other identified CSIs are specific for others clusters of sarbecoviruses including a clade consisting of bat SARS-CoVs (BM48-31/BGR/2008 and SARS_BtKY72). Structural mappings studies show that the identified CSIs are located within surface-exposed loops and form distinct patches on the surface of the spike protein. These surface loops/patches are predicted to interact with other host components and play important role in the biology/pathology of SARS-CoV-2 virus. Lastly, the CSIs specific for the SARS-CoV-2r clade provide novel means for development of new diagnostic and therapeutic targets for these viruses.

PeerJ ◽  
2021 ◽  
Vol 9 ◽  
pp. e12434
Author(s):  
Bijendra Khadka ◽  
Radhey S. Gupta

Both SARS-CoV-2 and SARS coronaviruses (CoVs) are members of the subgenus Sarbecovirus. To understand the origin of SARS-CoV-2, sequences for the spike and nucleocapsid proteins from sarbecoviruses were analyzed to identify molecular markers consisting of conserved inserts or deletions (termed CSIs) that are specific for either a particular clade of Sarbecovirus or are commonly shared by two or more clades of these viruses. Three novel CSIs in the N-terminal domain (NTD) of the spike protein S1-subunit (S1-NTD) are uniquely shared by SARS-CoV-2, Bat-CoV-RaTG13 and most pangolin CoVs (SARS-CoV-2r clade). Three other sarbecoviruses viz. bat-CoVZXC21, -CoVZC45 and -PrC31 (forming CoVZC/PrC31 clade), and a pangolin-CoV_MP789 also contain related CSIs in the same positions. In contrast to the S1-NTD, both SARS and SARS-CoV-2r viruses contain two large CSIs in the S1-C-terminal domain (S1-CTD) that are absent in the CoVZC/PrC31 clade. One of these CSIs, consisting of a 12 aa insert, is also present in the RShSTT clade (Cambodia-CoV strains). Sequence similarity studies show that the S1-NTD of SARS-CoV-2r viruses is most similar to the CoVZC/PrC31 clade, whereas their S1-CTD exhibits highest similarity to the RShSTT- (and the SARS-related) CoVs. Results from the shared presence of CSIs and sequence similarity studies on different CoV lineages support the inference that the SARS-CoV-2r cluster of viruses has originated by a genetic recombination between the S1-NTD of the CoVZC/PrC31 clade of CoVs and the S1-CTD of RShSTT/SARS viruses, respectively. We also present compelling evidence, based on the shared presence of CSIs and sequence similarity studies, that the pangolin-CoV_MP789, whose receptor-binding domain is most similar to the SARS-CoV-2 virus, has resulted from another independent recombination event involving the S1-NTD of the CoVZC/PrC31 CoVs and the S1-CTD of an unidentified SARS-CoV-2r related virus. The SARS-CoV-2 virus involved in this latter recombination event is postulated to be most similar to the SARS-CoV-2. Several other CSIs reported here are specific for other clusters of sarbecoviruses including a clade consisting of bat-SARS-CoVs (BM48-31/BGR/2008 and SARS_BtKY72). Structural mapping studies show that the identified CSIs form distinct loops/patches on the surface of the spike protein. It is hypothesized that these novel loops/patches on the spike protein, through their interactions with other host components, should play important roles in the biology/pathology of SARS-CoV-2 virus. Lastly, the CSIs specific for different clades of sarbecoviruses including SARS-CoV-2r clade provide novel means for the identification of these viruses and other potential applications.


Author(s):  
Christopher A. Beaudoin ◽  
Samir W. Hamaia ◽  
Christopher L.-H. Huang ◽  
Tom L. Blundell ◽  
Antony P. Jackson

The RGD motif in the Severe Acute Syndrome Coronavirus 2 (SARS-CoV-2) spike protein has been predicted to bind RGD-recognizing integrins. Recent studies have shown that the spike protein does, indeed, interact with αVβ3 and α5β1 integrins, both of which bind to RGD-containing ligands. However, computational studies have suggested that binding between the spike RGD motif and integrins is not favourable, even when unfolding occurs after conformational changes induced by binding to the canonical host entry receptor, angiotensin-converting enzyme 2 (ACE2). Furthermore, non-RGD-binding integrins, such as αx, have been suggested to interact with the SARS-CoV-2 spike protein. Other viral pathogens, such as rotaviruses, have been recorded to bind integrins in an RGD-independent manner to initiate host cell entry. Thus, in order to consider the potential for the SARS-CoV-2 spike protein to bind integrins independent of the RGD sequence, we investigate several factors related to the involvement of integrins in SARS-CoV-2 infection. First, we review changes in integrin expression during SARS-CoV-2 infection to identify which integrins might be of interest. Then, all known non-RGD integrin-binding motifs are collected and mapped to the spike protein receptor-binding domain and analyzed for their 3D availability. Several integrin-binding motifs are shown to exhibit high sequence similarity with solvent accessible regions of the spike receptor-binding domain. Comparisons of these motifs with other betacoronavirus spike proteins, such as SARS-CoV and RaTG13, reveal that some have recently evolved while others are more conserved throughout phylogenetically similar betacoronaviruses. Interestingly, all of the potential integrin-binding motifs, including the RGD sequence, are conserved in one of the known pangolin coronavirus strains. Of note, the most recently recorded mutations in the spike protein receptor-binding domain were found outside of the putative integrin-binding sequences, although several mutations formed inside and close to one motif, in particular, may potentially enhance binding. These data suggest that the SARS-CoV-2 spike protein may interact with integrins independent of the RGD sequence and may help further explain how SARS-CoV-2 and other viruses can evolve to bind to integrins.


Author(s):  
Feroza Begum ◽  
Debica Mukherjee ◽  
Dluya Thagriki ◽  
Sandeepan Das ◽  
Prem Prakash Tripathi ◽  
...  

AbstractIndia has recently started sequencing SARS-CoV2 genome from clinical isolates. Currently only few sequences are available from three states in India. Kerala was the first state to deposit complete sequence from two isolates followed by one from Gujarat. On April 27, 2020, the first five sequences from the state of West Bengal (Eastern India) were deposited on ‘a global initiative on sharing avian flu data’ (GISAID) platform. In this paper we have analysed the spike protein sequences from all these five isolates and also compared for their similarities or differences with other sequences reported in India and with isolates of Wuhan origin. We report one unique mutation at position 723 and the other at 1124 in the S2 domain of spike protein of the isolates from West Bengal only and one mutation downstream of the receptor binding domain at position 614 in S1 domain which was common with the sequence from Gujarat (a state of western part of India). Mutation in the S2 domain showed changes in the secondary structure of the spike protein at region of mutation. We also studied molecular dynamics using normal mode analyses and found that this mutation decreases the flexibility of S2 domain. Since both S1 and S2 are important in receptor binding followed by entry in the host cells, such mutations may define the affinity or avidity of receptor binding.


Author(s):  
Mayanka Awasthi ◽  
Sahil Gulati ◽  
Debi P. Sarkar ◽  
Swasti Tiwari ◽  
Suneel Kateriya ◽  
...  

Abstract COVID-19 novel coronavirus disease caused by SARS-CoV-2 causes severe lethal respiratory illness in humans and has recently developed into a worldwide pandemic. The lack of effective treatment strategy and vaccines against SARS-CoV-2 poses a threat to human health. Extremely high infection rate and multi-organ secondary infection within a short time period makes this virus more deadly and challenging for therapeutic interventions. Despite of high sequence similarity and utilization of common host-cell receptor, human angiotensin-converting enzyme-2 (ACE2) for virus entry, SARS-CoV-2 is much infectious than SARS-CoV. Structure-based sequence comparison of the N-terminal domain (NTD) of spike protein of MERS-CoV, SARS-CoV and SARS-CoV-2 illustrate three short stretches of amino acid motifs in SARS-CoV-2 , which appears to be the reminiscent of MERS-CoV sialoside binding pockets. These key differences with SARS-CoV and similarity with MERS-CoV, suggest an evolutionary adaptation of SARS-CoV-2 spike protein reciprocal interaction with host surface sialosides.


Author(s):  
Lalitha Guruprasad

<p>The human SARS-CoV-2 spike protein sequences from Asia, Africa, Europe, North America, South America and Oceania were analyzed by comparing with the reference SARS-CoV-2 protein sequence from Wuhan-Hu-1, China. Out of 10,333 spike protein sequences analyzed, 8,155 proteins comprised one or more mutations. A total of 9,654 mutations were observed that correspond to 400 distinct mutation sites. The receptor binding domain (RBD) which is involved in the interactions with human ACE-2 receptor and causes infection leading to the COVID-19 disease comprised 44 mutations that included residues within 3.2 Å interacting distance from the ACE-2 receptor. The mutations observed in the spike proteins are discussed in the context of their distribution according to the geographical locations, mutation sites, mutation types, distribution of the number of mutations at the mutation sites and mutations at the glycosylation sites. The density of mutations in different regions of the spike protein sequence and location of the mutations in protein three-dimensional structure corresponding to the RBD are discussed. The mutations identified in the present work are important considerations for antibody, vaccine and drug development.</p>


2021 ◽  
Author(s):  
Tzu-Jing Yang ◽  
Pei-Yu Yu ◽  
Yuan-Chih Chang ◽  
Ning-En Chang ◽  
Yu-Xi Tsai ◽  
...  

The surge of COVID-19 infection cases is spurred by emerging SARS-CoV-2 variants such as B.1.617. Here we report 38 cryo-EM structures, corresponding to the spike protein of the Beta (B.1.351), Gamma (P.1), Delta (B.1.617.2) and Kappa (B.1.617.1) variants in different functional states with and without its receptor, ACE2. Mutations on the N-terminal domain not only alter the conformation of the highly antigenic supersite of the Delta variant, but also remodel the glycan shield by deleting or adding N-glycans of the Delta and Gamma variants, respectively. Substantially enhanced ACE2 binding was observed for all variants, whose mutations on the receptor binding domain modulate the electrostatics of the binding interfaces. Despite their abilities to escape host immunity, all variants can be potently neutralized by three unique antibodies.


Science ◽  
2020 ◽  
Vol 369 (6504) ◽  
pp. 650-655 ◽  
Author(s):  
Xiangyang Chi ◽  
Renhong Yan ◽  
Jun Zhang ◽  
Guanying Zhang ◽  
Yuanyuan Zhang ◽  
...  

Developing therapeutics against severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) could be guided by the distribution of epitopes, not only on the receptor binding domain (RBD) of the Spike (S) protein but also across the full Spike (S) protein. We isolated and characterized monoclonal antibodies (mAbs) from 10 convalescent COVID-19 patients. Three mAbs showed neutralizing activities against authentic SARS-CoV-2. One mAb, named 4A8, exhibits high neutralization potency against both authentic and pseudotyped SARS-CoV-2 but does not bind the RBD. We defined the epitope of 4A8 as the N-terminal domain (NTD) of the S protein by determining with cryo–eletron microscopy its structure in complex with the S protein to an overall resolution of 3.1 angstroms and local resolution of 3.3 angstroms for the 4A8-NTD interface. This points to the NTD as a promising target for therapeutic mAbs against COVID-19.


Author(s):  
Xiangyang Chi ◽  
Renhong Yan ◽  
Jun Zhang ◽  
Guanying Zhang ◽  
Yuanyuan Zhang ◽  
...  

AbstractThe pandemic of coronavirus disease 2019 (COVID-19) caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) presents a global public health threat. Most research on therapeutics against SARS-CoV-2 focused on the receptor binding domain (RBD) of the Spike (S) protein, whereas the vulnerable epitopes and functional mechanism of non-RBD regions are poorly understood. Here we isolated and characterized monoclonal antibodies (mAbs) derived from convalescent COVID-19 patients. An mAb targeting the N-terminal domain (NTD) of the SARS-CoV-2 S protein, named 4A8, exhibits high neutralization potency against both authentic and pseudotyped SARS-CoV-2, although it does not block the interaction between angiotensin-converting enzyme 2 (ACE2) receptor and S protein. The cryo-EM structure of the SARS-CoV-2 S protein in complex with 4A8 has been determined to an overall resolution of 3.1 Angstrom and local resolution of 3.4 Angstrom for the 4A8-NTD interface, revealing detailed interactions between the NTD and 4A8. Our functional and structural characterizations discover a new vulnerable epitope of the S protein and identify promising neutralizing mAbs as potential clinical therapy for COVID-19.


F1000Research ◽  
2020 ◽  
Vol 9 ◽  
pp. 371 ◽  
Author(s):  
Feroza Begum ◽  
Debica Mukherjee ◽  
Dluya Thagriki ◽  
Sandeepan Das ◽  
Prem Prakash Tripathi ◽  
...  

India has recently started sequencing SARS-CoV2 genome from clinical isolates. Currently only few sequences are available from three states in India. Kerala was the first state to deposit complete sequence from two isolates followed by one from Gujarat. On April 27, 2020, the first five sequences from the state of West Bengal (Eastern India) were deposited on GISAID, a global initiative for sharing avian flu data. In this study, we have analysed the spike protein sequences from all five isolates and also compared their similarities or differences with other sequences reported in India and with isolates of Wuhan origin. We report one unique mutation at position 723 and another at 1124 in the S2 domain of spike protein of the isolates from West Bengal only.  There was one mutation downstream of the receptor binding domain at position 614 in S1 domain which was common with the sequence from Gujarat (a state of western India).  Mutation in the S2 domain showed changes in the secondary structure of the spike protein at region of the mutation. We also studied molecular dynamics using normal mode analyses and found that this mutation decreases the flexibility of S2 domain.  Since both S1 and S2 are important in receptor binding followed by entry in the host cells, such mutations may define the affinity or avidity of receptor binding.


2022 ◽  
Vol 8 (1) ◽  
pp. 84
Author(s):  
Marilia M. Knychala ◽  
Angela A. dos Santos ◽  
Leonardo G. Kretzer ◽  
Fernanda Gelsleichter ◽  
Maria José Leandro ◽  
...  

In previous work, we developed a Saccharomyces cerevisiae strain (DLG-K1) lacking the main monosaccharide transporters (hxt-null) and displaying high xylose reductase, xylitol dehydrogenase and xylulokinase activities. This strain proved to be a useful chassis strain to study new glucose/xylose transporters, as SsXUT1 from Scheffersomyces stipitis. Proteins with high amino acid sequence similarity (78–80%) to SsXUT1 were identified from Spathaspora passalidarum and Spathaspora arborariae genomes. The characterization of these putative transporter genes (SpXUT1 and SaXUT1, respectively) was performed in the same chassis strain. Surprisingly, the cloned genes could not restore the ability to grow in several monosaccharides tested (including glucose and xylose), but after being grown in maltose, the uptake of 14C-glucose and 14C-xylose was detected. While SsXUT1 lacks lysine residues with high ubiquitinylation potential in its N-terminal domain and displays only one in its C-terminal domain, both SpXUT1 and SaXUT1 transporters have several such residues in their C-terminal domains. A truncated version of SpXUT1 gene, deprived of the respective 3′-end, was cloned in DLG-K1 and allowed growth and fermentation in glucose or xylose. In another approach, two arrestins known to be involved in the ubiquitinylation and endocytosis of sugar transporters (ROD1 and ROG3) were knocked out, but only the rog3 mutant allowed a significant improvement of growth and fermentation in glucose when either of the XUT permeases were expressed. Therefore, for the efficient heterologous expression of monosaccharide (e.g., glucose/xylose) transporters in S. cerevisiae, we propose either the removal of lysines involved in ubiquitinylation and endocytosis or the use of chassis strains hampered in the specific mechanism of membrane protein turnover.


Sign in / Sign up

Export Citation Format

Share Document