A multi-pronged approach targeting SARS-CoV-2 proteins using ultra-large virtual screening

SARS-CoV-2-induced phosphorylation and its pharmacotherapy backed by artificial intelligence and machine learning

Aims: To investigate the role of phosphorylation in SARS-CoV-2 infection, potential therapeutic targets and its harmful genetic sequences. Materials & Methods: Data mining techniques were employed to identify upregulated kinases responsible for proteomic changes induced by SARS-CoV-2. Spike and nucleocapsid proteins' sequences were analyzed using predictive tools, including SNAP2, MutPred2, PhD-SNP, SNPs&Go, MetaSNP, Predict-SNP and PolyPhen-2. Missense variants were identified using ensemble-based algorithms and homology/structure-based models like SIFT, PROVEAN, Predict-SNP and MutPred-2. Results: Eight missense variants were identified in viral sequences. Four damaging variants were found, with SNPs&Go and PolyPhen-2. Promising therapeutic candidates, including gilteritinib, pictilisib, sorafenib, RO5126766 and omipalisib, were identified. Conclusion: This research offers insights into SARS-CoV-2 pathogenicity, highlighting potential treatments and harmful variants in viral proteins.

Virtual Screening of a Chemically Diverse "Superscaffold" Library Enables Ligand Discovery for a Key GPCR Target

The advent of ultra-large libraries of drug-like compounds has significantly broadened the possibilities in structure-based virtual screening, accelerating the discovery and optimization of high-quality lead chemotypes for diverse clinical targets. Compared to traditional high-throughput screening, which is constrained to libraries of approximately one million compounds, the ultra-large virtual screening approach offers substantial advantages in both cost and time efficiency. By expanding the chemical space with compounds synthesized from easily accessible and reproducible reactions and utilizing a large, diverse set of building blocks, we can enhance both the diversity and quality of the discovered lead chemotypes. In this study, we explore new chemical spaces using reactions of sulfur(VI) fluorides to create a combinatorial library consisting of several hundred million compounds. We screened this virtual library for cannabinoid type II receptor (CB2) antagonists using the high-resolution structure in conjunction with a rationally designed antagonist, AM10257. The top-predicted compounds were then synthesized and tested in vitro for CB2 binding and functional antagonism, achieving an experimentally validated hit rate of 55%. Our findings demonstrate the effectiveness of reliable reactions, such as sulfur fluoride exchange, in diversifying ultra-large chemical spaces and facilitate the discovery of new lead compounds for important biological targets.

AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays or AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold-Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

Large-Scale Pretraining Improves Sample Efficiency of Active Learning-Based Virtual Screening

Virtual screening of large compound libraries to identify potential hit candidates is one of the earliest steps in drug discovery. As the size of commercially available compound collections grows exponentially to the scale of billions, active learning and Bayesian optimization have recently been proven as effective methods of narrowing down the search space. An essential component of those methods is a surrogate machine learning model that predicts the desired properties of compounds. An accurate model can achieve high sample efficiency by finding hits with only a fraction of the entire library being virtually screened. In this study, we examined the performance of a pretrained transformer-based language model and graph neural network in a Bayesian optimization active learning framework. The best pretrained model identifies 58.97% of the top-50,000 compounds after screening only 0.6% of an ultralarge library containing 99.5 million compounds, improving 8% over the previous state-of-the-art baseline. Through extensive benchmarks, we show that the superior performance of pretrained models persists in both structure-based and ligand-based drug discovery. Pretrained models can serve as a boost to the accuracy and sample efficiency of active learning-based virtual screening.

Shape-Aware Synthon Search (SASS) for Virtual Screening of Synthon-Based Chemical Spaces

Virtual screening of large-scale chemical libraries has become increasingly useful for identifying high-quality candidates for drug discovery. While it is possible to exhaustively screen chemical spaces that number on the order of billions, indirect combinatorial approaches are needed to efficiently navigate larger, synthon-based virtual spaces. We describe Shape-Aware Synthon Search (SASS), a synthon-based virtual screening method that carries out shape similarity searches in the synthon space instead of the enumerated product space. SASS can replicate results from exhaustive searches in ultralarge, combinatorial spaces with high recall on a variety of query molecules while only scoring a small subspace of possible enumerated products, thereby significantly accelerating large-scale, shape-based virtual screening.

DockOpt: A Tool for Automatic Optimization of Docking Models

Molecular docking is a widely used technique for leveraging protein structure for ligand discovery, but it remains difficult to utilize due to limitations that have not been adequately addressed. Despite some progress toward automation, docking still requires expert guidance, hindering its adoption by a broader range of investigators. To make docking more accessible, we developed a new utility called DockOpt, which automates the creation, evaluation, and optimization of docking models prior to their deployment in large-scale prospective screens. DockOpt outperforms our previous automated pipeline across all 43 targets in the DUDE-Z benchmark data set, and the generated models for 84% of targets demonstrate sufficient enrichment to warrant their use in prospective screens, with normalized LogAUC values of at least 15%. DockOpt is available as part of the Python package Pydock3 included in the UCSF DOCK 3.8 distribution, which is available for free to academic researchers at https://dock.compbio.ucsf.edu and free for everyone upon registration at https://tldr.docking.org.

Docking and other computing tools in drug design against SARS-CoV-2

The use of computer simulation methods has become an indispensable component in identifying drugs against the SARS-CoV-2 coronavirus. There is a huge body of literature on application of molecular modelling to predict inhibitors against target proteins of SARS-CoV-2. To keep our review clear and readable, we limited ourselves primarily to works that use computational methods to find inhibitors and test the predicted compounds experimentally either in target protein assays or in cell culture with live SARS-CoV-2. Some works containing results of experimental discovery of corresponding inhibitors without using computer modelling are included as examples of a success. Also, some computational works without experimental confirmations are also included if they attract our attention either by simulation methods or by databases used. This review collects studies that use various molecular modelling methods: docking, molecular dynamics, quantum mechanics, machine learning, and others. Most of these studies are based on docking, and other methods are used mainly for post-processing to select the best compounds among those found through docking. Simulation methods are presented concisely, information is also provided on databases of organic compounds that can be useful for virtual screening, and the review itself is structured in accordance with coronavirus target proteins.

Medicinal chemistry strategies towards the development of non-covalent SARS-CoV-2 Mpro inhibitors

The main protease (Mpro) of SARS-CoV-2 is an attractive target in anti-COVID-19 therapy for its high conservation and major role in the virus life cycle. The covalent Mpro inhibitor nirmatrelvir (in combination with ritonavir, a pharmacokinetic enhancer) and the non-covalent inhibitor ensitrelvir have shown efficacy in clinical trials and have been approved for therapeutic use. Effective antiviral drugs are needed to fight the pandemic, while non-covalent Mpro inhibitors could be promising alternatives due to their high selectivity and favorable druggability. Numerous non-covalent Mpro inhibitors with desirable properties have been developed based on available crystal structures of Mpro. In this article, we describe medicinal chemistry strategies applied for the discovery and optimization of non-covalent Mpro inhibitors, followed by a general overview and critical analysis of the available information. Prospective viewpoints and insights into current strategies for the development of non-covalent Mpro inhibitors are also discussed.

