KinFragLib: Exploring the Kinase Inhibitor Space Using Subpocket-Focused Fragmentation and Recombination

Design, synthesis and bio-evaluation of 2,5-disubstituted thiazole derivatives for potential treatment of acute myeloid leukemia through targeting CDK9

CDK9 plays a vital role in cellular transcriptional regulation. Hyper-activation of CDK9 leads to the occurrence of various cancers including acute myeloid leukemia, thereby rendering CDK9 an attractive target for cancer treatment. Based on hit compound A4 with 2,5-disubstituted thiazole core identified through the SyntaLinker-Hybrid scheme that shows weak inhibitory activity against both CDK9 and MOLM-13 cells, we designed and synthesized 32 derivatives through structural modification. In vitro anti-proliferative test screened and confirmed that 14 compounds showed highly inhibitory activity against MOLM-13 cells with IC50 values in the micromolar range. Among them, compound 24 displayed the best antiproliferative activity against MOLM-13 cells with an IC50 value of 0.034 μM, which was comparable to the positive drug (BAY1251152, IC50 = 0.031 μM). In vitro kinase inhibition assay results demonstrated that compound 24 had considerable inhibitory activity against CDK9 with an IC50 value of 5.5 nM and a weak inhibitory activity on other CDKs. Further cellular mechanism assays revealed that 24 affected CDK9 signaling pathways, induced cellular apoptosis and arrested cell cycle in the G2/M phase. Finally, further studies of compound 24 about molecular docking, molecular dynamics simulations and ADMET prediction were investigated. Collectively, compound 24 deserves further structural optimization and development for the treatment of acute myeloid leukemia.

Identification of pyrazole scaffold inhibitors targeting cyclin-dependent kinase 8 for potential use in pulmonary fibrosis

Idiopathic Pulmonary Fibrosis (IPF) is a disease that includes inflammation and scarring of the lung tissues. Cyclin-dependent kinase 8 (CDK8) is a target of interest due to its role in inflammatory pathways. CDK8 can also modulate the TGF-β/Smad signaling associated with IPF. Herein, a structure-based virtual screening (SBVS) campaign led to the identification of three CDK8 inhibitors. Testing of candidate inhibitors in protein and cellular assays confirmed CDK8 inhibition, with the most potent inhibitor producing an IC50 value of 398.8 nM. Computational analysis identified pharmacological interactions that lead to CDK8 inhibition. No significant cytotoxicity was observed when the inhibitor was treated in vitro. Further results showed that the inhibitor can disrupt proteins associated with the epithelial-mesenchymal transition (EMT) and reduce cell migration. Additionally, the inhibitor can disrupt the TGF- β1/Smad signaling axis in the nucleus, potentially impacting the transcription of IPF related protein expression, when treated in cells at 5 µM. Comparisons to structures of known CDK8 inhibitors showed the identified inhibitor to be structurally novel. When tested against a panel of kinases at 1 µM, the most potent inhibitor demonstrated a favorable CDK8 selectivity profile. The identification of the CDK8 inhibitors in this study can be used in future drug design studies and as CDK8 probes to explore alternative therapeutics for IPF.

Designing Macrocyclic Kinase Inhibitors Using Macrocycle Scaffold Hopping with Reinforced Learning (Macro-Hop)

Macrocycles have gained significant attention in drug discovery, with over 70 macrocyclic compounds currently in clinical use. Despite this progress, the effective methods for designing macrocycles remain elusive. In this study, we present Macro-Hop, a reinforced learning framework designed to rapidly and comprehensively explore the macrocycle chemical space. Macro-Hop efficiently generates novel macrocyclic scaffolds that not only align with predefined physicochemical properties but also exhibit 3D structural similarities to a specified reference compound. As a proof of concept, we applied Macro-Hop to design a new series of macrocycle inhibitors targeting PDGFRαD842 V kinase. The representative compound L7 exhibited high potency against PDGFRαD842 V in both biochemical and cellular assays with IC50 values of 23.8 and 2.1 nM, respectively. L7 effectively inhibited clinically relevant secondary mutants PDGFRαD842 V/G680R (IC50 = 64.1 nM) and PDGFRαD842 V/T674I (IC50 = 27.6 nM), highlighting the rapid effectiveness of wet-leb validation with Macro-Hop.

