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

Exploration of structure-activity relationships for the SARS-CoV-2 macrodomain from shape-based fragment linking and active learning

The macrodomain of severe acute respiratory syndrome coronavirus 2 nonstructural protein 3 is required for viral pathogenesis and is an emerging antiviral target. We previously performed an x-ray crystallography-based fragment screen and found submicromolar inhibitors by fragment linking. However, these compounds had poor membrane permeability and liabilities that complicated optimization. Here, we developed a shape-based virtual screening pipeline-FrankenROCS. We screened the Enamine high-throughput collection of 2.1 million compounds, selecting 39 compounds for testing, with the most potent binding with a 130 μM median inhibitory concentration (IC50). We then paired FrankenROCS with an active learning algorithm (Thompson sampling) to efficiently search the Enamine REAL database of 22 billion molecules, testing 32 compounds with the most potent binding with a 220 μM IC50. Further optimization led to analogs with IC50 values better than 10 μM. This lead series has improved membrane permeability and is poised for optimization. FrankenROCS is a scalable method for fragment linking to exploit synthesis-on-demand libraries.

BitBIRCH Clustering Refinement Strategies

Chemical libraries are becoming not only increasingly bigger, but they are doing so at an accelerated pace. Keeping up with this explosion in chemical data demands more than just hardware upgrades; we need dramatically more efficient algorithms as well. We have been working in this direction, with the introduction of the instant similarity (iSIM) framework, which uses n-ary similarity to speed up the processing of very large sets. Recently, we showed how to use this technique to cluster billions of molecules with unprecedented efficiency through the BitBIRCH algorithm. In this Application Note, we present a package fully dedicated to expanding on the BitBIRCH method, including multiple options that give the user appreciable control over the tree structure, while dramatically improving the quality of the final partitions. Remarkably, this is achieved without compromising the efficiency of the original method. We also present new postprocessing tools that help dissect the clustering information, as well as ample examples showcasing the new functionalities. BitBIRCH is publicly available at: https://github.com/mqcomplab/bitbirch.

On the Difficulty to Rescore Hits from Ultralarge Docking Screens

Docking-based virtual screening tools customized to mine ultralarge chemical spaces are consistently reported to yield both higher hit rates and more potent ligands than that achieved by conventional docking of smaller million-sized compound libraries. This remarkable achievement is however counterbalanced by the absolute necessity to design an efficient postprocessing of the millions of potential virtual hits for selecting a few chemically diverse compounds for synthesis and biological evaluation. We here retrospectively analyzed ten successful ultralarge virtual screening hit lists that underwent in vitro binding assays, for binding affinity prediction using eight rescoring methods including simple empirical scoring functions, machine learning, molecular-mechanics and quantum-mechanics approaches. Although the best predictions usually rely on the most sophisticated methods, none of the tested rescoring methods could robustly distinguish known binders from inactive compounds, across all assays. Energy refinement of protein-ligand complexes, prior to rescoring, marginally helped molecular mechanics and quantum mechanics approaches but deteriorates predictions from empirical and machine learning scoring functions.

Computational Hit Finding: An Industry Perspective

Computational hit finding, particularly virtual screening, has been a mainstay of drug discovery campaigns for decades, providing a cost-efficient complement to wet experiments. Innovation in this space slowed considerably as these approaches converged around mature software programs and stock chemical libraries up to ∼10 million in size. Recently, however, powered by massive increases in computational power, the emergence of ultralarge make-on-demand virtual libraries, the development of large capacity neural networks, the expansion of the domain of applicability of free energy calculations, and advances in protein structure prediction, the virtual screening field is currently seeing major change. We present a guide from industry practitioners summarizing key aspects on the changing computational hit finding landscape including practical recommendations for building a performant end-to-end screening workflow, strategies to mitigate risk by avoiding common pitfalls, determining success criteria, and a brief discussion of emerging technologies likely to impact drug discovery in the near future.

Synthon-Based Strategies Exploiting Molecular Similarity and Protein-Ligand Interactions for Efficient Screening of Ultra-Large Chemical Libraries

The rapid expansion of ultralarge chemical libraries has revolutionized drug discovery, providing access to billions of compounds. However, this growth poses relevant challenges for traditional virtual screening (VS) methods. To address these limitations, synthon-based approaches have emerged as scalable alternatives, exploiting combinatorial chemistry principles to prioritize building blocks over enumerated molecules. In this work, we present exaScreen and exaDock, two novel synthon-based methodologies designed for ligand-based and structure-based VS, respectively. In the former case, synthon selection is guided by the 3D hydrophobic/philic distribution pattern in conjunction with a specific synthon alignment protocol based on a quadrupolar expansion over the atoms that participate in the linking bonds between fragments. On the other hand, accommodation to the binding site under a geometrically restrained docking of synthon-based hybrid compounds is used in the selection of the optimal synthon combinations. These strategies exhibit comparable performance to the search performed using fully enumerated libraries in identifying active compounds with significantly lower computational cost, offering computationally efficient strategies for VS in ultralarge chemical spaces.

BitBIRCH: efficient clustering of large molecular libraries

The widespread use of Machine Learning (ML) techniques in chemical applications has come with the pressing need to analyze extremely large molecular libraries. In particular, clustering remains one of the most common tools to dissect the chemical space. Unfortunately, most current approaches present unfavorable time and memory scaling, which makes them unsuitable to handle million- and billion-sized sets. Here, we propose to bypass these problems with a time- and memory-efficient clustering algorithm, BitBIRCH. This method uses a tree structure similar to the one found in the Balanced Iterative Reducing and Clustering using Hierarchies (BIRCH) algorithm to ensure O(N) time scaling. BitBIRCH leverages the instant similarity (iSIM) formalism to process binary fingerprints, allowing the use of Tanimoto similarity, and reducing memory requirements. Our tests show that BitBIRCH is already >1000 times faster than standard implementations of the Taylor-Butina clustering for libraries with 1 500 000 molecules. BitBIRCH increases efficiency without compromising the quality of the resulting clusters. We explore strategies to handle large sets, which we applied in the clustering of one billion molecules under 5 hours using a parallel/iterative BitBIRCH approximation.

AI/ML methodologies and the future-will they be successful in designing the next generation of new chemical entities?

