Simulation of the omicron variant of SARS-CoV-2 shows broad antibody escape, weakened ACE2 binding, and modest increase in furin binding

2021 ◽  
Author(s):  
Muhammad- Zaki Jawaid ◽  
Avinash Baidya ◽  
Rustin Mahboubi-Ardakani ◽  
Richard L Davis ◽  
Daniel L Cox
Keyword(s):  
Molecular Dynamics ◽  
Viral Load ◽  
Spike Protein ◽  
Receptor Protein ◽  
Furin Cleavage ◽  
Binding Strength ◽  
Large Numbers ◽  
Recent Emergence ◽  
Alpha Variant ◽  
Cleavage Domain

The recent emergence of the omicron variant of the SARS-CoV-2 virus with large numbers of mutations has raised concern about a potential new surge in infections. Here we use molecular dynamics to study the biophysics of the interface of the omicron spike protein binding to (i) the ACE2 receptor protein, (ii) antibodies from all known binding regions, and (iii) the furin binding domain. Our simulations suggest that while there is significant reduction of antibody binding strength corresponding to escape, the omicron spike pays a cost in terms of weaker receptor binding. The furin cleavage domain is the same or weaker binding than the alpha variant, suggesting less viral load and disease intensity than the extant delta variant.

scholarly journals Elucidating Interactions Between SARS-CoV-2 Trimeric Spike Protein and ACE2 Using Homology Modeling and Molecular Dynamics Simulations

2021 ◽  
Vol 8 ◽  
Author(s):  
Sugunadevi Sakkiah ◽  
Wenjing Guo ◽  
Bohu Pan ◽  
Zuowei Ji ◽  
Gokhan Yavas ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Homology Modeling ◽  
Molecular Dynamics Simulations ◽  
Host Cell ◽  
Tertiary Structure ◽  
Cell Receptor ◽  
Host Cells ◽  
Spike Protein ◽  
Receptor Protein ◽  
Dynamics Simulations

Severe Acute Respiratory Syndrome Coronavirus-2 (SARS-CoV-2) causes coronavirus disease 2019 (COVID-19). As of October 21, 2020, more than 41.4 million confirmed cases and 1.1 million deaths have been reported. Thus, it is immensely important to develop drugs and vaccines to combat COVID-19. The spike protein present on the outer surface of the virion plays a major role in viral infection by binding to receptor proteins present on the outer membrane of host cells, triggering membrane fusion and internalization, which enables release of viral ssRNA into the host cell. Understanding the interactions between the SARS-CoV-2 trimeric spike protein and its host cell receptor protein, angiotensin converting enzyme 2 (ACE2), is important for developing drugs and vaccines to prevent and treat COVID-19. Several crystal structures of partial and mutant SARS-CoV-2 spike proteins have been reported; however, an atomistic structure of the wild-type SARS-CoV-2 trimeric spike protein complexed with ACE2 is not yet available. Therefore, in our study, homology modeling was used to build the trimeric form of the spike protein complexed with human ACE2, followed by all-atom molecular dynamics simulations to elucidate interactions at the interface between the spike protein and ACE2. Molecular Mechanics Poisson-Boltzmann Surface Area (MMPBSA) and in silico alanine scanning were employed to characterize the interacting residues at the interface. Twenty interacting residues in the spike protein were identified that are likely to be responsible for tightly binding to ACE2, of which five residues (Val445, Thr478, Gly485, Phe490, and Ser494) were not reported in the crystal structure of the truncated spike protein receptor binding domain (RBD) complexed with ACE2. These data indicate that the interactions between ACE2 and the tertiary structure of the full-length spike protein trimer are different from those between ACE2 and the truncated monomer of the spike protein RBD. These findings could facilitate the development of drugs and vaccines to prevent SARS-CoV-2 infection and combat COVID-19.

scholarly journals SARS-Cov-2 Spike binding to ACE2 is stronger and longer ranged with glycans

2021 ◽  
Author(s):  
Yihan Huang ◽  
Bradley S Harris ◽  
Shiaki A Minami ◽  
Seongwon Jung ◽  
Priya Shah ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Hydrogen Bonds ◽  
Electrostatic Interactions ◽  
Spike Protein ◽  
Binding Strength ◽  
Lennard Jones ◽  
Binding Domains ◽  
Interaction Mode ◽  
Dynamics Simulations ◽  
New Interactions

Highly detailed steered Molecular Dynamics simulations are performed on differently glycosylated receptor binding domains of the SARS-CoV-2 spike protein. The binding strength and the binding range increases with glycosylation. The interaction energy rises very quickly with pulling the proteins apart and only slowly drops at larger distances. We see a catch slip type behavior where interactions during pulling break and are taken over by new interactions forming. The dominant interaction mode are hydrogen bonds but Lennard-Jones and electrostatic interactions are relevant as well.

scholarly journals Computational study of the furin cleavage domain of SARS-CoV-2: delta binds strongest of extant variants

