Prediction analysis of porcine AXL protein as a potential receptor for SARS-CoV-2.

Dear editor,

Recent published study reported that angiotensin-converting enzyme 2 (ACE2) of various mammals may contribute to cross-species transmission of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). In fact, there is strong evidence that SARS-CoV-2 is a zoonotic spillover from bats into human populations through bridging hosts. Therefore, ongoing investigation on possible intermediate hosts would facilitate understanding the mechanisms of cross-species transmission of SARS-CoV-2 with potential benefits for the control of the coronavirus disease 2019 (COVID-19).

Notably, the receptor is the crucial host restrict factor in SARS-CoV-2 infection. ACE2 is proved to serve as a primary factor for the cellular entry of SARS-CoV-2, which may determine the host tropism. , However, the low-abundance expression of ACE2 in human lung and trachea indicated that additional cellular mediators may exist that promote cellular uptake of SARS-CoV-2. Recently research demonstrated that the tyrosine-protein kinase receptor UFO (AXL) is a novel candidate receptor, especially in the respiratory system. Since pigs as important livestock species and the huge size of their global population, it is of great significance to reveal the putative susceptibility of swine to SARS-CoV-2. Although previous data suggested that porcine ACE2 could bind SARS-CoV-2 spike (S), the interaction between porcine AXL and S that predict infection capacity is still unclear. Here, the amino acid (aa) sequences of AXLs in different species containing the residues interfacing with the N-terminal domain (NTD) of S were aligned to screen the critical sites by ESPript 3.0 (https://espript.ibcp.fr/ESPript/cgi-bin/ESPript.cgi). According the structure of human AXL-S NTD complex, P57–58, E59, H61, I68, E70, E85, F113, and G115 are located in the interface and interact with S NTD. As illustrated in Fig. 1 A, the key residues of porcine AXL were consistent with those of mink and ferret, which are efficiently infected with SARS-CoV-2, and showed weaker affinity characteristics than human AXL. Among these sites, T61H and G85E substitutions may reduce the binding affinities, while V68I and E115G mutations would not disrupt the formation of hydrogen bonds. Specifically, G116H substitution existing in hamster, mouse, mink, ferret, and pig AXLs result in a new interaction on the interface, which may have little effect on interfacial binding properties. The distribution of electrostatic potential also confirmed that the more identical pattern in AXL-S interface between pig and ferret, which appears to be significantly different from mouse (Fig. 1B). In consequence, the local structural difference of AXLs among different species suggested that porcine AXL may exhibit the similar binding affinity with S NTD to those of mink and ferret.

Fig. 1

Alignment and surface potential analysis of crucial amino acids in AXL proteins. (A) Comparative analysis of the residues of AXLs at the interface binding to spike of SARS-CoV-2 from human (GenBank accession no. NP_068713), rhesus (GenBank accession no. EHH30062.1), hamster (GenBank accession no. XP_040590606), mouse (GenBank accession no. NP_033491), ferret (GenBank accession no. XP_004776133), mink (GenBank accession no. XP_044113292.1), and pig (GenBank accession no. NP_001121930). The ALX residues at position 68 and 115 are marked in cyan triangles, position 61, 85, 113, and 116 are labeled with purple triangles. (B) Surface potential diagram of interface zone of AXLs. The structural superposition of the AXL region 27–128 from human (green, PDB code 4RA0), mouse (cyan, Uniprot entry Q0093), ferret (gray), and pig (wheat) is in the center. The AXL structures of pig and ferret were modeled using the homology models of human AXL (PDB code 4RA0 and 5VXZ, respectively) as the templates by SWISS-MODEL (https://swissmodel.expasy.org/). The six key differential residues of AXL interacting with spike protein of SARS-CoV-2 are represented by yellow sticks in the structural superposition, the details of which are further circled by dash line in the surface electrostatic potential maps of each AXL. The electrostatic potential color range is -/+5.

Interestingly, the experimental infection results in two independent studies did not support that pig is susceptible to SARS-CoV-2 infection, , in contrast to previous ACE2 results and the above AXL analysis associated with the spike-binding peculiarity. , Nevertheless, several porcine cell lines are identified to be permissive for SARS-CoV-2 infection. Owing to the limited pig breeds and numbers, age and health condition should been taken into account in future study, which may combine the pre-existing illnesses to develop typical COVID-19 symptoms. In addition, SARS-CoV-2 possesses the largest RNA genome, and evolves quickly to generate genetic diversity. The newly emerging variants omicron and delta have stronger ACE2 binding affinity and higher contiguous than these predecessors, with the mild clinical manifestations. To address the impact of the certain changes in S NTD of different SARS-CoV-2 strains, comparative analysis of aa sequences was performed. The result in Fig. 2 uncovered that the key residues K147, K150, W152, R246, S247, P251, and S256 in the AXL-S NTD complexes were unaltered, while multiple residue deletions and substitutes occurred in the same positions between two variants. More importantly, whether these aa differences would have an influence on the ability of SARS-CoV-2 to cause COVID-19 related diseases in pigs is urgently needed to be explored. Once the SARS-CoV-2 variants enhance the infectivity and transmission in pigs, the outcome for evolution of zoonotic potential would become more unpredictable and complicated when these variants are introduced into swine herds.

Fig. 2

Sequence alignment of spike protein from three SARS-CoV-2 strains including Wuhan-Hu-1 (GenBank accession no.YP_009724390), delta variant (GenBank accession no. QYM88683), and omicron variant (GenBank accession no. UFS23237). The important residues in NTD of S interacting with AXL are marked with green stars.

In summary, alignment of AXL from different species and S NTD of three representative SARS-CoV-2 strains would be conducive to quest for the potential recipient hosts in zoonotic transmission. Considering pig as the dominant livestock and the idea human biomedical model, a specific focus on interaction between porcine AXL and S of SARS-CoV-2 variants may also help prevent the risk of next pandemic and food supply stability.

Declaration of Competing Interest

The authors declare no conflict of interest.