Ultra-large library docking for discovering new chemotypes

Dockey: a modern integrated tool for large-scale molecular docking and virtual screening

Molecular docking is a structure-based and computer-aided drug design approach that plays a pivotal role in drug discovery and pharmaceutical research. AutoDock is the most widely used molecular docking tool for study of protein-ligand interactions and virtual screening. Although many tools have been developed to streamline and automate the AutoDock docking pipeline, some of them still use outdated graphical user interfaces and have not been updated for a long time. Meanwhile, some of them lack cross-platform compatibility and evaluation metrics for screening lead compound candidates. To overcome these limitations, we have developed Dockey, a flexible and intuitive graphical interface tool with seamless integration of several useful tools, which implements a complete docking pipeline covering molecular sanitization, molecular preparation, paralleled docking execution, interaction detection and conformation visualization. Specifically, Dockey can detect the non-covalent interactions between small molecules and proteins and perform cross-docking between multiple receptors and ligands. It has the capacity to automatically dock thousands of ligands to multiple receptors and analyze the corresponding docking results in parallel. All the generated data will be kept in a project file that can be shared between any systems and computers with the pre-installation of Dockey. We anticipate that these unique characteristics will make it attractive for researchers to conduct large-scale molecular docking without complicated operations, particularly for beginners. Dockey is implemented in Python and freely available at https://github.com/lmdu/dockey.

Pharmacological characterization of novel small molecule agonists and antagonists for the orphan receptor GPR139

The orphan G protein-coupled receptor GPR139 is predominantly expressed in the central nervous system and has attracted considerable interest as a therapeutic target. However, the biological role of this receptor remains somewhat elusive, in part due to the lack of quality pharmacological tools to investigate GPR139 function. In an effort to understand GPR139 signaling and to identify improved compounds, in this study we performed virtual screening and analog searches, in combination with multiple pharmacological assays. We characterized GPR139-dependent signaling using previously published reference agonists in Ca²⁺ mobilization and inositol monophosphate accumulation assays, as well as a novel real-time GPR139 internalization assay. For the four reference agonists tested, the rank order of potency was conserved across signaling and internalization assays: JNJ-63533054 > Compound 1a » Takeda > AC4 > DL43, consistent with previously reported values. We noted an increased efficacy of JNJ-63533054-mediated inositol monophosphate signaling and internalization, relative to Compound 1a. We then performed virtual screening for GPR139 agonist and antagonist compounds that were screened and validated in GPR139 functional assays. We identified four GPR139 agonists that were active in all assays, with similar or reduced potency relative to known compounds. Likewise, compound analogs selected based on GPR139 agonist and antagonist substructure searches behaved similarly to their parent compounds. Thus, we have characterized GPR139 signaling for multiple new ligands using G protein-dependent assays and a new real-time internalization assay. These data add to the GPR139 tool compound repertoire, which could be optimized in future medical chemistry campaigns.

Target-driven machine learning-enabled virtual screening (TAME-VS) platform for early-stage hit identification

High-throughput screening (HTS) methods enable the empirical evaluation of a large scale of compounds and can be augmented by virtual screening (VS) techniques to save time and money by using potential active compounds for experimental testing. Structure-based and ligand-based virtual screening approaches have been extensively studied and applied in drug discovery practice with proven outcomes in advancing candidate molecules. However, the experimental data required for VS are expensive, and hit identification in an effective and efficient manner is particularly challenging during early-stage drug discovery for novel protein targets. Herein, we present our TArget-driven Machine learning-Enabled VS (TAME-VS) platform, which leverages existing chemical databases of bioactive molecules to modularly facilitate hit finding. Our methodology enables bespoke hit identification campaigns through a user-defined protein target. The input target ID is used to perform a homology-based target expansion, followed by compound retrieval from a large compilation of molecules with experimentally validated activity. Compounds are subsequently vectorized and adopted for machine learning (ML) model training. These machine learning models are deployed to perform model-based inferential virtual screening, and compounds are nominated based on predicted activity. Our platform was retrospectively validated across ten diverse protein targets and demonstrated clear predictive power. The implemented methodology provides a flexible and efficient approach that is accessible to a wide range of users. The TAME-VS platform is publicly available at https://github.com/bymgood/Target-driven-ML-enabled-VS to facilitate early-stage hit identification.

The European Federation for Medicinal Chemistry and Chemical Biology (EFMC) Best Practice Initiative: Hit Generation

Hit generation is a crucial step in drug discovery that will determine the speed and chance of success of identifying drug candidates. Many strategies are now available to identify chemical starting points, or hits, and each biological target warrants a tailored approach. In this set of best practices, we detail the essential approaches for target centric hit generation and the opportunities and challenges they come with. We then provide guidance on how to validate hits to ensure medicinal chemistry is only performed on compounds and scaffolds that engage the target of interest and have the desired mode of action. Finally, we discuss the design of integrated hit generation strategies that combine several approaches to maximize the chance of identifying high quality starting points to ensure a successful drug discovery campaign.

Targeting in silico GPCR conformations with ultra-large library screening for hit discovery

The use of deep machine learning (ML) in protein structure prediction has made it possible to easily access a large number of annotated conformations that can potentially compensate for missing experimental structures in structure-based drug discovery (SBDD). However, it is still unclear whether the accuracy of these predicted conformations is sufficient for screening chemical compounds that will effectively interact with a protein target for pharmacological purposes. In this opinion article, we examine the potential benefits and limitations of using state-annotated conformations for ultra-large library screening (ULLS) in light of the growing size of ultra-large libraries (ULLs). We believe that targeting different conformational states of common drug targets like G-protein-coupled receptors (GPCRs), which can regulate human physiology by switching between different conformations, can offer multiple advantages.

Developing novel antifungals: lessons from G protein-coupled receptors

