Skip to main content

ORIGINAL RESEARCH article

Front. Genet., 23 October 2020
Sec. Computational Genomics
This article is part of the Research Topic Coronavirus Disease (COVID-19): Molecular Mechanisms, Translational Approaches and Therapeutics View all 118 articles

Predicting the Animal Susceptibility and Therapeutic Drugs to SARS-CoV-2 Based on Spike Glycoprotein Combined With ACE2

\r\nMin ShenMin Shen1Chao LiuChao Liu1Run XuRun Xu1Zijing RuanZijing Ruan2Shiying ZhaoShiying Zhao1Huidong ZhangHuidong Zhang3Wen WangWen Wang3Xinhe HuangXinhe Huang2Li YangLi Yang1Yong TangYong Tang1Tai Yang*Tai Yang4*Xu Jia*Xu Jia1*
  • 1Non-coding RNA and Drug Discovery Key Laboratory of Sichuan Province, Chengdu Medical College, Chengdu, China
  • 2School of Life Sciences and Engineering, Southwest Jiaotong University, Chengdu, China
  • 3West China School of Public Health, Sichuan University, Chengdu, China
  • 4School of Pharmacy, Chengdu Medical College, Chengdu, China

Recently, a few animals have been frequently reported to have been diagnosed with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Whether they are SARS-CoV-2 intermediate hosts is worthy of great attention. The interaction of SARS-CoV-2 spike protein and its acceptor protein ACE2 is an important issue in determining viral host range and cross-species infection, while the binding capacity of Spike protein to ACE2 of different species is unknown. Here, we used the atomic structure model of SARS-CoV-2 and human ACE2 to assess the receptor utilization capacity of ACE2s from 10 kinds of animals. Results show that chimpanzees, domestic cats and cattles are more susceptible to infection by SARS-CoV-2. Cats in particular, such as pet cats and stray cats, interact very closely with humans, implying the necessity to carefully evaluate the risk of cats during the current COVID-19 pandemic. Furthermore, based on ACE2(cats)-SARS-CoV-2-RBD model, through high-throughput screening methods using a pool of 30,000 small molecules, eight compounds were selected for binding free energy calculations. All the eight compounds can effectively interfere with the binding of ACE2 and Spike protein, especially Nelfinavir, providing drug candidates for the treatment and prevention of SARS-CoV-2, suggesting further assessment of the anti-SARS-CoV-2 activity of these compounds in cell culture. Although we only reported the results of the simulation, and more laboratory and epidemiological investigation are required. Like cats are a risk factor, we can further detect SARS-CoV-2 according to the susceptibility of different animals, find the potential host of infection, and completely cut off the living space of the virus. Especially, cats could be a choice of animal model for screening antiviral drugs or vaccine candidates against SARS-CoV-2.

Introduction

By June 14, 2020, the COVID-19 outbreak in China would have been well under control. SARS-CoV-2 was suddenly detected in Xinfadi seafood market in Beijing, which makes it more urgent to assess the susceptibility of animals in close contact with humans. SARS-CoV-2 has spread all over the world, and COVID-19 has posed a serious global public health emergency (Zhu et al., 2020). On one hand, the homology of SARS-CoV-2 and bat SARS-like coronavirus (bat-sl-covzc45) is higher than 85% (Boni et al., 2020; Lu et al., 2020), on the other hand, the virus is unlikely to spread directly from bats to humans due to the lack of direct contact, meaning there is an unknown intermediate host. Although bats are likely reservoir hosts for SARS-CoV-2, the identity of any intermediate host that might have facilitated transfer to humans is unknown.

