Knowledge based prediction of ligand binding modes and rational inhibitor design for kinase drug discovery

Discovery of a potent, selective, and orally available EGFR C797S mutant inhibitor (DS06652923) with in vivo antitumor activity

The C797S mutation is one of the major factors behind resistance to the third-generation epidermal growth factor receptor (EGFR) tyrosine kinase inhibitors (TKIs). Herein, we describe the discovery of DS06652923, a novel, potent, and orally available EGFR-triple-mutant inhibitor. Through scaffold hopping from the previously reported nicotinamide derivative, a novel biaryl scaffold was obtained. The potency was successfully enhanced by the introduction of basic substituents based on analysis of the docking study results. In addition, the difluoromethoxy group on the pyrazole ring improved the kinase selectivity by inducing steric clash with the other kinases. The most optimized compound, DS06652923, achieved tumor regression in the Ba/F3 allograft model upon its oral administration.

Optimization of Novel Pyrido-pyridazinone Derivatives as FER Tyrosine Kinase Inhibitors, Leading to the Potent DS08701581

Previously, we reported the new pyrido-pyridazinone template as a feline sarcoma-related (FER) tyrosine kinase inhibitor. Representative compound 1 (DS21360717) showed strong enzyme inhibitory activity (IC50 = 0.5 nM), however, its antitumor effect was insufficient, probably due to poor solubility and resultant low bioavailability (BA). In addition, the kinase selectivity was inadequate, which may result in certain safety risks. Here, we focused on derivatization of the unoptimized C-5 position to obtain promising FER inhibitors possessing strong antitumor effects and improved selectivity, referring to their X-ray crystal structure and the docking model with FES proto-oncogene tyrosine kinase as an FER surrogate. While establishing the synthetic route of the pyrido-pyridazinone scaffold, we obtained a desired compound via our derivatization. Our optimized compound 17c (DS08701581) showed the highest class cell-free and cell activities in this template, good oral BA, and improved kinase selectivity, resulting in significant tumor growth inhibition in the Ba/F3-FER tumor model without body weight loss.

Design, synthesis, and evaluation of dual EGFR/AURKB inhibitors as anticancer agents for non-small cell lung cancer

The epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKIs) are first-line agents for mutant EGFR-positive (mEGFR+) NSCLC. However, secondary resistant mutations develop following therapy that prevent EGFR-TKI binding. The EGFR-TKIs are rendered ineffective in NSCLC expressing EGFR resistant mutations (rmEGFR+). Mutations in Kirsten rat sarcoma virus protein (mKRAS) support persistent signaling downstream of EGFR regardless of EGFR-TKI earlier in the signaling cascade. The EGFR-TKIs are ineffective in mKRAS+ NSCLC. Thus, newer anticancer agents are needed for rmEGFR+ and mKRAS+ NSCLC. Aurora kinase B (AURKB) is a mitosis related kinase that is overexpressed in NSCLC and supports cancer cell proliferation and survival. Literature reports have suggested that AURKB inhibitors if given concurrently with an EGFR-TKI could overcome EGFR-TKI resistance in mKRAS+ NSCLC and rmEGFR + NSCLC, and showed improved anticancer effects compared to current single-targeted EGFR-TKIs. Molecular modeling was used to identify similarities between the kinase pockets of EGFR and AURKB. An overlap was observed for the inactive conformation of EGFR and the active conformation of AURKB. Compounds 3-7 were synthesized as dual EGFR/AURKB inhibitors for mKRAS+ and rmEGFR+ NSCLC. Compounds 5, 6 and 7 were identified as dual EGFR/AURKB inhibitors. Compound 5 demonstrated modest micromolar inhibition of rmEGFR+ NSCLC. All investigated compounds showed moderate inhibition of mKRAS+ NSCLC cells. Compound 7 demonstrated single-digit micromolar inhibition of mKRAS+ NSCLC.

Poor Generalization by Current Deep Learning Models for Predicting Binding Affinities of Kinase Inhibitors

The extreme surge of interest over the past decade surrounding the use of neural networks has inspired many groups to deploy them for predicting binding affinities of drug-like molecules to their receptors. A model that can accurately make such predictions has the potential to screen large chemical libraries and help streamline the drug discovery process. However, despite reports of models that accurately predict quantitative inhibition using protein kinase sequences and inhibitors' SMILES strings, it is still unclear whether these models can generalize to previously unseen data. Here, we build a Convolutional Neural Network (CNN) analogous to those previously reported and evaluate the model over four datasets commonly used for inhibitor/kinase predictions. We find that the model performs comparably to those previously reported, provided that the individual data points are randomly split between the training set and the test set. However, model performance is dramatically deteriorated when all data for a given inhibitor is placed together in the same training/testing fold, implying that information leakage underlies the models' performance. Through comparison to simple models in which the SMILES strings are tokenized, or in which test set predictions are simply copied from the closest training set data points, we demonstrate that there is essentially no generalization whatsoever in this model. In other words, the model has not learned anything about molecular interactions, and does not provide any benefit over much simpler and more transparent models. These observations strongly point to the need for richer structure-based encodings, to obtain useful prospective predictions of not-yet-synthesized candidate inhibitors.

Integrating a phenotypic screening with a structural simplification strategy to identify 4-phenoxy-quinoline derivatives to potently disrupt the mitotic localization of Aurora kinase B