Up to 1.5 million people die yearly from fungal disease, but the repertoire of antifungal drug classes is minimal and the incidence of drug resistance is rising rapidly. This dilemma was recently declared by the World Health Organization as a global health emergency, but the discovery of new antifungal drug classes remains excruciatingly slow. This process could be accelerated by focusing on novel targets, such as G protein-coupled receptor (GPCR)-like proteins, that have a high likelihood of being druggable and have well-defined biology and roles in disease. We discuss recent successes in understanding the biology of virulence and in structure determination of yeast GPCRs, and highlight new approaches that might pay significant dividends in the urgent search for novel antifungal drugs.

GB-score: Minimally designed machine learning scoring function based on distance-weighted interatomic contact features

In recent years, thanks to advances in computer hardware and dataset availability, data-driven approaches (like machine learning) have become one of the essential parts of the drug design framework to accelerate drug discovery procedures. Constructing a new scoring function, a function that can predict the binding score for a generated protein-ligand pose during docking procedure or a crystal complex, based on machine and deep learning has become an active research area in computer-aided drug design. GB-Score is a state-of-the-art machine learning-based scoring function that utilizes distance-weighted interatomic contact features, PDBbind-v2019 general set, and Gradient Boosting Trees algorithm to the binding affinity prediction. The distance-weighted interatomic contact featurization method used the distance between different ligand and protein atom types for numerical representation of the protein-ligand complex. GB-Score attains Pearson's correlation 0.862 and RMSE 1.190 on the CASF-2016 benchmark test in the scoring power metric. GB-Score's codes are freely available on the web at https://github.com/miladrayka/GB_Score.

Advances in the discovery of new chemotypes through ultra-large library docking

INTRODUCTION

The size and complexity of virtual screening libraries in drug discovery have skyrocketed in recent years, reaching up to multiple billions of accessible compounds. However, virtual screening of such ultra-large libraries poses several challenges associated with preparing the libraries, sampling, and pre-selection of suitable compounds. The utilization of artificial intelligence (AI)-assisted screening approaches, such as deep learning, poses a promising countermeasure to deal with this rapidly expanding chemical space. For example, various AI-driven methods were recently successfully used to identify novel small molecule inhibitors of the SARS-CoV-2 main protease (M^pro).

AREAS COVERED

This review focuses on presenting various kinds of virtual screening methods suitable for dealing with ultra-large libraries. Challenges associated with these computational methodologies are discussed, and recent advances are highlighted in the example of the discovery of novel M^pro inhibitors targeting the SARS-CoV-2 virus.

EXPERT OPINION

With the rapid expansion of the virtual chemical space, the methodologies for docking and screening such quantities of molecules need to keep pace. Employment of AI-driven screening compounds has already been shown to be effective in a range from a few thousand to multiple billion compounds, furthered by de novo generation of drug-like molecules without human interference.

Quantum Mechanical-Cluster Approach to Solve the Bioisosteric Replacement Problem in Drug Design

Bioisosteres are molecules that differ in substituents but still have very similar shapes. Bioisosteric replacements are ubiquitous in modern drug design, where they are used to alter metabolism, change bioavailability, or modify activity of the lead compound. Prediction of relative affinities of bioisosteres with computational methods is a long-standing task; however, the very shape closeness makes bioisosteric substitutions almost intractable for computational methods, which use standard force fields. Here, we design a quantum mechanical (QM)-cluster approach based on the GFN2-xTB semi-empirical quantum-chemical method and apply it to a set of H → F bioisosteric replacements. The proposed methodology enables advanced prediction of biological activity change upon bioisosteric substitution of -H with -F, with the standard deviation of 0.60 kcal/mol, surpassing the ChemPLP scoring function (0.83 kcal/mol), and making QM-based ΔΔG estimation comparable to ∼0.42 kcal/mol standard deviation of in vitro experiment. The speed of the method and lack of tunable parameters makes it affordable in current drug research.

ZINC-22─A Free Multi-Billion-Scale Database of Tangible Compounds for Ligand Discovery

Purchasable chemical space has grown rapidly into the tens of billions of molecules, providing unprecedented opportunities for ligand discovery but straining the tools that might exploit these molecules at scale. We have therefore developed ZINC-22, a database of commercially accessible small molecules derived from multi-billion-scale make-on-demand libraries. The new database and tools enable analog searching in this vast new space via a facile GUI, CartBlanche, drawing on similarity methods that scale sublinearly in the number of molecules. The new library also uses data organization methods, enabling rapid lookup of molecules and their physical properties, including conformations, partial atomic charges, c Log P values, and solvation energies, all crucial for molecule docking, which had become slow with older database organizations in previous versions of ZINC. As the libraries have continued to grow, we have been interested in finding whether molecular diversity has suffered, for instance, because certain scaffolds have come to dominate via easy analoging. This has not occurred thus far, and chemical diversity continues to grow with database size, with a log increase in Bemis-Murcko scaffolds for every two-log unit increase in database size. Most new scaffolds come from compounds with the highest heavy atom count. Finally, we consider the implications for databases like ZINC as the libraries grow toward and beyond the trillion-molecule range. ZINC is freely available to everyone and may be accessed at cartblanche22.docking.org, via Globus, and in the Amazon AWS and Oracle OCI clouds.

'Chemistry at the speed of sound': automated 1536-well nanoscale synthesis of 16 scaffolds in parallel

Screening of large and diverse libraries is the 'bread and butter' in the first phase of the discovery of novel drugs. However, maintenance and periodic renewal of high-quality large compound collections pose considerable logistic, environmental and monetary problems. Here, we exercise an alternative, the 'on-the-fly' synthesis of large and diverse libraries on a nanoscale in a highly automated fashion. For the first time, we show the feasibility of the synthesis of a large library based on 16 different chemistries in parallel on several 384-well plates using the acoustic dispensing ejection (ADE) technology platform. In contrast to combinatorial chemistry, we produced 16 scaffolds at the same time and in a sparse matrix fashion, and each compound was produced by a random combination of diverse large building blocks. The synthesis, analytics, resynthesis of selected compounds, and chemoinformatic analysis of the library are described. The advantages of the herein described automated nanoscale synthesis approach include great library diversity, absence of library storage logistics, superior economics, speed of synthesis by automation, increased safety, and hence sustainable chemistry.