Cheminformatics and chemical databases are essential to drug discovery. However, machine learning (ML) and artificial intelligence (AI) methodologies are changing the way in which chemical data is used. How will the use of chemical data change in drug discovery moving forward? How do the new ML methods in molecular property prediction, hit and lead and target identification and structure prediction differ and compare with previous computational methods? Will new ML methodologies improve chemical diversity in ligand design, and offer computational enhancements. There are still many advantages to physics based methods and they offer something lacking in ML/ AI based methods. Additionally, ML training methods often give the best results when experimental assay measurements are fed back into the model. Often modeling and experimental methods are not diametrically opposed but offer the greatest advantage when used complementary.

Rapid traversal of vast chemical space using machine learning-guided docking screens

The accelerating growth of make-on-demand chemical libraries provides unprecedented opportunities to identify starting points for drug discovery with virtual screening. However, these multi-billion-scale libraries are challenging to screen, even for the fastest structure-based docking methods. Here we explore a strategy that combines machine learning and molecular docking to enable rapid virtual screening of databases containing billions of compounds. In our workflow, a classification algorithm is trained to identify top-scoring compounds based on molecular docking of 1 million compounds to the target protein. The conformal prediction framework is then used to make selections from the multi-billion-scale library, reducing the number of compounds to be scored by docking. The CatBoost classifier showed an optimal balance between speed and accuracy and was used to adapt the workflow for screens of ultralarge libraries. Application to a library of 3.5 billion compounds demonstrated that our protocol can reduce the computational cost of structure-based virtual screening by more than 1,000-fold. Experimental testing of predictions identified ligands of G protein-coupled receptors and demonstrated that our approach enables discovery of compounds with multi-target activity tailored for therapeutic effect.

BitBIRCH Clustering Refinement Strategies

Chemical libraries are becoming not only increasingly bigger, but they are doing so at an accelerated pace. Keeping up with this explosion in chemical data demands more than just hardware upgrades, we need dramatically more efficient algorithms as well. We have been working in this direction, with the introduction of the iSIM framework, which uses n-ary similarity to speed up the processing of very large sets. Recently, we showed how to use this technique to cluster billions of molecules with unprecedented efficiency through the BitBIRCH algorithm. In this Application Note we present a package fully-dedicated to expanding on the BitBIRCH method, including multiple options that give the user appreciable control over the tree structure, while dramatically improving the quality of the final partitions. Remarkably, this is achieved without compromising the efficiency of the original method. We also present new post-processing tools that help dissect the clustering information, as well as ample examples showcasing the new functionalities. BitBIRCH is publicly available at: https://github.com/mqcomplab/bitbirch.

The Need for Continuing Blinded Pose- and Activity Prediction Benchmarks

Computational tools for structure-based drug design (SBDD) are widely used in drug discovery and can provide valuable insights to advance projects in an efficient and cost-effective manner. However, despite the importance of SBDD to the field, the underlying methodologies and techniques have many limitations. In particular, binding pose and activity predictions (P-AP) are still not consistently reliable. We strongly believe that a limiting factor is the lack of a widely accepted and established community benchmarking process that independently assesses the performance and drives the development of methods, similar to the CASP benchmarking challenge for protein structure prediction. Here, we provide an overview of P-AP, unblinded benchmarking data sets, and blinded benchmarking initiatives (concluded and ongoing) and offer a perspective on learnings and the future of the field. To accelerate a breakthrough on the development of novel P-AP methods, it is necessary for the community to establish and support a long-term benchmark challenge that provides nonbiased training/test/validation sets, a systematic independent validation, and a forum for scientific discussions.

Discovery of Dual ROCK1/2 Inhibitors from Nocardiopsis sp. under Metal Stress

Rho-associated protein kinase (ROCK) inhibitors are promising therapeutic agents for reducing elevated intraocular pressure in patients with glaucoma. We explored new ROCK inhibitors derived from bioactive metabolites produced by microbes, specifically cryptic metabolites from Nocardiopsis sp. MCY7, using a liquid chromatography-mass spectrometry-based chemical analysis approach integrated with metal stress-driven isolation. This strategy led to the identification of two previously undescribed linear peptides, nocarnickelamides A and B (1 and 2), and an unreported cittilin derivative, cittilin C (3). The planar structures of 1-3 were elucidated using UV spectroscopy, high-resolution mass spectrometry, and nuclear magnetic resonance. The absolute configurations of 1 and 2 were assigned using the advanced Marfey's method. Biological assays demonstrated that nocarnickelamides (1 and 2) exhibited dual inhibitory activity against ROCK1 (IC50 29.8 and 14.9 μM, respectively) and ROCK2 (IC50 27.0 and 21.9 μM, respectively), with molecular simulations suggesting binding to the ATP-binding site. In human trabecular meshwork cells, 2 significantly inhibited the activation of ROCK-regulated cytoskeletal contraction markers such as the myosin light chain. Nocarnickelamide B (2) is a novel dual ROCK1/2 inhibitor and a potential pharmacophore for designing new therapeutic agents to reduce intraocular pressure in glaucoma.

Discovery of TRPV4-Targeting Small Molecules with Anti-Influenza Effects Through Machine Learning and Experimental Validation

Transient receptor potential vanilloid 4 (TRPV4) is a calcium-permeable cation channel critical for maintaining intracellular Ca2+ homeostasis and is essential in regulating immune responses, metabolic processes, and signal transduction. Recent studies have shown that TRPV4 activation enhances influenza A virus infection, promoting viral replication and transmission. However, there has been limited exploration of antiviral drugs targeting the TRPV4 channel. In this study, we developed the first machine learning model specifically designed to predict TRPV4 inhibitory small molecules, providing a novel approach for rapidly identifying repurposed drugs with potential antiviral effects. Our approach integrated machine learning, virtual screening, data analysis, and experimental validation to efficiently screen and evaluate candidate molecules. For high-throughput virtual screening, we employed computational methods to screen open-source molecular databases targeting the TRPV4 receptor protein. The virtual screening results were ranked based on predicted scores from our optimized model and binding energy, allowing us to prioritize potential inhibitors. Fifteen small-molecule drugs were selected for further in vitro and in vivo antiviral testing against influenza. Notably, glecaprevir and everolimus demonstrated significant inhibitory effects on the influenza virus, markedly improving survival rates in influenza-infected mice (protection rates of 80% and 100%, respectively). We also validated the mechanisms by which these drugs interact with the TRPV4 channel. In summary, our study presents the first predictive model for identifying TRPV4 inhibitors, underscoring TRPV4 inhibition as a promising strategy for antiviral drug development against influenza. This pioneering approach lays the groundwork for future clinical research targeting the TRPV4 channel in antiviral therapies.