AI-assisted structural consensus-proteome prediction of human monkeypox viruses isolated within a year after the 2022 multi-country outbreak

IMPORTANCE

The 2022 outbreak of the monkeypox virus already involves, by April 2023, 110 countries with 86,956 confirmed cases and 119 deaths. Understanding an emerging disease on a molecular level is essential to study infection processes and eventually guide drug discovery at an early stage. To support this, we provide the so far most comprehensive structural proteome of the monkeypox virus, which includes 210 structural models, each computed with three state-of-the-art structure prediction methods. Instead of building on a single-genome sequence, we generated our models from a consensus of 3,713 high-quality genome sequences sampled from patients within 1 year of the outbreak. Therefore, we present an average structural proteome of the currently isolated viruses, including mutational analyses with a special focus on drug-binding sites. Continuing dynamic mutation monitoring within the structural proteome presented here is essential to timely predict possible physiological changes in the evolving virus.

Unsupervised deep learning for molecular dynamics simulations: a novel analysis of protein-ligand interactions in SARS-CoV-2 Mpro

Molecular dynamics (MD) simulations, which are central to drug discovery, offer detailed insights into protein-ligand interactions. However, analyzing large MD datasets remains a challenge. Current machine-learning solutions are predominantly supervised and have data labelling and standardisation issues. In this study, we adopted an unsupervised deep-learning framework, previously benchmarked for rigid proteins, to study the more flexible SARS-CoV-2 main protease (Mpro). We ran MD simulations of Mpro with various ligands and refined the data by focusing on binding-site residues and time frames in stable protein conformations. The optimal descriptor chosen was the distance between the residues and the center of the binding pocket. Using this approach, a local dynamic ensemble was generated and fed into our neural network to compute Wasserstein distances across system pairs, revealing ligand-induced conformational differences in Mpro. Dimensionality reduction yielded an embedding map that correlated ligand-induced dynamics and binding affinity. Notably, the high-affinity compounds showed pronounced effects on the protein's conformations. We also identified the key residues that contributed to these differences. Our findings emphasize the potential of combining unsupervised deep learning with MD simulations to extract valuable information and accelerate drug discovery.

A practical guide to machine-learning scoring for structure-based virtual screening

Structure-based virtual screening (SBVS) via docking has been used to discover active molecules for a range of therapeutic targets. Chemical and protein data sets that contain integrated bioactivity information have increased both in number and in size. Artificial intelligence and, more concretely, its machine-learning (ML) branch, including deep learning, have effectively exploited these data sets to build scoring functions (SFs) for SBVS against targets with an atomic-resolution 3D model (e.g., generated by X-ray crystallography or predicted by AlphaFold2). Often outperforming their generic and non-ML counterparts, target-specific ML-based SFs represent the state of the art for SBVS. Here, we present a comprehensive and user-friendly protocol to build and rigorously evaluate these new SFs for SBVS. This protocol is organized into four sections: (i) using a public benchmark of a given target to evaluate an existing generic SF; (ii) preparing experimental data for a target from public repositories; (iii) partitioning data into a training set and a test set for subsequent target-specific ML modeling; and (iv) generating and evaluating target-specific ML SFs by using the prepared training-test partitions. All necessary code and input/output data related to three example targets (acetylcholinesterase, HMG-CoA reductase, and peroxisome proliferator-activated receptor-α) are available at https://github.com/vktrannguyen/MLSF-protocol , can be run by using a single computer within 1 week and make use of easily accessible software/programs (e.g., Smina, CNN-Score, RF-Score-VS and DeepCoy) and web resources. Our aim is to provide practical guidance on how to augment training data to enhance SBVS performance, how to identify the most suitable supervised learning algorithm for a data set, and how to build an SF with the highest likelihood of discovering target-active molecules within a given compound library.

Identification of motif-based interactions between SARS-CoV-2 protein domains and human peptide ligands pinpoint antiviral targets

The virus life cycle depends on host-virus protein-protein interactions, which often involve a disordered protein region binding to a folded protein domain. Here, we used proteomic peptide phage display (ProP-PD) to identify peptides from the intrinsically disordered regions of the human proteome that bind to folded protein domains encoded by the SARS-CoV-2 genome. Eleven folded domains of SARS-CoV-2 proteins were found to bind 281 peptides from human proteins, and affinities of 31 interactions involving eight SARS-CoV-2 protein domains were determined (KD ∼ 7-300 μM). Key specificity residues of the peptides were established for six of the interactions. Two of the peptides, binding Nsp9 and Nsp16, respectively, inhibited viral replication. Our findings demonstrate how high-throughput peptide binding screens simultaneously identify potential host-virus interactions and peptides with antiviral properties. Furthermore, the high number of low-affinity interactions suggest that overexpression of viral proteins during infection may perturb multiple cellular pathways.

AI-guided pipeline for protein-protein interaction drug discovery identifies a SARS-CoV-2 inhibitor

Protein-protein interactions (PPIs) offer great opportunities to expand the druggable proteome and therapeutically tackle various diseases, but remain challenging targets for drug discovery. Here, we provide a comprehensive pipeline that combines experimental and computational tools to identify and validate PPI targets and perform early-stage drug discovery. We have developed a machine learning approach that prioritizes interactions by analyzing quantitative data from binary PPI assays and AlphaFold-Multimer predictions. Using the quantitative assay LuTHy together with our machine learning algorithm, we identified high-confidence interactions among SARS-CoV-2 proteins for which we predicted three-dimensional structures using AlphaFold Multimer. We employed VirtualFlow to target the contact interface of the NSP10-NSP16 SARS-CoV-2 methyltransferase complex by ultra-large virtual drug screening. Thereby, we identified a compound that binds to NSP10 and inhibits its interaction with NSP16, while also disrupting the methyltransferase activity of the complex, and SARS-CoV-2 replication. Overall, this pipeline will help to prioritize PPI targets to accelerate the discovery of early-stage drug candidates targeting protein complexes and pathways.

Computational approaches streamlining drug discovery