Advancing Targeted Protein Degradation via Multiomics Profiling and Artificial Intelligence

Only around 20% of the human proteome is considered to be druggable with small-molecule antagonists. This leaves some of the most compelling therapeutic targets outside the reach of ligand discovery. The concept of targeted protein degradation (TPD) promises to overcome some of these limitations. In brief, TPD is dependent on small molecules that induce the proximity between a protein of interest (POI) and an E3 ubiquitin ligase, causing ubiquitination and degradation of the POI. In this perspective, we want to reflect on current challenges in the field, and discuss how advances in multiomics profiling, artificial intelligence, and machine learning (AI/ML) will be vital in overcoming them. The presented roadmap is discussed in the context of small-molecule degraders but is equally applicable for other emerging proximity-inducing modalities.

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.

Surely you are joking, Mr Docking!

In the wake of recent COVID-19 pandemics scientists around the world rushed to deliver numerous CADD (Computer-Aided Drug Discovery) methods and tools that could be reliably used to discover novel drug candidates against the SARS-CoV-2 virus. With that, there emerged a trend of a significant democratization of CADD that contributed to the rapid development of various COVID-19 drug candidates currently undergoing different stages of validation. On the other hand, this democratization also inadvertently led to the surge rapidly performed molecular docking studies to nominate multiple scores of novel drug candidates supported by computational arguments only. Albeit driven by best intentions, most of such studies also did not follow best practices in the field that require experience and expertise learned through multiple rigorously designed benchmarking studies and rigorous experimental validation. In this Viewpoint we reflect on recent disbalance between small number of rigorous and comprehensive studies and the proliferation of purely computational studies enabled by the ease of docking software availability. We further elaborate on the hyped oversale of CADD methods' ability to rapidly yield viable drug candidates and reiterate the critical importance of rigor and adherence to the best practices of CADD in view of recent emergence of AI and Big Data in the field.

Docking is not enough: 17-trifluoromethylphenyl trinor PGF2α is only a very weak ligand of neurokinin-1 receptor

17-trifluoromethylphenyl trinor prostaglandin F_2α (17-CF₃PTPGF_2α) was reported recently to exhibit in vitro and in vivo anticancer activity. Based solely on the results of in silico molecular docking, it was claimed that this compound is NK1 receptor (NK1R) antagonist and that its activity is through this receptor. In this contribution we show that 17-CF₃PTPGF_2α is only a very weak NK1R ligand (IC₅₀ > 200 μM). In connection with that we discuss the issue of this compound's molecular target. Finally, we briefly narrate on the proper use of molecular docking in biomedical research.

Computer-Aided Drug Design towards New Psychotropic and Neurological Drugs

Central nervous system (CNS) disorders are a therapeutic area in drug discovery where demand for new treatments greatly exceeds approved treatment options. This is complicated by the high failure rate in late-stage clinical trials, resulting in exorbitant costs associated with bringing new CNS drugs to market. Computer-aided drug design (CADD) techniques minimise the time and cost burdens associated with drug research and development by ensuring an advantageous starting point for pre-clinical and clinical assessments. The key elements of CADD are divided into ligand-based and structure-based methods. Ligand-based methods encompass techniques including pharmacophore modelling and quantitative structure activity relationships (QSARs), which use the relationship between biological activity and chemical structure to ascertain suitable lead molecules. In contrast, structure-based methods use information about the binding site architecture from an established protein structure to select suitable molecules for further investigation. In recent years, deep learning techniques have been applied in drug design and present an exciting addition to CADD workflows. Despite the difficulties associated with CNS drug discovery, advances towards new pharmaceutical treatments continue to be made, and CADD has supported these findings. This review explores various CADD techniques and discusses applications in CNS drug discovery from 2018 to November 2022.

Dopamine-induced arrestin recruitment and desensitization of the dopamine D4 receptor is regulated by G protein-coupled receptor kinase-2

The dopamine D₄ receptor (D₄R) is expressed in the retina, prefrontal cortex, and autonomic nervous system and has been implicated in attention deficit hyperactivity disorder (ADHD), substance use disorders, and erectile dysfunction. D₄R has also been investigated as a target for antipsychotics due to its high affinity for clozapine. As opposed to the closely related dopamine D₂ receptor (D₂R), dopamine-induced arrestin recruitment and desensitization at the D₄R have not been studied in detail. Indeed, some earlier investigations could not detect arrestin recruitment and desensitization of this receptor upon its activation by agonist. Here, we used a novel nanoluciferase complementation assay to study dopamine-induced recruitment of β-arrestin2 (βarr2; also known as arrestin3) and G protein-coupled receptor kinase-2 (GRK2) to the D₄R in HEK293T cells. We also studied desensitization of D₄R-evoked G protein-coupled inward rectifier potassium (GIRK; also known as Kir3) current responses in Xenopus oocytes. Furthermore, the effect of coexpression of GRK2 on βarr2 recruitment and GIRK response desensitization was examined. The results suggest that coexpression of GRK2 enhanced the potency of dopamine to induce βarr2 recruitment to the D₄R and accelerated the rate of desensitization of D₄R-evoked GIRK responses. The present study reveals new details about the regulation of arrestin recruitment to the D₄R and thus increases our understanding of the signaling and desensitization of this receptor.

Repurposing amine and carboxylic acid building blocks with an automatable esterification reaction

New methodologies to unite amines and carboxylic acids that complement the popular amide coupling can significantly expand accessible chemical space if they yield products distinct from the classic R-NHC(O)-R' amide arrangement. Here we have developed an amine-acid esterification reaction based on pyridinium salt activation of amine C-N bonds to create products of type R-OC(O)-R' upon reaction with alkyl and aryl carboxylic acids. The protocol is robust and facile as demonstrated by automation on open-source robotics.