Active Learning to Select the Most Suitable Reagents and One-Step Organic Chemistry Reactions for Prioritizing Target-Specific Hits from Ultralarge Chemical Spaces

Designing chemically novel and synthesizable ligands from the largest possible chemical space is a major issue in modern drug discovery to identify early hits that are easily amenable to medicinal chemistry optimization. Starting from the sole three-dimensional structure of a protein binding site, we herewith describe a fully automated active learning protocol to propose the commercial chemical reagents and one-step organic chemistry reactions necessary to enumerate target-specific primary hits from ultralarge chemical spaces. When applied in different scenarios (single transform and multiple transforms) addressing chemical spaces of various sizes (from 670 million to 4.5 billion compounds), the method was able to recover up to 98% of virtual hits discovered by an exhaustive docking-based approach while scanning only 5% of the full chemical space. It is therefore applicable to the structure-based screening of trillion-sized chemical spaces at a very high throughput with minimal computational resources.

Ultra-Large Virtual Screening: Definition, Recent Advances, and Challenges in Drug Design

Virtual screening (VS) in drug design employs computational methodologies to systematically rank molecules from a virtual compound library based on predicted features related to their biological activities or chemical properties. The recent expansion in commercially accessible compound libraries and the advancements in artificial intelligence (AI) and computational power - including enhanced central processing units (CPUs), graphics processing units (GPUs), high-performance computing (HPC), and cloud computing - have significantly expanded our capacity to screen libraries containing over 109 molecules. Herein, we review the concept of ultra-large virtual screening (ULVS), focusing on the various algorithms and methodologies employed for virtual screening at this scale. In this context, we present the software utilized, applications, and results of different approaches, such as brute force docking, reaction-based docking approaches, machine learning (ML) strategies applied to docking or other VS methods, and similarity/pharmacophore search-based techniques. These examples represent a paradigm shift in the drug discovery process, demonstrating not only the feasibility of billion-scale compound screening but also their potential to identify hit candidates and increase the structural diversity of novel compounds with biological activities.

Introducing SpaceGA: A Search Tool to Accelerate Large Virtual Screenings of Combinatorial Libraries

The growth of make-on-demand libraries in recent years has provided completely new possibilities for virtual screening for discovering new hit compounds with specific and favorable properties. However, since these libraries now contain billions of compounds, screening them using traditional methods such as molecular docking has become challenging and requires substantial computational resources. Thus, to take real advantage of the new possibilities introduced by the make-on-demand libraries, different methods have been proposed to accelerate the screening process and prioritize molecules for evaluation. Here, we introduce SpaceGA, a genetic algorithm that leverages the rapid similarity search tool SpaceLight (Bellmann, L.; Penner, P.; Rarey, M. Topological similarity search in large combinatorial fragment spaces. J. Chem. Inf. Model. 2021, 61, 238-251). to constrain the optimization process to accessible compounds within desired combinatorial libraries. As shown herein, SpaceGA is able to efficiently identify molecules with desired properties from trillions of synthesizable compounds by enumerating and evaluating only a small fraction of them. On this basis, SpaceGA represents a promising new tool for accelerating and simplifying virtual screens of ultralarge combinatorial databases.

Navigating Ultralarge Virtual Chemical Spaces with Product-of-Experts Chemical Language Models

Ultralarge virtual chemical spaces have emerged as a valuable resource for drug discovery, providing access to billions of make-on-demand compounds with high synthetic success rates. Chemical language models can potentially accelerate the exploration of these vast spaces through direct compound generation. However, existing models are not designed to navigate specific virtual chemical spaces and often overlook synthetic accessibility. To address this gap, we introduce product-of-experts (PoE) chemical language models, a modular and scalable approach to navigating ultralarge virtual chemical spaces. This method allows for controlled compound generation within a desired chemical space by combining a prior model pretrained on the target space with expert and anti-expert models fine-tuned using external property-specific data sets. We demonstrate that the PoE chemical language model can generate compounds with desirable properties, such as those that favorably dock to dopamine receptor D2 (DRD2) and are predicted to cross the blood-brain barrier (BBB), while ensuring that the majority of generated compounds are present within the target chemical space. Our results highlight the potential of chemical language models for navigating ultralarge virtual chemical spaces, and we anticipate that this study will motivate further research in this direction. The source code and data are freely available at https://github.com/shuyana/poeclm.

Recent advances from computer-aided drug design to artificial intelligence drug design

Computer-aided drug design (CADD), a cornerstone of modern drug discovery, can predict how a molecular structure relates to its activity and interacts with its target using structure-based and ligand-based methods. Fueled by ever-increasing data availability and continuous model optimization, artificial intelligence drug design (AIDD), as an enhanced iteration of CADD, has thrived in the past decade. AIDD demonstrates unprecedented opportunities in protein folding, property prediction, and molecular generation. It can also facilitate target identification, high-throughput screening (HTS), and synthetic route prediction. With AIDD involved, the process of drug discovery is greatly accelerated. Notably, AIDD offers the potential to explore uncharted territories of chemical space beyond current knowledge. In this perspective, we began by briefly outlining the main workflows and components of CADD. Then through showcasing exemplary cases driven by AIDD in recent years, we describe the evolving role of artificial intelligence (AI) in drug discovery from three distinct stages, that is, chemical library screening, linker generation, and de novo molecular generation. In this process, we attempted to draw comparisons between the features of CADD and AIDD.

Identifying Artifacts from Large Library Docking

While large library docking has discovered potent ligands for multiple targets, as the libraries have grown the hit lists can become dominated by rare artifacts that cheat our scoring functions. Here, we investigate rescoring top-ranked docked molecules with orthogonal methods to identify these artifacts, exploring implicit solvent models and absolute binding free energy perturbation as cross-filters. In retrospective studies, this approach deprioritized high-ranking nonbinders for nine targets while leaving true ligands relatively unaffected. We tested the method prospectively against hits from docking against AmpC β-lactamase. We prioritized 128 high-ranking molecules for synthesis and testing, a mixture of 39 molecules flagged as likely cheaters and 89 that were plausible inhibitors. None of the predicted cheating compounds inhibited AmpC detectably, while 57% of the 89 plausible compounds did so. As our libraries continue to grow, deprioritizing docking artifacts by rescoring with orthogonal methods may find wide use.

