Dataset Information

Molecular dynamic simulation analysis of SARS-CoV-2 spike mutations and evaluation of ACE2 from pets and wild animals for infection risk.

ABSTRACT: Coronavirus Disease 2019 (COVID-19) is an ongoing global health emergency that has caused tremendous stress and loss of life worldwide. The viral spike glycoprotein is a critical molecule mediating transmission of SARS-CoV-2 by interacting with human ACE2. However, through the course of the pandemics, there has not been a thorough analysis of the spike protein mutations, and on how these mutants influence the transmission of SARS-CoV-2. Besides, cases of SARS-CoV-2 infection among pets and wild animals have been reported, so the susceptibility of these animals requires great attention to investigate, as they may also link to the renewed question of a possible intermediate host for SARS-CoV-2 before it was transmitted to humans. With over 226,000 SARS-CoV-2 sequences obtained, we found 1573 missense mutations in the spike gene, and 226 of them were within the receptor-binding domain (RBD) region that directly interacts with human ACE2. Modeling the interactions between SARS-CoV-2 spike mutants and ACE2 molecules showed that most of the 74 missense mutations in the RBD region of the interaction interface had little impact on spike binding to ACE2, whereas several within the spike RBD increased the binding affinity toward human ACE2 thus making the virus likely more contagious. On the other hand, modeling the interactions between animal ACE2 molecules and SARS-CoV-2 spike revealed that many pets and wild animals' ACE2 had a variable binding ability. Particularly, ACE2 of bamboo rat had stronger binding to SARS-CoV-2 spike protein, whereas that of mole, vole, Mus pahari, palm civet, and pangolin had a weaker binding compared to human ACE2. Our results provide structural insights into the impact on interactions of the SARS-CoV-2 spike mutants to human ACE2, and shed light on SARS-CoV-2 transmission in pets and wild animals, and possible clues to the intermediate host(s) for SARS-CoV-2.

SUBMITTER: Chen P

PROVIDER: S-EPMC8634692 | biostudies-literature |

REPOSITORIES: biostudies-literature

ACCESS DATA

Similar Datasets

Project description:The emergence of novel coronavirus mutants is a main factor behind the deterioration of the epidemic situation. Further studies into the pathogenicity of these mutants are thus urgently needed. Binding of the spinous protein receptor binding domain (RBD) of SARS-CoV-2 to the angiotensin-converting enzyme 2 (ACE2) receptor was shown to initiate coronavirus entry into host cells and lead to their infection. The receptor-binding motif (RBM, 438-506) is a region that directly interacts with ACE2 receptor in the RBD and plays a crucial role in determining affinity. To unravel how mutations in the non-RBM regions impact the interaction between RBD and ACE2, we selected three non-RBM mutant systems (N354D, D364Y, and V367F) from the documented clinical cases, and the Q498A mutant system located in the RBM region served as the control. Molecular dynamics simulation was conducted on the mutant systems and the wild-type (WT) system, and verified experiments also performed. Non-RBM mutations have been shown not only to change conformation of the RBM region but also to significantly influence its hydrogen bonding and hydrophobic interactions. In particular, the D364Y and V367F systems showed a higher affinity for ACE2 owing to their electrostatic interactions and polar solvation energy changes. In addition, although the binding free energy at this point increased after the mutation of N354D, the conformation of the random coil (Pro384-Asp389) was looser than that of other systems, and the combined effect weakened the binding free energy between RBD and ACE2. Interestingly, we also found a random coil (Ala475-Gly485). This random coil is very sensitive to mutations, and both types of mutations increase the binding free energy of residues in this region. We found that the binding loop (Tyr495-Tyr505) in the RBD domain strongly binds to Lys353, an important residue of the ACE2 domain previously identified. The binding free energy of the non-RBM mutant group at the binding loop had positive and negative changes, and these changes were more obvious than that of the Q498A system. The results of this study elucidate the effect of non-RBM mutation on ACE2-RBD binding, and provide new insights for SARS-CoV-2 mutation research.

Project description:The SARS-CoV-2 claimed millions of lives, globally. Occurring from Wuhan (wild type) in December, 2019, it constantly mutated to Omicron (B.1.1.529), the predecessor to Delta. Omicron having ~ 32 spike mutations has variable infectivity-multiplicity-immuno-invasive properties. Understanding of its mutational effect on ACE2-binding/disease severity and developing preventive/therapeutic strategies are important. The binding affinities of Wuhan/Delta/Omicron spikes (PDB/GISAID/SWISS-MODEL) were docked (HADDOCK2.4) with ACE2 and compared by competitive-docking (PRODIGY). The protein structural stability was verified by kinetic-data/Ramachandran-plot (Zlab/UMassMedBioinfo). After several trials, a 59 amino acid (453ARG-510VAL) peptide-cut (Expasy-server) of the wild-type spike RBD with some desired mutants (THR500SER/THR500GLY/THR500ALA/THR500CYS) was blindly/competitively docked (PyMOL-V2.2.2) to block the Omicron-ACE2 binding. We examined molecular dynamic simulation (iMOD-server, with 9000 cycles/300 k-heating/1 atm pressure for system equilibration for 50 ns-run) of ACE2 and two CUTs with different SARS-CoV-2 variants. The binding-affinity of Omicron-ACE2 is slightly higher than the rest two in competitive docking setup. During individual (1:1) docking, Omicron showed little higher than wild type but much weaker binding affinity than Delta. Competitive docking suggests ten H-bonding (1.3-2.4 Å) with highly favorable energy values/Van-der-Walls-force/Haddock score for more stable-binding of Omicron-RBD with ACE2. Blind docking of different CUTs (wild/mutants) and Omicron to ACE2 completely rejected the Omicron-RBD from ACE2-target. The best blocking/binding affinity of -16.4 and -13 kcal/mole were observed in the case of THR500SER and THR500GLY, respectively, with multiple H-bonding 1.9-2.2 Å. These are supported by the MD-simulation results. So, the spike binding affinities were Delta > Omicron > wild in 1:1 docking with ACE2. Considering the wild type is non-existing nowadays, Omicron showed less ACE2 binding properties. The 59 cut of spike-RBD and its mutant THR500SER/THR500GLY may be further screened as universal blockers of this virus.Supplementary informationThe online version contains supplementary material available at 10.1007/s11224-022-02022-x.

Dataset Information

Molecular dynamic simulation analysis of SARS-CoV-2 spike mutations and evaluation of ACE2 from pets and wild animals for infection risk.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets