key: cord-1025487-z05sks20 authors: Sirakov, I.; Bakalov, D.; Popova, R.; Mitov, I. title: Analysis of Host Cell Receptor GRP78 for Potential Natural Reservoirs of SARS-CoV-2 date: 2020-09-03 journal: J Epidemiol Glob Health DOI: 10.2991/jegh.k.200806.001 sha: ce1d7aeb55e2b75e380cbe56cbf5894f1fe0ae6c doc_id: 1025487 cord_uid: z05sks20 nan Since 31 December 2019, over 9 million cases of SARS-CoV-2 have been reported on five continents. During the 2003 SARS outbreak, the reported burden of asymptomatic infections was 7.5% among healthcare workers, and 13% among SARS-CoV positive cases [1] . In the current pandemic (COVID-19), the percentage of asymptomatic carriers is even higher up to 56.5% [2] , posing a serious challenge to controlling the spread. A key factor in the infectious process is the interaction between the viral surface receptor(s) and the host-cell receptor(s). A study using combined molecular docking and structural bioinformatics has demonstrated that the SARS-CoV-2 spike protein can interact with the Glucose Regulating Protein 78 (GRP78) on the surface of human respiratory cells [3] . The same approach could shed light on the possible interactions between SARS-CoV-2 spike protein and GRP78 in different animal species. We hypothesize that some animals that come into close contact with humans might become a reservoir for the disease. To explore this hypothesis, we studied the structure of GRP78 protein in different species and compared it with that in humans. Identifying any potential ecological reservoirs of the SARS-CoV-2 virus could help predict possible future outbreaks. In our research, we retrieved the amino-acid sequences of the GRP78 proteins in 84 mammals, amphibians and birds from the National Center for Biotechnology Information (NCBI) database ( Figure 1 ). For data analysis, we applied statistical Maximum Likelihood (ML) analysis with MEGA X software, Cn3Dincluding the online version in NCBI, and Phyre2 [4] . In the phylogenetic analysis, we used human DnaJ homolog subfamily C member 3 protein (NCBI accession number Q13217.1) as an outgroup; and the differences between bovine and human GRP78 protein, determined by Yang et al. [5] , as an ingroup reference. The human GRP78 deposited by Yang et al. [5] under sequence number P11021.1. Multiple sequence alignment and phylogenetic analysis identified several taxa ( Figure 1 ). The species with the highest similarity to the human GRP78 protein were apes from the genera Macaca, Nomascus, Pongo and Piliocolobus. This similarity was also present in the 3D visualization of the proteins. In humans, GRP78 consists of two domains: a nucleotide-binding domain and a Substrate-binding Domain (SBD). SBD has two subdomains: SBDa, a 10-kDa subdomain, and SBDb, a 25-kDa subdomain. The active site of the protein consists of amino acid residues located in SBD: L 1,2 (I426, T428, V429, V432, T434) and L 3,4 (F451, S452, V461, and I463), which interact with the substrate through Van der Waals forces [5] . Ibrahim et al. [3] described the interaction between the spike protein of SARS-CoV-2 and cell-surface GRP78 SBDb with the corresponding amino-acid residues [3] . While analyzing the structure of GRP78 in the phylogenetically similar species, using the data reported by Ibrahim et al. [3] in their Table 1 , we found complete congruence with some of the amino acid residues in human GRP78. Other positions did not have a corresponding match in the human reference sequences of GRP78, e.g.: residues forming hydrogen bonds in positions G489 (P), V453 (T), V490 (A), G489 (P), G454 (A), We associate the lack of matches in these positions with possible technical errors, due to the immense amount of data analyzed by the authors; else the role of GRP78 receptors would not be as significant as thought. We found differences in SBD (TASDNQP → TASDTQP) in Macaca fuscata. However, this modification in L 3,4 does not have a significant effect on the polypeptide-binding pocket of the protein [6] . Our analysis suggested that domestic animals are not likely to act as an epidemiological factor due to significant differences in the structure of their GRP78 protein. The species that stood out as potential SARS-CoV-2 reservoirs belong to the so-called Old World Мonkey (OWM) family Cercopithecidae, with a broad distribution in Asia (Macaca mulatta), Southeast Asia (Macaca fascicularis), Japan (M. fuscata), and Africa -Tanzania, Uganda, Democratic Republic of Congo, Rwanda, Burundi (Piliocolobus tephrosceles). Potential reservoirs other than OWM are Nomascus leucogenys, found in Vietnam and Laos, and Pongo abelii, from Sumatra (superfamily Hominoidea). The animals from the group of OWM are one of the important reservoirs for Lentiviruses and have played an essential role in the evolution of HIV [7] . In superfamily Hominoidea, there are reports about interspecies transmissions between man and some closely-related species (Pan, Gorilla, Pongo and Hylobates), e.g. evidence of HBV of gibbon and human origin in chimpanzees [8] . Based on these data and our results, we hypothesize that these species could act as potential reservoirs for SARS-CoV-2. Interestingly, in our analysis, Manis javanica did not demonstrate a significant risk of becoming a reservoir for SARS-CoV-2 despite the reported link between Pangolin-CoV and SARS-CoV-2 [9] . This could be because Li et al. [9] show the probable path of crossspecies transmission of the virus from bat or pangolin and its adaptation to the human body. After the virus has adapted to the new host, it can become a single-host pathogen infecting humans only or it can also infect species other than the ones the virus has originated from. An example of the latter evolutionary scenario is the human measles virus, which has resulted from a cattle-to-human host-jump of rinderpest virus [10] or HBV. In our study, we discuss the possibility of transmission of SARS-CoV-2 from humans to primates to humans. Our findings based on GRP78 protein analysis corroborated with published reports support the hypothesis that some primate species -but not likely domesticated animals -could serve as a reservoir for SARS-CoV-2. Surveillance of those potential reservoirs for SARS-COV-2 can help prevent future outbreaks. The authors declare they have no conflicts of interest. This research did not receive any specific grant from funding agencies in the public, commercial, or not-for-profit sectors. Asymptomatic SARS coronavirus infection among healthcare workers Asymptomatic and presymptomatic SARS-CoV-2 infections in residents of a long-term care skilled nursing facility COVID-19 spike-host cell receptor GRP78 binding site prediction The Phyre2 web portal for protein modeling, prediction and analysis Close and allosteric opening of the polypeptide-binding site in a human Hsp70 chaperone BiP Conformation transitions of the polypeptide-binding pocket support an active substrate release from Hsp70s Evolutionary and functional analysis of old world primate TRIM5 reveals the ancient emergence of primate lentiviruses and convergent evolution targeting a conserved capsid interface Hepatitis B virus lineages in mammalian hosts: potential for bidirectional cross-species transmission Evolutionary history, potential intermediate animal host, and cross-species analyses of SARS-CoV-2 Mapping the evolutionary trajectories of morbilliviruses: what, where and whither We would like to acknowledge Prof. T. Strateva (Department "Medical Microbiology", Medical University -Sofia) for her invaluable advice and guidance while writing the manuscript.