We combined a mechanism-informed phenotypic screening (MIPS) assay with a structural simplification strategy to guide the discovery of compounds that disrupt the localization of the mitotic regulator, Aurora kinase B (AURKB), rather than inhibiting its catalytic activity. An initial hit 4-(4-methylthiophen-2-yl)-N-(4-(quinolin-4-yloxy)phenyl)phthalazin-1-amine was identified after screening an in-house library of small molecules and phenocopied the loss of function mutations in AURKB without inhibiting its catalytic activity. We isolated this hit compound activity to its 4-phenoxy-quinoline moiety. The fragment was further optimized into a class of new chemical entities that potently disrupt the mitotic localization of AURKB at low nanomolar concentrations and consequently elicit severe growth inhibition in diverse human cancer cell lines. A lead compound, N-(3-methoxy-5-(6-methoxyquinolin-4-yl)oxy)phenyl)acetamide possessed desirable pharmacokinetic properties such as AUC_0-∞: 227.15 [ng∙h/mL/(mg/kg)]; C_max: 3378.52 ng/mL T_1/2: 3.52 h; and F%: 42 % and produced the AURKB-inhibitory phenotypes in a mouse xenograft model. A lead compound is a powerful tool for interrogating the regulation of AURKB and has the potential to be further developed as a first-in-class oncology therapeutic.

SGC-CAMKK2-1: A Chemical Probe for CAMKK2

The serine/threonine protein kinase calcium/calmodulin-dependent protein kinase kinase 2 (CAMKK2) plays critical roles in a range of biological processes. Despite its importance, only a handful of inhibitors of CAMKK2 have been disclosed. Having a selective small molecule tool to interrogate this kinase will help demonstrate that CAMKK2 inhibition can be therapeutically beneficial. Herein, we disclose SGC-CAMKK2-1, a selective chemical probe that targets CAMKK2.

Evolutionary divergence in the conformational landscapes of tyrosine vs serine/threonine kinases

Inactive conformations of protein kinase catalytic domains where the DFG motif has a "DFG-out" orientation and the activation loop is folded present a druggable binding pocket that is targeted by FDA-approved 'type-II inhibitors' in the treatment of cancers. Tyrosine kinases (TKs) typically show strong binding affinity with a wide spectrum of type-II inhibitors while serine/threonine kinases (STKs) usually bind more weakly which we suggest here is due to differences in the folded to extended conformational equilibrium of the activation loop between TKs vs. STKs. To investigate this, we use sequence covariation analysis with a Potts Hamiltonian statistical energy model to guide absolute binding free-energy molecular dynamics simulations of 74 protein-ligand complexes. Using the calculated binding free energies together with experimental values, we estimated free-energy costs for the large-scale (~17-20 Å) conformational change of the activation loop by an indirect approach, circumventing the very challenging problem of simulating the conformational change directly. We also used the Potts statistical potential to thread large sequence ensembles over active and inactive kinase states. The structure-based and sequence-based analyses are consistent; together they suggest TKs evolved to have free-energy penalties for the classical 'folded activation loop' DFG-out conformation relative to the active conformation, that is, on average, 4-6 kcal/mol smaller than the corresponding values for STKs. Potts statistical energy analysis suggests a molecular basis for this observation, wherein the activation loops of TKs are more weakly 'anchored' against the catalytic loop motif in the active conformation and form more stable substrate-mimicking interactions in the inactive conformation. These results provide insights into the molecular basis for the divergent functional properties of TKs and STKs, and have pharmacological implications for the target selectivity of type-II inhibitors.

Knowledge-driven design and optimization of potent symmetric anticancer molecules: A case study on PKM2 activators

BACKGROUND

Pyruvate kinase M2 (PKM2) is preferentially expressed as a low-activity dimer over the active tetramer in proliferating tumor cells, resulting in metabolic reprogramming to achieve high energy requirements and nutrient uptake. This leads to a shift from the normal glycolytic pathway causing tumor cells to proliferate uncontrollably. This study utilizes knowledge-based drug discovery to determine the critical features from experimentally known PKM2 activators and design compounds that would significantly confer a stable structural and functional edge over the known compounds which are still at the preclinical stage.

METHODS

Conscientious molecular modeling studies were carried out and critical structural features were identified and validated from the knowledge of experimentally known PKM2 activators to confer high-binding affinities. A virtual library of 200 palindromic and non-palindromic activators was designed based on these identified critical features to target a distinct activator binding-site. This binding would favor specific dimer-dimer association and subsequent protein tetramerization. The resultant compounds strongly correlated with identified structural features and binding affinities which further strengthened our findings. The designed activators were then subjected to pharmacokinetic profiling and toxicity prediction, followed by free-binding energy calculations and MD simulations.

RESULTS

All the virtually designed activators comprising the identified critical features were observed to confer high-binding affinities ranging from -9.1 to -15.0 kcal/mol to the receptor protein. The designed activators also demonstrated optimum pharmacokinetic and toxicity profiles.

CONCLUSION

The best activators selected for MD simulations studies were conclusively observed to stabilize the required tetrameric conformation suggesting that these activators could potentially target PKM2 tetramerization that might restore the normal glycolytic pathway and suppress tumor progression.

Structural features of the protein kinase domain and targeted binding by small-molecule inhibitors

Protein kinases are key components in cellular signaling pathways as they carry out the phosphorylation of proteins, primarily on Ser, Thr, and Tyr residues. The catalytic activity of protein kinases is regulated, and they can be thought of as molecular switches that are controlled through protein-protein interactions and post-translational modifications. Protein kinases exhibit diverse structural mechanisms of regulation and have been fascinating subjects for structural biologists from the first crystal structure of a protein kinase over 30 years ago, to recent insights into kinase assemblies enabled by the breakthroughs in cryo-EM. Protein kinases are high-priority targets for drug discovery in oncology and other disease settings, and kinase inhibitors have transformed the outcomes of specific groups of patients. Most kinase inhibitors are ATP competitive, deriving potency by occupying the deep hydrophobic pocket at the heart of the kinase domain. Selectivity of inhibitors depends on exploiting differences between the amino acids that line the ATP site and exploring the surrounding pockets that are present in inactive states of the kinase. More recently, allosteric pockets outside the ATP site are being targeted to achieve high selectivity and to overcome resistance to current therapeutics. Here, we review the key regulatory features of the protein kinase family, describe the different types of kinase inhibitors, and highlight examples where the understanding of kinase regulatory mechanisms has gone hand in hand with the development of inhibitors.

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.

Hinge Binder Scaffold Hopping Identifies Potent Calcium/Calmodulin-Dependent Protein Kinase Kinase 2 (CAMKK2) Inhibitor Chemotypes

CAMKK2 is a serine/threonine kinase and an activator of AMPK whose dysregulation is linked with multiple diseases. Unfortunately, STO-609, the tool inhibitor commonly used to probe CAMKK2 signaling, has limitations. To identify promising scaffolds as starting points for the development of high-quality CAMKK2 chemical probes, we utilized a hinge-binding scaffold hopping strategy to design new CAMKK2 inhibitors. Starting from the potent but promiscuous disubstituted 7-azaindole GSK650934, a total of 32 compounds, composed of single-ring, 5,6-, and 6,6-fused heteroaromatic cores, were synthesized. The compound set was specifically designed to probe interactions with the kinase hinge-binding residues. Compared to GSK650394 and STO-609, 13 compounds displayed similar or better CAMKK2 inhibitory potency in vitro, while compounds 13g and 45 had improved selectivity for CAMKK2 across the kinome. Our systematic survey of hinge-binding chemotypes identified several potent and selective inhibitors of CAMKK2 to serve as starting points for medicinal chemistry programs.

Engineering Selectivity for Reduced Toxicity of Bacterial Kinase Inhibitors Using Structure-Guided Medicinal Chemistry

Tuberculosis is a major global public health concern, and new drugs are needed to combat both the typical form and the increasingly common drug-resistant form of this disease. The essential tuberculosis kinase PknB is an attractive drug development target because of its central importance in several critical signaling cascades. A major hurdle in kinase inhibitor development is the reduction of toxicity due to nonspecific kinase activity in host cells. Here a novel class of PknB inhibitors was developed from hit aminopyrimidine 1 (GW779439X), which was originally designed for human CDK4 but failed to progress clinically because of high toxicity and low specificity. Replacing the pyrazolopyridazine headgroup of the original hit with substituted pyridine or phenyl headgroups resulted in a reduction of Cdk activity and a 3-fold improvement in specificity over the human kinome while maintaining PknB activity. This also resulted in improved microbiological activity and reduced toxicity in THP-1 cells and zebrafish.

Molecular simulation studies to reveal the binding mechanisms of shikonin derivatives inhibiting VEGFR-2 kinase

Traditional vascular endothelial growth factor receptor 2 (VEGFR-2) inhibitors can manage angiogenesis; however, severe toxicity and resistance limit their long-term applications in clinical therapy. Shikonin (SHK) and its derivatives could be promising to inhibit the VEGFR-2 mediated angiogenesis, as they are reported to bind in the catalytic kinase domain with low affinity. However, the detailed molecular insights and binding dynamics of these natural inhibitors are unknown, which is crucial for potential SHK based lead design. Therefore, the present study employed molecular modeling and simulations techniques to get insight into the binding behaviors of SHK and its two derivates, β-hydroxyisovalerylshikonin (β-HIVS) and acetylshikonin (ACS). Here the intermolecular interactions between protein and ligands were studied by induced fit docking approach, which were further evaluated by treating QM/MM (quantum mechanics/molecular mechanics) and molecular dynamics (MD) simulation. The result showed that the naphthazarin ring of the SHK derivates is vital for strong binding to the catalytic domain; however, the binding stability can be modulated by the side chain modification. Because of having electrostatic potential, this ring makes essential interactions with the DFG (Asp1046 and Phe1047) motif and also allows interacting with the allosteric binding site. Taken together, the studies will advance our knowledge and scope for the development of new selective VEGFR-2 inhibitors based on SHK and its analogs.

Concise, gram-scale synthesis of furo[2,3-b]pyridines with functional handles for chemoselective cross-coupling

A concise 4-step synthesis of furo[2,3-b]pyridines, with handles in the 3- and 5-positions for palladium mediated cross-coupling reactions, is described. The synthetic route has been optimized, with only one step requiring purification by column chromatography. The route is amenable to scale-up, and was successfully executed on a multi-gram scale. Furopyridines are of growing interest in medicinal chemistry, and this route should enable easy access to the core for structure-activity relationship (SAR) studies.

MovableType Software for Fast Free Energy-Based Virtual Screening: Protocol Development, Deployment, Validation, and Assessment