Recently, pet dogs and cats have been tested to be positive for SARS-CoV-2 infection (https://finance.sina.com.cn/7x24/2020-03-19/doc-iimxxsth0356081.shtml; https://news.sina.com.cn/w/2020-03-28/doc-iimxyqwa3732342.shtml). On March 5, 2020, samples from the mouth and nose of a pet dog in Hong Kong showed SARS-CoV-2 positive but no symptom of COVID-19. The dog’s owner is a COVID-19 patient who was diagnosed in Hong Kong on February 25, 2020. Experts believe the dog was probably infected by humans, and the dog’s low level of infection with a novel coronavirus could be the first case in the world (https://finance.sina.com.cn/7x24/2020-03-19/doc-iimxxsth0356081.shtml). On March 27, 2020, Belgian public health authorities announced a pet cat in liege had been diagnosed with COVID-19. Researchers found SARS-CoV-2 in the feces of the cat, which was suffering from breathing difficulties. Belgian health officials stated subsequently the cat was infected with a novel coronavirus from its owner, and this is the first confirmed case of a pet infection in Europe (https://news.sina.com.cn/w/2020-03-28/doc-iimxyqwa3732342.shtml). Back to 2003, studies about SARS have shown that ferrets (Mustela furo) and domestic cats (Felis domesticus) were susceptible to infection by SARS-CoV and that they can efficiently transmit the virus to previously uninfected animals that are housed with them (Martina et al., 2003). Meanwhile, one laboratory study reported that cats could be infected with SARS-CoV-2 (Shi et al., 2020), and another reported that cats in Wuhan had been infected with SARS-CoV-2 during the outbreak (Zhang et al., 2020).

The specificity of the interaction between virus and its receptor is an important issue in regulating both the cross-species and human-to-human transmissions infection. SARS-CoV-2 invades cells mainly through Spike protein recognition of host cell receptor angiotensin-converting enzyme 2 (ACE2), such binding triggers a cascade events leading to the fusion between cell and viral membranes for cell entry (Ortega et al., 2020; Wan et al., 2020). ACE2 is expressed in most mammals, but not all ACE2 can be utilized by SARS-CoV-2 as the receptor (Boni et al., 2020). The utilization of ACE2 by SARS-CoV-2 can rapidly screen and narrow the range of intermediate hosts of SARS-CoV-2. Currently, it is not clear which mammals are involved in the evolution of SARS-CoV-2 and which animals are infected by SARS-CoV-2.

Like SARS-CoV, the evolutionary rate of SARS-CoV evolution were analyzed by modeling with genetic method, the results showed that SARS-CoV may spread in the bats and transmitted to other hosts such as pangolins, and have a tendency to spread across species (O’Brien et al., 2006; Le Poder, 2011; Boni et al., 2020). Objects or environments that humans and animals come into contact with can be hidden carriers of the virus. There are no clinical drugs or vaccines for SARS-CoV-2, and most experts agree that the epidemic is only under control if SARS-CoV-2 is contained in poor countries. Controlling the source of SARS-CoV-2 will therefore be a protracted battle.

In this study, we use the reported crystal structure of Spike protein binding to human ACE2 to construct 3D models of 10 kinds of animals by homologous modeling, among which 7 kinds of animals are easy to contact with humans in daily life, and the other 3 kinds of animals are commonly used model animals in the laboratory. According to the algorithm, the amino acid positions, types and sizes of the interactions between proteins were calculated, and the binding free energy was finally calculated. At the same time, high-throughput screening was used to screen out drugs that might interfere with the combination of SARS-CoV-2-RBD and ACE2 of cats. Our results may help to screen out SARS-CoV-2 intermediate hosts and to figure out the transmission model of SARS-CoV-2 and finally to control the COVID-19 disease.

Materials and Methods

ACE2 Sequences Alignment

The amino acid sequences of the extracellular domains of these eleven ACE2 (residues Ser19-Asp615): Homo sapiens (human), Pan troglodytes (chimpanze), Macaca mulatta (Rhesus monkey), Felis catusc (domestic cat), Equus caballus (horse), Oryctolagus cuniculus (rabbit), Canis lupus familiaris (dog), Sus scrofa (pig), Ovis aries (sheep), Bos taurus (cattle), Mus musculus (house mouse) were downloaded from GenBank (Supplementary File S1). Protein sequence alignments were done using Clustal Omega.

Protein–Protein Docking

Homology model of the target ACE2 protein was built by modeller9.18 using crystal structure of human ACE2 (PDB: 1R42) as template. Hundred independent structures were constructed and the one with best DOPE score was chosen for further energy minimization in Amber18 using ff14SB force field. Rosetta3.7 was used to perform protein-protein docking to get the ACE2-SARS-CoV-2-RBD complex, and default values are used for parameter setting. The PDB of the SARS-CoV-2 spike protein is 6VSB (Lan et al., 2020).

Molecular Dynamics Simulation

The ACE2-SARS-CoV-2-RBD complex was immersed in an octahedron box of TIP3P water that was extended by 10 Å from the solute. In the simulation system, the amount of Na + is exactly equal to the amount of negative charge, so as to neutralize the system charge and make the total charge zero. Amber ff14SB force field was used to parameterize the protein. 10,000 steps of minimization with constraints (10 kcal/mol/Å2) on heavy atoms of complex, including 5,000 steps of steepest descent minimization and 5,000 steps of conjugate gradient minimization, was used to optimize each system. Then each system was heated to 300 K within 50 ps followed by 50 ps equilibration in NPT ensemble. Finally, 100 ns MD simulation on each system at 300 K was performed. The minimization, heating and equilibrium are performed with sander program in Amber18. The 100 ns production run was performed with pmemd.cuda. Based on the 100 ns MD simulation trajectory, binding free energy (ΔG) between spike and viral ACE2 receptor was calculated with MM/GBSA method (The MMGBsa.py module in the Amber program was used for the MM/GBSA calculations)according to the following equation (Genheden and Ryde, 2015):ΔG = ΔHTΔS = ΔEele + ΔEVDW + ΔGgb + ΔGnpTΔS. Where ΔEele and ΔEVDW refer to electrostatic and van der Waals energy terms, respectively. ΔGgb and ΔGnp refer to polar and non-polar solvation free energies, respectively. Conformational entropy (TΔS) was not calculated for saving time. Besides, the ligands were compared based on the same target, so it is reasonable to ignore the entropy.σ=1NI=1N(xi-u)2, SD was calculated according this equation where xi represents the binding free energy in each frame.

Construction of Small Molecular Ligands

Each sub-library (FDA, investigational-only, world-not-FDA1) was downloaded from the zinc database. The 2D structure of the compound was then converted into the corresponding 3D coordinates using the Babel server2. Then the model was converted to pdbqt format by prepare_receptor4.py script (From the AutodockTool package)with assigning atomic types and atomic charges. All rotatable bonds in the molecule are set to be flexible for flexible docking. Vina1.1.2 was used for molecular docking.

Binding Free Energy Calculation Between ACE2(cats)-SARS-CoV-2-RBD and Small Molecular Ligand

The ACE2-SARS-CoV-2-RBD complex structure of cats that obtained from Section “Protein–Protein Docking” was used for docking with small molecular ligand to calculate binding free energy. The minimum conformation of the docking was used as the initial position of the drug molecule.

Each simulation system was immersed in a cubic box of TIP3P water with 10 Å distance from the solute. The Na+ or Cl was applied to neutralize the system. General Amber force field (GAFF) 15 and Amber ff14SB force field were used to parameterize the ligand and protein respectively. 10,000 steps of minimization with constraints (10 kcal/mol/Å2) on heavy atoms of complex, including 5,000 steps of steepest descent minimization and 5,000 steps of conjugate gradient minimization, was used to optimize each system. Then each system was heated to 300 K within 0.2 ns followed by 0.1 ns equilibration in NPT ensemble. Finally, 5 ns MD simulation on each system at 300 K was performed. The minimization, heating and equilibrium are performed with sander program in Amber18. The 5 ns production run was performed with pmemd.cuda. Based on the 5 ns MD simulation trajectory, binding free energy (ΔG) was calculated with MM/GBSA method according to the following equation: ΔGcal = ΔH-TΔS = ΔEvdw + ΔEele + ΔGgb + ΔGnp-TΔS, where ΔEele and ΔEvdw refer to electrostatic and van der Waals energy terms respectively. ΔGgb and ΔGnp refer to polar and non-polar solvation free energies respectively. Conformational entropy (TΔS) was not calculated for saving time. Besides, the ligands were compared based on the same target, so it is reasonable to ignore the entropy.

Results and Discussion

Comparison of Amino Acid Sequences of ACE2 Among 10 Kinds of Animals

In order to analyze the possibility of SARS-CoV-2 infection on mammals that humans may come into contact with in daily life, we collected the amino acid sequence of ACE2 of pet cats, dogs, livestock cattle, horses, sheep and pigs, as well as laboratory animals Rhesus monkey, chimpanzees, mice and rabbits in Genbank. We aligned the amino acid sequence of the extracellular domains of these eleven ACE2, the results are shown in Table 1. Compared with the ACE2 of human, the homology of amino acid sequence was between 95.0%∼67.4%. It’s worth noting that, as a commonly used model animal in the laboratory, mice had the lowest homology. That is to say, the N-terminal peptidase domain of ACE2 in these animals is very similar to that of the known human hosts of infection. Chimpanzees and Rhesus monkey seem to be the most susceptible to SARS-CoV-2, while mice is the least sensitive.

TABLE 1
www.frontiersin.org

Table 1. Comparison of amino acids homology of ACE2 between human and other organisms.

Binding Affinity (ΔG) Values of the Interaction Between Spike and Viral ACE2 Receptor

To further analyze the possibility of infection in these animals, we determined the crystal structure of the SARS-CoV-2 spike receptor-binding domain (RBD) bound with the cell receptor ACE2, calculated the binding free energy. The predicted results are shown in Table 2, chimpanzees have the highest binding affinity, even higher than human, while gradually decreases in order of cats, cattle, Rhesus monkey, dogs, pigs, horses, sheep, mice, and rabbits. Although they belong to different species, they all have the binding affinity of the interaction. The binding affinity of cats and chimpanzees are very similar to human, while rabbits and mice are the lowest. Higher affinity values might be related to the dynamic of infection and the rapid spread observed for this virus. These data suggested that the higher binding affinity of RBD of coronavirus to ACE2 will confer the virus higher infectivity and pathogenicity (Ortega et al., 2020; Xu et al., 2020).

TABLE 2
www.frontiersin.org

Table 2. Binding affinity (ΔG) predicted values for the interaction between spike and viral ACE2 receptor.

Computer modeling of interaction between SARS-CoV-2 RBD and ACE2 has identified some residues potentially involved in the actual interaction. Spike protein contacts with the helical structure of 19-83aa of human ACE2 and the folding structure which in the 347-358aa region of human ACE2, producing intermolecular interactions. The results are shown in Table 3 and Figure 1. Structural analysis revealed a total of 11 residues of the SARS-CoV-2 RBD contact 13 residues of the human ACE2, and there are 15 hydrogen bonds at the SARS-CoV-2 RBD/ACE2 interface. It is found that cat and human ACE2 are very similar when it binds to RBD, a total of 11 residues of the RBD contact 12 residues of the ACE2, 14 hydrogen bonds at the SARS-CoV-2 RBD/ACE2 interface. Structural analysis revealed that most of the cats-critical ACE2 binding residues in ACE2-SARS-CoV-2-RBD had highly conserved or side chain characteristics similar to those in humans (Table 3 and Figure 2). Taken together, these results show that the SARS-CoV-2 RBD/ACE2 of cats interfaces share substantial similarity in the number of interacting residues, and hydrophilic interaction networks. In the study of SARS, the infection experiments on many kinds of animals show that ferrets and domestic cats can be used as potential animal infection models such as vaccine and drug screening (Martina et al., 2003). This is consistent with the results of recent virus infection experiments: ferrets and cats have effective replication ability. Virus RNA was found in the nose, soft jaw, tonsil and small intestine of ferrets, but no virus was detected in other organs, which proved that SARS-CoV-2 was only replicated in the upper respiratory tract (Chen, 2020). A research in human body have been reported that: ACE2 receptor is highly expressed in human nasal cells, but not detected in lung cells (Wu et al., 2020). Meanwhile, the expression of ACE2 in the small intestine is high, which is consistent with the recently reported gastrointestinal tract as a potential route of SARS-CoV-2 infection (Lu et al., 2020; Wang et al., 2020). Compared with human, ACE2 and SARS-CoV-2 binding affinity of mice are very low, similar interacting residuals were fewer, and hydrophobic interaction networks are weaker. Therefore, SARS-CoV-2 may replicate inefficiently in mice and rats, ruling them out as animal models to test vaccine or antiviral drugs candidates against SARS-CoV-2.

TABLE 3
www.frontiersin.org

Table 3. Potential interaction between the S protein receptor binding domain (RBD) and ACE2 of a variety of animals.

FIGURE 1
www.frontiersin.org

Figure 1. Comparison of amino acids homology of ACE2 in the helical structure and the folding between human and other organisms. (a) The helical structure of 19–83 aa. (b) The folding structure of 347–358 aa. The different amino acids are shown in yellow.

FIGURE 2
www.frontiersin.org

Figure 2. Overall structure of SARS-CoV-2 RBD bound with ACE2. (A–F) Are respectively human, dogs, chimpanzees, cattles, cats, Rhesus monkey. ACE2 is colored light coffee. SARS-CoV-2 RBD core is colored cyan. Amino acid interaction sites are also shown.

Docking Results of 7496 Drugs Against ACE2(Cats)-SARS-CoV-2-RBD Model

The 7,496 drugs obtained from the zinc database were screened for molecular docking. Among them, 57 of the 2,100 compounds approved by the FDA have a docking score better than −8.0 kcal/mol. There were 4,264 compounds approved by regulatory agencies other than the FDA. Among them, there were 124 compounds with the docking score better than −8.0 kcal/mol. 67 of the 1132 compounds in clinical trials have a docking score better than −8.0 kcal/mol. Further, eight top compounds showed the docking score in a range of −8.1 to −9.9 kcal/mol were selected from docking results of homology model (Table 4).

TABLE 4
www.frontiersin.org

Table 4. Eight drugs selected from the ACE2(cats)-SARS-CoV-2-RBD model.

Docking Results of Eight Drugs Against ACE2(Cats)-SARS-CoV-2-RBD Model

So far, only Remdesivir has been approved for the treatment of COVID-19 in Japan, and in the world no corresponding vaccine has been approved. For these confirmed cases, they need antiviral drugs or other drugs to treat SARS-CoV-2 infection. More potentiality compounds need to be tested. In the presence of ACE2-RBD complex formation, that is to say, when symptoms of infection may already be present, the interference of small molecules with ACE2(cats)-SARS-CoV-2-RBD complex can be assessed by assessing the likelihood of small molecules entering the lumen of the complex. In order to reduce the cost of compounds screening, eight compounds were selected based on virtual screening and docking scores and its showed as Table 5.

TABLE 5
www.frontiersin.org

Table 5. The calculated binding energies of ligand to ACE2(cats)-SARS-CoV-2-RBD.

Nelfinavir is a clinically important antiviral drug, which can inhibit the production of human immunodeficiency virus (HIV) and it can inhibits the Herpes Simples Virus 1(HSV-1) (Kalu et al., 2014; Gantt et al., 2015). Our docking results show that five hydrogen bonds, including GLU-30, ASN-33, PRO-391 and LYS-417, are maintained on the binding of nelfinavir with the ACE2 and Spike glycoprotein complex. Among them, GLU-30 and nelfinavir maintain two hydrogen bonds. In the previous studies, some scientists have demonstrated the antiviral effect of nefinavir using Vero cell lines infected with SARS CoV (Yamamoto et al., 2004; Hsieh et al., 2010).

Saquinavir is the first protease inhibitor used to treat patients with HIV infection (Figgitt and Plosker, 2000). One hydrogen bonds involving THR-27 and van der Waals energy maintained upon the binding of saquinavir with ACE2 and Spike glycoprotein complex interface. Previous studies shown that saquinavir could inhaibit the RNA-dependent RNA polymerase (RdRP) activity (Beck et al., 2020).

Tegobuvir and Setrobuvir, Non-Nucleoside Analog inhibitors of the Hepatitis C virus(HCV) (Zeuzem et al., 2012; Mallalieu et al., 2014; Wyles et al., 2014). Tegobuvir, the hydrogen bond involving ARG-403 maintained upon binding of tegobuvir with ACE2 and Spike glycoprotein complex interface, with additionally van der Waals energy, but for Setrobuvir, it was mainly combined with ACE2 and Spike glycoprotein complex interface through van der Waals potential energy.

Our simulation results show that these four antiviral drugs can enter the interface of ACE2 and Spike glycoprotein complex, and they all have antiviral activity and can treat viral infections (Figure 3). Therefore, we speculate that these four drugs have potential activity for the treatment of COVID-19. But this conclusion needs to be needs to be further verified by in vivo experiments.

FIGURE 3
www.frontiersin.org

Figure 3. Binding model of antiviral drugs with ACE2 and Spike glycoprotein complex interface. Antiviral drug binding model with ACE2 and Spike glycoprotein complex interface. (A,C,E,G) are the binding models of ACE2 and Spike glycoprotein complexes with antiviral drugs, Olysio (cyan) and related residues (off-white), (B,D,F,H) are compound Olysio (blue). The interaction-green model is located at the interface of the Spike glycoprotein complex (between coffee and white surface). The yellow dotted line indicates the interaction distance (A). (A,C,E,G) or (B,D,F,H) are nelfinavir, saquinavir, tegobuvir, and setrobuvir, respectively.

Dihydroergotamine, Dihydroergotoxine and Ergotamine are the synthetic drugs developed in the 20th century for treating migraine (Metra et al., 1995; Villamil-Hernandez et al., 2014; Connelly et al., 2020) (Figure 4). Our docking results show that one hydrogen bonds involving LYS-26 and van der Waals energy maintained upon the binding of DHE with the interface of ACE2 and Spike glycoprotein complex. The interface between DDE and ACE2 and Spike glycoprotein complex maintains a hydrogen bond involving GLU-37 and Van der Waals energy. As for ergotamine, when ergotamine binds to the ACE2 and Spike glycoprotein complex, three hydrogen bonds involving GLU-23, GLU-30 and LYS-417 are maintained. Our simulation results further confirmed that DHE, DDE and ERG can bind to the interface of ACE2 and spike glycoprotein complex.

FIGURE 4
www.frontiersin.org

Figure 4. The binding model of non-antiviral drugs with the interface of ACE2 and Spike glycoprotein complexes. Antiviral drug binding model with ACE2 and Spike glycoprotein complexes interface. (A,C,E,G) are the binding models of ACE2 and Spike glycoprotein complexes with antiviral drugs, Olysio (cyan) and related residues (off-white), (B,D,F,H) are compound Olysio (blue). The interaction-green model is located at the interface of the Spike glycoprotein complex (between coffee and white surface). The yellow dotted line indicates the interaction distance (A). (A,C,E,G) or (B,D,F,H) are dihydroergotamine, dihydroergotoxine, ergotamine, fiduxosin, respectively.

Fiduxosin, one of the α1-Adrenocepter Antagonists, is used to treat the lower urinary tract symptoms (LUTS) (Brune et al., 2002; Hancock et al., 2002) (Figure 4). Our ducking results showed that one hydrogen bonds involving THR-421 and van der Waals energy maintained upon the binding of fiduxosin with the interface of ACE2 and Spike glycoprotein complex. Thus making them as candidates for further in vitro evaluation of anti-SARS-CoV-2 activity.

Binding Free Energy Calculated by MM/GBSA

Through the simulation trajectory of 100 ns molecular dynamics simulations, we calculated the binding free energy of five drugs by MM/GBSA methods. The bonding free energy for the interface of ACE2 and spike glycoprotein of eight compounds were calculated. Nelfinavir has the strongest binding free energy, suggesting it can be tested their anti-SARS-CoV-2 infection in vitro (Table 5).

Conclusion

Based on the potential interaction between S protein and mammalian ACE2, it was speculated that SARS-CoV-2 preserved the ability to infect many mammals including chimpanzees, cats, cattle and Rhesus monkey. Cats’ binding free energy, as well as key amino acids, are highly similar to humans, meaning they could serve as model animals for developing vaccines and drugs, implicating these animal species as possible intermediate hosts or animal models for SARS-CoV-2 infections. About 60 percent of people infected with the virus are asymptomatic carriers and bats are a special case of SARS asymptomatic carriers, so these animals can carry the virus and infect people (Case Study Shows Asymptomatic Transmission of COVID-19 in China; Kupferschmidt, 2020). Modeling with genetic method, according to the HcoV-OC43 and MERS-CoV the evolutionary rate of SARS-CoV evolution are analyzed, the results showed that SARS-CoV may spread in the bats is transmitted to other hosts such as pangolins, dozens of years later it was the third time in the 17 years of coronary virus outbreak, will probably have a virus across species boundaries (O’Brien et al., 2006; Le Poder, 2011; Boni et al., 2020). Considering the widespread of stray cats, wildlife markets and stock farms in Wuhan, it was not strange that these animals could serve as potential intermediate hosts of SARS-CoV-2. Therefore, pets exposed to the patient should be screened for SARS-CoV-2. During the process of epidemic prevention, we should prevent the predict possible of zoonosis event or cross-infection in the future. However, these are still preliminary results predicted by sequence analysis, and more laboratory and epidemiological investigation are required.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

MS, YT, and XJ designed the study. MS, CL, and RX performed most studies, but a few were carried out by ZR, SZ, HZ, WW, XH, and LY. All authors contributed thoughts and advice. MS and CL analyzed and interpreted the data. MS wrote the text, and the other authors contributed to the final text presentation. All authors have approved the submission.

Funding

This work was supported by grants from the National Natural Science Foundation of China (Nos. 31870135 and 31600116), the “1000 Talent Plan” of Sichuan Province (Nos. 980 and 1060), and the Research Fund of Chengdu Medical College (CYZ18-03).

Conflict of Interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Acknowledgments

This manuscript has been released as a pre-print at: https://www.preprints.org/manuscript/202004.0116/v1 (Shen et al., 2020).

Footnotes

  1. ^ http://zinc.docking.org/
  2. ^ http://openbabel.sf.net

References

Beck, B. R., Shin, B., Choi, Y., Park, S., and Kang, K. (2020). Predicting commercially available antiviral drugs that may act on the novel coronavirus (SARS-CoV-2) through a drug-target interaction deep learning model. Comput. Struct. Biotechnol. J. 18, 784-790. doi: 10.1016/j.csbj.2020.03.025

PubMed Abstract | CrossRef Full Text | Google Scholar

Boni, M. F., Lemey, P., Jiang, X., Tsan-Yuk Lam, T., Perry, B., Castoe, T., et al. (2020). Evolutionary origins of the SARS-CoV-2 sarbecovirus lineage responsible for the COVID-19 pandemic. bioRxiv [Preprint]. doi: 10.1101/2020.03.30.015008

CrossRef Full Text | Google Scholar

Brune, M. E., Katwala, S. P., Milicic, I., Witte, D. G., Kerwin, J. F. Jr., Meyer, M. D., et al. (2002). Effect of fiduxosin, an antagonist selective for alpha(1A)- and alpha(1D)-adrenoceptors, on intraurethral and arterial pressure responses in conscious dogs. J. Pharmacol. Exp. Ther. 300, 487–494. Case Study Shows Asymptomatic Transmission of COVID-19 in China. doi: 10.1124/jpet.300.2.487

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, H. (2020). Susceptibility of ferrets, cats, dogs, and different domestic animals to SARS-coronavirus-2. bioRxiv [Preprint]. doi: 10.1101/2020.03.30.015347

CrossRef Full Text | Google Scholar

Connelly, M., Sekhon, S., Stephens, D., Boorigie, M., and Bickel, J. (2020). Enhancing outpatient dihydroergotamine infusion with interdisciplinary care to treat refractory pediatric migraine: preliminary outcomes from the comprehensive aggressive migraine protocol (”CAMP”). Headache 60, 101–109. doi: 10.1111/head.13685

PubMed Abstract | CrossRef Full Text | Google Scholar

Figgitt, D. P., and Plosker, G. L. (2000). Saquinavir soft-gel capsule: an updated review of its use in the management of HIV infection. Drugs 60, 481–516. doi: 10.2165/00003495-200060020-00016

PubMed Abstract | CrossRef Full Text | Google Scholar

Gantt, S., Gachelet, E., Carlsson, J., Barcy, S., Casper, C., and Lagunoff, M. (2015). Nelfinavir impairs glycosylation of herpes simplex virus 1 envelope proteins and blocks virus maturation. Adv. Virol. 2015:687162. doi: 10.1155/2015/687162

PubMed Abstract | CrossRef Full Text | Google Scholar

Genheden, S., and Ryde, U. (2015). The MM/PBSA and MM/GBSA methods to estimate ligand-binding affinities. Expert Opin. Drug Discov. 10, 449–461. doi: 10.1517/17460441.2015.1032936

PubMed Abstract | CrossRef Full Text | Google Scholar

Hancock, A. A., Buckner, S. A., Brune, M. E., Esbenshade, T. A., Ireland, L. M., Katwala, S., et al. (2002). Preclinical pharmacology of fiduxosin, a novel alpha(1)-adrenoceptor antagonist with uroselective properties. J. Pharmacol. Exp. Ther. 300, 478–486. doi: 10.1124/jpet.300.2.478

PubMed Abstract | CrossRef Full Text | Google Scholar

Hsieh, L. E., Lin, C. N., Su, B. L., Jan, T. R., Chen, C. M., Wang, C. H., et al. (2010). Synergistic antiviral effect of Galanthus nivalis agglutinin and nelfinavir against feline coronavirus. Antiv. Res. 88, 25–30. doi: 10.1016/j.antiviral.2010.06.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kalu, N. N., Desai, P. J., Shirley, C. M., Gibson, W., Dennis, P. A., and Ambinder, R. F. (2014). Nelfinavir inhibits maturation and export of herpes simplex virus 1. J. Virol. 88, 5455–5461. doi: 10.1128/JVI.03790-13

PubMed Abstract | CrossRef Full Text | Google Scholar

Kupferschmidt, K. (2020). Study claiming new coronavirus can be transmitted by people without symptoms was flawed. Science doi: 10.1126/science.abb1524

CrossRef Full Text | Google Scholar

Lan, J., Ge, J., Yu, J., Shan, S., Zhou, H., Fan, S., et al. (2020). Crystal structure of the 2019-nCoV spike receptor-binding domain bound with the ACE2 receptor. bioRxiv [Preprint]. doi: 10.1101/2020.02.19.956235

CrossRef Full Text | Google Scholar

Le Poder, S. (2011). Feline and canine coronaviruses: common genetic and pathobiological features. Adv. Virol. 2011, 1–11. doi: 10.1155/2011/609465

PubMed Abstract | CrossRef Full Text | Google Scholar

Lu, R., Zhao, X., Li, J., Niu, P., Yang, B., Wu, H., et al. (2020). Genomic characterisation and epidemiology of 2019 novel coronavirus: implications for virus origins and receptor binding. Lancet 395, 565–574. doi: 10.1016/S0140-6736(20)30251-8

CrossRef Full Text | Google Scholar

Mallalieu, N. L., Rahimy, M. H., Crowley, C. A., Appleman, J. R., Smith, P. F., and Freddo, J. L. (2014). Pharmacokinetics and pharmacodynamics of setrobuvir, an orally administered hepatitis C virus non-nucleoside analogue inhibitor. Clin. Ther. 36, 2047.e3–2063.e3. doi: 10.1016/j.clinthera.2014.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Martina, B. E., Haagmans, B. L., Kuiken, T., Fouchier, R. A., Rimmelzwaan, G. F., Van Amerongen, G., et al. (2003). Virology: SARS virus infection of cats and ferrets. Nature 425:915. doi: 10.1038/425915a

PubMed Abstract | CrossRef Full Text | Google Scholar

Metra, M., Missale, C., Spano, P. F., and Cas, L. D. (1995). Dopaminergic drugs in congestive heart failure: hemodynamic and neuroendocrine responses to ibopamine, dopamine, and dihydroergotoxine. J. Cardiovasc. Pharmacol. 25, 732–740. doi: 10.1097/00005344-199505000-00008

PubMed Abstract | CrossRef Full Text | Google Scholar

O’Brien, S. J., Troyer, J. L., Roelke, M., Marker, L., and Pecon-Slattery, J. (2006). Plagues and adaptation: lessons from the Felidae models for SARS and AIDS. Biol. Conserv. 131, 255–267. doi: 10.1016/j.biocon.2006.05.001

PubMed Abstract | CrossRef Full Text | Google Scholar

Ortega, J. T., Serrano, M. L., Pujol, F. H., and Rangel, H. R. (2020). Role of changes in SARS-CoV-2 spike protein in the interaction with the human ACE2 receptor: an in silico analysis. EXCLI J. 19, 410–417.

Google Scholar

Shen, M., Liu, C., Xu, R., Ruan, Z., Zhao, S., Zhang, H., et al. (2020). SARS-CoV-2 infection of Cats and Dogs? Preprints 2020, 2020040116.

Google Scholar

Shi, J., Wen, Z., Zhong, G., Yang, H., Wang, C., Huang, B., et al. (2020). Susceptibility of ferrets, cats, dogs, and other domesticated animals to SARS–coronavirus 2. Science 368, 1016–1020. doi: 10.1126/science.abb7015

PubMed Abstract | CrossRef Full Text | Google Scholar

Villamil-Hernandez, M. T., Alcantara-Vazquez, O., Sanchez-Lopez, A., Gutierrez-Lara, E. J., and Centurion, D. (2014). Pharmacological evidence that 5-HT1A/1B/1D, alpha2-adrenoceptors and D2-like receptors mediate ergotamine-induced inhibition of the vasopressor sympathetic outflow in pithed rats. Eur. J. Pharmacol. 740, 512–521. doi: 10.1016/j.ejphar.2014.06.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Wan, Y., Shang, J., Graham, R., Baric, R. S., and Li, F. (2020). Receptor recognition by the novel coronavirus from wuhan: an analysis based on decade-long structural studies of SARS coronavirus. J. Virol. 94:e00127-20. doi: 10.1128/JVI.00127-20

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, J., Zhao, S., Liu, M., Zhao, Z., Xu, Y., Wang, P., et al. (2020). ACE2 expression by colonic epithelial cells is associated with viral infection, immunity and energy metabolism. medRxiv [Preprint]. doi: 10.1101/2020.02.05.20020545

CrossRef Full Text | Google Scholar

Wu, C., Zheng, S., Chen, Y., and Zheng, M. (2020). Single-cell RNA expression profiling of ACE2, the putative receptor of Wuhan 2019-nCoV, in the nasal tissue. medRxiv [Preprint]. doi: 10.1101/2020.02.11.20022228v2

CrossRef Full Text | Google Scholar

Wyles, D. L., Rodriguez-Torres, M., Lawitz, E., Shiffman, M. L., Pol, S., Herring, R. W., et al. (2014). All-oral combination of ledipasvir, vedroprevir, tegobuvir, and ribavirin in treatment-naive patients with genotype 1 HCV infection. Hepatology 60, 56–64. doi: 10.1002/hep.27053

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, X., Chen, P., Wang, J., Feng, J., Zhou, H., Li, X., et al. (2020). Evolution of the novel coronavirus from the ongoing Wuhan outbreak and modeling of its spike protein for risk of human transmission. Sci. China Life Sci. 63, 457–460. doi: 10.1007/s11427-020-1637-5

PubMed Abstract | CrossRef Full Text | Google Scholar

Yamamoto, N., Yang, R., Yoshinaka, Y., Amari, S., Nakano, T., Cinatl, J., et al. (2004). HIV protease inhibitor nelfinavir inhibits replication of SARS-associated coronavirus. Biochem. Biophys. Res. Commun. 318, 719–725. doi: 10.1016/j.bbrc.2004.04.083

PubMed Abstract | CrossRef Full Text | Google Scholar

Zeuzem, S., Buggisch, P., Agarwal, K., Marcellin, P., Sereni, D., Klinker, H., et al. (2012). The protease inhibitor, GS-9256, and non-nucleoside polymerase inhibitor tegobuvir alone, with ribavirin, or pegylated interferon plus ribavirin in hepatitis C. Hepatology 55, 749–758. doi: 10.1002/hep.24744

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, Q., Zhang, H., Huang, K., Yang, Y., Hui, X., Gao, J., et al. (2020). SARS-CoV-2 neutralizing serum antibodies in cats: a serological investigation. bioRxiv [Preprint]. doi: 10.1101/2020.04.01.021196

CrossRef Full Text | Google Scholar

Zhu, N., Zhang, D., Wang, W., Li, X., Yang, B., Song, J., et al. (2020). A Novel Coronavirus from Patients with Pneumonia in China, 2019. New Engl. J. Med. 382, 727–733. doi: 10.1056/NEJMoa2001017

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: SARS-CoV-2, angiotensin-converting enzyme 2 (ACE2), nelfinavir (NFV), binding affinity [K(b)], cross-species transmission

Citation: Shen M, Liu C, Xu R, Ruan Z, Zhao S, Zhang H, Wang W, Huang X, Yang L, Tang Y, Yang T and Jia X (2020) Predicting the Animal Susceptibility and Therapeutic Drugs to SARS-CoV-2 Based on Spike Glycoprotein Combined With ACE2. Front. Genet. 11:575012. doi: 10.3389/fgene.2020.575012

Received: 22 June 2020; Accepted: 30 September 2020;
Published: 23 October 2020.

Edited by:

Xian-Tao Zeng, Wuhan University, China

Reviewed by:

Andras Szilagyi, Hungarian Academy of Sciences (MTA), Hungary
Fuhai Li, Washington University in St. Louis, United States

Copyright © 2020 Shen, Liu, Xu, Ruan, Zhao, Zhang, Wang, Huang, Yang, Tang, Yang and Jia. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Tai Yang, taiyang@cmc.edu.cn; Xu Jia, jiaxu@cmc.edu.cn

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.