Computer-aided drug discovery has been around for decades, although the past few years have seen a tectonic shift towards embracing computational technologies in both academia and pharma. This shift is largely defined by the flood of data on ligand properties and binding to therapeutic targets and their 3D structures, abundant computing capacities and the advent of on-demand virtual libraries of drug-like small molecules in their billions. Taking full advantage of these resources requires fast computational methods for effective ligand screening. This includes structure-based virtual screening of gigascale chemical spaces, further facilitated by fast iterative screening approaches. Highly synergistic are developments in deep learning predictions of ligand properties and target activities in lieu of receptor structure. Here we review recent advances in ligand discovery technologies, their potential for reshaping the whole process of drug discovery and development, as well as the challenges they encounter. We also discuss how the rapid identification of highly diverse, potent, target-selective and drug-like ligands to protein targets can democratize the drug discovery process, presenting new opportunities for the cost-effective development of safer and more effective small-molecule treatments.

SARS-CoV2 billion-compound docking

This dataset contains ligand conformations and docking scores for 1.4 billion molecules docked against 6 structural targets from SARS-CoV2, representing 5 unique proteins: MPro, NSP15, PLPro, RDRP, and the Spike protein. Docking was carried out using the AutoDock-GPU platform on the Summit supercomputer and Google Cloud. The docking procedure employed the Solis Wets search method to generate 20 independent ligand binding poses per compound. Each compound geometry was scored using the AutoDock free energy estimate, and rescored using RFScore v3 and DUD-E machine-learned rescoring models. Input protein structures are included, suitable for use by AutoDock-GPU and other docking programs. As the result of an exceptionally large docking campaign, this dataset represents a valuable resource for discovering trends across small molecule and protein binding sites, training AI models, and comparing to inhibitor compounds targeting SARS-CoV-2. The work also gives an example of how to organize and process data from ultra-large docking screens.

AI-Accelerated Design of Targeted Covalent Inhibitors for SARS-CoV-2

Direct-acting antivirals for the treatment of the COVID-19 pandemic caused by the SARS-CoV-2 virus are needed to complement vaccination efforts. Given the ongoing emergence of new variants, automated experimentation, and active learning based fast workflows for antiviral lead discovery remain critical to our ability to address the pandemic's evolution in a timely manner. While several such pipelines have been introduced to discover candidates with noncovalent interactions with the main protease (Mpro), here we developed a closed-loop artificial intelligence pipeline to design electrophilic warhead-based covalent candidates. This work introduces a deep learning-assisted automated computational workflow to introduce linkers and an electrophilic "warhead" to design covalent candidates and incorporates cutting-edge experimental techniques for validation. Using this process, promising candidates in the library were screened, and several potential hits were identified and tested experimentally using native mass spectrometry and fluorescence resonance energy transfer (FRET)-based screening assays. We identified four chloroacetamide-based covalent inhibitors of Mpro with micromolar affinities (KI of 5.27 μM) using our pipeline. Experimentally resolved binding modes for each compound were determined using room-temperature X-ray crystallography, which is consistent with the predicted poses. The induced conformational changes based on molecular dynamics simulations further suggest that the dynamics may be an important factor to further improve selectivity, thereby effectively lowering KI and reducing toxicity. These results demonstrate the utility of our modular and data-driven approach for potent and selective covalent inhibitor discovery and provide a platform to apply it to other emerging targets.

Lessons Learnt from COVID-19: Computational Strategies for Facing Present and Future Pandemics

Since its outbreak in December 2019, the COVID-19 pandemic has caused the death of more than 6.5 million people around the world. The high transmissibility of its causative agent, the SARS-CoV-2 virus, coupled with its potentially lethal outcome, provoked a profound global economic and social crisis. The urgency of finding suitable pharmacological tools to tame the pandemic shed light on the ever-increasing importance of computer simulations in rationalizing and speeding up the design of new drugs, further stressing the need for developing quick and reliable methods to identify novel active molecules and characterize their mechanism of action. In the present work, we aim at providing the reader with a general overview of the COVID-19 pandemic, discussing the hallmarks in its management, from the initial attempts at drug repurposing to the commercialization of Paxlovid, the first orally available COVID-19 drug. Furthermore, we analyze and discuss the role of computer-aided drug discovery (CADD) techniques, especially those that fall in the structure-based drug design (SBDD) category, in facing present and future pandemics, by showcasing several successful examples of drug discovery campaigns where commonly used methods such as docking and molecular dynamics have been employed in the rational design of effective therapeutic entities against COVID-19.

AI-accelerated protein-ligand docking for SARS-CoV-2 is 100-fold faster with no significant change in detection

Protein-ligand docking is a computational method for identifying drug leads. The method is capable of narrowing a vast library of compounds down to a tractable size for downstream simulation or experimental testing and is widely used in drug discovery. While there has been progress in accelerating scoring of compounds with artificial intelligence, few works have bridged these successes back to the virtual screening community in terms of utility and forward-looking development. We demonstrate the power of high-speed ML models by scoring 1 billion molecules in under a day (50 k predictions per GPU seconds). We showcase a workflow for docking utilizing surrogate AI-based models as a pre-filter to a standard docking workflow. Our workflow is ten times faster at screening a library of compounds than the standard technique, with an error rate less than 0.01% of detecting the underlying best scoring 0.1% of compounds. Our analysis of the speedup explains that another order of magnitude speedup must come from model accuracy rather than computing speed. In order to drive another order of magnitude of acceleration, we share a benchmark dataset consisting of 200 million 3D complex structures and 2D structure scores across a consistent set of 13 million "in-stock" molecules over 15 receptors, or binding sites, across the SARS-CoV-2 proteome. We believe this is strong evidence for the community to begin focusing on improving the accuracy of surrogate models to improve the ability to screen massive compound libraries 100 × or even 1000 × faster than current techniques and reduce missing top hits. The technique outlined aims to be a fast drop-in replacement for docking for screening billion-scale molecular libraries.

Speed vs Accuracy: Effect on Ligand Pose Accuracy of Varying Box Size and Exhaustiveness in AutoDock Vina

Structure-based virtual high-throughput screening involves docking chemical libraries to targets of interest. A parameter pertinent to the accuracy of the resulting pose is the root mean square deviation (RMSD) from a known crystallographic structure, i. e., the 'docking power'. Here, using a popular algorithm, Autodock Vina, as a model program, we evaluate the effects of varying two common docking parameters: the box size (the size of docking search space) and the exhaustiveness of the global search (the number of independent runs starting from random ligand conformations) on the RMSD from the PDBbind v2017 refined dataset of experimental protein-ligand complexes. Although it is clear that exhaustiveness is an important parameter, there is wide variation in the values used, with variation between 1 and >100. We, therefore, evaluated a combination of cubic boxes of different sizes and five exhaustiveness values (1, 8, 25, 50, 75, 100) within the range of those commonly adopted. The results show that the default exhaustiveness value of 8 performs well overall for most box sizes. In contrast, for all box sizes, but particularly for large boxes, an exhaustiveness value of 1 led to significantly higher median RMSD (mRMSD) values. The docking power was slightly improved with an exhaustiveness of 25, but the mRMSD changes little with values higher than 25. Therefore, although low exhaustiveness is computationally faster, the results are more likely to be far from reality, and, conversely, values >25 led to little improvement at the expense of computational resources. Overall, we recommend users to use at least the default exhaustiveness value of 8 for virtual screening calculations.