Recent advances in computational prediction of molecular properties in food chemistry

The combination of food chemistry and computational simulation has brought many impacts to food research, moving from experimental chemistry to computer chemistry. This paper will systematically review in detail the important role played by computational simulations in the development of the molecular structure of food, mainly from the atomic, molecular, and multicomponent dimension. It will also discuss how different computational chemistry models can be constructed and analyzed to obtain reliable conclusions. From the calculation principle to case analysis, this paper focuses on the selection and application of quantum mechanics, molecular mechanics and coarse-grained molecular dynamics in food chemistry research. Finally, experiments and computations of food chemistry are compared and summarized to obtain the best balance between them. The above review and outlook will provide an important reference for the intersection of food chemistry and computational chemistry, and is expected to provide innovative thinking for structural research in food chemistry.

Caught between a ROCK and a hard place: current challenges in structure-based drug design

The discipline of structure-based drug design (SBDD) is several decades old and it is tempting to think that the proliferation of experimental structures for many drug targets might make computer-aided drug design (CADD) straightforward. However, this is far from true. In this review, we illustrate some of the challenges that CADD scientists face every day in their work, even now. We use Rho-associated protein kinase (ROCK), and public domain structures and data, as an example to illustrate some of the challenges we have experienced during our project targeting this protein. We hope that this will help to prevent unrealistic expectations of what CADD can accomplish and to educate non-CADD scientists regarding the challenges still facing their CADD colleagues.

Structure-Based Drug Design with a Deep Hierarchical Generative Model

Recently, the remarkable growth of available crystal structure data and libraries of commercially available or readily synthesizable molecules have unlocked previously inaccessible regions of chemical space for drug development. Paired with improvements in virtual ligand screening methods, these expanded libraries are having a notable impact on early drug design efforts. Yet screening-based methods still face scalability limits, due to computational constraints and the sheer scale of drug-like space. Machine learning approaches are overcoming these limitations by learning the fundamental intra- and intermolecular relationships in drug-target systems from existing data. Here, we introduce DrugHIVE, a deep hierarchical variational autoencoder that outperforms state-of-the-art autoregressive and diffusion-based methods in both speed and performance on common generative benchmarks. DrugHIVE's hierarchical design enables improved control over molecular generation. Its capabilities include dramatically increasing virtual screening efficiency and accelerating a wide range of common drug design tasks, including de novo generation, molecular optimization, scaffold hopping, linker design, and high-throughput pattern replacement. Our highly scalable method can even be applied to receptors with high-confidence AlphaFold-predicted structures, extending the ability to generate high-quality drug-like molecules to a majority of the unsolved human proteome.

Extensive exploration of structure activity relationships for the SARS-CoV-2 macrodomain from shape-based fragment merging and active learning

The macrodomain contained in the SARS-CoV-2 non-structural protein 3 (NSP3) is required for viral pathogenesis and lethality. Inhibitors that block the macrodomain could be a new therapeutic strategy for viral suppression. We previously performed a large-scale X-ray crystallography-based fragment screen and discovered a sub-micromolar inhibitor by fragment linking. However, this carboxylic acid-containing lead had poor membrane permeability and other liabilities that made optimization difficult. Here, we developed a shape-based virtual screening pipeline - FrankenROCS - to identify new macrodomain inhibitors using fragment X-ray crystal structures. We used FrankenROCS to exhaustively screen the Enamine high-throughput screening (HTS) collection of 2.1 million compounds and selected 39 compounds for testing, with the most potent compound having an IC50 value equal to 130 μM. We then paired FrankenROCS with an active learning algorithm (Thompson sampling) to efficiently search the Enamine REAL database of 22 billion molecules, testing 32 compounds with the most potent having an IC50 equal to 220 μM. Further optimization led to analogs with IC50 values better than 10 μM, with X-ray crystal structures revealing diverse binding modes despite conserved chemical features. These analogs represent a new lead series with improved membrane permeability that is poised for optimization. In addition, the collection of 137 X-ray crystal structures with associated binding data will serve as a resource for the development of structure-based drug discovery methods. FrankenROCS may be a scalable method for fragment linking to exploit ever-growing synthesis-on-demand libraries.

Efficient clustering of large molecular libraries

The widespread use of Machine Learning (ML) techniques in chemical applications has come with the pressing need to analyze extremely large molecular libraries. In particular, clustering remains one of the most common tools to dissect the chemical space. Unfortunately, most current approaches present unfavorable time and memory scaling, which makes them unsuitable to handle million- and billion-sized sets. Here, we propose to bypass these problems with a time- and memory-efficient clustering algorithm, BitBIRCH. This method uses a tree structure similar to the one found in the Balanced Iterative Reducing and Clustering using Hierarchies (BIRCH) algorithm to ensure O N time scaling. BitBIRCH leverages the instant similarity (iSIM) formalism to process binary fingerprints, allowing the use of Tanimoto similarity, and reducing memory requirements. Our tests show that BitBIRCH is already > 1,000 times faster than standard implementations of the Taylor-Butina clustering for libraries with 1,500,000 molecules. BitBIRCH increases efficiency without compromising the quality of the resulting clusters. We explore strategies to handle large sets, which we applied in the clustering of one billion molecules under 5 hours using a parallel/iterative BitBIRCH approximation.

Structure-based virtual screening of vast chemical space as a starting point for drug discovery

Structure-based virtual screening aims to find molecules forming favorable interactions with a biological macromolecule using computational models of complexes. The recent surge of commercially available chemical space provides the opportunity to search for ligands of therapeutic targets among billions of compounds. This review offers a compact overview of structure-based virtual screens of vast chemical spaces, highlighting successful applications in early drug discovery for therapeutically important targets such as G protein-coupled receptors and viral enzymes. Emphasis is placed on strategies to explore ultra-large chemical libraries and synergies with emerging machine learning techniques. The current opportunities and future challenges of virtual screening are discussed, indicating that this approach will play an important role in the next-generation drug discovery pipeline.