2022 ◽  
Author(s):  
Muhammad Zaki Jawaid ◽  
Avinash Baidya ◽  
Sofia Jakovcevic ◽  
Jacob Lusk ◽  
Rustin Mahboubi-Ardakani ◽  
...  
Keyword(s):  
Hydrogen Bond ◽  
Hydrogen Bonds ◽  
Computational Study ◽  
Binding Free Energy ◽  
Energy Estimates ◽  
Strong Positive Correlation ◽  
Statistical Accuracy ◽  
Furin Cleavage ◽  
Binding Strength ◽  
Cleavage Domain

We demonstrate that AlphaFold and AlphaFold Multimer, implemented within the ColabFold suite, can accurately predict the structures of the furin enzyme with known six residue inhibitory peptides. Noting the similarity of the peptide inhibitors to polybasic furin cleavage domain insertion region of the SARS-CoV-2, which begins at P681, we implement this approach to study the wild type furin cleavage domain for the virus and several mutants. We introduce mutations in silico for alpha, omicron, and delta variants, for several sequences which have been rarely observed, for sequences which have not yet been observed, for other coronaviruses (NL63, OC43, HUK1a, HUK1b, MERS, and 229E), and for the H5N1 flu. We show that interfacial hydrogen bonds between the furin cleavage domain and furin are a good measure of binding strength that correlate well with endpoint binding free energy estimates, and conclude that among all candidate viral sequences studied, delta is near the very top binding strength within statistical accuracy. However, the binding strength of several rare sequences match delta within statistical accuracy. We find that the furin S1 pocket is optimized for binding arginine as opposed to lysine. This residue, typically at sequence position five, contains the most hydrogen bonds to the furin, and hydrogen bond count for just this residue shows a strong positive correlation with the overall hydrogen bond count . We demonstrate that the root mean square backbone C-alpha fluctuation of the first residue in the furin cleavage domain has a strong negative correlation with the interfacial hydrogen bond count. We show by considering the variation with the number of basic residues that the maximum mean number of interfacial hydrogen bonds expected is 15.7 at 4 basic residues.

scholarly journals The Functional Consequences of the Novel Ribosomal Pausing Site in SARS-CoV-2 Spike Glycoprotein RNA

2021 ◽  
Vol 22 (12) ◽  
pp. 6490
Author(s):  
Olga A. Postnikova ◽  
Sheetal Uppal ◽  
Weiliang Huang ◽  
Maureen A. Kane ◽  
Rafael Villasmil ◽  
...  
Keyword(s):  
Amino Acid ◽  
Insertion Sequence ◽  
Experimental Studies ◽  
Spike Protein ◽  
Spike Glycoprotein ◽  
Furin Cleavage ◽  
New Host ◽  
S Protein ◽  
Furin Cleavage Site ◽  
Ribosome Pausing

The SARS-CoV-2 Spike glycoprotein (S protein) acquired a unique new 4 amino acid -PRRA- insertion sequence at amino acid residues (aa) 681–684 that forms a new furin cleavage site in S protein as well as several new glycosylation sites. We studied various statistical properties of the -PRRA- insertion at the RNA level (CCUCGGCGGGCA). The nucleotide composition and codon usage of this sequence are different from the rest of the SARS-CoV-2 genome. One of such features is two tandem CGG codons, although the CGG codon is the rarest codon in the SARS-CoV-2 genome. This suggests that the insertion sequence could cause ribosome pausing as the result of these rare codons. Due to population variants, the Nextstrain divergence measure of the CCU codon is extremely large. We cannot exclude that this divergence might affect host immune responses/effectiveness of SARS-CoV-2 vaccines, possibilities awaiting further investigation. Our experimental studies show that the expression level of original RNA sequence “wildtype” spike protein is much lower than for codon-optimized spike protein in all studied cell lines. Interestingly, the original spike sequence produces a higher titer of pseudoviral particles and a higher level of infection. Further mutagenesis experiments suggest that this dual-effect insert, comprised of a combination of overlapping translation pausing and furin sites, has allowed SARS-CoV-2 to infect its new host (human) more readily. This underlines the importance of ribosome pausing to allow efficient regulation of protein expression and also of cotranslational subdomain folding.

Viral and Host Attributes Underlying the Origins of Zoonotic Coronaviruses in Bats

2021 ◽  
Author(s):  
Alison E Stout ◽  
Qinghua Guo ◽  
Jean K Millet ◽  
Gary R Whittaker
Keyword(s):  
Viral Pathogenesis ◽  
Spike Protein ◽  
Future Research ◽  
Receptor Binding Domain ◽  
Viral Reservoirs ◽  
Furin Cleavage ◽  
Proteolytic Activation ◽  
Furin Cleavage Site ◽  
Unique Aspect ◽  
Genome Level