A review on structural, non-structural, and accessory proteins of SARS-CoV-2: Highlighting drug target sites

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the causative agent of COVID-19, is a highly transmittable and pathogenic human coronavirus that first emerged in China in December 2019. The unprecedented outbreak of SARS-CoV-2 devastated human health within a short time leading to a global public health emergency. A detailed understanding of the viral proteins including their structural characteristics and virulence mechanism on human health is very crucial for developing vaccines and therapeutics. To date, over 1800 structures of non-structural, structural, and accessory proteins of SARS-CoV-2 are determined by cryo-electron microscopy, X-ray crystallography, and NMR spectroscopy. Designing therapeutics to target the viral proteins has several benefits since they could be highly specific against the virus while maintaining minimal detrimental effects on humans. However, for ongoing and future research on SARS-CoV-2, summarizing all the viral proteins and their detailed structural information is crucial. In this review, we compile comprehensive information on viral structural, non-structural, and accessory proteins structures with their binding and catalytic sites, different domain and motifs, and potential drug target sites to assist chemists, biologists, and clinicians finding necessary details for fundamental and therapeutic research.

Valproate-coenzyme A conjugate blocks opening of receptor binding domains in the spike trimer of SARS-CoV-2 through an allosteric mechanism

The receptor-binding domains (RBDs) of the SARS-CoV-2 spike trimer exhibit "up" and "down" conformations often targeted by neutralizing antibodies. Only in the "up" configuration can RBDs bind to the ACE2 receptor of the host cell and initiate the process of viral multiplication. Here, we identify a lead compound (3-oxo-valproate-coenzyme A conjugate or Val-CoA) that stabilizes the spike trimer with RBDs in the down conformation. Val-CoA interacts with three R408 residues, one from each RBD, which significantly reduces the inter-subunit R408-R408 distance by ∼ 13 Å and closes the central pore formed by the three RBDs. Experimental evidence is presented that R408 is part of a triggering mechanism that controls the prefusion to postfusion state transition of the spike trimer. By stabilizing the RBDs in the down configuration, this and other related compounds can likely attenuate viral transmission. The reported findings for binding of Val-CoA to the spike trimer suggest a new approach for the design of allosteric antiviral drugs that do not have to compete for specific virus-receptor interactions but instead hinder the conformational motion of viral membrane proteins essential for interaction with the host cell. Here, we introduce an approach to target the spike protein by identifying lead compounds that stabilize the RBDs in the trimeric "down" configuration. When these compounds trimerize monomeric RBD immunogens as co-immunogens, they could also induce new types of non-ACE2 blocking antibodies that prevent local cell-to-cell transmission of the virus, providing a novel approach for inhibition of SARS-CoV-2.

Recent changes in the mutational dynamics of the SARS-CoV-2 main protease substantiate the danger of emerging resistance to antiviral drugs

Introduction

The current coronavirus pandemic is being combated worldwide by nontherapeutic measures and massive vaccination programs. Nevertheless, therapeutic options such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) main-protease (Mpro) inhibitors are essential due to the ongoing evolution toward escape from natural or induced immunity. While antiviral strategies are vulnerable to the effects of viral mutation, the relatively conserved Mpro makes an attractive drug target: Nirmatrelvir, an antiviral targeting its active site, has been authorized for conditional or emergency use in several countries since December 2021, and a number of other inhibitors are under clinical evaluation. We analyzed recent SARS-CoV-2 genomic data, since early detection of potential resistances supports a timely counteraction in drug development and deployment, and discovered accelerated mutational dynamics of Mpro since early December 2021.

Methods

We performed a comparative analysis of 10.5 million SARS-CoV-2 genome sequences available by June 2022 at GISAID to the NCBI reference genome sequence NC_045512.2. Amino-acid exchanges within high-quality regions in 69,878 unique Mpro sequences were identified and time- and in-depth sequence analyses including a structural representation of mutational dynamics were performed using in-house software.

Results

The analysis showed a significant recent event of mutational dynamics in Mpro. We report a remarkable increase in mutational variability in an eight-residue long consecutive region (R188-G195) near the active site since December 2021.

Discussion

The increased mutational variability in close proximity to an antiviral-drug binding site as described herein may suggest the onset of the development of antiviral resistance. This emerging diversity urgently needs to be further monitored and considered in ongoing drug development and lead optimization.

Therapeutic potential of compounds targeting SARS-CoV-2 helicase

The economical and societal impact of COVID-19 has made the development of vaccines and drugs to combat SARS-CoV-2 infection a priority. While the SARS-CoV-2 spike protein has been widely explored as a drug target, the SARS-CoV-2 helicase (nsp13) does not have any approved medication. The helicase shares 99.8% similarity with its SARS-CoV-1 homolog and was shown to be essential for viral replication. This review summarizes and builds on existing research on inhibitors of SARS-CoV-1 and SARS-CoV-2 helicases. Our analysis on the toxicity and specificity of these compounds, set the road going forward for the repurposing of existing drugs and the development of new SARS-CoV-2 helicase inhibitors.

Discovering new potential inhibitors to SARS-CoV-2 RNA dependent RNA polymerase (RdRp) using high throughput virtual screening and molecular dynamics simulations

RNA dependent RNA polymerase (RdRp), is an essential in the RNA replication within the life cycle of the severely acute respiratory coronavirus-2 (SARS-CoV-2), causing the deadly respiratory induced sickness COVID-19. Remdesivir is a prodrug that has seen some success in inhibiting this enzyme, however there is still the pressing need for effective alternatives. In this study, we present the discovery of four non-nucleoside small molecules that bind favorably to SARS-CoV-2 RdRp over the active form of the popular drug remdesivir (RTP) and adenosine triphosphate (ATP) by utilizing high-throughput virtual screening (HTVS) against the vast ZINC compound database coupled with extensive molecular dynamics (MD) simulations. After post-trajectory analysis, we found that the simulations of complexes containing both ATP and RTP remained stable for the duration of their trajectories. Additionally, it was revealed that the phosphate tail of RTP was stabilized by both the positive amino acid pocket and magnesium ions near the entry channel of RdRp which includes residues K551, R553, R555 and K621. It was also found that residues D623, D760, and N691 further stabilized the ribose portion of RTP with U10 on the template RNA strand forming hydrogen pairs with the adenosine motif. Using these models of RdRp, we employed them to screen the ZINC database of ~ 17 million molecules. Using docking and drug properties scoring, we narrowed down our selection to fourteen candidates. These were subjected to 200 ns simulations each underwent free energy calculations. We identified four hit compounds from the ZINC database that have similar binding poses to RTP while possessing lower overall binding free energies, with ZINC097971592 having a binding free energy two times lower than RTP.