For decades, the complicated energy surfaces found in macromolecular protein:ligand structures, which require large amounts of computational time and resources for energy state sampling, have been an inherent obstacle to fast, routine free energy estimation in industrial drug discovery efforts. Beginning in 2013, the Merz research group addressed this cost with the introduction of a novel sampling methodology termed “Movable Type” (MT). Using numerical integration methods, the MT method reduces the computational expense for energy state sampling by independently calculating each atomic partition function from an initial molecular conformation in order to estimate the molecular free energy using ensembles of the atomic partition functions. In this work, we report a software package, the DivCon Discovery Suite with the MovableType module from QuantumBio Inc., that performs this MT free energy estimation protocol in a fast, fully encapsulated manner. We discuss the computational procedures and improvements to the original work, and we detail the corresponding settings for this software package. Finally, we introduce two validation benchmarks to evaluate the overall robustness of the method against a broad range of protein:ligand structural cases. With these publicly available benchmarks, we show that the method can use a variety of input types and parameters and exhibits comparable predictability whether the method is presented with “expensive” X-ray structures or “inexpensively docked” theoretical models. We also explore some next steps for the method. The MovableType software is available at http://www.quantumbioinc.com/

Development of Potent PfCLK3 Inhibitors Based on TCMDC-135051 as a New Class of Antimalarials

The protein kinase PfCLK3 plays a critical role in the regulation of malarial parasite RNA splicing and is essential for the survival of blood stage Plasmodium falciparum. We recently validated PfCLK3 as a drug target in malaria that offers prophylactic, transmission blocking, and curative potential. Herein, we describe the synthesis of our initial hit TCMDC-135051 (1) and efforts to establish a structure-activity relationship with a 7-azaindole-based series. A total of 14 analogues were assessed in a time-resolved fluorescence energy transfer assay against the full-length recombinant protein kinase PfCLK3, and 11 analogues were further assessed in asexual 3D7 (chloroquine-sensitive) strains of P. falciparum parasites. SAR relating to rings A and B was established. These data together with analysis of activity against parasites collected from patients in the field suggest that TCMDC-135051 (1) is a promising lead compound for the development of new antimalarials with a novel mechanism of action targeting PfCLK3.

Discovery and development of ASK1 inhibitors

Aberrant activation of mitogen-activated protein kinases (MAPKs) like c-Jun N-terminal kinase (JNK) and p38 is an event involved in the pathophysiology of numerous human diseases. The apoptosis signal-regulating kinase 1 (ASK1) is an upstream target that gets activated only under pathological conditions and as such is a promising target for therapeutic intervention. In the first part of this review the molecular mechanisms leading to ASK1 activation and regulation will be described as well as the evidences supporting a pathogenic role for ASK1 in human disease. In the second part, an update on drug discovery efforts towards the discovery and development of ASK1-targeting therapies will be provided.

Unveiling the Structural Insights into the Selective Inhibition of Protein Kinase D1

BACKGROUND

Although protein kinase D1 (PKD1) has been proved to be an efficient target for anticancer drug development, lack of structural details and substrate binding mechanisms are the main obstacles for the development of selective inhibitors with therapeutic benefits.

OBJECTIVE

The present study described the in silico dynamics behaviors of PKD1 in binding with selective and non-selective inhibitors and revealed the critical binding site residues for the selective kinase inhibition.

METHODS

Here, the three dimensional model of PKD1 was initially constructed by homology modeling along with binding site characterization to explore the non-conserved residues. Subsequently, two known inhibitors were docked to the catalytic site and the detailed ligand binding mechanisms and post binding dyanmics were investigated by molecular dynamics simulation and binding free energy calculations.

RESULTS

According to the binding site analysis, PKD1 serves several non-conserved residues in the G-loop, hinge and catalytic subunits. Among them, the residues including Leu662, His663, and Asp665 from hinge region made polar interactions with selective PKD1 inhibitor in docking simulation, which were further validated by the molecular dynamics simulation. Both inhibitors strongly influenced the structural dynamics of PKD1 and their computed binding free energies were in accordance with experimental bioactivity data.

CONCLUSION

The identified non-conserved residues likely to play critical role on molecular reorganization and inhibitor selectivity. Taken together, this study explained the molecular basis of PKD1 specific inhibition, which may help to design new selective inhibitors for better therapies to overcome cancer and PKD1 dysregulated disorders.

In Depth Analysis of Kinase Cross Screening Data to Identify CAMKK2 Inhibitory Scaffolds

The calcium/calmodulin-dependent protein kinase kinase 2 (CAMKK2) activates CAMK1, CAMK4, AMPK, and AKT, leading to numerous physiological responses. The deregulation of CAMKK2 is linked to several diseases, suggesting the utility of CAMKK2 inhibitors for oncological, metabolic and inflammatory indications. In this work, we demonstrate that STO-609, frequently described as a selective inhibitor for CAMKK2, potently inhibits a significant number of other kinases. Through an analysis of literature and public databases, we have identified other potent CAMKK2 inhibitors and verified their activities in differential scanning fluorimetry and enzyme inhibition assays. These inhibitors are potential starting points for the development of selective CAMKK2 inhibitors and will lead to tools that delineate the roles of this kinase in disease biology.

Solvent-Free C-3 Coupling of Azaindoles with Cyclic Imines