Active Learning Guided Drug Design Lead Optimization Based on Relative Binding Free Energy Modeling

In silico identification of potent protein inhibitors commonly requires prediction of a ligand binding free energy (BFE). Thermodynamics integration (TI) based on molecular dynamics (MD) simulations is a BFE calculation method capable of acquiring accurate BFE, but it is computationally expensive and time-consuming. In this work, we have developed an efficient automated workflow for identifying compounds with the lowest BFE among thousands of congeneric ligands, which requires only hundreds of TI calculations. Automated machine learning (AutoML) orchestrated by active learning (AL) in an AL-AutoML workflow allows unbiased and efficient search for a small set of best-performing molecules. We have applied this workflow to select inhibitors of the SARS-CoV-2 papain-like protease and were able to find 133 compounds with improved binding affinity, including 16 compounds with better than 100-fold binding affinity improvement. We obtained a hit rate that outperforms that expected of traditional expert medicinal chemist-guided campaigns. Thus, we demonstrate that the combination of AL and AutoML with free energy simulations provides at least 20× speedup relative to the naïve brute force approaches.

ChemFlow─From 2D Chemical Libraries to Protein-Ligand Binding Free Energies

The accurate prediction of protein-ligand binding affinities is a fundamental problem for the rational design of new drug entities. Current computational approaches are either too expensive or inaccurate to be effectively used in virtual high-throughput screening campaigns. In addition, the most sophisticated methods, e.g., those based on configurational sampling by molecular dynamics, require significant pre- and postprocessing to provide a final ranking, which hinders straightforward applications by nonexpert users. We present a novel computational platform named ChemFlow to bridge the gap between 2D chemical libraries and estimated protein-ligand binding affinities. The software is designed to prepare a library of compounds provided in SMILES or SDF format, dock them into the protein binding site, and rescore the poses by simplified free energy calculations. Using a data set of 626 protein-ligand complexes and GPU computing, we demonstrate that ChemFlow provides relative binding free energies with an RMSE < 2 kcal/mol at a rate of 1000 ligands per day on a midsize computer cluster. The software is publicly available at https://github.com/IFMlab/ChemFlow.

Iterative computational design and crystallographic screening identifies potent inhibitors targeting the Nsp3 macrodomain of SARS-CoV-2

The nonstructural protein 3 (NSP3) of the severe acute respiratory syndrome-coronavirus-2 (SARS-CoV-2) contains a conserved macrodomain enzyme (Mac1) that is critical for pathogenesis and lethality. While small-molecule inhibitors of Mac1 have great therapeutic potential, at the outset of the COVID-19 pandemic, there were no well-validated inhibitors for this protein nor, indeed, the macrodomain enzyme family, making this target a pharmacological orphan. Here, we report the structure-based discovery and development of several different chemical scaffolds exhibiting low- to sub-micromolar affinity for Mac1 through iterations of computer-aided design, structural characterization by ultra-high-resolution protein crystallography, and binding evaluation. Potent scaffolds were designed with in silico fragment linkage and by ultra-large library docking of over 450 million molecules. Both techniques leverage the computational exploration of tangible chemical space and are applicable to other pharmacological orphans. Overall, 160 ligands in 119 different scaffolds were discovered, and 153 Mac1-ligand complex crystal structures were determined, typically to 1 Å resolution or better. Our analyses discovered selective and cell-permeable molecules, unexpected ligand-mediated conformational changes within the active site, and key inhibitor motifs that will template future drug development against Mac1.

Predicting the Likelihood of Molecules to Act as Modulators of Protein-Protein Interactions

Targeting protein-protein interactions (PPIs) by small molecule modulators (iPPIs) is an attractive strategy for drug therapy, and some iPPIs have already been introduced into the clinic. Blocking PPIs is however considered to be a more difficult task than inhibiting enzymes or antagonizing receptor activity. In this paper, we examine whether it is possible to predict the likelihood of molecules to act as iPPIs. Using our in-house iterative stochastic elimination (ISE) algorithm, we constructed two classification models that successfully distinguish between iPPIs from the iPPI-DB database and decoy molecules from either the Enamine HTS collection (ISE 1) or the ZINC database (ISE 2). External test sets of iPPIs taken from the TIMBAL database and decoys from Enamine HTS or ZINC were screened by the models: the area under the curve for the receiver operating characteristic curve was 0.85-0.89, and the Enrichment Factor increased from an initial 1 to as much as 66 for ISE 1 and 57 for ISE 2. Screening of the Enamine HTS and ZINC data sets through both models results in a library of ∼1.3 million molecules that pass either one of the models. This library is enriched with iPPI candidates that are structurally different from known iPPIs, and thus, it is useful for target-specific screenings and should accelerate the discovery of iPPI drug candidates. The entire library is available in Table S6.

Prioritizing Small Sets of Molecules for Synthesis through in-silico Tools: A Comparison of Common Ranking Methods

Prioritizing molecules for synthesis is a key role of computational methods within medicinal chemistry. Multiple tools exist for ranking molecules, from the cheap and popular molecular docking methods to more computationally expensive molecular-dynamics (MD)-based methods. It is often questioned whether the accuracy of the more rigorous methods justifies the higher computational cost and associated calculation time. Here, we compared the performance on ranking the binding of small molecules for seven scoring functions from five docking programs, one end-point method (MM/GBSA), and two MD-based free energy methods (PMX, FEP+). We investigated 16 pharmaceutically relevant targets with a total of 423 known binders. The performance of docking methods for ligand ranking was strongly system dependent. We observed that MD-based methods predominantly outperformed docking algorithms and MM/GBSA calculations. Based on our results, we recommend the application of MD-based free energy methods for prioritization of molecules for synthesis in lead optimization, whenever feasible.

Biological Mechanism-based Neurology and Psychiatry: A BACE1/2 and Downstream Pathway Model

