key: cord-0884863-4y1q8oqc authors: Soté, William O.; Franca, Eduardo F.; Hora, Aline S.; Comar, Moacyr title: A computational study of the interface interaction between SARS‐CoV‐2 RBD and ACE2 from human, cat, dog, and ferret date: 2021-07-22 journal: Transbound Emerg Dis DOI: 10.1111/tbed.14234 sha: d424041edc85a5d83da803bc1c00d3d5fc64ab6c doc_id: 884863 cord_uid: 4y1q8oqc The total impact of the worldwide COVID‐19 pandemic is still emerging, changing all relationships as a result, including those with pet animals. In the infection process, the use of angiotensin‐converting enzyme 2 (ACE2) as a cellular receptor to the spike protein of the new coronavirus is a fundamental step. In this sense, understanding which residue plays what role in the interaction between SARS‐CoV‐2 spike glycoprotein and ACE2 from cats, dogs, and ferrets is an important guide for helping to choose which animal model can be used to study the pathology of COVID‐19, and if there are differences between these interactions and those occurring in the human system. To help answer these questions, we performed classical molecular dynamics simulations to evaluate, from an atomistic point of view, the interactions in these systems. Our results show that there are significant differences in the interacting residues between the systems from different animal species, and the role of ACE2 key residues are different in each system, and can assist in the search for different inhibitors for each animal. The emergence of a new coronavirus (Morens et al., 2020; Wu et al., 2020; Zhou et al., 2020; Zhu et al., 2020) with high infection capacity clearly shows that there is an increasing need to study many aspects of this pathogen. A key point during infection is cellular attachment of the virus, mediated by interaction between the SARS-CoV-2 receptor binding domain (RBD), which plays a pivotal role in host selectivity (Li, 2015; Li et al., 2003) , and host angiotensin-converting enzyme 2 (ACE2) receptor. The origin of the new coronavirus has been demonstrated to be a species of bat (Cui et al., 2019; Li et al., 2006) , but recent reports suggest that animals which share some specific residues in ACE2 could be hosts of SARS-CoV-2, including cats, dogs, and ferrets (Luan et al., 2020) . Hence, considering the large population and economic importance of companion animals (Gapper, 2019) , and the difference in their ACE2 sequences (Figure 1 ), it is important to understand what the role of these ACE2 key residues is, in the interaction process between SARS-CoV-2 RBD and ACE2, because the participation of companion animals in the epidemiological chain of transmission remains in debate (Hora, 2020; Mallapaty, 2020; Shi et al., 2020; Sit et al., 2020) . In this sense, an evaluation of this process from an atomistic point of view can provide important information about the role of these specific residues; and the use of in silico strategies to simulate the interaction of proteins is ubiquitous, mainly using molecular docking and molecular dynamics simulations, since the advancement of both hardware and software has allowed for the study of bigger and larger systems (Perilla & Schulten, 2017 This is, as far we know, the first study of its kind. The structure used for the human system was SARS-CoV-2 RBD complexed with human receptor ACE2 (PDB: 6LZG). For cat (Uniprot: Q56H28), dog (Uniprot: J9P7Y2), and ferret (Uniprot: Q2WG88) systems, the FASTA sequences of ACE2 receptors were used in the Swiss Model webserver (Waterhouse et al., 2018) to build models which were equilibrated by a 1 ns simulation (details in section 2.2). The quality of all the models was confirmed with the help of Ramachandran plots (Fig- ure S1 in the Supporting Information), using PROCHECK online server (Laskowski et al., 1993) . The amino acid sequence of the SARS-CoV-2 spike obtained from human (Uniprot: P0DTC2) was compared with those obtained from dogs (GenBank: QIT08256.1; QIT08292.1) and cats (GenBank; QOF07648.1; QLG96797.1) naturally infected with SARS-CoV-2, to assess possible differences in RBD by visual comparison after alignment using the BioEdit 7.0.5.3 software (Hall, 1999) . No SARS-CoV-2 spike sequences from naturally infected ferrets were found in the public database. In RBD, differences between amino acids from human, dog, and cat were not found ( Figure the number of trials for rigid body minimization was five; the number of structures for rigid body docking (it0) was 1000, and the semi-flexible residues were automatically defined based on intermolecular contacts that were less than 5 Å. The fraction of common contacts (FCC) clustering parameters used were 0.6 for cut-off and 4 for minimum cluster size. The docking results are given in Table 1 . The Protein Data Bank (PDB) format of the docked systems, described in Section 2.1, were placed in a rectangular cuboid simulation box and were solvated by filling the box with Three Site Transferrable Intermolecular Potential (TIP3P) water model (Jorgensen et al., 1983) molecules. Sodium and chloride ions were used to neutralize the total charge of the macromolecules. Molecular dynamics simulations were carried out using GROMACS 2016.4 (Abraham et al., 2015) and OPLS-AA force field (Jorgensen & Tirado-Rives, 1988) . All systems were simulated in triplicate to ensure reliability in the simulation and the obtained results ( Figure S3 -S6 in the Supporting Information confirmed reproducibility). Each system was then energy minimized using 50,000 steps of the steepest descent method, and a convergence criterion with a maximum force of 1000 kJ mol −1 nm −1 , followed by two steps of solvent equilibration with all non-hydrogen protein atoms constrained by a force constant of 1000 kJ mol −1 nm −2 . The first equilibration consisted of a molecular dynamics simulation performed by 125 ps in an isochoric-isothermal (NVT, number of particles, volume, and temperature) ensemble at 310 K, maintained using velocity-rescale coupling method (Bussi et al., 2007) , with a coupling time constant of 0.1 ps; the second equilibration was performed in an isobaric-isothermal (NPT, number of particles, pressure, and temperature) ensemble at 1.0 bar, isotropically applied and maintained by the Berendsen weak coupling method (Berendsen et al., 1984) , with a coupling time constant of 0.5 ps and compressibility of 4.5 × 10 −5 bar −1 . Each production step of 50 ns used a time step of 1 fs with no position restraints; the temperature was maintained at 310 K using a Nosé-Hoover thermostat (Hoover, 1985; Nosé, 1984 ) with a coupling time constant of 0.2 ps. Pressure was maintained at 1.0 bar using a Parrinello-Rahman barostat (Nosé & Klein, 1983; Parrinello & Rahman, 1981 ) with a coupling time constant of 0.5 ps. All bonds' stretching and bending motions were constrained using LINCS algorithm (Hess et al., 1997) . Longrange electrostatics were evaluated using particle-mesh Ewald (PME) method, and a 1.0 nm cut-off was considered for short-range interactions. All analyses were performed with GROMACS and VMD software (Humphrey et al., 1996) . Using the described methodology, the following results were obtained to the structural and interaction behavior for the studied systems. All simulations were done in triplicate of 50 ns each, totalizing 600 ns of simulation, and the first analysis, shown in Figure (Li et al., 2005; Wan et al., 2020) , the average fluctuation for human system residues is 10 times lower than that of the other systems (Table 2) . This behavior, and the average fluctuation of the residues for the human system, suggest that human interface of interaction is more stable, which can also be confirmed by the formation and maintenance of hydrogen bonds along all simulation time ( Figure 3 ). Human ACE2 presents five key residues (K31, E35, D38, M82, and K353) while cats, dogs, and ferrets present three (K31, E35, and K353), and an important task is to know what interactions these residues can make, indicating possible sites of inhibition. In order to identify the interactions of these residues, we used the structures obtained from the clustering of each replica. The interaction of all residues for the systems are presented in Table S1 in the Supporting Information, and the principal interactions are in Figures S11-S14 in the Supporting Information. 3.2.1 Human system Considering the human system, in all replicas, the key residues of hACE2 had interactions with RBD residues, in agreement with the results of Wan et al. (2020) and Lan et al. (2020) . Furthermore, two key residues, E35 and M82, each interacted with only one RBD residue, Q493 and F486, respectively. All other residues had more than one interaction, but the RBD residues were basically the same in all replicas; the only difference was for K31, that had one interaction in Replica Considering this system, the interactions presented the same pattern that was noted in the other systems: residues from RBD and dACE2 formed more than one interaction. Two of the three key residues, one a hot spot, K31, were found having interactions in all replicas. The interaction pattern shows some residues from RBD that form more than one interaction, such as R466, K462, and L517, for example. From dACE2, E328 and K67 were found to have more than one interaction as well. In the cat system, among all interacting residues in the RBD, the residues E57, N330, and F72 had several interactions. However, as can be observed in Table S1 in the Supporting Information, the number of interaction residues change between replicas, suggesting that the interaction interface is moving and changing the interacting residues. Considering the key residues in cACE2, there are three of them, K31, E35, and K353. However, K31 only interacts in replica number 3, which is different from the dog system, where K31 (dACE2) interacts in all three simulations. The residue K458 from RBD interacts with E35 (cACE2) and K353 (cACE2) , showing differences in interaction pattern between this system and the dog's system. In the ferret system, two fACE2 key residues were found: K31 (a hot spot) and E35. The residue F486 from a ferret's RBD seems to be a key residue, because it forms several interactions with different residues in all replicas. The interactions of the two key residues, K31 (fACE2) :V445 (RBD) and K31 (fACE2) :G498 (RBD) , as well as E35 (fACE2) :G446 (RBD) , E35 (fACE2) :G447 (RBD) , and E35 (fACE2) :N448 (RBD) , show that the pattern of interaction in all three systems is different. Along with the analysis of the interacting residues, we extracted from simulations, 10 most persistent hydrogen bonding and among them we showed for how long the key residues, if they appear, maintained their interactions (Figure 4 ). For the human system, almost all hACE2 key residues are present, showing an occupancy of above 33%. Residues D38 (hACE2) and K353 (hACE2) interacted with the same RBD residue Q498, but only D38 (hACE2) interacted with Y499 (RBD) . The key residues E35 (hACE2) , D38 (hACE2) , and K353 (hACE2) presented the most stable interactions, Hydrogen bonds occupancy during the simulation time through the replicas. (Yellow for cat system, blue for dog system, green for ferret system, and orange for human system) forming a group that can contribute to stable interface interaction. The only residue that was not found doing any interaction was T82 (hACE2) . For the dog system, E35 (dACE2) was the most common interacting residue, as observed in the cat system. However, the occupancy was far below than what was observed in the cat interface, maintaining in all simulations below 50%. This residue interacts in a groove formed by S470 (RBD) and T471 (RBD) residues. T471 (RBD) interacted with another dACE2 key residue, E30 (dACE2) , that showed an occupancy of almost 40% in the first replica, but did not appear in the other two. Considering the cat system, hydrogen bond occupancy shows that only the E35 (cACE2) key residue is found in them, and presents a maximum occupancy of 71.73% as an acceptor of hydrogen bonding in one replica, and the donor was K458 (RBD) , that was the most common interacting residue to E35 (cACE2) , along with G476 (RBD) , although the occupancy of the latter residue was lower than 50% of the simulation time. Another interaction for the cat system involving E35 (cACE2) residue was with Y473 (RBD) . The K353 (cACE2) residue had interactions with N460 (RBD) and K458 (RBD) . K31 (cACE2) residue did not have any interactions. For the ferret system, the only fACE2 key residue to figure among the first 10 with higher occupancy, was E35 (fACE2) , showing a maximum value of 90.51% when interacting with G446 (RBD) . It also interacted with N448 (RBD) and G447 (RBD) during 79.20% and 67.86% of the simulated time, respectively. Another characteristic of the ferret system was the remarkable difference between replicas, suggesting that this interface is not well established. In this work, homology modeling, molecular docking, and molecular dynamics simulations, involving SARS-CoV-2 RBD and ACE2 from human (hACE2), cat (cACE2), dog (dACE2), and ferret (fACE2) systems, were used in order to evaluate and compare the interaction interface in each system. The interactions and hydrogen bonds in the human system were similar to that found in a recent report (Ali & Vijayan, 2020) , suggesting that the interaction interface does not fluctuate as much as those of the pets studied here. Analyses using ACE2 amino acids homology between animal species and humans have been published (Luan et al., 2020; Zhai et al., 2020) as a tool to predict which animal can be infected by SARS-CoV-2. However, these analyses only show resemblance between the residue sequences, but nothing can be inferred about the dynamic behavior of the systems regarding how they interact and should not be con-sidered to assess cellular permissiveness to infection, since from the atomistic simulation, the interaction between the RBD and the ACE2 is much more complex than a simple amino acids homology analysis, as has been shown. Furthermore, in these systems, key residues are considered crucial regions of interaction in the process of SARS-CoV and SARS-CoV-2 host cell infection (Li et al., 2005; Wan et al., 2020) , and divergences among the key residues in the human (K31, E35, D38, M82, and K353) and animal (K31, E35, and K353) systems were found, which leads to a clue that the interaction between the virus and the cell receptor in dog, cat, and ferret could behave differently. The COVID-19 pandemic caused by SARS-CoV-2 infections has created an urgent need for treatment and vaccine research, which require testing in appropriate animal models of the disease. An ideal animal model of a human disease is one in which the pattern of infection is similar to that which occurs in humans (Cleary et al., 2020) . Our results show that human interface of interaction is more stable, and there are differences in their interaction modes and binding pattern compared to pet's systems that suggest that the use of dog, cat, or ferret as experimental models will not reflect what happens during infection in humans, so results obtained using in vivo experiments with these species may not be applicable to addressing COVID-19. Additionally, the pathogenicity and high susceptibility to SARS-CoV and SARS-CoV-2 were reproduced only in transgenic mice that express human ACE2 (hACE2), and not in wild-type mice (mACE2), concluding that hACE2 was essential for SARS-CoV-2 infection and replication in mice (Bao et al., 2020; McCray et al., 2007; Netland et al., 2008; Sun et al., 2020) , which by analogy can be extended to the data found here for dog, cat, and ferret. Studies using SARS-CoV experimental infection in transgenic dog, cat, or ferret model expressing hACE2 were not found. Dogs and cats do not suffer from the presence of SARS-CoV-2, showing an absence of persistent infection and clinical signs (Hosie et al, 2020; Sit et al., 2020) , which can be suggested, from an atomistic point of view, by the differences in interactions observed in our results. Furthermore, SARS-CoV-2 led to acute bronchiolitis in experimentally infected ferrets; however, fatalities were not observed (Kim et al., 2020) , and reports of natural infections in this species were not found; high-affinity virus receptor interaction might be one of the crucial factors that determines the virulence of this pathogen in the host (Sarkar & Guha, 2020) . Ferrets are experimental models of respiratory viruses often used in laboratories, and are also kept as pets with some popularity in certain countries. Under natural conditions, 29 ferrets in prolonged and direct contact with two humans with COVID-19 proved that they are resistant to infection, with no detection of viral RNA (RT-qPCR), and no RBD-specific seroconversion (ELISA) (Sawatzki et al., 2021) . Beyond the question of these species as animal models, the possibility of detection of SARS-CoV-2 in companion animals has generated concerns about abandonment (Hora, 2020; Huang et al., 2020) . Cat and dog populations are extremely large in several countries, and the increase of these animals circulating in the streets could result in serious consequences to public health, in addition to questions related to animal welfare. Results obtained in this work suggest a very low-level participation of these animals in pandemic maintenance, which is sup-ported by reports stating that natural infections are rarely observed in dogs and cats (fewer than 200 cases), and absent in ferrets. Considering the world population of dogs and cats kept as pets, the number of infected animals already reported is irrelevant in the pandemic context. In this work, molecular dynamics simulations were used to study the interaction interface between SARS-CoV-2-RBD and ACE2 from humans, cats, dogs, and ferrets in order to identify commonalities across species. The analysis of the structural results suggest that the interaction interface of the human system is more stable, corroborated by differences in the fluctuation values of the key residues and the behavior of the hydrogen bonds along the simulation time. The differences in the number of interactions made by key residues during the replicas, between the systems, suggest that the infection process in humans is more effective. Additionally, the hydrogen bonds formed over time present high differences in the occupancy time during the simulation, showing that the human RBD-hACE2 interface can be considered more stable among the studied systems. Furthermore, the few reports of natural occurrence of SARS-CoV-2 in pet animals worldwide suggest that infection in them occurs rarely, and the results discussed here can help evaluate the reasons for these numbers. In addition, our data suggest that the use of dogs, cats, and ferrets as experimental models for treatment and vaccines for SARS-CoV-2 should be considered cautiously, if not actively discouraged. Gromacs: High performance molecular simulations through multi-level parallelism from laptops to supercomputers. SoftwareX, 1-2, 19-25 Dynamics of the ACE2-SARS-CoV-2/SARS-CoV spike protein interface reveal unique mechanisms Fighting COVID-19 using molecular dynamics simulations The pathogenicity of SARS-CoV-2 in hACE2 transgenic mice Molecular docking study of potential phytochemicals and their effects on the complex of SARS-CoV2 spike protein and human ACE2 Molecular dynamics with coupling to an external bath Canonical sampling through velocity rescaling Animal models of mechanisms of SARS-CoV-2 infection and COVID-19 pathology Origin and evolution of pathogenic coronaviruses Pets are making Wall Street wag its tail. Financial Times BioEdit: A user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT LINCS: A linear constraint solver for molecular simulations Canonical dynamics: Equilibrium phase-space distributions Coronavirus: A veterinary perspective SARS-Coronavirus (CoV)-2 and cats. European Advisory Board on Cat Diseases COVID-19 pandemic: Stop panic abandonment of household pets VMD: Visual molecular dynamics Comparison of simple potential functions for simulating liquid water The OPLS potential functions for proteins. Energy minimizations for crystals of cyclic peptides and crambin Infection and rapid transmission of SARS-CoV-2 in ferrets Structure of the SARS-CoV-2 spike receptor-binding domain bound to the ACE2 receptor PROCHECK: A program to check the stereochemical quality of protein structures Receptor recognition mechanisms of coronaviruses: A decade of structural studies Structural biology: Structure of SARS coronavirus spike receptor-binding domain complexed with receptor Angiotensin-converting enzyme 2 is a functional receptor for the SARS coronavirus Animal origins of the severe acute respiratory syndrome coronavirus: Insight from ACE2-S-protein interactions Spike protein recognition of mammalian ACE2 predicts the host range and an optimized ACE2 for SARS-CoV-2 infection Coronavirus can infect cats -dogs, not so much Lethal infection of K18-hACE2 mice infected with severe acute respiratory syndrome coronavirus Escaping Pandora's box -Another novel coronavirus Severe acute respiratory syndrome coronavirus infection causes neuronal death in the absence of encephalitis in mice transgenic for human ACE2 Constant pressure molecular dynamics for molecular systems A unified formulation of the constant temperature molecular dynamics methods Polymorphic transitions in single crystals: A new molecular dynamics method Physical properties of the HIV-1 capsid from all-atom molecular dynamics simulations Identification of a potential SARS-CoV2 inhibitor via molecular dynamics simulations and amino acid decomposition analysis Infectivity, virulence, pathogenicity, hostpathogen interactions of SARS and SARS-CoV-2 in experimental animals: A systematic review Host barriers to SARS-CoV-2 demonstrated by ferrets in a high-exposure domestic setting Susceptibility of ferrets, cats, dogs, and other domesticated animals to SARS-coronavirus 2 Infection of dogs with SARS-CoV-2 A mouse model of SARS-CoV-2 infection and pathogenesis The HADDOCK2.2 web server: User-friendly integrative modeling of biomolecular complexes Receptor recognition by the novel coronavirus from Wuhan: An analysis based on decade-long structural studies of SARS coronavirus SWISS-MODEL: Homology modelling of protein structures and complexes A new coronavirus associated with human respiratory disease in China Comparison of severe acute respiratory syndrome coronavirus 2 spike protein binding to ACE2 receptors from human, pets, farm animals, and putative intermediate hosts A pneumonia outbreak associated with a new coronavirus of probable bat origin & China Novel Coronavirus Investigating and Research Team A computational study of the interface interaction between SARS-CoV-2 RBD and ACE2 from human, cat, dog, and ferret The authors want to thank Cesar Simenes who allowed the simulations to be done in the Orion Computer at UFU-Ituiutaba. The authors confirm that there is no conflict of interest. The authors declare that ethics statements are not applicable for this work. There were no questionnaires or sample collections. All used information was retrieved from published articles. Data sharing not applicable-the article describes entirely theoretical research.