Different applications and differentiated libraries for crystallographic fragment screening

Macromolecular X-ray crystallography allows detection and characterisation of the binding of small, low-affinity chemical fragments. Here we review the utility of fragment screening for drug discovery, its potential for use in discovery science, as well as some of the distinct types of fragments that have been compiled into libraries.

Identification of a Potent CDK8 Inhibitor Using Structure-Based Virtual Screening

Pulmonary fibrosis is excessive scarring of the lung tissues. Transforming growth factor-beta (TGF-β) has been implicated in pulmonary fibrosis due to its ability to induce the epithelial-to-mesenchymal transition (EMT) and promote epithelial cell migration. Cyclin-dependent kinase 8 (CDK8) can mediate the TGF-β signaling pathways and could function as an alternative therapeutic target for treating pulmonary fibrosis. Here, we performed a structure-based virtual screening campaign to identify CDK8 inhibitors from a library of 1.6 million compounds. The screening process ended with the identification of a novel CDK8 inhibitor, P162-0948 (IC50: 50.4 nM). An interaction analysis highlighted important CDK8-ligand interactions that support its binding and inhibitory activity. Testing against a panel of 60 different kinases demonstrated P162-0948 selectivity toward CDK8. Crucially, the inhibitor was found to be structurally novel when compared to known CDK8 inhibitors. Testing in A549 human alveolar epithelial cell lines showed that the P162-0948 can reduce cell migration and protein expression of EMT-related proteins. When P162-0948 was treated in cells at 5 μM, phosphorylation of Smad in the nucleus was reduced, which suggests disruption of the TGF-β/Smad signaling pathway. The identification of P162-0948 shows that it is not only potent, but its structural novelty can inform future design studies for potential therapeutics targeting pulmonary fibrosis.

Subpocket Similarity-Based Hit Identification for Challenging Targets: Application to the WDR Domain of LRRK2

We herewith applied a priori a generic hit identification method (POEM) for difficult targets of known three-dimensional structure, relying on the simple knowledge of physicochemical and topological properties of a user-selected cavity. Searching for local similarity to a set of fragment-bound protein microenvironments of known structure, a point cloud registration algorithm is first applied to align known subpockets to the target cavity. The resulting alignment then permits us to directly pose the corresponding seed fragments in a target cavity space not typically amenable to classical docking approaches. Last, linking potentially connectable atoms by a deep generative linker enables full ligand enumeration. When applied to the WD40 repeat (WDR) central cavity of leucine-rich repeat kinase 2 (LRRK2), an unprecedented binding site, POEM was able to quickly propose 94 potential hits, five of which were subsequently confirmed to bind in vitro to LRRK2-WDR.

Guided Docking as a Data Generation Approach Facilitates Structure-Based Machine Learning on Kinases

Drug discovery pipelines nowadays rely on machine learning models to explore and evaluate large chemical spaces. While including 3D structural information is considered beneficial, structural models are hindered by the availability of protein-ligand complex structures. Exemplified for kinase drug discovery, we address this issue by generating kinase-ligand complex data using template docking for the kinase compound subset of available ChEMBL assay data. To evaluate the benefit of the created complex data, we use it to train a structure-based E(3)-invariant graph neural network. Our evaluation shows that binding affinities can be predicted with significantly higher precision by models that take synthetic binding poses into account compared to ligand- or drug-target interaction models alone.

Constructing Innovative Covalent and Noncovalent Compound Libraries: Insights from 3D Protein-Ligand Interactions