In oncology, comprehensive omics and functional enrichment studies have led to an extensive profiling of (epi)genetic and neurobiological alterations that can be mapped onto a single tumor's clinical phenotype and divergent clinical phenotypes expressing common pathophysiological pathways. Consequently, molecular pathway-based therapeutic interventions for different cancer typologies, namely tumor type- and site-agnostic treatments, have been developed, encouraging the real-world implementation of a paradigm shift in medicine. Given the breakthrough nature of the new-generation translational research and drug development in oncology, there is an increasing rationale to transfertilize this blueprint to other medical fields, including psychiatry and neurology. In order to illustrate the emerging paradigm shift in neuroscience, we provide a state-of-the-art review of translational studies on the β-site amyloid precursor protein cleaving enzyme (BACE) and its most studied downstream effector, neuregulin, which are molecular orchestrators of distinct biological pathways involved in several neurological and psychiatric diseases. This body of data aligns with the evidence of a shared genetic/biological architecture among Alzheimer's disease, schizoaffective disorder, and autism spectrum disorders. To facilitate a forward-looking discussion about a potential first step towards the adoption of biological pathway-based, clinical symptom-agnostic, categorization models in clinical neurology and psychiatry for precision medicine solutions, we engage in a speculative intellectual exercise gravitating around BACE-related science, which is used as a paradigmatic case here. We draw a perspective whereby pathway-based therapeutic strategies could be catalyzed by highthroughput techniques embedded in systems-scaled biology, neuroscience, and pharmacology approaches that will help overcome the constraints of traditional descriptive clinical symptom and syndrome-focused constructs in neurology and psychiatry.

Discovery of a cryptic pocket in the AI-predicted structure of PPM1D phosphatase explains the binding site and potency of its allosteric inhibitors

Virtual screening is a widely used tool for drug discovery, but its predictive power can vary dramatically depending on how much structural data is available. In the best case, crystal structures of a ligand-bound protein can help find more potent ligands. However, virtual screens tend to be less predictive when only ligand-free crystal structures are available, and even less predictive if a homology model or other predicted structure must be used. Here, we explore the possibility that this situation can be improved by better accounting for protein dynamics, as simulations started from a single structure have a reasonable chance of sampling nearby structures that are more compatible with ligand binding. As a specific example, we consider the cancer drug target PPM1D/Wip1 phosphatase, a protein that lacks crystal structures. High-throughput screens have led to the discovery of several allosteric inhibitors of PPM1D, but their binding mode remains unknown. To enable further drug discovery efforts, we assessed the predictive power of an AlphaFold-predicted structure of PPM1D and a Markov state model (MSM) built from molecular dynamics simulations initiated from that structure. Our simulations reveal a cryptic pocket at the interface between two important structural elements, the flap and hinge regions. Using deep learning to predict the pose quality of each docked compound for the active site and cryptic pocket suggests that the inhibitors strongly prefer binding to the cryptic pocket, consistent with their allosteric effect. The predicted affinities for the dynamically uncovered cryptic pocket also recapitulate the relative potencies of the compounds (τ_b = 0.70) better than the predicted affinities for the static AlphaFold-predicted structure (τ_b = 0.42). Taken together, these results suggest that targeting the cryptic pocket is a good strategy for drugging PPM1D and, more generally, that conformations selected from simulation can improve virtual screening when limited structural data is available.

Machine Learning and Computational Chemistry for the Endocannabinoid System

Computational methods in medicinal chemistry facilitate drug discovery and design. In particular, machine learning methodologies have recently gained increasing attention. This chapter provides a structured overview of the current state of computational chemistry and its applications for the interrogation of the endocannabinoid system (ECS), highlighting methods in structure-based drug design, virtual screening, ligand-based quantitative structure-activity relationship (QSAR) modeling, and de novo molecular design. We emphasize emerging methods in machine learning and anticipate a forecast of future opportunities of computational medicinal chemistry for the ECS.

Overlap of On-demand Ultra-large Combinatorial Spaces with On-the-shelf Drug-like Libraries

On-demand combinatorial spaces are shifting paradigms in early drug discovery, by considerably increasing the searchable chemical space to several billions of compounds while securing their synthetic accessibility. We here systematically compared the on-the-shelf available drug-like chemical space (9 million compounds) to three on-demand ultra-large (ODUL) combinatorial fragment spaces (REAL, CHEMriya, GalaXi) covering 32 billion of readily accessible molecules. Surprisingly, only one space (REAL) intersects almost entirely the currently available drug-like space, suggesting that it is the only ODUL widely suitable for in-stock hit expansion. Of course, expanding a preliminary ODUL hit in the same chemical space is the best possible strategy to rapidly generate structure-activity relationships. All three spaces remain well suited to early hit finding initiatives since they all provide numerous unique scaffolds that are not described by on-the shelf collections.

Development of Optimal Virtual Screening Strategies to Identify Novel Toll-Like Receptor Ligands Using the DockBox Suite

Toll-like receptors (TLRs) represent attractive targets for developing modulators for the treatment of many pathologies, including inflammation, cancer, and autoimmune diseases. Here, we describe a protocol based on the DockBox package that enables to set up and perform structure-based virtual screening in order to increase the chance of identifying novel TLR ligands from chemical libraries.

Computer aided drug design in the development of proteolysis targeting chimeras

Proteolysis targeting chimeras represent a class of drug molecules with a number of attractive properties, most notably a potential to work for targets that, so far, have been in-accessible for conventional small molecule inhibitors. Due to their different mechanism of action, and physico-chemical properties, many of the methods that have been designed and applied for computer aided design of traditional small molecule drugs are not applicable for proteolysis targeting chimeras. Here we review recent developments in this field focusing on three aspects: de-novo linker-design, estimation of absorption for beyond-rule-of-5 compounds, and the generation and ranking of ternary complex structures. In spite of this field still being young, we find that a good number of models and algorithms are available, with the potential to assist the design of such compounds in-silico, and accelerate applied pharmaceutical research.