By direct coupling 7-azaindole and cyclic imines, such as 3,4-dihydroisoquinoline, 6,7-dihydrothieno[3,2-c]pyridine, 3,4-dihydro-β-carboline, and 4,5-dihydro-3H-benz[c]azepine, new 3-substituted 7-azaindole derivatives have been synthesized. The reaction was extended to 4-azaindoles and 6-azaindoles, as electron-rich aromatic compounds. The lowest reactivity was observed in the case of C-3 substitution of 5-azaindole. In this case, the aza-Friedel-Crafts reaction took place by using 10 mol % of p-toluenesulfonic acid (p-TSA) as the catalyst. The role of the acid catalyst can be explained by the different pKa values of the azaindoles. All reactions were performed in solvent-free conditions by using both classical heating and microwave irradiation. In all cases, microwave heating proved to be more convenient to synthesize new C-3-substituted azaindole derivatives.

Structure-based identification of potent VEGFR-2 inhibitors from in vivo metabolites of a herbal ingredient

Vascular endothelial growth factor receptor-2 (VEGFR-2) is one of the regulatory elements of angiogenesis that is expressed highly in various diseases and is also essential for solid tumor growth. The present study was aimed at identifying potent inhibitors of VEGFR-2 by considering herbal secondary metabolites; as natural molecules are less toxic than synthetic derivatives. A structure-based virtual screening protocol consisting of molecular docking, MM-GBSA and ADME/T analysis was initially used to screen a library of in vivo metabolites of the herbal ingredient. Using a fixed cutoff value, four potent virtual hits were identified from molecular docking, ADME/T and binding affinity calculations, which were considered further for molecular dynamics (MD) simulation to broadly describe the binding mechanisms to VEGFR-2. The results suggested that these molecules have high affinity for the catalytic region of VEGFR-2, and form strong hydrophobic and polar interactions with the amino acids involved in the binding site of ATP and linker regions of the catalytic site. Subsequently, the stability of the docked complexes and binding mechanisms were evaluated by MD simulations, and the energy of binding was calculated through MM-PBSA analysis. The results uncovered two virtual hits, designated ZINC14762520 and ZINC36470466, as VEGFR-2 inhibitors, and suggested that they bind to kinase domain in an ATP-competitive manner. These virtual hits will offer a suitable starting point for the further design of their various analogs, allowing a rational search for more effective inhibitors in the future. Graphical abstract.

Inhibitory-κB Kinase (IKK) α and Nuclear Factor-κB (NFκB)-Inducing Kinase (NIK) as Anti-Cancer Drug Targets

The cellular kinases inhibitory-κB kinase (IKK) α and Nuclear Factor-κB (NF-κB)-inducing kinase (NIK) are well recognised as key central regulators and drivers of the non-canonical NF-κB cascade and as such dictate the initiation and development of defined transcriptional responses associated with the liberation of p52-RelB and p52-p52 NF-κB dimer complexes. Whilst these kinases and downstream NF-κB complexes transduce pro-inflammatory and growth stimulating signals that contribute to major cellular processes, they also play a key role in the pathogenesis of a number of inflammatory-based conditions and diverse cancer types, which for the latter may be a result of background mutational status. IKKα and NIK, therefore, represent attractive targets for pharmacological intervention. Here, specifically in the cancer setting, we reflect on the potential pathophysiological role(s) of each of these kinases, their associated downstream signalling outcomes and the stimulatory and mutational mechanisms leading to their increased activation. We also consider the downstream coordination of transcriptional events and phenotypic outcomes illustrative of key cancer 'Hallmarks' that are now increasingly perceived to be due to the coordinated recruitment of both NF-κB-dependent as well as NF-κB⁻independent signalling. Furthermore, as these kinases regulate the transition from hormone-dependent to hormone-independent growth in defined tumour subsets, potential tumour reactivation and major cytokine and chemokine species that may have significant bearing upon tumour-stromal communication and tumour microenvironment it reiterates their potential to be drug targets. Therefore, with the emergence of small molecule kinase inhibitors targeting each of these kinases, we consider medicinal chemistry efforts to date and those evolving that may contribute to the development of viable pharmacological intervention strategies to target a variety of tumour types.

Rational approach to highly potent and selective apoptosis signal-regulating kinase 1 (ASK1) inhibitors

Many diseases are believed to be driven by pathological levels of reactive oxygen species (ROS) and oxidative stress has long been recognized as a driver for inflammatory disorders. Apoptosis signal-regulating kinase 1 (ASK1) has been reported to be activated by intracellular ROS and its inhibition leads to a down regulation of p38-and JNK-dependent signaling. Consequently, ASK1 inhibitors may have the potential to treat clinically important inflammatory pathologies including renal, pulmonary and liver diseases. Analysis of the ASK1 ATP-binding site suggested that Gln756, an amino acid that rarely occurs at the GK+2 position, offered opportunities for achieving kinase selectivity for ASK1 which was applied to the design of a parallel medicinal chemistry library that afforded inhibitors of ASK1 with nanomolar potency and excellent kinome selectivity. A focused optimization strategy utilizing structure-based design resulted in the identification of ASK1 inhibitors with low nanomolar potency in a cellular assay, high selectivity when tested against kinase and broad pharmacology screening panels, and attractive physicochemical properties. The compounds we describe are attractive tool compounds to inform the therapeutic potential of ASK1 inhibition.

2D-QSAR study, molecular docking, and molecular dynamics simulation studies of interaction mechanism between inhibitors and transforming growth factor-beta receptor I (ALK5)