Noncovalent interactions between small-molecule drugs and protein targets assume a pivotal role in drug design. Moreover, the design of covalent inhibitors, forming covalent bonds with amino acid residues, requires rational reactivity for their covalent warheads, presenting a key challenge as well. Understanding the intricacies of these interactions provides a more comprehensive understanding of molecular binding mechanisms, thereby guiding the rational design of potent inhibitors. In this study, we adopted the fragment-based drug design approach, introducing a novel methodology to extract noncovalent and covalent fragments according to distinct three-dimensional (3D) interaction modes from noncovalent and covalent compound libraries. Additionally, we systematically replaced existing ligands with rational fragment substitutions, based on the spatial orientation of fragments in 3D space. Furthermore, we adopted a molecular generation approach to create innovative covalent inhibitors. This process resulted in the recombination of a noncovalent compound library and several covalent compound libraries, constructed by two commonly encountered covalent amino acids: cysteine and serine. We utilized noncovalent ligands in KLIFS and covalent ligands in CovBinderInPDB as examples to recombine noncovalent and covalent libraries. These recombined compound libraries cover a substantial portion of the chemical space present in the original compound libraries and exhibit superior performance in terms of molecular scaffold diversity compared to the original compound libraries and other 11 commercial libraries. We also recombined BTK-focused libraries, and 23 compounds within our libraries have been validated by former researchers to possess potential biological activity. The establishment of these compound libraries provides valuable resources for virtual screening of covalent and noncovalent drugs targeting similar molecular targets.

TWN-FS method: A novel fragment screening method for drug discovery

Fragment-based drug discovery (FBDD) is a well-established and effective method for generating diverse and novel hits in drug design. Kinases are suitable targets for FBDD due to their well-defined structure. Water molecules contribute to structure and function of proteins and also influence the environment within the binding pocket. Water molecules form a variety of hydrogen-bonded cyclic water-ring networks, collectively known as topological water networks (TWNs). Analyzing the TWNs in protein binding sites can provide valuable insights into potential locations and shapes for fragments within the binding site. Here, we introduce TWN-based fragment screening (TWN-FS) method, a novel screening method that suggests fragments through grouped TWN analysis within the protein binding site. We used this method to screen known CDK2, CHK1, IGF1R and ERBB4 inhibitors. Our findings suggest that TWN-FS method has the potential to effectively screen fragments. The TWN-FS method package is available on GitHub at https://github.com/pkj0421/TWN-FS.

A novel strategy for designing the magic shotguns for distantly related target pairs

Due to its promising capacity in improving drug efficacy, polypharmacology has emerged to be a new theme in the drug discovery of complex disease. In the process of novel multi-target drugs (MTDs) discovery, in silico strategies come to be quite essential for the advantage of high throughput and low cost. However, current researchers mostly aim at typical closely related target pairs. Because of the intricate pathogenesis networks of complex diseases, many distantly related targets are found to play crucial role in synergistic treatment. Therefore, an innovational method to develop drugs which could simultaneously target distantly related target pairs is of utmost importance. At the same time, reducing the false discovery rate in the design of MTDs remains to be the daunting technological difficulty. In this research, effective small molecule clustering in the positive dataset, together with a putative negative dataset generation strategy, was adopted in the process of model constructions. Through comprehensive assessment on 10 target pairs with hierarchical similarity-levels, the proposed strategy turned out to reduce the false discovery rate successfully. Constructed model types with much smaller numbers of inhibitor molecules gained considerable yields and showed better false-hit controllability than before. To further evaluate the generalization ability, an in-depth assessment of high-throughput virtual screening on ChEMBL database was conducted. As a result, this novel strategy could hierarchically improve the enrichment factors for each target pair (especially for those distantly related/unrelated target pairs), corresponding to target pair similarity-levels.

Target-Focused Library Design by Pocket-Applied Computer Vision and Fragment Deep Generative Linking