In silico and in vitro studies of potential inhibitors against Dengue viral protein NS5 Methyl Transferase from Ginseng and Notoginseng

Background and aim

Dengue is a potentially deadly tropical infectious disease transmitted by mosquito vector Aedes aegypti with no antiviral drug available to date. Dengue NS5 protein is crucial for viral replication and is the most conserved among all four Dengue serotypes, making it an attractive drug target. Both Ginseng and Notoginseng extracts and isolates have been shown to be effective against various viral infections yet against Dengue Virus is understudied. We aim to identify potential inhibitors against Dengue NS5 Methyl transferase from small molecular compounds found in Ginseng and Notoginseng.

Experimental procedure

A molecular docking model of Dengue NS5 Methyl transferase (MTase) domain was tested with decoys and then used to screen 91 small molecular compounds found in Ginseng and Notoginseng followed by Molecular dynamics simulations and the per-residue free energy decompositions based on molecular mechanics/Poisson-Boltzmann (generalised Born) surface area (MM/PB(GB)SA) calculations of the hit. ADME predictions and drug-likeness analyses were discussed to evaluate the viability of the hit as a drug candidate. To confirm our findings, in vitro studies of antiviral activities against RNA and a E protein synthesis and cell toxicity were carried out.

Results and conclusion

The virtual screening resulted in Isoquercitrin as a single hit. Further analyses of the Isoquercitrin-MTase complex show that Isoquercitrin can reside within both of the NS5 Methyl Transferase active sites; the AdoMet binding site and the RNA capping site. The Isoquercitrin is safe for consumption and accessible on multikilogram scale. In vitro studies showed that Isoquercitrin can inhibit Dengue virus by reducing viral RNA and viral protein synthesis with low toxicity to cells (CC50 > 20 μM). Our work provides evidence that Isoquercitrin can serve as an inhibitor of Dengue NS5 protein at the Methyl Transferase domain, further supporting its role as an anti-DENV agent.

A Comprehensive Survey of Prospective Structure-Based Virtual Screening for Early Drug Discovery in the Past Fifteen Years

Structure-based virtual screening (SBVS), also known as molecular docking, has been increasingly applied to discover small-molecule ligands based on the protein structures in the early stage of drug discovery. In this review, we comprehensively surveyed the prospective applications of molecular docking judged by solid experimental validations in the literature over the past fifteen years. Herein, we systematically analyzed the novelty of the targets and the docking hits, practical protocols of docking screening, and the following experimental validations. Among the 419 case studies we reviewed, most virtual screenings were carried out on widely studied targets, and only 22% were on less-explored new targets. Regarding docking software, GLIDE is the most popular one used in molecular docking, while the DOCK 3 series showed a strong capacity for large-scale virtual screening. Besides, the majority of identified hits are promising in structural novelty and one-quarter of the hits showed better potency than 1 μM, indicating that the primary advantage of SBVS is to discover new chemotypes rather than highly potent compounds. Furthermore, in most studies, only in vitro bioassays were carried out to validate the docking hits, which might limit the further characterization and development of the identified active compounds. Finally, several successful stories of SBVS with extensive experimental validations have been highlighted, which provide unique insights into future SBVS drug discovery campaigns.

The era of high-quality chemical probes

Small-molecule chemical probes are among the most important tools to study the function of proteins in cells and organisms. Regrettably, the use of weak and non-selective small molecules has generated an abundance of erroneous conclusions in the scientific literature. More recently, minimal criteria have been outlined for investigational compounds, encouraging the selection and use of high-quality chemical probes. Here, we briefly recall the milestones and key initiatives that have paved the way to this new era, illustrate examples of recent high-quality chemical probes and provide our perspective on future challenges and developments.

Challenges for chemistry in Ukraine after the war: Ukrainian science requires rebuilding and support

The unprovoked Russian invasion has created considerable challenges for Ukrainian science. In this article, we discuss actions needed to support and rebuild Ukrainian science and educational systems. The proposed actions take into account past Ukrainian scientific achievements including developments in organic chemistry.

Magnet for the Needle in Haystack: "Crystal Structure First" Fragment Hits Unlock Active Chemical Matter Using Targeted Exploration of Vast Chemical Spaces

Fragment-based drug discovery (FBDD) has successfully led to approved therapeutics for challenging and "undruggable" targets. In the context of FBDD, we introduce a novel, multidisciplinary method to identify active molecules from purchasable chemical space. Starting from four small-molecule fragment complexes of protein kinase A (PKA), a template-based docking screen using Enamine's multibillion REAL Space was performed. A total of 93 molecules out of 106 selected compounds were successfully synthesized. Forty compounds were active in at least one validation assay with the most active follow-up having a 13,500-fold gain in affinity. Crystal structures for six of the most promising binders were rapidly obtained, verifying the binding mode. The overall success rate for this novel fragment-to-hit approach was 40%, accomplished in only 9 weeks. The results challenge the established fragment prescreening paradigm since the standard industrial filters for fragment hit identification in a thermal shift assay would have missed the initial fragments.

Targeting protein conformations with small molecules to control protein complexes

Dynamic protein complexes function in all cellular processes, from signaling to transcription, using distinct conformations that regulate their activity. Conformational switching of proteins can turn on or off their activity through protein-protein interactions, catalytic function, cellular localization, or membrane interaction. Recent advances in structural, computational, and chemical methodologies have enabled the discovery of small-molecule activators and inhibitors of conformationally dynamic proteins by using a more rational design than a serendipitous screening approach. Here, we discuss such recent examples, focusing on the mechanism of protein conformational switching and its regulation by small molecules. We emphasize the rational approaches to control protein oligomerization with small molecules that offer exciting opportunities for investigation of novel biological mechanisms and drug discovery.

MILCDock: Machine Learning Enhanced Consensus Docking for Virtual Screening in Drug Discovery