With a presumed origin in bats, the COVID-19 pandemic has been a major source of morbidity and mortality in the humanpopulation, and the causative agent, SARS-CoV-2, aligns most closely at the genome level with the bat coronavirusesRaBtCoV4991/RaTG13 and RmYN02. The ability of bats to provide reservoirs of numerous viruses in addition to coronavirusesremains an active area of research. Unique aspects of the physiology of the chiropteran immune system may contributeto the ability of bats to serve as viral reservoirs. The coronavirus spike protein plays important roles in viral pathogenesis and the immune response. Although much attention has focused on the spike receptor-binding domain, a unique aspect of SARS-CoV-2 as compared with its closest relatives is the presence of a furin cleavage site in the S1–S2 region of the spike protein. Proteolytic activation is likely an important feature that allows SARS-CoV-2—and other coronaviruses—to overcome the species barriers and thus cause human disease. The diversity of bat species limits the ability to draw broad conclusions about viral pathogenesis, but comparisons across species and with reference to humans and other susceptible mammals may guide future research in this regard.

Compared with SARS-CoV2 wild type's spike protein, the SARS-CoV2 omicron's receptor binding motif has adopted a more SARS-CoV1 and/or bat/civet-like structure.

2021 ◽  
Author(s):  
Michael O. Glocker ◽  
Kwabena F. M. Opuni ◽  
Hans-Juergen Thiesen
Keyword(s):  
Free Energy ◽  
Receptor Binding ◽  
Free Energy Calculations ◽  
Viral Fitness ◽  
Spike Protein ◽  
Receptor Protein ◽  
Binding Motif ◽  
Energy Calculations ◽  
Selection Processes ◽  
Protein Receptor

Our study focuses on free energy calculations of SARS-Cov2 spike protein receptor binding motives (RBMs) from wild type and variants-of-concern with particular emphasis on currently emerging SARS- CoV2 omicron variants of concern (VOC). Our computational free energy analysis underlines the occurrence of positive selection processes that specify omicron host adaption and bring changes on the molecular level into context with clinically relevant observations. Our free energy calculations studies regarding the interaction of omicron's RBM with human ACE2 shows weaker binding to ACE2 than alpha's, delta's, or wild type's RBM. Thus, less virus is predicted to be generated in time per infected cell. Our mutant analyses predict with focus on omicron variants a reduced spike-protein binding to ACE2--receptor protein possibly enhancing viral fitness / transmissibility and resulting in a delayed induction of danger signals as trade-off. Finally, more virus is produced but less per cell accompanied with delayed Covid-19 immunogenicity and pathogenicity. Regarding the latter, more virus is assumed to be required to initiate inflammatory immune responses.

ReaxFF-based molecular dynamics simulation of the interaction between plasma ROS and the receptor-binding domain of the spike protein in the capsid protein of SARS-CoV-2

2021 ◽  
Author(s):  
Huichao Wang ◽  
Tong Zhao ◽  
Shuhui Yang ◽  
Liang Zou ◽  
Xiaolong Wang ◽  
...  
Keyword(s):  
Molecular Dynamics ◽  
Amino Acid ◽  
Capsid Protein ◽  
Receptor Binding ◽  
Hydrogen Abstraction ◽  
Binding Domain ◽  
Spike Protein ◽  
Receptor Binding Domain ◽  
Amino Acid Residues ◽  
Global Epidemic

Abstract Under the severe situation of the current global epidemic, researchers have been working hard to find a reliable way to suppress the infection of the virus and prevent the spread of the epidemic. Studies have shown that the recognition and binding of human angiotensin-converting enzyme 2 (ACE2) by the receptor-binding domain (BRD) of spike protein on the surface of SARS-CoV-2 is a crucial step for SARS-CoV-2 to invade human receptor cells, and blocking this process can inhibit the virus from invading human normal cells. Plasma treatment can disrupt the structure of the RBD and effectively block the binding process. However, the mechanism by which plasma blocks the recognition and binding between the two is not clear. In this study, reaction process between reactive oxygen species (ROS) in plasma and the molecular model of RBD was simulated using a reactive molecular dynamics method. The results showed that the destruction of RBD molecule by ROS was triggered by hydrogen abstraction reactions. O and OH abstracted H atoms from RBD, while the H atoms of H2O2 and HO2 were abstracted by RBD. The hydrogen abstraction resulted in the breakage of C-H, N-H, O-H and C=O bonds and the formation of C=C, C=N bonds. The addition reaction of OH increased the number of O-H bonds and caused the formation of C-O, N-O and O-H bonds. The dissociation of N-H bonds led to the destruction of the original structure of peptide bonds and amino acid residues, change the type of amino acid residues, and caused the conversion of N-C and N=C, C=O and C-O. The simulation partially elucidated the microscopic mechanism of the interaction between ROS in plasma and the capsid protein of SARS-CoV-2, providing theoretical support for the control of SARS-CoV-2 infection by plasma, a contribution to overcoming the global epidemic problem.

Sign in / Sign up

Export Citation Format

Share Document