Transforming growth factor type 1 receptor (ALK5) is kinase associated with a wide variety of pathological processes, and inhibition of ALK5 is a good strategy to treat many kinds of cancer and fibrotic diseases. Recently, a series of compounds have been synthesized as ALK5 inhibitors. However, the study of their selectivity against other potential targets remains elusive. In this research, a data-set of ALK5 inhibitors were collected and studied based on the combination of 2D-QSAR, molecular docking and molecular dynamics simulation. The quality of QSAR models were assessed statistically by F, R², and R²_ADJ, proved to be credible. The cross-validations for the models (q²_LOO = 0.571 and 0.629, respectively) showed their robustness, while the external validations (r²_test = 0.703 and 0.764, respectively) showed their predictive power. Besides, the predicted binding free energy results calculated by MM/GBSA method were in accordance with the experimental data, and the van der Waals energy term was the factor that had the most significant impact on ligand binding. What is more, several important residues were found to significantly affect the binding affinity. Finally, based on our analyses above, a proposed series of molecules were designed.

In Silico Identification of a Novel Hinge-Binding Scaffold for Kinase Inhibitor Discovery

To explore novel kinase hinge-binding scaffolds, we carried out structure-based virtual screening against p38α MAPK as a model system. With the assistance of developed kinase-specific structural filters, we identify a novel lead compound that selectively inhibits a panel of kinases with threonine as the gatekeeper residue, including BTK and LCK. These kinases play important roles in lymphocyte activation, which encouraged us to design novel kinase inhibitors as drug candidates for ameliorating inflammatory diseases and cancers. Therefore, we chemically modified our substituted triazole-class lead compound to improve the binding affinity and selectivity via a "minimal decoration" strategy, which resulted in potent and selective kinase inhibitors against LCK (18 nM) and BTK (8 nM). Subsequent crystallographic experiments validated our design. These rationally designed compounds exhibit potent on-target inhibition against BTK in B cells or LCK in T cells, respectively. Our work demonstrates that structure-based virtual screening can be applied to facilitate the development of novel chemical entities in crowded chemical space in the field of kinase inhibitor discovery.

Inhibitory Kappa B Kinase α (IKKα) Inhibitors That Recapitulate Their Selectivity in Cells against Isoform-Related Biomarkers

IKKβ plays a central role in the canonical NF-kB pathway, which has been extensively characterized. The role of IKKα in the noncanonical NF-kB pathway, and indeed in the canonical pathway as a complex with IKKβ, is less well understood. One major reason for this is the absence of chemical tools designed as selective inhibitors for IKKα over IKKβ. Herein, we report for the first time a series of novel, potent, and selective inhibitors of IKKα. We demonstrate effective target engagement and selectivity with IKKα in U2OS cells through inhibition of IKKα-driven p100 phosphorylation in the noncanonical NF-kB pathway without affecting IKKβ-dependent IKappa-Bα loss in the canonical pathway. These compounds represent the first chemical tools that can be used to further characterize the role of IKKα in cellular signaling, to dissect this from IKKβ and to validate it in its own right as a target in inflammatory diseases.

Discovery of novel furanone derivatives as potent Cdc7 kinase inhibitors

Cdc7 is a serine-threonine kinase and plays a conserved and important role in DNA replication, and it has been recognized as a potential anticancer target. Herein, we report the design, synthesis and structure-activity relationship of novel furanone derivatives as Cdc7 kinase inhibitors. Compound 13 was identified as a strong inhibitor of Cdc7 with an IC₅₀ value of 0.6 nM in the presence of 1 mM ATP and showed excellent kinase selectivity. In addition, it exhibited slow off-rate characteristics, which may offer advantages over known Cdc7 inhibitors in its potential to yield prolonged inhibitory effects in vivo. Compound 13 potently inhibited Cdc7 activity in cancer cells, and effectively induced cell death.

2-Methyl-4/5-nitroimidazole derivatives potentiated against sexually transmitted Trichomonas: Design, synthesis, biology and 3D-QSAR study

Trichomoniasis is the most prevalent, non-viral sexually transmitted diseases (STD) caused by amitochondriate protozoan Trichomonas vaginalis. Increased resistance of T. vaginalis to the marketed drug Metronidazole necessitates the development of newer chemical entities. A library of sixty 2-methyl-4/5-nitroimidazole derivatives was synthesized via nucleophilic ring opening reaction of epoxide and the efficacies against drug-susceptible and -resistant Trichomonas vaginalis were evaluated. All the molecules except two were found to be active against both susceptible and resistant strains with MICs ranging 8.55-336.70 μM and 28.80-1445.08 μM, respectively. Most of the compounds were remarkably more effective than the standard Metronidazole. This study analyzes the in vitro and in vivo activities of the new 5-nitroimidazoles, which were found to be safe against human cervical HeLa cells with good selectivity index. The exploration of SAR by the synthesis of four different prototypes and 3D-QSAR study has shown the importance of prototype 1 over other prototypes.

Discovery and Rational Design of Pteridin-7(8H)-one-Based Inhibitors Targeting FMS-like Tyrosine Kinase 3 (FLT3) and Its Mutants