Molecular docking tools are regularly used to computationally identify new molecules in virtual screening for drug discovery. However, docking tools suffer from inaccurate scoring functions with widely varying performance on different proteins. To enable more accurate ranking of active over inactive ligands in virtual screening, we created a machine learning consensus docking tool, MILCDock, that uses predictions from five traditional molecular docking tools to predict the probability a ligand binds to a protein. MILCDock was trained and tested on data from both the DUD-E and LIT-PCBA docking datasets and shows improved performance over traditional molecular docking tools and other consensus docking methods on the DUD-E dataset. LIT-PCBA targets proved to be difficult for all methods tested. We also find that DUD-E data, although biased, can be effective in training machine learning tools if care is taken to avoid DUD-E's biases during training.

Computer simulation of molecular recognition in biomolecular system: from in silico screening to generalized ensembles

Prediction of ligand-receptor complex structure is important in both the basic science and the industry such as drug discovery. We report various computation molecular docking methods: fundamental in silico (virtual) screening, ensemble docking, enhanced sampling (generalized ensemble) methods, and other methods to improve the accuracy of the complex structure. We explain not only the merits of these methods but also their limits of application and discuss some interaction terms which are not considered in the in silico methods. In silico screening and ensemble docking are useful when one focuses on obtaining the native complex structure (the most thermodynamically stable complex). Generalized ensemble method provides a free-energy landscape, which shows the distribution of the most stable complex structure and semi-stable ones in a conformational space. Also, barriers separating those stable structures are identified. A researcher should select one of the methods according to the research aim and depending on complexity of the molecular system to be studied.

Assessment of the Generalization Abilities of Machine-Learning Scoring Functions for Structure-Based Virtual Screening

In structure-based virtual screening (SBVS), it is critical that scoring functions capture protein-ligand atomic interactions. By focusing on the local domains of ligand binding pockets, a standardized pocket Pfam-based clustering (Pfam-cluster) approach was developed to assess the cross-target generalization ability of machine-learning scoring functions (MLSFs). Subsequently, 12 typical MLSFs were evaluated using random cross-validation (Random-CV), protein sequence similarity-based cross-validation (Seq-CV), and pocket Pfam-based cross-validation (Pfam-CV) methods. Surprisingly, all of the tested models showed decreased performances from Random-CV to Seq-CV to Pfam-CV experiments, not showing satisfactory generalization capacity. Our interpretable analysis suggested that the predictions on novel targets by MLSFs were dependent on buried solvent-accessible surface area (SASA)-related features of complex structures, with greater predicted binding affinities on complexes owning larger protein-ligand interfaces. By combining buried SASA-related features with target-specific patterns that were only shared among structurally similar compounds in the same cluster, the random forest (RF)-Score attained a good performance in the Random-CV test. Based on these findings, we strongly advise assessing the generalization ability of MLSFs with the Pfam-cluster approach and being cautious with the features learned by MLSFs.

Construction of a synthetic methodology-based library and its application in identifying a GIT/PIX protein-protein interaction inhibitor