Identifying Artifacts from Large Library Docking

While large library docking has discovered potent ligands for multiple targets, as the libraries have grown, the very top of the hit-lists can become populated with artifacts that cheat our scoring functions. Though these cheating molecules are rare, they become ever-more dominant with library growth. Here, we investigate rescoring top-ranked molecules from docking screens with orthogonal methods to identify these artifacts, exploring implicit solvent models and absolute binding free energy perturbation (AB-FEP) as cross-filters. In retrospective studies, this approach deprioritized high-ranking non-binders for nine targets while leaving true ligands relatively unaffected. We tested the method prospectively against results from large library docking AmpC β-lactamase. From the very top of the docking hit lists, we prioritized 128 molecules for synthesis and experimental testing, a mixture of 39 molecules that rescoring flagged as likely cheaters and another 89 that were plausible true actives. None of the 39 predicted cheating compounds inhibited AmpC up to 200 μ M in enzyme assays, while 57% of the 89 plausible true actives did do so, with 19 of them inhibiting the enzyme with apparentK i values better than 50 μ M . As our libraries continue to grow, a strategy of catching docking artifacts by rescoring with orthogonal methods may find wide use in the field.

Emerging structure-based computational methods to screen the exploding accessible chemical space

Structure-based virtual screening can be a valuable approach to computationally select hit candidates based on their predicted interaction with a protein of interest. The recent explosion in the size of chemical libraries increases the chances of hitting high-quality compounds during virtual screening exercises but also poses new challenges as the number of chemically accessible molecules grows faster than the computing power necessary to screen them. We review here two novel approaches rapidly gaining in popularity to address this problem: machine learning-accelerated and synthon-based library screening. We summarize the results from seminal proof-of-concept studies, highlight the latest developments, and discuss limitations and future directions.

Correlation of protein binding pocket properties with hits' chemistries used in generation of ultra-large virtual libraries

Although the size of virtual libraries of synthesizable compounds is growing rapidly, we are still enumerating only tiny fractions of the drug-like chemical universe. Our capability to mine these newly generated libraries also lags their growth. That is why fragment-based approaches that utilize on-demand virtual combinatorial libraries are gaining popularity in drug discovery. These à la carte libraries utilize synthetic blocks found to be effective binders in parts of target protein pockets and a variety of reliable chemistries to connect them. There is, however, no data on the potential impact of the chemistries used for making on-demand libraries on the hit rates during virtual screening. There are also no rules to guide in the selection of these synthetic methods for production of custom libraries. We have used the SAVI (Synthetically Accessible Virtual Inventory) library, constructed using 53 reliable reaction types (transforms), to evaluate the impact of these chemistries on docking hit rates for 40 well-characterized protein pockets. The data shows that the virtual hit rates differ significantly for different chemistries with cross coupling reactions such as Sonogashira, Suzuki-Miyaura, Hiyama and Liebeskind-Srogl coupling producing the highest hit rates. Virtual hit rates appear to depend not only on the property of the formed chemical bond but also on the diversity of available building blocks and the scope of the reaction. The data identifies reactions that deserve wider use through increasing the number of corresponding building blocks and suggests the reactions that are more effective for pockets with certain physical and hydrogen bond-forming properties.

A virtual library of small molecules mimicking dipeptides

Virtual combinatorial libraries are prevalent in drug discovery due to improvements in the prediction of synthetic reactions that can be performed. This has gone hand in hand with the development of virtual screening capabilities to effectively screen the large chemical spaces spanned by exhaustive enumeration of reaction products. In this study, we generated a small-molecule dipeptide mimic library to target proteins binding small peptides. The library was created based on the general idea of peptide synthesis, that is, amino acid mimics were reacted in silico to form the dipeptide mimics, yielding 2,036,819 unique compounds. After docking calculations, two compounds from the library were synthesized and tested against WD repeat-containing protein 5 (WDR5) and histamine receptors H1-H4 to evaluate whether these molecules are viable in assays. The compounds showed the highest potency at the histamine H3 receptor, with Ki values in the two-digit micromolar range.

Regression-Based Active Learning for Accessible Acceleration of Ultra-Large Library Docking

Structure-based drug discovery is a process for both hit finding and optimization that relies on a validated three-dimensional model of a target biomolecule, used to rationalize the structure-function relationship for this particular target. An ultralarge virtual screening approach has emerged recently for rapid discovery of high-affinity hit compounds, but it requires substantial computational resources. This study shows that active learning with simple linear regression models can accelerate virtual screening, retrieving up to 90% of the top-1% of the docking hit list after docking just 10% of the ligands. The results demonstrate that it is unnecessary to use complex models, such as deep learning approaches, to predict the imprecise results of ligand docking with a low sampling depth. Furthermore, we explore active learning meta-parameters and find that constant batch size models with a simple ensembling method provide the best ligand retrieval rate. Finally, our approach is validated on the ultralarge size virtual screening data set, retrieving 70% of the top-0.05% of ligands after screening only 2% of the library. Altogether, this work provides a computationally accessible approach for accelerated virtual screening that can serve as a blueprint for the future design of low-compute agents for exploration of the chemical space via large-scale accelerated docking. With recent breakthroughs in protein structure prediction, this method can significantly increase accessibility for the academic community and aid in the rapid discovery of high-affinity hit compounds for various targets.

On the hunt for metalloenzyme inhibitors: Investigating the presence of metal-coordinating compounds in screening libraries and chemical spaces

Metalloenzymes play vital roles in various biological processes, requiring the search for inhibitors to develop treatment options for diverse diseases. While compound library screening is a conventional approach, the exploration of virtual chemical spaces housing trillions of compounds has emerged as an alternative strategy. In this study, we investigated the suitability of selected screening libraries and chemical spaces for discovering inhibitors of metalloenzymes featuring common ions (Mg2+, Mn2+, and Zn2+). First, metal-coordinating groups from ligands interacting with ions in the Protein Data Bank were extracted. Subsequently, the prevalence of these groups in two focused screening libraries (Life Chemicals' chelator library, comprising 6,428 compounds, and Otava's chelator fragment library, with 1,784 fragments) as well as two chemical spaces (GalaXi and REAL space, containing billions of virtual products) was investigated. In total, 1,223 metal-coordinating groups were identified, with about a quarter of these groups found within the examined libraries and spaces. Our results indicate that these can serve as valuable starting points for drug discovery targeting metalloenzymes. In addition, this study suggests ways to improve libraries and spaces for better success in finding potential inhibitors for metalloenzymes.