Chemical space docking enables large-scale structure-based virtual screening to discover ROCK1 kinase inhibitors

With the ever-increasing number of synthesis-on-demand compounds for drug lead discovery, there is a great need for efficient search technologies. We present the successful application of a virtual screening method that combines two advances: (1) it avoids full library enumeration (2) products are evaluated by molecular docking, leveraging protein structural information. Crucially, these advances enable a structure-based technique that can efficiently explore libraries with billions of molecules and beyond. We apply this method to identify inhibitors of ROCK1 from almost one billion commercially available compounds. Out of 69 purchased compounds, 27 (39%) have K_i values < 10 µM. X-ray structures of two leads confirm their docked poses. This approach to docking scales roughly with the number of reagents that span a chemical space and is therefore multiple orders of magnitude faster than traditional docking.

Direct Interaction of Coronavirus Nonstructural Protein 3 with Melanoma Differentiation-Associated Gene 5 Modulates Type I Interferon Response during Coronavirus Infection

Coronavirus nonstructural protein 3 (nsp3) is a multi-functional protein, playing a critical role in viral replication and in regulating host antiviral innate immunity. In this study, we demonstrate that nsp3 from severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and avian coronavirus infectious bronchitis virus (IBV) directly interacts with melanoma differentiation-associated gene 5 (MDA5), rendering an inhibitory effect on the MDA5-mediated type I interferon (IFN) response. By the co-expression of MDA5 with wild-type and truncated nsp3 constructs, at least three interacting regions mapped to the papain-like protease (PLpro) domain and two other domains located at the N- and C-terminal regions were identified in SARS-CoV-2 nsp3. Furthermore, by introducing point mutations to the catalytic triad, the deubiquitylation activity of the PLpro domain from both SARS-CoV-2 and IBV nsp3 was shown to be responsible for the suppression of the MDA5-mediated type I IFN response. It was also demonstrated that both MDA5 and nsp3 were able to interact with ubiquitin and ubiquitinated proteins, contributing to the interaction between the two proteins. This study confirms the antagonistic role of nsp3 in the MDA5-mediated type I IFN signaling, highlighting the complex interaction between a multi-functional viral protein and the innate immune response.

High-resolution structures of the SARS-CoV-2 N7-methyltransferase inform therapeutic development

Emergence of SARS-CoV-2 coronavirus has led to millions of deaths globally. We present three high-resolution crystal structures of the SARS-CoV-2 nsp14 N7-methyltransferase core bound to S-adenosylmethionine (1.62 Å), S-adenosylhomocysteine (1.55 Å) and sinefungin (1.41 Å). We identify features of the methyltransferase core that are crucial for the development of antivirals and show SAH as the best scaffold for the design of antivirals against SARS-CoV-2 and other pathogenic coronaviruses.

Structural bioinformatics analysis of SARS-CoV-2 variants reveals higher hACE2 receptor binding affinity for Omicron B.1.1.529 spike RBD compared to wild type reference

To date, more than 263 million people have been infected with SARS-CoV-2 during the COVID-19 pandemic. In many countries, the global spread occurred in multiple pandemic waves characterized by the emergence of new SARS-CoV-2 variants. Here we report a sequence and structural-bioinformatics analysis to estimate the effects of amino acid substitutions on the affinity of the SARS-CoV-2 spike receptor binding domain (RBD) to the human receptor hACE2. This is done through qualitative electrostatics and hydrophobicity analysis as well as molecular dynamics simulations used to develop a high-precision empirical scoring function (ESF) closely related to the linear interaction energy method and calibrated on a large set of experimental binding energies. For the latest variant of concern (VOC), B.1.1.529 Omicron, our Halo difference point cloud studies reveal the largest impact on the RBD binding interface compared to all other VOC. Moreover, according to our ESF model, Omicron achieves a much higher ACE2 binding affinity than the wild type and, in particular, the highest among all VOCs except Alpha and thus requires special attention and monitoring.

When Virtual Screening Yields Inactive Drugs: Dealing with False Theoretical Friends

The search of antivirals against SARS-CoV-2 in available libraries of compounds was initiated as soon as WHO announced that the coronavirus outbreak became a pandemic. That pivotal task has been conducted by both experimental groups in wet-labs as well as by theoretical chemists in supercomputing centers. The combination of biochemical and clinical intuitions yields first to remdesivir, a broad-spectrum antiviral that remains as the standard solution for the treatment of severe cases, while paxlovid, molnupiravir and fluvoxamine have been recently proposed as oral alternatives. Unfortunately, the intensive publication of standard virtual screening (VS) simulations might be not the best strategy to increase that short list of antivirals. This contribution joins theory and biological assays to rescore massive VS. Our goal is to critically assess pros and cons of using molecular models for drug repurposing.

Interaction of panduratin A and derivatives with the SARS-CoV-2 main protease (mpro): a molecular docking study

Panduratin A (Pa-A) is a prenylated cyclohexenyl chalcone isolated from the rhizomes of the medicinal and culinary plant Boesenbergia rotunda (L.) Mansf., commonly called fingerroots. Both an ethanolic plant extract and Pa-A have shown a marked antiviral activity against the severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2), responsible for the COVID-19 pandemic disease. Pa-A functions as a protease inhibitor inhibiting infection of human cells by the virus. We have modeled the interaction of Pa-A, and 26 panduratin analogues with the main protease (M^pro) of SARS-CoV-2 using molecular docking. The natural product 4-hydroxypanduratin showed a higher M^pro binding capacity than Pa-A and isopanduratin A. The interaction with M^Pro of all known panduratin derivatives (Pa-A to Pa-Y) have been compared, together with more than 60 reference products. Three compounds emerged as potential robust M^Pro binders: Pa-R, Pa-V, Pa-S, with a binding capacity significantly higher than 4-OH-Pa-A and Pa-A. The empirical energy of interaction (ΔE) calculated with the best compound in the panduratin series, Pa-R bound to M^pro, surpassed that measured with the top reference protease inhibitors such a ruprintrivir, lufotrelvir, and glecaprevir. Structure-binding relationships are discussed. Compounds with a flavanone moiety (PA-R/S) are the best binders, better than those with a chromene unit (Pa-F/G). The extended molecules (such as Pa-V) exhibit good M^pro binding, but the dimeric compound Pa-Y is too long and protrudes outside the binding cavity. The work provides novel ideas to guide the design of new molecules interacting with M^pro.Communicated by Ramaswamy H. Sarma.

Kinetic Characterization of SARS-CoV-2 nsp13 ATPase Activity and Discovery of Small-Molecule Inhibitors

