SARS-CoV-2 omicron spike simulations: broad antibody escape, weakened ACE2 binding, and modest furin cleavage

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 BA1 and BA2 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 a significant reduction of antibody (Ab) binding strength corresponding to escape, the omicron spikes pay a cost in terms of weaker receptor binding as measured by interfacial hydrogen bonds (H-bond). The furin cleavage domain (FCD) is the same or weaker binding than the delta variant, suggesting lower fusogenicity resulting in less viral load and disease intensity than the delta variant. IMPORTANCE The BA1 and BA2 and closely related BA2.12.2 and BA.5 omicron variants of SARS-CoV-2 dominate the current global infection landscape. Given the high number of mutations, particularly those which will lead to antibody escape, it is important to establish accurate methods that can guide developing health policy responses that identify at a fundamental level whether omicron and its variants are more threatening than its predecessors, especially delta. The importance of our work is to demonstrate that simple in silico simulations can predict biochemical binding details of the omicron spike protein that have epidemiological consequences, especially for binding to the cells and for fusing the viral membrane with the cells. In each case, we predicted weaker binding of the omicron spike, which agreed with subsequent experimental results. Future virology experiments will be needed to test these predictions further.

Using machine learning to detect coronaviruses potentially infectious to humans

Establishing the host range for novel viruses remains a challenge. Here, we address the challenge of identifying non-human animal coronaviruses that may infect humans by creating an artificial neural network model that learns from spike protein sequences of alpha and beta coronaviruses and their binding annotation to their host receptor. The proposed method produces a human-Binding Potential (h-BiP) score that distinguishes, with high accuracy, the binding potential among coronaviruses. Three viruses, previously unknown to bind human receptors, were identified: Bat coronavirus BtCoV/133/2005 and Pipistrellus abramus bat coronavirus HKU5-related (both MERS related viruses), and Rhinolophus affinis coronavirus isolate LYRa3 (a SARS related virus). We further analyze the binding properties of BtCoV/133/2005 and LYRa3 using molecular dynamics. To test whether this model can be used for surveillance of novel coronaviruses, we re-trained the model on a set that excludes SARS-CoV-2 and all viral sequences released after the SARS-CoV-2 was published. The results predict the binding of SARS-CoV-2 with a human receptor, indicating that machine learning methods are an excellent tool for the prediction of host expansion events.