Protein Structure-Based Organic Chemistry-Driven Ligand Design from Ultralarge Chemical Spaces

Ultralarge chemical spaces describing several billion compounds are revolutionizing hit identification in early drug discovery. Because of their size, such chemical spaces cannot be fully enumerated and require ad-hoc computational tools to navigate them and pick potentially interesting hits. We here propose a structure-based approach to ultralarge chemical space screening in which commercial chemical reagents are first docked to the target of interest and then directly connected according to organic chemistry and topological rules, to enumerate drug-like compounds under three-dimensional constraints of the target. When applied to bespoke chemical spaces of different sizes and chemical complexity targeting two receptors of pharmaceutical interest (estrogen β receptor, dopamine D3 receptor), the computational method was able to quickly enumerate hits that were either known ligands (or very close analogs) of targeted receptors as well as chemically novel candidates that could be experimentally confirmed by in vitro binding assays. The proposed approach is generic, can be applied to any docking algorithm, and requires few computational resources to prioritize easily synthesizable hits from billion-sized chemical spaces.

DrugHIVE: Target-specific spatial drug design and optimization with a hierarchical generative model

Rapid advancement in the computational methods of structure-based drug design has led to their widespread adoption as key tools in the early drug development process. Recently, the remarkable growth of available crystal structure data and libraries of commercially available or readily synthesizable molecules have unlocked previously inaccessible regions of chemical space for drug development. Paired with improvements in virtual ligand screening methods, these expanded libraries are having a significant impact on the success of early drug design efforts. However, screening-based methods are limited in their scalability due to computational limits and the sheer scale of drug-like space. An approach within the quickly evolving field of artificial intelligence (AI), deep generative modeling, is extending the reach of molecular design beyond classical methods by learning the fundamental intra- and inter-molecular relationships in drug-target systems from existing data. In this work we introduce DrugHIVE, a deep hierarchical structure-based generative model that enables fine-grained control over molecular generation. Our model outperforms state of the art autoregressive and diffusion-based methods on common benchmarks and in speed of generation. Here, we demonstrate DrugHIVEs capacity to accelerate a wide range of common drug design tasks such as de novo generation, molecular optimization, scaffold hopping, linker design, and high throughput pattern replacement. Our method is highly scalable and can be applied to high confidence AlphaFold predicted receptors, extending our ability to generate high quality drug-like molecules to a majority of the unsolved human proteome.

SpaceGrow: efficient shape-based virtual screening of billion-sized combinatorial fragment spaces

The growing size of make-on-demand chemical libraries is posing new challenges to cheminformatics. These ultra-large chemical libraries became too large for exhaustive enumeration. Using a combinatorial approach instead, the resource requirement scales approximately with the number of synthons instead of the number of molecules. This gives access to billions or trillions of compounds as so-called chemical spaces with moderate hardware and in a reasonable time frame. While extremely performant ligand-based 2D methods exist in this context, 3D methods still largely rely on exhaustive enumeration and therefore fail to apply. Here, we present SpaceGrow: a novel shape-based 3D approach for ligand-based virtual screening of billions of compounds within hours on a single CPU. Compared to a conventional superposition tool, SpaceGrow shows comparable pose reproduction capacity based on RMSD and superior ranking performance while being orders of magnitude faster. Result assessment of two differently sized subsets of the eXplore space reveals a higher probability of finding superior results in larger spaces highlighting the potential of searching in ultra-large spaces. Furthermore, the application of SpaceGrow in a drug discovery workflow was investigated in four examples involving G protein-coupled receptors (GPCRs) with the aim to identify compounds with similar binding capabilities and molecular novelty.

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.

Thompson Sampling─An Efficient Method for Searching Ultralarge Synthesis on Demand Databases

Over the last five years, virtual screening of ultralarge synthesis on-demand libraries has emerged as a powerful tool for hit identification in drug discovery programs. As these libraries have grown to tens of billions of molecules, we have reached a point where it is no longer cost-effective to screen every molecule virtually. To address these challenges, several groups have developed heuristic search methods to rapidly identify the best molecules on a virtual screen. This article describes the application of Thompson sampling (TS), an active learning approach that streamlines the virtual screening of large combinatorial libraries by performing a probabilistic search in the reagent space, thereby never requiring the full enumeration of the library. TS is a general technique that can be applied to various virtual screening modalities, including 2D and 3D similarity search, docking, and application of machine-learning models. In an illustrative example, we show that TS can identify more than half of the top 100 molecules from a docking-based virtual screen of 335 million molecules by evaluating 1% of the data set.

Considerations Around Structure-Based Drug Discovery for KRAS Using DOCK

Molecular docking is a popular computational tool in drug discovery. Leveraging structural information, docking software predicts binding poses of small molecules to cavities on the surfaces of proteins. Virtual screening for ligand discovery is a useful application of docking software. In this chapter, using the enigmatic KRAS protein as an example system, we endeavor to teach the reader about best practices for performing molecular docking with UCSF DOCK. We discuss methods for virtual screening and docking molecules on KRAS. We present the following six points to optimize our docking setup for prosecuting a virtual screen: protein structure choice, pocket selection, optimization of the scoring function, modification of sampling spheres and sampling procedures, choosing an appropriate portion of chemical space to dock, and the choice of which top scoring molecules to pick for purchase.

Molecular Dynamics as a Tool for Virtual Ligand Screening