SARS-CoV-2 non-structural protein 13 (nsp13) is a highly conserved helicase and RNA 5'-triphosphatase. It uses the energy derived from the hydrolysis of nucleoside triphosphates for directional movement along the nucleic acids and promotes the unwinding of double-stranded nucleic acids. Nsp13 is essential for replication and propagation of all human and non-human coronaviruses. Combined with its defined nucleotide binding site and druggability, nsp13 is one of the most promising candidates for the development of pan-coronavirus therapeutics. Here, we report the development and optimization of bioluminescence assays for kinetic characterization of nsp13 ATPase activity in the presence and absence of single-stranded DNA. Screening of a library of 5000 small molecules in the presence of single-stranded DNA resulted in the discovery of six nsp13 small-molecule inhibitors with IC50 values ranging from 6 ± 0.5 to 50 ± 6 μM. In addition to providing validated methods for high-throughput screening of nsp13 in drug discovery campaigns, the reproducible screening hits we present here could potentially be chemistry starting points toward the development of more potent and selective nsp13 inhibitors, enabling the discovery of antiviral therapeutics.

Emerging frontiers in virtual drug discovery: From quantum mechanical methods to deep learning approaches

Virtual screening-based approaches to discover initial hit and lead compounds have the potential to reduce both the cost and time of early drug discovery stages, as well as to find inhibitors for even challenging target sites such as protein-protein interfaces. Here in this review, we provide an overview of the progress that has been made in virtual screening methodology and technology on multiple fronts in recent years. The advent of ultra-large virtual screens, in which hundreds of millions to billions of compounds are screened, has proven to be a powerful approach to discover highly potent hit compounds. However, these developments are just the tip of the iceberg, with new technologies and methods emerging to propel the field forward. Examples include novel machine-learning approaches, which can reduce the computational costs of virtual screening dramatically, while progress in quantum-mechanical approaches can increase the accuracy of predictions of various small molecule properties.

The Main Protease of SARS-CoV-2 as a Target for Phytochemicals against Coronavirus

In late December 2019, the first cases of COVID-19 emerged as an outbreak in Wuhan, China that later spread vastly around the world, evolving into a pandemic and one of the worst global health crises in modern history. The causative agent was identified as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Although several vaccines were authorized for emergency use, constantly emerging new viral mutants and limited treatment options for COVID-19 drastically highlighted the need for developing an efficient treatment for this disease. One of the most important viral components to target for this purpose is the main protease of the coronavirus (Mpro). This enzyme is an excellent target for a potential drug, as it is essential for viral replication and has no closely related homologues in humans, making its inhibitors unlikely to be toxic. Our review describes a variety of approaches that could be applied in search of potential inhibitors among plant-derived compounds, including virtual in silico screening (a data-driven approach), which could be structure-based or fragment-guided, the classical approach of high-throughput screening, and antiviral activity cell-based assays. We will focus on several classes of compounds reported to be potential inhibitors of Mpro, including phenols and polyphenols, alkaloids, and terpenoids.

The SARS-CoV-2 helicase as a target for antiviral therapy: Identification of potential small molecule inhibitors by in silico modelling

Although vaccines that provide protection against severe illness from coronavirus disease (COVID-19) have been made available, emerging variant strains of severe acute respiratory syndrome 2 coronavirus 2 (SARS-CoV-2) are of concern. A different research direction involves investigation of antiviral therapeutics. In addition to structural proteins, the SARS-CoV-2 non-structural proteins are of interest and this includes the helicase (nsp13). In this study, an initial screen of 300 ligands was performed to identify potential inhibitors of the SARS-CoV-2 nsp13 examining the nucleoside triphosphatase site (NTPase activity) as the target region. The antiviral activity of polyphenols has been previously reported in the literature and as a result, the phenolic compounds and fatty acids from the OliveNet™ library were utilised. Synthetic compounds with antimicrobial and anti-inflammatory properties were also selected. The structures of the SARS-CoV and MERS-CoV helicases, as well as the human RECQ-like DNA helicase, DHX9 helicase, PcrA helicase, hepatitis C NS3 helicase, and mouse Dna2 nuclease-helicase were used for comparison. As expected, sequence and structural homology between the various species was evident. A number of broad-spectrum and well-known inhibitors interacted with the NTPase active site highlighting the need to potentially identify more specific inhibitors for SARS-CoV-2. Acetylcysteine, clavulanic acid and homovanillic acid were identified as potential lead compounds for the SARS-CoV-2 helicase. Molecular dynamics simulations were performed with the leads bound to the SARS-CoV-2 helicase for 200 ns in triplicate, with favourable binding free energies to the NTPase site. Given their availability, further exploration of their potential inhibitory activity could be considered.

Antiviral Drug Discovery for the Treatment of COVID-19 Infections

The coronavirus disease 2019 (COVID-19) pandemic is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), a recently emerged human coronavirus. COVID-19 vaccines have proven to be successful in protecting the vaccinated from infection, reducing the severity of disease, and deterring the transmission of infection. However, COVID-19 vaccination faces many challenges, such as the decline in vaccine-induced immunity over time, and the decrease in potency against some SARS-CoV-2 variants including the recently emerged Omicron variant, resulting in breakthrough infections. The challenges that COVID-19 vaccination is facing highlight the importance of the discovery of antivirals to serve as another means to tackle the pandemic. To date, neutralizing antibodies that block viral entry by targeting the viral spike protein make up the largest class of antivirals that has received US FDA emergency use authorization (EUA) for COVID-19 treatment. In addition to the spike protein, other key targets for the discovery of direct-acting antivirals include viral enzymes that are essential for SARS-CoV-2 replication, such as RNA-dependent RNA polymerase and proteases, as judged by US FDA approval for remdesivir, and EUA for Paxlovid (nirmatrelvir + ritonavir) for treating COVID-19 infections. This review presents an overview of the current status and future direction of antiviral drug discovery for treating SARS-CoV-2 infections, covering important antiviral targets such as the viral spike protein, non-structural protein (nsp) 3 papain-like protease, nsp5 main protease, and the nsp12/nsp7/nsp8 RNA-dependent RNA polymerase complex.

Opportunities and Challenges in Targeting the Proofreading Activity of SARS-CoV-2 Polymerase Complex

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) is the cause of the COVID-19 pandemic. While the development of vaccines and the emergence of antiviral therapeutics is promising, alternative strategies to combat COVID-19 (and potential future pandemics) remain an unmet need. Coronaviruses feature a unique mechanism that may present opportunities for therapeutic intervention: the RNA polymerase complex of coronaviruses is distinct in its ability to proofread and remove mismatched nucleotides during genome replication and transcription. The proofreading activity has been linked to the exonuclease (ExoN) activity of non-structural protein 14 (NSP14). Here, we review the role of NSP14, and other NSPs, in SARS-CoV-2 replication and describe the assays that have been developed to assess the ExoN function. We also review the nucleoside analogs and non-nucleoside inhibitors known to interfere with the proofreading activity of NSP14. Although not yet validated, the potential use of non-nucleoside proofreading inhibitors in combination with chain-terminating nucleosides may be a promising avenue for the development of anti-CoV agents.