We here describe a computational approach (POEM: Pocket Oriented Elaboration of Molecules) to drive the generation of target-focused libraries while taking advantage of all publicly available structural information on protein-ligand complexes. A collection of 31 384 PDB-derived images with key shapes and pharmacophoric properties, describing fragment-bound microenvironments, is first aligned to the query target cavity by a computer vision method. The fragments of the most similar PDB subpockets are then directly positioned in the query cavity using the corresponding image transformation matrices. Lastly, suitable connectable atoms of oriented fragment pairs are linked by a deep generative model to yield fully connected molecules. POEM was applied to generate a library of 1.5 million potential cyclin-dependent kinase 8 inhibitors. By synthesizing and testing as few as 43 compounds, a few nanomolar inhibitors were quickly obtained with limited resources in just two iterative cycles.

KiSSim: Predicting Off-Targets from Structural Similarities in the Kinome

Protein kinases are among the most important drug targets because their dysregulation can cause cancer, inflammatory and degenerative diseases, and many more. Developing selective inhibitors is challenging due to the highly conserved binding sites across the roughly 500 human kinases. Thus, detecting subtle similarities on a structural level can help explain and predict off-targets among the kinase family. Here, we present the kinase-focused, subpocket-enhanced KiSSim fingerprint (Kinase Structural Similarity). The fingerprint builds on the KLIFS pocket definition, composed of 85 residues aligned across all available protein kinase structures, which enables residue-by-residue comparison without a computationally expensive alignment. The residues' physicochemical and spatial properties are encoded within their structural context including key subpockets at the hinge region, the DFG motif, and the front pocket. Since structure was found to contain information complementary to sequence, we used the fingerprint to calculate all-against-all similarities within the structurally covered kinome. We could identify off-targets that are unexpected if solely considering the sequence-based kinome tree grouping; for example, Erlobinib's known kinase off-targets SLK and LOK show high similarities to the key target EGFR (TK group), although belonging to the STE group. KiSSim reflects profiling data better or at least as well as other approaches such as KLIFS pocket sequence identity, KLIFS interaction fingerprints (IFPs), or SiteAlign. To rationalize observed (dis)similarities, the fingerprint values can be visualized in 3D by coloring structures with residue and feature resolution. We believe that the KiSSim fingerprint is a valuable addition to the kinase research toolbox to guide off-target and polypharmacology prediction. The method is distributed as an open-source Python package on GitHub and as a conda package: https://github.com/volkamerlab/kissim.

ChemInformatics Model Explorer (CIME): exploratory analysis of chemical model explanations

The introduction of machine learning to small molecule research- an inherently multidisciplinary field in which chemists and data scientists combine their expertise and collaborate - has been vital to making screening processes more efficient. In recent years, numerous models that predict pharmacokinetic properties or bioactivity have been published, and these are used on a daily basis by chemists to make decisions and prioritize ideas. The emerging field of explainable artificial intelligence is opening up new possibilities for understanding the reasoning that underlies a model. In small molecule research, this means relating contributions of substructures of compounds to their predicted properties, which in turn also allows the areas of the compounds that have the greatest influence on the outcome to be identified. However, there is no interactive visualization tool that facilitates such interdisciplinary collaborations towards interpretability of machine learning models for small molecules. To fill this gap, we present CIME (ChemInformatics Model Explorer), an interactive web-based system that allows users to inspect chemical data sets, visualize model explanations, compare interpretability techniques, and explore subgroups of compounds. The tool is model-agnostic and can be run on a server or a workstation.

Fragment-to-Lead Medicinal Chemistry Publications in 2020

Fragment-based drug discovery (FBDD) continues to evolve and make an impact in the pharmaceutical sciences. We summarize successful fragment-to-lead studies that were published in 2020. Having systematically analyzed annual scientific outputs since 2015, we discuss trends and best practices in terms of fragment libraries, target proteins, screening technologies, hit-optimization strategies, and the properties of hit fragments and the leads resulting from them. As well as the tabulated Fragment-to-Lead (F2L) programs, our 2020 literature review identifies several trends and innovations that promise to further increase the success of FBDD. These include developing structurally novel screening fragments, improving fragment-screening technologies, using new computer-aided design and virtual screening approaches, and combining FBDD with other innovative drug-discovery technologies.