Rational drug design is essential for new drugs to emerge, especially when the structure of a target protein or nucleic acid is known. To that purpose, high-throughput virtual ligand screening campaigns aim at discovering computationally new binding molecules or fragments to modulate particular biomolecular interactions or biological activities, related to a disease process. The structure-based virtual ligand screening process primarily relies on docking methods which allow predicting the binding of a molecule to a biological target structure with a correct conformation and the best possible affinity. The docking method itself is not sufficient as it suffers from several and crucial limitations (lack of full protein flexibility information, no solvation and ion effects, poor scoring functions, and unreliable molecular affinity estimation).At the interface of computer techniques and drug discovery, molecular dynamics (MD) allows introducing protein flexibility before or after a docking protocol, refining the structure of protein-drug complexes in the presence of water, ions, and even in membrane-like environments, describing more precisely the temporal evolution of the biological complex and ranking these complexes with more accurate binding energy calculations. In this chapter, we describe the up-to-date MD, which plays the role of supporting tools in the virtual ligand screening (VS) process.Without a doubt, using docking in combination with MD is an attractive approach in structure-based drug discovery protocols nowadays. It has proved its efficiency through many examples in the literature and is a powerful method to significantly reduce the amount of required wet experimentations (Tarcsay et al, J Chem Inf Model 53:2990-2999, 2013; Barakat et al, PLoS One 7:e51329, 2012; De Vivo et al, J Med Chem 59:4035-4061, 2016; Durrant, McCammon, BMC Biol 9:71-79, 2011; Galeazzi, Curr Comput Aided Drug Des 5:225-240, 2009; Hospital et al, Adv Appl Bioinforma Chem 8:37-47, 2015; Jiang et al, Molecules 20:12769-12786, 2015; Kundu et al, J Mol Graph Model 61:160-174, 2015; Mirza et al, J Mol Graph Model 66:99-107, 2016; Moroy et al, Future Med Chem 7:2317-2331, 2015; Naresh et al, J Mol Graph Model 61:272-280, 2015; Nichols et al, J Chem Inf Model 51:1439-1446, 2011; Nichols et al, Methods Mol Biol 819:93-103, 2012; Okimoto et al, PLoS Comput Biol 5:e1000528, 2009; Rodriguez-Bussey et al, Biopolymers 105:35-42, 2016; Sliwoski et al, Pharmacol Rev 66:334-395, 2014).

Deep Learning with Geometry-Enhanced Molecular Representation for Augmentation of Large-Scale Docking-Based Virtual Screening

Structure-based virtual screening has been a crucial tool in drug discovery for decades. However, as the chemical space expands, the existing structure-based virtual screening techniques based on molecular docking and scoring struggle to handle billion-entry ultralarge libraries due to the high computational cost. To address this challenge, people have resorted to machine learning techniques to enhance structure-based virtual screening for efficiently exploring the vast chemical space. In those cases, compounds are usually treated as sequential strings or two-dimensional topology graphs, limiting their ability to incorporate three-dimensional structural information for downstream tasks. We herein propose a novel deep learning protocol, GEM-Screen, which utilizes the geometry-enhanced molecular representation of the compounds docking to a specific target and is trained on docking scores of a small fraction of a library through an active learning strategy to approximate the docking outcome for yet nontraining entries. This protocol is applied to virtual screening campaigns against the AmpC and D4 targets, demonstrating that GEM-Screen enriches more than 90% of the hit scaffolds for AmpC in the top 4% of model predictions and more than 80% of the hit scaffolds for D4 in the same top-ranking size of library. GEM-Screen can be used in conjunction with traditional docking programs for docking of only the top-ranked compounds to avoid the exhaustive docking of the whole library, thus allowing for discovering top-scoring compounds from billion-entry libraries in a rapid yet accurate fashion.

"You've got the Body I've got the Brains" - Could the current AI-based tools replace the human ingenuity for designing new drug candidates?

The emergence of artificial intelligence (AI) tools has transformed the landscape of drug discovery, providing unprecedented speed, efficiency, and cost-effectiveness in the search for new therapeutics. From target identification to drug formulation and delivery, AI-driven algorithms have revolutionized various aspects of medicinal chemistry, significantly accelerating the drug design process. Despite the transformative power of AI, this perspective article emphasizes the limitations of AI tools in drug discovery, requiring inventive skills of medicinal chemists. However, the article highlighted that there is a need for a harmonious integration of AI-based tools and human expertise in drug discovery. Such a synergistic approach promises to lead to groundbreaking therapies that address unmet medical needs and benefit humankind. As the world evolves technologically, the question remains: When will AI tools effectively design and develop drugs? The answer may lie in the seamless collaboration between AI and human researchers, unlocking transformative therapies that combat diseases effectively.

Back in Person: Frontiers in Medicinal Chemistry 2023

The Frontiers in Medicinal Chemistry (FiMC) is the largest international Medicinal Chemistry conference in the German speaking area and took place from April 3rd to 5th 2023 in Vienna (Austria). Fortunately, after being cancelled in 2020 and two years (2021-2022) of entirely virtual meetings, due to the COVID-19 pandemic, the FiMC could be held in a face-to-face format again. Organized by the Division of Medicinal Chemistry of the German Chemical Society (GDCh), the Division of Pharmaceutical and Medicinal Chemistry of the German Pharmaceutical Society (DPhG), together with the Division of Medicinal Chemistry of the Austrian Chemical Society (GÖCH), the Austrian Pharmaceutical Society (ÖPhG), and a local organization committee from the University of Vienna headed by Thierry Langer, the meeting brought together 260 participants from 21 countries. The program included 38 lectures by leading scientists from industry and academia as well as early career investigators. Moreover, 102 posters were presented in two highly interactive poster sessions.

Machine Learning-Boosted Docking Enables the Efficient Structure-Based Virtual Screening of Giga-Scale Enumerated Chemical Libraries

The emergence of ultra-large screening libraries, filled to the brim with billions of readily available compounds, poses a growing challenge for docking-based virtual screening. Machine learning (ML)-boosted strategies like the tool HASTEN combine rapid ML prediction with the brute-force docking of small fractions of such libraries to increase screening throughput and take on giga-scale libraries. In our case study of an anti-bacterial chaperone and an anti-viral kinase, we first generated a brute-force docking baseline for 1.56 billion compounds in the Enamine REAL lead-like library with the fast Glide high-throughput virtual screening protocol. With HASTEN, we observed robust recall of 90% of the true 1000 top-scoring virtual hits in both targets when docking only 1% of the entire library. This reduction of the required docking experiments by 99% significantly shortens the screening time. In the kinase target, the employment of a hydrogen bonding constraint resulted in a major proportion of unsuccessful docking attempts and hampered ML predictions. We demonstrate the optimization potential in the treatment of failed compounds when performing ML-boosted screening and benchmark and showcase HASTEN as a fast and robust tool in a growing arsenal of approaches to unlock the chemical space covered by giga-scale screening libraries for everyday drug discovery campaigns.