Chemical biology and medicinal chemistry of RNA methyltransferases

RNA methyltransferases (MTases) are ubiquitous enzymes whose hitherto low profile in medicinal chemistry, contrasts with the surging interest in RNA methylation, the arguably most important aspect of the new field of epitranscriptomics. As MTases become validated as drug targets in all major fields of biomedicine, the development of small molecule compounds as tools and inhibitors is picking up considerable momentum, in academia as well as in biotech. Here we discuss the development of small molecules for two related aspects of chemical biology. Firstly, derivates of the ubiquitous cofactor S-adenosyl-l-methionine (SAM) are being developed as bioconjugation tools for targeted transfer of functional groups and labels to increasingly visible targets. Secondly, SAM-derived compounds are being investigated for their ability to act as inhibitors of RNA MTases. Drug development is moving from derivatives of cosubstrates towards higher generation compounds that may address allosteric sites in addition to the catalytic centre. Progress in assay development and screening techniques from medicinal chemistry have led to recent breakthroughs, e.g. in addressing human enzymes targeted for their role in cancer. Spurred by the current pandemic, new inhibitors against coronaviral MTases have emerged at a spectacular rate, including a repurposed drug which is now in clinical trial.

GROMACS in the Cloud: A Global Supercomputer to Speed Up Alchemical Drug Design

We assess costs and efficiency of state-of-the-art high-performance cloud computing and compare the results to traditional on-premises compute clusters. Our use case is atomistic simulations carried out with the GROMACS molecular dynamics (MD) toolkit with a particular focus on alchemical protein–ligand binding free energy calculations. We set up a compute cluster in the Amazon Web Services (AWS) cloud that incorporates various different instances with Intel, AMD, and ARM CPUs, some with GPU acceleration. Using representative biomolecular simulation systems, we benchmark how GROMACS performs on individual instances and across multiple instances. Thereby we assess which instances deliver the highest performance and which are the most cost-efficient ones for our use case. We find that, in terms of total costs, including hardware, personnel, room, energy, and cooling, producing MD trajectories in the cloud can be about as cost-efficient as an on-premises cluster given that optimal cloud instances are chosen. Further, we find that high-throughput ligand-screening can be accelerated dramatically by using global cloud resources. For a ligand screening study consisting of 19 872 independent simulations or ∼200 μs of combined simulation trajectory, we made use of diverse hardware available in the cloud at the time of the study. The computations scaled-up to reach peak performance using more than 4 000 instances, 140 000 cores, and 3 000 GPUs simultaneously. Our simulation ensemble finished in about 2 days in the cloud, while weeks would be required to complete the task on a typical on-premises cluster consisting of several hundred nodes.

Non-covalent SARS-CoV-2 Mpro inhibitors developed from in silico screen hits

M^pro, the main protease of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is essential for the viral life cycle. Accordingly, several groups have performed in silico screens to identify M^pro inhibitors that might be used to treat SARS-CoV-2 infections. We selected more than five hundred compounds from the top-ranking hits of two very large in silico screens for on-demand synthesis. We then examined whether these compounds could bind to M^pro and inhibit its protease activity. Two interesting chemotypes were identified, which were further evaluated by characterizing an additional five hundred synthesis on-demand analogues. The compounds of the first chemotype denatured M^pro and were considered not useful for further development. The compounds of the second chemotype bound to and enhanced the melting temperature of M^pro. The most active compound from this chemotype inhibited M^pro in vitro with an IC₅₀ value of 1 μM and suppressed replication of the SARS-CoV-2 virus in tissue culture cells. Its mode of binding to M^pro was determined by X-ray crystallography, revealing that it is a non-covalent inhibitor. We propose that the inhibitors described here could form the basis for medicinal chemistry efforts that could lead to the development of clinically relevant inhibitors.

Ultralarge Virtual Screening Identifies SARS-CoV-2 Main Protease Inhibitors with Broad-Spectrum Activity against Coronaviruses

Drugs targeting SARS-CoV-2 could have saved millions of lives during the COVID-19 pandemic, and it is now crucial to develop inhibitors of coronavirus replication in preparation for future outbreaks. We explored two virtual screening strategies to find inhibitors of the SARS-CoV-2 main protease in ultralarge chemical libraries. First, structure-based docking was used to screen a diverse library of 235 million virtual compounds against the active site. One hundred top-ranked compounds were tested in binding and enzymatic assays. Second, a fragment discovered by crystallographic screening was optimized guided by docking of millions of elaborated molecules and experimental testing of 93 compounds. Three inhibitors were identified in the first library screen, and five of the selected fragment elaborations showed inhibitory effects. Crystal structures of target-inhibitor complexes confirmed docking predictions and guided hit-to-lead optimization, resulting in a noncovalent main protease inhibitor with nanomolar affinity, a promising in vitro pharmacokinetic profile, and broad-spectrum antiviral effect in infected cells.

High-Throughput Virtual Screening and Validation of a SARS-CoV-2 Main Protease Noncovalent Inhibitor

Despite the recent availability of vaccines against the acute respiratory syndrome coronavirus 2 (SARS-CoV-2), the search for inhibitory therapeutic agents has assumed importance especially in the context of emerging new viral variants. In this paper, we describe the discovery of a novel noncovalent small-molecule inhibitor, MCULE-5948770040, that binds to and inhibits the SARS-Cov-2 main protease (Mpro) by employing a scalable high-throughput virtual screening (HTVS) framework and a targeted compound library of over 6.5 million molecules that could be readily ordered and purchased. Our HTVS framework leverages the U.S. supercomputing infrastructure achieving nearly 91% resource utilization and nearly 126 million docking calculations per hour. Downstream biochemical assays validate this Mpro inhibitor with an inhibition constant (Ki) of 2.9 μM (95% CI 2.2, 4.0). Furthermore, using room-temperature X-ray crystallography, we show that MCULE-5948770040 binds to a cleft in the primary binding site of Mpro forming stable hydrogen bond and hydrophobic interactions. We then used multiple μs-time scale molecular dynamics (MD) simulations and machine learning (ML) techniques to elucidate how the bound ligand alters the conformational states accessed by Mpro, involving motions both proximal and distal to the binding site. Together, our results demonstrate how MCULE-5948770040 inhibits Mpro and offers a springboard for further therapeutic design.

RNA helicases required for viral propagation in humans