FLT3 has been validated as a therapeutic target for the treatment of acute myeloid leukemia (AML). In this paper, we describe for the first time, pteridin-7(8H)-one as a scaffold for potent FLT3 inhibitors derived from structural optimizations on irreversible EGFR inhibitors. The representative inhibitor (31) demonstrates single-digit nanomolar inhibition against FLT3 and subnanomolar KD for drug-resistance FLT3 mutants. In profiling of the in vitro tumor cell lines, it shows good selectivity against AML cells harboring FLT3-ITD mutations over other leukemia and solid tumor cell lines. The mechanism of action study illustrates that pteridin-7(8H)-one derivatives suppress the phosphorylation of FLT3 and its downstream pathways, thereby inducing G0/G1 cell cycle arrest and apoptosis in AML cells. In in vivo studies, 31 significantly suppresses the tumor growth in MV4-11 xenograft model. Overall, we provide a structurally distinct chemical scaffold with which to develop FLT3 mutants-selective inhibitors for AML treatment.

5-Hydroxy-7-azaindolin-2-one, a novel hybrid of pyridinol and sunitinib: design, synthesis and cytotoxicity against cancer cells

Synthesis of a series of hybrid compounds of pyridinol and sunitinib and their cytotoxicity against human cancer cell lines and improved safety windows are described.

Identification and optimisation of a series of N-(4-anilino-2-pyridyl)acetamide activin receptor-like kinase 1 (ALK1) inhibitors

A probe compound is described to inhibit ALK1, a target of interest for oncology due to its connection to angiogenesis.

Miltirone Induces G2/M Cell Cycle Arrest and Apoptosis in CCRF-CEM Acute Lymphoblastic Leukemia Cells

Miltirone (1) is a diterpene quinone extracted from a well-known Chinese traditional herb (Salvia miltiorrhiza). We investigated the cytotoxic effects of miltirone toward sensitive and multidrug-resistant acute lymphoblastic leukemia cell lines. Miltirone inhibited multidrug-resistant P-glycoprotein (P-gp)-overexpressing CEM/ADR5000 cells better than drug-sensitive CCRF-CEM wild-type cells, a phenomenon termed collateral sensitivity. Flow cytometric analyses revealed that miltirone induced G2/M arrest and apoptosis. Furthermore, miltirone stimulated reactive oxygen species (ROS) generation and mitochondrial membrane potential (MMP) disruption, which in turn induced DNA damage and activation of caspases and poly ADP-ribose polymerase (PARP). Downregulation of CCNB1 (cyclin B1) and CDC2 mRNA and upregulation of CDKN1A (p21) mRNA were in accord with miltirone-induced G2/M arrest. Moreover, miltirone decreased cell adherence to fibronectin. Molecular docking revealed that miltirone bound to the ATP-binding site of IKK-β. In conclusion, miltirone was collateral sensitive in multidrug-resistant P-gp-overexpressing cells, induced G2/M arrest, and triggered apoptosis via ROS-generated breakdown of MMP and DNA damage. Therefore, miltirone may be a promising candidate for cancer chemotherapy.

Development of small molecules targeting the pseudokinase Her3

Her3 is a member of the human epidermal growth factor receptor (EGFR) tyrosine kinase family, and it is often either overexpressed or deregulated in many types of human cancer. Her3 has not been the subject of small-molecule inhibitor development because it is a pseudokinase and does not possess appreciable kinase activity. We recently reported on the development of the first selective irreversible Her3 ligand (TX1-85-1) that forms a covalent bond with cysteine 721 which is unique to Her3 among all kinases. We also developed a bi-functional compound (TX2-121-1) containing a hydrophobic adamantane moiety and the same warhead of TX1-85-1 that is capable of inhibiting Her3-dependent signaling and growth. Here we report on the structure-based medicinal chemistry effort that resulted in the discovery of these two compounds.

The ABC of protein kinase conformations

Due to their involvement in human diseases, protein kinases are an important therapeutic target class. Conformation is a key concept for understanding how functional activity, inhibition and sequence are linked. We assemble and annotate the mammalian structural kinome from the Protein Data Bank on the basis of a universal residue nomenclature. We identify a torsion angle around the Gly of the DFG-motif whose sharp distribution profile corresponds to three eclipsed conformations. This allows the definition a small set of clusters whose distribution shows a bias for the active conformation. A common rationale links the active and inactive state: stabilization of the active conformation, as well as inactivation by displacement of helix-αC or the DFG-motif is governed by the interaction between helix-αC and the DFG motif. In particular, the conformation of the DFG-motif is tightly correlated with the propensity of helix-αC displacement. Our analysis reveals detailed mechanisms for the displacement of helix-αC and the DFG and improves our understanding of the role of individual residues. By pooling conformations from the whole structural kinome, the energetic contributions of sequence and extrinsic factors can be estimated in free energy analyses. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases.

Can structural features of kinase receptors provide clues on selectivity and inhibition? A molecular modeling study

Cancer is a complex disease resulting from the uncontrolled proliferation of cell signaling events. Protein kinases have been identified as central molecules that participate overwhelmingly in oncogenic events, thus becoming key targets for anticancer drugs. A majority of studies converged on the idea that ligand-binding pockets of kinases retain clues to the inhibiting abilities and cross-reacting tendencies of inhibitor drugs. Even though these ideas are critical for drug discovery, validating them using experiments is not only difficult, but also in some cases infeasible. To overcome these limitations and to test these ideas at the molecular level, we present here the results of receptor-focused in-silico docking of nine marketed drugs to 19 different wild-type and mutated kinases chosen from a wide range of families. This investigation highlights the need for using relevant models to explain the correct inhibition trends and the results are used to make predictions that might be able to influence future experiments. Our simulation studies are able to correctly predict the primary targets for each drug studied in majority of cases and our results agree with the existing findings. Our study shows that the conformations a given receptor acquires during kinase activation, and their micro-environment, defines the ligand partners. Type II drugs display high compatibility and selectivity for DFG-out kinase conformations. On the other hand Type I drugs are less selective and show binding preferences for both the open and closed forms of selected kinases. Using this receptor-focused approach, it is possible to capture the observed fold change in binding affinities between the wild-type and disease-centric mutations in ABL kinase for Imatinib and the second-generation ABL drugs. The effects of mutation are also investigated for two other systems, EGFR and B-Raf. Finally, by including pathway information in the design it is possible to model kinase inhibitors with potentially fewer side-effects.