Evaluating Scalable Supervised Learning for Synthesize-on-Demand Chemical Libraries

Traditional small-molecule drug discovery is a time-consuming and costly endeavor. High-throughput chemical screening can only assess a tiny fraction of drug-like chemical space. The strong predictive power of modern machine-learning methods for virtual chemical screening enables training models on known active and inactive compounds and extrapolating to much larger chemical libraries. However, there has been limited experimental validation of these methods in practical applications on large commercially available or synthesize-on-demand chemical libraries. Through a prospective evaluation with the bacterial protein-protein interaction PriA-SSB, we demonstrate that ligand-based virtual screening can identify many active compounds in large commercial libraries. We use cross-validation to compare different types of supervised learning models and select a random forest (RF) classifier as the best model for this target. When predicting the activity of more than 8 million compounds from Aldrich Market Select, the RF substantially outperforms a naïve baseline based on chemical structure similarity. 48% of the RF's 701 selected compounds are active. The RF model easily scales to score one billion compounds from the synthesize-on-demand Enamine REAL database. We tested 68 chemically diverse top predictions from Enamine REAL and observed 31 hits (46%), including one with an IC50 value of 1.3 μM.

Recent Developments in Ultralarge and Structure-Based Virtual Screening Approaches

Drug development is a wide scientific field that faces many challenges these days. Among them are extremely high development costs, long development times, and a small number of new drugs that are approved each year. New and innovative technologies are needed to solve these problems that make the drug discovery process of small molecules more time and cost efficient, and that allow previously undruggable receptor classes to be targeted, such as protein-protein interactions. Structure-based virtual screenings (SBVSs) have become a leading contender in this context. In this review, we give an introduction to the foundations of SBVSs and survey their progress in the past few years with a focus on ultralarge virtual screenings (ULVSs). We outline key principles of SBVSs, recent success stories, new screening techniques, available deep learning-based docking methods, and promising future research directions. ULVSs have an enormous potential for the development of new small-molecule drugs and are already starting to transform early-stage drug discovery.

Structure-Based Discovery of Receptor Activator of Nuclear Factor-κB Ligand (RANKL)-Induced Osteoclastogenesis Inhibitors

Receptor activator of nuclear factor-κB ligand (RANKL) has been actively pursued as a therapeutic target for osteoporosis, given that RANKL is the master mediator of bone resorption as it promotes osteoclast differentiation, activity and survival. We employed a structure-based virtual screening approach comprising two stages of experimental evaluation and identified 11 commercially available compounds that displayed dose-dependent inhibition of osteoclastogenesis. Their inhibitory effects were quantified through TRAP activity at the low micromolar range (IC₅₀ < 5 μΜ), but more importantly, 3 compounds displayed very low toxicity (LC₅₀ > 100 μΜ). We also assessed the potential of an N-(1-aryl-1H-indol-5-yl)aryl-sulfonamide scaffold that was based on the structure of a hit compound, through synthesis of 30 derivatives. Their evaluation revealed 4 additional hits that inhibited osteoclastogenesis at low micromolar concentrations; however, cellular toxicity concerns preclude their further development. Taken together with the structure-activity relationships provided by the hit compounds, our study revealed potent inhibitors of RANKL-induced osteoclastogenesis of high therapeutic index, which bear diverse scaffolds that can be employed in hit-to-lead optimization for the development of therapeutics against osteolytic diseases.

Relevance of the Trillion-Sized Chemical Space "eXplore" as a Source for Drug Discovery

Within the past two decades, virtual combinatorial compound collections, so-called chemical spaces, became an important molecule source for pharmaceutical research all over the world. The emergence of compound vendor chemical spaces with rapidly growing numbers of molecules raises questions about their application suitability and the quality of the content. Here, we examine the composition of the recently published and, so far, biggest chemical space, "eXplore", which comprises approximately 2.8 trillion virtual product molecules. The utility of eXplore to retrieve interesting chemistry around approved drugs and common Bemis Murcko scaffolds has been assessed with several methods (FTrees, SpaceLight, SpaceMACS). Further, the overlap between several vendor chemical spaces and a physicochemical property distribution analysis has been performed. Despite the straightforward chemical reactions underlying its setup, eXplore is demonstrated to provide relevant and, most importantly, easily accessible molecules for drug discovery campaigns.

Navigating large chemical spaces in early-phase drug discovery

The size of actionable chemical spaces is surging, owing to a variety of novel techniques, both computational and experimental. As a consequence, novel molecular matter is now at our fingertips that cannot and should not be neglected in early-phase drug discovery. Huge, combinatorial, make-on-demand chemical spaces with high probability of synthetic success rise exponentially in content, generative machine learning models go hand in hand with synthesis prediction, and DNA-encoded libraries offer new ways of hit structure discovery. These technologies enable to search for new chemical matter in a much broader and deeper manner with less effort and fewer financial resources. These transformational developments require new cheminformatics approaches to make huge chemical spaces searchable and analyzable with low resources, and with as little energy consumption as possible. Substantial progress has been made in the past years with respect to computation as well as organic synthesis. First examples of bioactive compounds resulting from the successful use of these novel technologies demonstrate their power to contribute to tomorrow's drug discovery programs. This article gives a compact overview of the state-of-the-art.

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.

Chemical Space Virtual Screening against Hard-to-Drug RNA Methyltransferases DNMT2 and NSUN6

Targeting RNA methyltransferases with small molecules as inhibitors or tool compounds is an emerging field of interest in epitranscriptomics and medicinal chemistry. For two challenging RNA methyltransferases that introduce the 5-methylcytosine (m⁵C) modification in different tRNAs, namely DNMT2 and NSUN6, an ultra-large commercially available chemical space was virtually screened by physicochemical property filtering, molecular docking, and clustering to identify new ligands for those enzymes. Novel chemotypes binding to DNMT2 and NSUN6 with affinities down to K_D,app = 37 µM and K_D,app = 12 µM, respectively, were identified using a microscale thermophoresis (MST) binding assay. These compounds represent the first molecules with a distinct structure from the cofactor SAM and have the potential to be developed into activity-based probes for these enzymes. Additionally, the challenges and strategies of chemical space docking screens with special emphasis on library focusing and diversification are discussed.