Search

The COVID-19 pandemic continues to pose a substantial threat to human lives and is likely to do so for years to come. Despite the availability of vaccines, searching for efficient small-molecule drugs that are widely available, including in low- and middle-income countries, is an ongoing challenge. In this work, we report the results of an open science community effort, the "Billion molecules against COVID-19 challenge", to identify small-molecule inhibitors against SARS-CoV-2 or relevant human receptors. Participating teams used a wide variety of computational methods to screen a minimum of 1 billion virtual molecules against 6 protein targets. Overall, 31 teams participated, and they suggested a total of 639,024 molecules, which were subsequently ranked to find 'consensus compounds'. The organizing team coordinated with various contract research organizations (CROs) and collaborating institutions to synthesize and test 878 compounds for biological activity against proteases (Nsp5, Nsp3, TMPRSS2), nucleocapsid N, RdRP (only the Nsp12 domain), and (alpha) spike protein S. Overall, 27 compounds with weak inhibition/binding were experimentally identified by binding-, cleavage-, and/or viral suppression assays and are presented here. Open science approaches such as the one presented here contribute to the knowledge base of future drug discovery efforts in finding better SARS-CoV-2 treatments., (© 2023 The Authors. Molecular Informatics published by Wiley-VCH GmbH.)

Multiple severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants continue to evolve carrying flexible amino acid substitutions in the spike protein's receptor binding domain (RBD). These substitutions modify the binding of the SARS-CoV-2 to human angiotensin-converting enzyme 2 (hACE2) receptor and have been implicated in altered host fitness, transmissibility, and efficacy against antibody therapeutics and vaccines. Reliably predicting the binding strength of SARS-CoV-2 variants RBD to hACE2 receptor and neutralizing antibodies (NAbs) can help assessing their fitness, and rapid deployment of effective antibody therapeutics, respectively. Here, we introduced a two-step computational framework with 3-fold validation that first identified dissociation constant as a reliable predictor of binding affinity in hetero- dimeric and trimeric protein complexes. The second step implements dissociation constant as descriptor of the binding strengths of SARS-CoV-2 variants RBD to hACE2 and NAbs. Then, we examined several variants of concerns (VOCs) such as Alpha, Beta, Gamma, Delta, and Omicron and demonstrated that these VOCs RBD bind to the hACE2 with enhanced affinity. Furthermore, the binding affinity of Omicron variant's RBD was reduced with majority of the RBD-directed NAbs, which is highly consistent with the experimental neutralization data. By studying the atomic contacts between RBD and NAbs, we revealed the molecular footprints of four NAbs (GH-12, P2B-1A1, Asarnow_3D11, and C118)-that may likely neutralize the recently emerged Omicron variant-facilitating enhanced binding affinity. Finally, our findings suggest a computational pathway that could aid researchers identify a range of current NAbs that may be effective against emerging SARS-CoV-2 variants., (© 2023 Wiley Periodicals LLC.)

The ongoing evolution of severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) has resulted in the recent emergence of a highly divergent variant of concern (VOC) defined as Omicron or B.1.1.529. This VOC is of particular concern because it has the potential to evade most therapeutic antibodies and has undergone a sustained genetic evolution, resulting in the emergence of five distinct sub-lineages. However, the evolutionary dynamics of the initially identified Omicron BA.1 and BA.2 sub-lineages remain poorly understood. Herein, we combined Bayesian phylogenetic analysis, mutational profiling, and selection pressure analysis to track the virus's genetic changes that drive the early evolutionary dynamics of the Omicron. Based on the Omicron dataset chosen for the improved temporal signals and sampled globally between November 2021 and January 2022, the most recent common ancestor (tMRCA) and substitution rates for BA.1 were estimated to be that of 18 September 2021 (95% highest posterior density (HPD), 4 August-22 October 2021) and 1.435 × 10 -3 (95% HPD  =  1.021 × 10 -3 - 1.869 × 10 -3 ) substitution/site/year, respectively, whereas 3 November 2021 (95% highest posterior density (HPD) 26 September-28 November 2021) and 1.074 × 10 -3 (95% HPD  =  6.444 × 10 -4 - 1.586 × 10 -3 ) substitution/site/year were estimated for the BA.2 sub-lineage. The findings of this study suggest that the Omicron BA.1 and BA.2 sub-lineages originated independently and evolved over time. Furthermore, we identified multiple sites in the spike protein undergoing continued diversifying selection that may alter the neutralization profile of BA.1. This study sheds light on the ongoing global genomic surveillance and Bayesian molecular dating analyses to better understand the evolutionary dynamics of the virus and, as a result, mitigate the impact of emerging variants on public health.

Understanding the origin of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been a highly debatable and unresolved issue for scientific communities all over the world. Understanding the mechanism of virus entry to the host cells is crucial to deciphering the susceptibility profiles of animal species to SARS-CoV-2. The interaction of SARS-CoV-2 ligands (receptor-binding domain on spike protein) with its host cell receptor, angiotensin-converting enzyme 2 (ACE2), is a critical determinant of host range and cross-species transmission. In this study, we developed and implemented a rigorous computational approach for predicting binding affinity between 299 ACE2 orthologs from diverse vertebrate species and the SARS-CoV-2 spike protein. The findings show that the SARS-CoV-2 spike protein can bind to a wide range of vertebrate species carrying evolutionary divergent ACE2, implying a broad host range at the virus entry level, which may contribute to cross-species transmission and further viral evolution. Furthermore, the current study facilitated the identification of genetic determinants that may differentiate susceptible from resistant host species based on the conservation of ACE2-spike protein interacting residues in vertebrate host species known to facilitate SARS-CoV-2 infection; however, these genetic determinants warrant in vivo experimental confirmation. The molecular interactions associated with varied binding affinity of distinct ACE2 isoforms in a specific bat species were identified using protein structure analysis, implying the existence of diversified bat species' susceptibility to SARS-CoV-2. The current study's findings highlight the importance of intensive surveillance programmes aimed at identifying susceptible hosts, especially those with the potential to transmit zoonotic pathogens, in order to prevent future outbreaks., (Copyright © 2022 Elsevier Inc. All rights reserved.)