In recent years, the flourishing of synthetic methodology studies has provided concise access to numerous molecules with new chemical space. These compounds form a large library with unique scaffolds, but their application in hit discovery is not systematically evaluated. In this work, we establish a synthetic methodology-based compound library (SMBL), integrated with compounds obtained from our synthetic researches, as well as their virtual derivatives in significantly larger scale. We screen the library and identify small-molecule inhibitors to interrupt the protein-protein interaction (PPI) of GIT1/β-Pix complex, an unrevealed target involved in gastric cancer metastasis. The inhibitor 14-5-18 with a spiro[bicyclo[2.2.1]heptane-2,3'-indolin]-2'-one scaffold, considerably retards gastric cancer metastasis in vitro and in vivo. Since the PPI targets are considered undruggable as they are hard to target, the successful application illustrates the structural specificity of SMBL, demonstrating its potential to be utilized as compound source for more challenging targets.

Systematic Investigation of Docking Failures in Large-Scale Structure-Based Virtual Screening

In recent years, large-scale structure-based virtual screening has attracted increasing levels of interest for identification of novel compounds corresponding to potential drug targets. It is critical to understand the strengths and weaknesses of docking algorithms to increase the success rate in practical applications. Here, we systematically investigated the docking successes and failures of two representative docking programs: UCSF DOCK 3.7 and AutoDock Vina. DOCK 3.7 performed better in early enrichment on the Directory of Useful Decoys: Enhanced (DUD-E) data set, although both docking methods were roughly comparable in overall enrichment performance. DOCK 3.7 also showed superior computational efficiency. Intriguingly, the Vina scoring function showed a bias toward compounds with higher molecular weights. Both the tested docking approaches yielded incorrectly predicted ligand binding poses caused by the limitations of torsion sampling. Based on a careful analysis of docking results from six representative cases, we propose the reasons underlying docking failures; furthermore, we provide a few solutions, representing practical guidance for large-scale virtual screening campaigns and future docking algorithm development.

The neural basis of psychedelic action

Psychedelics are serotonin 2A receptor agonists that can lead to profound changes in perception, cognition and mood. In this review, we focus on the basic neurobiology underlying the action of psychedelic drugs. We first discuss chemistry, highlighting the diversity of psychoactive molecules and the principles that govern their potency and pharmacokinetics. We describe the roles of serotonin receptors and their downstream molecular signaling pathways, emphasizing key elements for drug discovery. We consider the impact of psychedelics on neuronal spiking dynamics in several cortical and subcortical regions, along with transcriptional changes and sustained effects on structural plasticity. Finally, we summarize neuroimaging results that pinpoint effects on association cortices and thalamocortical functional connectivity, which inform current theories of psychedelic action. By synthesizing knowledge across the chemical, molecular, neuronal, and network levels, we hope to provide an integrative perspective on the neural mechanisms responsible for the acute and enduring effects of psychedelics on behavior.

Olfactory Receptors as an Emerging Chemical Sensing Scaffold

Chemical biosensors are an increasingly ubiquitous part of our lives. Beyond enzyme-coupled assays, recent synthetic biology advances now allow us to hijack more complex biosensing systems to respond to difficult to detect analytes, such as chemical small molecules. Here, we briefly overview recent advances in the biosensing of small molecules, including nucleic acid aptamers, allosteric transcription factors, and two-component systems. We then look more closely at a recently developed chemical sensing system, G protein-coupled receptor (GPCR)-based sensors. Finally, we consider the chemical sensing capabilities of the largest GPCR subfamily, olfactory receptors (ORs). We examine ORs' role in nature, their potential as a biomedical target, and their ability to detect compounds not amenable for detection using other biological scaffolds. We conclude by evaluating the current challenges, opportunities, and future applications of GPCR- and OR-based sensors.

One class classification for the detection of β2 adrenergic receptor agonists using single-ligand dynamic interaction data

G protein-coupled receptors are involved in many biological processes, relaying the extracellular signal inside the cell. Signaling is regulated by the interactions between receptors and their ligands, it can be stimulated by agonists, or inhibited by antagonists or inverse agonists. The development of a new drug targeting a member of this family requires to take into account the pharmacological profile of the designed ligands in order to elicit the desired response. The structure-based virtual screening of chemical libraries may prioritize a specific class of ligands by combining docking results and ligand binding information provided by crystallographic structures. The performance of the method depends on the relevance of the structural data, in particular the conformation of the targeted site, the binding mode of the reference ligand, and the approach used to compare the interactions formed by the docked ligand with those formed by the reference ligand in the crystallographic structure. Here, we propose a new method based on the conformational dynamics of a single protein–ligand reference complex to improve the biased selection of ligands with specific pharmacological properties in a structure-based virtual screening exercise. Interactions patterns between a reference agonist and the receptor, here exemplified on the β2 adrenergic receptor, were extracted from molecular dynamics simulations of the agonist/receptor complex and encoded in graphs used to train a one-class machine learning classifier. Different conditions were tested: low to high affinity agonists, varying simulation duration, considering or ignoring hydrophobic contacts, and tuning of the classifier parametrization. The best models applied to post-process raw data from retrospective virtual screening obtained by docking of test libraries effectively filtered out irrelevant poses, discarding inactive and non-agonist ligands while identifying agonists. Taken together, our results suggest that consistency of the binding mode during the simulation is a key to the success of the method.

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.

Optimizing interactions to protein binding sites by integrating docking-scoring strategies into generative AI methods

The identification and optimization of promising lead molecules is essential for drug discovery. Recently, artificial intelligence (AI) based generative methods provided complementary approaches for generating molecules under specific design constraints of relevance in drug design. The goal of our study is to incorporate protein 3D information directly into generative design by flexible docking plus an adapted protein-ligand scoring function, thereby moving towards automated structure-based design. First, the protein-ligand scoring function RFXscore integrating individual scoring terms, ligand descriptors, and combined terms was derived using the PDBbind database and internal data. Next, design results for different workflows are compared to solely ligand-based reward schemes. Our newly proposed, optimal workflow for structure-based generative design is shown to produce promising results, especially for those exploration scenarios, where diverse structures fitting to a protein binding site are requested. Best results are obtained using docking followed by RFXscore, while, depending on the exact application scenario, it was also found useful to combine this approach with other metrics that bias structure generation into "drug-like" chemical space, such as target-activity machine learning models, respectively.

Bespoke library docking for 5-HT2A receptor agonists with antidepressant activity

There is considerable interest in screening ultralarge chemical libraries for ligand discovery, both empirically and computationally1-4. Efforts have focused on readily synthesizable molecules, inevitably leaving many chemotypes unexplored. Here we investigate structure-based docking of a bespoke virtual library of tetrahydropyridines-a scaffold that is poorly sampled by a general billion-molecule virtual library but is well suited to many aminergic G-protein-coupled receptors. Using three inputs, each with diverse available derivatives, a one pot C-H alkenylation, electrocyclization and reduction provides the tetrahydropyridine core with up to six sites of derivatization5-7. Docking a virtual library of 75 million tetrahydropyridines against a model of the serotonin 5-HT2A receptor (5-HT2AR) led to the synthesis and testing of 17 initial molecules. Four of these molecules had low-micromolar activities against either the 5-HT2A or the 5-HT2B receptors. Structure-based optimization led to the 5-HT2AR agonists (R)-69 and (R)-70, with half-maximal effective concentration values of 41 nM and 110 nM, respectively, and unusual signalling kinetics that differ from psychedelic 5-HT2AR agonists. Cryo-electron microscopy structural analysis confirmed the predicted binding mode to 5-HT2AR. The favourable physical properties of these new agonists conferred high brain permeability, enabling mouse behavioural assays. Notably, neither had psychedelic activity, in contrast to classic 5-HT2AR agonists, whereas both had potent antidepressant activity in mouse models and had the same efficacy as antidepressants such as fluoxetine at as low as 1/40th of the dose. Prospects for using bespoke virtual libraries to sample pharmacologically relevant chemical space will be considered.

Structure-based discovery of nonopioid analgesics acting through the α2A-adrenergic receptor

Because nonopioid analgesics are much sought after, we computationally docked more than 301 million virtual molecules against a validated pain target, the α2A-adrenergic receptor (α2AAR), seeking new α2AAR agonists chemotypes that lack the sedation conferred by known α2AAR drugs, such as dexmedetomidine. We identified 17 ligands with potencies as low as 12 nanomolar, many with partial agonism and preferential Gi and Go signaling. Experimental structures of α2AAR complexed with two of these agonists confirmed the docking predictions and templated further optimization. Several compounds, including the initial docking hit '9087 [mean effective concentration (EC50) of 52 nanomolar] and two analogs, '7075 and PS75 (EC50 4.1 and 4.8 nanomolar), exerted on-target analgesic activity in multiple in vivo pain models without sedation. These newly discovered agonists are interesting as therapeutic leads that lack the liabilities of opioids and the sedation of dexmedetomidine.