Fragment-based design of kinase inhibitors: a practical guide

Fragment-based drug design has become an important strategy for drug design and development over the last decade. It has been used with particular success in the development of kinase inhibitors, which are one of the most widely explored classes of drug targets today. The application of fragment-based methods to discovering and optimizing kinase inhibitors can be a complicated and daunting task; however, a general process has emerged that has been highly fruitful. Here a practical outline of the fragment process used in kinase inhibitor design and development is laid out with specific examples. A guide to the overall process from initial discovery through fragment screening, including the difficulties in detection, to the computational methods available for use in optimization of the discovered fragments is reported.

Identification of new hit scaffolds by INPHARMA-guided virtual screening

INPHARMA, an NMR-based method that determines ligand binding poses, allows finding novel scaffolds as inhibitors of a target protein.

Specificity rendering 'hot-spots' for aurora kinase inhibitor design: the role of non-covalent interactions and conformational transitions

The present study examines the conformational transitions occurring among the major structural motifs of Aurora kinase (AK) concomitant with the DFG-flip and deciphers the role of non-covalent interactions in rendering specificity. Multiple sequence alignment, docking and structural analysis of a repertoire of 56 crystal structures of AK from Protein Data Bank (PDB) has been carried out. The crystal structures were systematically categorized based on the conformational disposition of the DFG-loop [in (DI) 42, out (DO) 5 and out-up (DOU) 9], G-loop [extended (GE) 53 and folded (GF) 3] and αC-helix [in (CI) 42 and out (CO) 14]. The overlapping subsets on categorization show the inter-dependency among structural motifs. Therefore, the four distinct possibilities a) 2W1C (DI, CI, GE) b) 3E5A (DI, CI, GF) c) 3DJ6 (DI, CO, GF) d) 3UNZ (DOU, CO, GF) along with their co-crystals and apo-forms were subjected to molecular dynamics simulations of 40 ns each to evaluate the variations of individual residues and their impact on forming interactions. The non-covalent interactions formed by the 157 AK co-crystals with different regions of the binding site were initially studied with the docked complexes and structure interaction fingerprints. The frequency of the most prominent interactions was gauged in the AK inhibitors from PDB and the four representative conformations during 40 ns. Based on this study, seven major non-covalent interactions and their complementary sites in AK capable of rendering specificity have been prioritized for the design of different classes of inhibitors.

c-Src binds to the cancer drug Ruxolitinib with an active conformation

The cancer drug Ruxolitinib is a potent janus kinase inhibitor approved for the treatment of the myeloproliferative neoplasms. In addition, Ruxolitinib has weak inhibitory activity against a panel of other kinases, including Src kinase. There is no structural information of Ruxolitinib binding to any kinase. In this paper, we determined the crystal structure of c-Src kinase domain in complex of Ruxolitinib at a resolution of 2.26 Å. C-Src kinase domain adopts the DFG-in active conformation upon Ruxolitinib binding, indicating Ruxolitinib is a type I inhibitor for c-Src. Ruxolitinib forms two hydrogen bonds with Met341, a water-mediated hydrogen bond with Thr338, and a number of van der Waals contacts with c-Src. Ruxolitinib was then docked into the ligand-binding pocket of a previously solved JAK1 structure. From the docking result, Ruxolitinib also binds JAK1 as a type I inhibitor, with more interactions and a higher shape complementarity with the ligand-binding pocket of JAK1 compared to that of c-Src. Since Ruxolitinib is a relatively small inhibitor and there is sizeable cavity between Ruxolitinib and c-Src ligand-binding pocket, we propose to modify Ruxolitinib to develop more potent inhibitors to c-Src.

Bioisosteric replacement of an acylureido moiety attached to an indolin-2-one scaffold with a malonamido or a 2/4-pyridinoylamido moiety produces a selectively potent Aurora-B inhibitor

Bioisosteric replacement of acylureido moiety in 6-acylureido-3-pyrrolylmethylidene-2-oxoindoline derivatives resulted in a series of malonamido derivatives with indolin-2-one scaffold (11-14). Further conformational restrictions of the malonamido moiety led to 2-oxo-1,2-dihydropyridine (21-25) or a 4-oxo-1,4-dihydropyridine derivatives (31-36). 4-Oxo-1,4-dihydropyridine derivatives were more potent Aurora B inhibitors than their 2-oxo-1,2-dihydropyridine counterparts and demonstrated cytotoxicities against A549 and HepG2 cells in the submicromolar range. In A549 cells, 31h decreased phosphorylation of histone H3, triggered polyploidy, induced expression of pro-apoptotic Fas and FasL with subsequent activation of caspase 8, resulting into apoptosis. In a Huh7-xenograft mouse model, 31h demonstrated potent in vivo efficacy with a daily dose of 5 mg/kg.