RNA viruses cause many routine illnesses, such as the common cold and the flu. Recently, more deadly diseases have emerged from this family of viruses. The hepatitis C virus has had a devastating impact worldwide. Despite the cures developed in the U.S. and Europe, economically disadvantaged countries remain afflicted by HCV infection due to the high cost of these medications. More recently, COVID-19 has swept across the world, killing millions and disrupting economies and lifestyles; the virus responsible for this pandemic is a coronavirus. Our understanding of HCV and SARS CoV-2 replication is still in its infancy. Helicases play a critical role in the replication, transcription and translation of viruses. These key enzymes need extensive study not only as an essential player in the viral lifecycle, but also as targets for antiviral therapeutics. In this review, we highlight the current knowledge for RNA helicases of high importance to human health.

Discovery of Zafirlukast as a novel SARS-CoV-2 helicase inhibitor using in silico modelling and a FRET-based assay

The coronavirus helicase is an essential enzyme required for viral replication/transcription pathways. Structural studies revealed a sulphate moiety that interacts with key residues within the nucleotide-binding site of the helicase. Compounds with a sulphoxide or a sulphone moiety could interfere with these interactions and consequently inhibit the enzyme. The molecular operating environment (MOE) was used to dock 189 sulphoxide and sulphone-containing FDA-approved compounds to the nucleotide-binding site. Zafirlukast, a leukotriene receptor antagonist used to treat chronic asthma, achieved the lowest docking score at -8.75 kcals/mol. The inhibitory effect of the compounds on the SARS-CoV-2 helicase dsDNA unwinding activity was tested by a FRET-based assay. Zafirlukast was the only compound to inhibit the enzyme (IC50 = 16.3 µM). The treatment of Vero E6 cells with 25 µM zafirlukast prior to SARS-CoV-2 infection decreased the cytopathic effects of SARS-CoV-2 significantly. These results suggest that zafirlukast alleviates SARS-CoV-2 pathogenicity by inhibiting the viral helicase and impairing the viral replication/transcription pathway. Zafirlukast could be clinically developed as a new antiviral treatment for SARS-CoV-2 and other coronavirus diseases. This discovery is based on molecular modelling, in vitro inhibition of the SARS-CoV helicase activity and cell-based SARS-CoV-2 viral replication.

Investigation of Thiocarbamates as Potential Inhibitors of the SARS-CoV-2 Mpro

In the present study we tested, using the microscale thermophoresis technique, a small library of thionocarbamates, thiolocarbamates, sulfide and disulfide as potential lead compounds for SARS-CoV-2 Mpro drug design. The successfully identified binder is a representative of the thionocarbamates group with a high potential for future modifications aiming for higher affinity and solubility. The experimental analysis was extended by computational studies that show insufficient accuracy of the simplest and widely applied approaches and underline the necessity of applying more advanced methods to properly evaluate the affinity of potential SARS-CoV-2 Mpro binders.

Efficient Exploration of Chemical Space with Docking and Deep Learning

With the advent of make-on-demand commercial libraries, the number of purchasable compounds available for virtual screening and assay has grown explosively in recent years, with several libraries eclipsing one billion compounds. Today's screening libraries are larger and more diverse, enabling the discovery of more-potent hit compounds and unlocking new areas of chemical space, represented by new core scaffolds. Applying physics-based in silico screening methods in an exhaustive manner, where every molecule in the library must be enumerated and evaluated independently, is increasingly cost-prohibitive. Here, we introduce a protocol for machine learning-enhanced molecular docking based on active learning to dramatically increase throughput over traditional docking. We leverage a novel selection protocol that strikes a balance between two objectives: (1) identifying the best scoring compounds and (2) exploring a large region of chemical space, demonstrating superior performance compared to a purely greedy approach. Together with automated redocking of the top compounds, this method captures almost all the high scoring scaffolds in the library found by exhaustive docking. This protocol is applied to our recent virtual screening campaigns against the D4 and AMPC targets that produced dozens of highly potent, novel inhibitors, and a blind test against the MT1 target. Our protocol recovers more than 80% of the experimentally confirmed hits with a 14-fold reduction in compute cost, and more than 90% of the hit scaffolds in the top 5% of model predictions, preserving the diversity of the experimentally confirmed hit compounds.

Biomolecular modeling thrives in the age of technology

The biomolecular modeling field has flourished since its early days in the 1970s due to the rapid adaptation and tailoring of state-of-the-art technology. The resulting dramatic increase in size and timespan of biomolecular simulations has outpaced Moore's law. Here, we discuss the role of knowledge-based versus physics-based methods and hardware versus software advances in propelling the field forward. This rapid adaptation and outreach suggests a bright future for modeling, where theory, experimentation and simulation define three pillars needed to address future scientific and biomedical challenges.

Resources and computational strategies to advance small molecule SARS-CoV-2 discovery: Lessons from the pandemic and preparing for future health crises

There is an urgent need to identify new therapies that prevent SARS-CoV-2 infection and improve the outcome of COVID-19 patients. This pandemic has thus spurred intensive research in most scientific areas and in a short period of time, several vaccines have been developed. But, while the race to find vaccines for COVID-19 has dominated the headlines, other types of therapeutic agents are being developed. In this mini-review, we report several databases and online tools that could assist the discovery of anti-SARS-CoV-2 small chemical compounds and peptides. We then give examples of studies that combined in silico and in vitro screening, either for drug repositioning purposes or to search for novel bioactive compounds. Finally, we question the overall lack of discussion and plan observed in academic research in many countries during this crisis and suggest that there is room for improvement.

Serine 477 plays a crucial role in the interaction of the SARS-CoV-2 spike protein with the human receptor ACE2

Since the worldwide outbreak of the infectious disease COVID-19, several studies have been published to understand the structural mechanism of the novel coronavirus SARS-CoV-2. During the infection process, the SARS-CoV-2 spike (S) protein plays a crucial role in the receptor recognition and cell membrane fusion process by interacting with the human angiotensin-converting enzyme 2 (hACE2) receptor. However, new variants of these spike proteins emerge as the virus passes through the disease reservoir. This poses a major challenge for designing a potent antigen for an effective immune response against the spike protein. Through a normal mode analysis (NMA) we identified the highly flexible region in the receptor binding domain (RBD) of SARS-CoV-2, starting from residue 475 up to residue 485. Structurally, the position S477 shows the highest flexibility among them. At the same time, S477 is hitherto the most frequently exchanged amino acid residue in the RBDs of SARS-CoV-2 mutants. Therefore, using MD simulations, we have investigated the role of S477 and its two frequent mutations (S477G and S477N) at the RBD during the binding to hACE2. We found that the amino acid exchanges S477G and S477N strengthen the binding of the SARS-COV-2 spike with the hACE2 receptor.