Efficient Hit-to-Lead Searching of Kinase Inhibitor Chemical Space via Computational Fragment Merging

In early-stage drug discovery, the hit-to-lead optimization (or "hit expansion") stage entails starting from a newly identified active compound and improving its potency or other properties. Traditionally, this process relies on synthesizing and evaluating a series of analogues to build up structure-activity relationships. Here, we describe a computational strategy focused on kinase inhibitors, intended to expedite the process of identifying analogues with improved potency. Our protocol begins from an inhibitor of the target kinase and generalizes the synthetic route used to access it. By searching for commercially available replacements for the individual building blocks used to make the parent inhibitor, we compile an enumerated library of compounds that can be accessed using the same chemical transformations; these huge libraries can exceed many millions─or billions─of compounds. Because the resulting libraries are much too large for explicit virtual screening, we instead consider alternate approaches to identify the top-scoring compounds. We find that contributions from individual substituents are well described by a pairwise additivity approximation, provided that the corresponding fragments position their shared core in precisely the same way relative to the binding site. This key insight allows us to determine which fragments are suitable for merging into single new compounds and which are not. Further, the use of pairwise approximation allows interaction energies to be assigned to each compound in the library without the need for any further structure-based modeling: interaction energies instead can be reliably estimated from the energies of the component fragments, and the reduced computational requirements allow for flexible energy minimizations that allow the kinase to respond to each substitution. We demonstrate this protocol using libraries built from six representative kinase inhibitors drawn from the literature, which target five different kinases: CDK9, CHK1, CDK2, EGFR^T790M, and ACK1. In each example, the enumerated library includes additional analogues reported by the original study to have activity, and these analogues are successfully prioritized within the library. We envision that the insights from this work can facilitate the rapid assembly and screening of increasingly large libraries for focused hit-to-lead optimization. To enable adoption of these methods and to encourage further analyses, we disseminate the computational tools needed to deploy this protocol.

Kinase Inhibitor Scaffold Hopping with Deep Learning Approaches

The protein kinase family contains many promising drug targets. Many kinase inhibitors target the ATP-binding pocket, leading to approved drugs in past decades. Scaffold hopping is an effective approach for drug design. The kinase ATP-binding pocket is highly conserved, crossing the whole kinase family. This provides an opportunity to develop a scaffold hopping approach to explore diversified scaffolds among various kinase inhibitors. In this work, we report the SyntaLinker-Hybrid scheme for kinase inhibitor scaffold hopping. With this scheme, we replace molecular fragments bound at the conserved kinase hinge region with deep generative models. Thus, we are able to generate new kinase-inhibitor-like structures hybridizing the privileged fragments against the hinge region. We demonstrate that this scheme allows generation of kinase-inhibitor-like molecules with novel scaffolds, while retaining the binding features of existing kinase inhibitors. This work can be employed in lead identification against kinase targets.

KLIFS: an overhaul after the first 5 years of supporting kinase research

Kinases are a prime target of drug development efforts with >60 drug approvals in the past two decades. Due to the research into this protein family, a wealth of data has been accumulated that keeps on growing. KLIFS-Kinase-Ligand Interaction Fingerprints and Structures-is a structural database focusing on how kinase inhibitors interact with their targets. The aim of KLIFS is to support (structure-based) kinase research through the systematic collection, annotation, and processing of kinase structures. Now, 5 years after releasing the initial KLIFS website, the database has undergone a complete overhaul with a new website, new logo, and new functionalities. In this article, we start by looking back at how KLIFS has been used by the research community, followed by a description of the renewed KLIFS, and conclude with showcasing the functionalities of KLIFS. Major changes include the integration of approved drugs and inhibitors in clinical trials, extension of the coverage to atypical kinases, and a RESTful API for programmatic access. KLIFS is available at the new domain https://klifs.net.