Structure-activity relationship study of 4-(thiazol-5-yl)benzoic acid derivatives as potent protein kinase CK2 inhibitors

Enhanced inhibitory activity of compounds containing purine scaffolds compared to protein kinase CK2α considering crystalline water

We recently reported novel purine-based CK2α inhibitors using the solvent ordering-based method as virtual screening. Among these, the X-ray crystal structure of a complex with CK2α was determined. The results showed that the crystalline water molecules observed in many previously reported complex structures of CK2α and its inhibitors had been eliminated. We then proposed a structure-based drug design. Since the removal of water molecules would be detrimental to inhibitor binding, new groups of compounds were designed by changing the position of the carboxy group located at the point where a water molecule would be present so as not to eliminate it. Compounds with (E)-2-carboxyethenyl and 3-carboxyphenyl substituted at the 2-position on the purine scaffold showed much higher inhibitory potency than 4-carboxyphenyl derivatives. Furthermore, in the presence of a 4-fluorophenyl group at the 9-position on the purine scaffold, the inhibitory activity of the 3-carboxyphenyl derivative against CK2α was 0.18 μM, a 167-fold improvement compared to the 4-carboxyphenyl derivative. The strategy of leaving crystalline water can significantly increase inhibitory activity.

Structure-guided discovery of adenosine triphosphate-competitive casein kinase 2 inhibitors

Casein kinase 2 (CK2) is a ubiquitous, highly pleiotropic serine-threonine kinase. CK2 has been identified as a potential drug target for the treatment of cancer and related disorders. Several adenosine triphosphate-competitive CK2 inhibitors have been identified and have progressed at different levels of clinical trials. This review presents details of CK2 protein, structural insights into adenosine triphosphate binding pocket, current clinical trial candidates and their analogues. Further, it includes the emerging structure-based drug design approaches, chemistry, structure-activity relationship and biological screening of potent and selective CK2 inhibitors. The authors tabulated the details of CK2 co-crystal structures because these co-crystal structures facilitated the structure-guided discovery of CK2 inhibitors. The narrow hinge pocket compared with related kinases provides useful insights into the discovery of CK2 inhibitors.

The pyridazine heterocycle in molecular recognition and drug discovery

The pyridazine ring is endowed with unique physicochemical properties, characterized by weak basicity, a high dipole moment that subtends π-π stacking interactions and robust, dual hydrogen-bonding capacity that can be of importance in drug-target interactions. These properties contribute to unique applications in molecular recognition while the inherent polarity, low cytochrome P450 inhibitory effects and potential to reduce interaction of a molecule with the cardiac hERG potassium channel add additional value in drug discovery and development. The recent approvals of the gonadotropin-releasing hormone receptor antagonist relugolix (24) and the allosteric tyrosine kinase 2 inhibitor deucravacitinib (25) represent the first examples of FDA-approved drugs that incorporate a pyridazine ring. In this review, the properties of the pyridazine ring are summarized in comparison to the other azines and its potential in drug discovery is illustrated through vignettes that explore applications that take advantage of the inherent physicochemical properties as an approach to solving challenges associated with candidate optimization.

Graphical Abstract

System truncation accelerates binding affinity calculations with the fragment molecular orbital method: A benchmark study

The fragment molecular orbital (FMO) method is a fast quantum-mechanical method that divides systems into pieces of fragments and performs ab initio calculations. The system truncation enables further speed improvement. In this article, we systematically study the effects of system truncations on binding affinity calculations obtained with FMO in combination with either the polarizable continuum model (FMO/PCM) or in combination with the Møller-Plesset method (FMO-MP2). We have used five protein complexes with ligands of several charged states. The calculated binding energies of the size variants of the truncated system, including only a restricted number of atoms around the ligand, are compared to the energy obtained from a full system. The result shows that the systems could be truncated to a radius of 8 Å from neutral ligands within an error of 0.7 kcal/mol, and 12 Å from charged ligands within an error of 1.1 kcal/mol for calculating the binding energy in solution.

Design, Synthesis and Structure-Activity Relationship Studies of Protein Kinase CK2 Inhibitors Containing a Purine Scaffold

Protein kinase CK2 (CK2) is involved in the suppression of gene expression, protein synthesis, cell proliferation, and apoptosis, thus making it a target protein for the development of therapeutics toward cancer, nephritis, and coronavirus disease 2019. Using the solvent dipole ordering-based method for virtual screening, we identified and designed new candidate CK2α inhibitors containing purine scaffolds. Virtual docking experiments supported by experimental structure-activity relationship studies identified the importance of the 4-carboxyphenyl group at the 2-position, a carboxamide group at the 6-position, and an electron-rich phenyl group at the 9-position of the purine scaffold. Docking studies based on the crystal structures of CK2α and inhibitor (PDBID: 5B0X) successfully predicted the binding mode of 4-(6-carbamoyl-8-oxo-9-phenyl-8,9-dihydro-7H-purin-2-yl) benzoic acid (11), and the results were used to design stronger small molecule targets for CK2α inhibition. Interaction energy analysis suggested that 11 bound around the hinge region without the water molecule (W1) near Trp176 and Glu81 that is frequently reported in crystal structures of CK2α inhibitor complexes. X-ray crystallographic data for 11 bound to CK2α was in very good agreement with the docking experiments, and consistent with activity. From the structure-activity relationship (SAR) studies presented here, 4-(6-Carbamoyl-9-(4-(dimethylamino)phenyl)-8-oxo-8,9-dihydro-7H-purin-2-yl) benzoic acid (12) was identified as an improved active purine-based CK2α inhibitor with an IC₅₀ of 4.3 µM. These active compounds with an unusual binding mode are expected to inspire new CK2α inhibitors and the development of therapeutics targeting CK2 inhibition.

Bivalent binding mode of an amino-pyrazole inhibitor indicates the potentials for CK2α1-selective inhibitors

Casein kinase 2 (CK2) is a vital protein kinase that consists of two catalytic subunits (CK2α1 and/or CK2α2) and two regulatory subunits (CK2β). CK2α1 is a drug target for nephritis and cancers, while CK2α2 is a serious off-target because its inhibition causes testicular toxicity. High similarity between the isozymes CK2α1 and CK2α2 make it difficult to design CK2α1-specific inhibitors. Herein, the crystal structures of CK2α1 and CK2α2 complexed with a 3-amino-pyrazole inhibitor revealed the remarkable differences in the protein-inhibitor interaction modes. This inhibitor bound to the ATP binding sites of both isozymes in apparently distinct orientations. In addition, another molecule of this inhibitor bound to CK2α1, but not to CK2α2, at the CK2β protein-protein interface. Binding energy calculations and biochemical experiments suggested that this inhibitor possesses the conventional ATP-competitive characteristics with moderate allosteric function in a molecular glue mechanism. These results will assist the potential design of potent and selective CK2α1 inhibitors.

Design, synthesis and biological evaluation of chromone derivatives as novel protein kinase CK2 inhibitors

Protein kinase CK2 is a potential target for the discovery of anticancer drugs. Flavonoids are reported to be effective CK2 inhibitors. Herein, based on structural trimming of flavonoids, a series of chromone-2-aminothiazole derivatives (1a-d, 2a-g, 4a-j, 5a-k) were designed and synthesized by hybridizing the chromone skeleton with 2-aminothiazole scaffold. Among these compounds, compound 5i was the most effective CK2 inhibitor (IC₅₀ = 0.08 μM) and possessed potent anti-proliferative activity against HL-60 tumor cells (IC₅₀ = 0.25 μM). Cellular thermal shift assay (CESTA) confirmed that 5i directly bound to the CK2, and the possible binding mode of 5i toward CK2 was also simulated. Further studies showed that 5i induced the apoptosis of HL-60 cells and arrested the cell cycle. Finally, western-blot analysis showed that 5i could inhibit the downstream of CK2, including α-catenin/Akt pathway and PARP/Survivin pathway.

Proposed Allosteric Inhibitors Bind to the ATP Site of CK2α

CK2α is a ubiquitous, well-studied kinase that is a target for small-molecule inhibition, for treatment of cancers. While many different classes of adenosine 5'-triphosphate (ATP)-competitive inhibitors have been described for CK2α, they tend to suffer from significant off-target activity and new approaches are needed. A series of inhibitors of CK2α has recently been described as allosteric, acting at a previously unidentified binding site. Given the similarity of these inhibitors to known ATP-competitive inhibitors, we have investigated them further. In our thorough structural and biophysical analyses, we have found no evidence that these inhibitors bind to the proposed allosteric site. Rather, we report crystal structures, competitive isothermal titration calorimetry (ITC) and NMR, hydrogen-deuterium exchange (HDX) mass spectrometry, and chemoinformatic analyses that all point to these compounds binding in the ATP pocket. Comparisons of our results and experimental approach with the data presented in the original report suggest that the primary reason for the disparity is nonspecific inhibition by aggregation.

Pyridazine as a privileged structure: An updated review on anticancer activity of pyridazine containing bioactive molecules

Identification of potent anticancer agents with high selectivity and low toxicity remains on the way to human health. Pyridazine featuring advantageous physicochemical properties and antitumor potential usually is regarded as a central core in numerous anticancer derivatives. There are several approved pyridazine-based drugs in the market and analogues currently going through different clinical phases or registration statuses, suggesting pyridazine as a promising drug-like scaffold. The current review is intended to provide a comprehensive and updated overview of pyridazine derivatives as potential anticancer agents. In particular, we focused on their structure-activity relationship (SAR) studies, design strategies, binding modes and biological activities in the hope of offering novel insights for further rational design of more active and less toxic anticancer drugs.

Optimization of pyrazolo[1,5-a]pyrimidines lead to the identification of a highly selective casein kinase 2 inhibitor

Casein kinase 2 (CK2) is a constitutively expressed serine/threonine kinase that has a large diversity of cellular substrates. Thus, CK2 has been associated with a plethora of regulatory functions and dysregulation of CK2 has been linked to disease development in particular to cancer. The broad implications in disease pathology makes CK2 an attractive target. To date, the most advanced CK2 inhibitor is silmitasertib, which has been investigated in clinical trials for treatment of various cancers, albeit several off-targets for silmitasertib have been described. To ascertain the role of CK2 inhibition in cancer, other disease and normal physiology the development of a selective CK2 inhibitor would be highly desirable. In this study we explored the pyrazolo [1,5-a]pyrimidine hinge-binding moiety for the development of selective CK2 inhibitors. Optimization of this scaffold, which included macrocyclization, led to IC20 (31) a compound that displayed high in vitro potency for CK2 (K_D = 12 nM) and exclusive selectivity for CK2. X-ray analysis revealed a canonical type-I binding mode for IC20 (31). However, the polar carboxylic acid moiety that is shared by many CK2 inhibitors including silmitasertib was required for potency but limits the cellular activity of IC20 (31) and the cellular IC₅₀ dropped to the low micromolar range. In summary, IC20 (31) represents a highly selective and potent inhibitor of CK2, which can be used as a tool compound to study CK2 biology and potential new applications for the treatment of diseases.

Structural insights for producing CK2α1-specific inhibitors

Casein kinase 2 catalytic subunit (CK2α) is classified into two subtypes CK2α1 and CK2α2. CK2α1 is a drug discovery target, whereas CK2α2 is an off-target of CK2α1 inhibitors. High amino acid sequence homology between these subtypes hampers efforts to produce ATP competitive inhibitors that are highly selective to CK2α1. Hematein was identified previously as a non-ATP-competitive inhibitor for CK2α1, whereas this compound acts as an ATP competitive CK2α2 inhibitor. Crystal structures of CK2α1 and CK2α2 in complex with hematein revealed distinct binding features that provide structural insights for producing CK2α1-selective inhibitors.

2-Aminothiazole Derivatives as Selective Allosteric Modulators of the Protein Kinase CK2. 1. Identification of an Allosteric Binding Site

CK2 is a ubiquitous Ser/Thr protein kinase involved in the control of various signaling pathways and is known to be constitutively active. In the present study, we identified aryl 2-aminothiazoles as a novel class of CK2 inhibitors, which displayed a non-ATP-competitive mode of action and stabilized an inactive conformation of CK2 in solution. Enzyme kinetics studies, STD NMR, circular dichroism spectroscopy, and native mass spectrometry experiments demonstrated that the compounds bind in an allosteric pocket outside the ATP-binding site. Our data, combined with molecular docking studies, strongly suggested that this new binding site was located at the interface between the αC helix and the flexible glycine-rich loop. A first hit optimization led to compound 7, exhibiting an IC₅₀ of 3.4 μM against purified CK2α in combination with a favorable selectivity profile. Thus, we identified a novel class of CK2 inhibitors targeting an allosteric pocket, offering great potential for further optimization into anticancer drugs.

Thiazole- and selenazole-comprising high-affinity inhibitors possess bright microsecond-scale photoluminescence in complex with protein kinase CK2

A previously disclosed protein kinase (PK) CK2-selective inhibitor 4-(2-amino-1,3-thiazol-5-yl)benzoic acid (ATB) and its selenium-containing counterpart (ASB) revealed remarkable room temperature phosphorescence when bound to the ATP pocket of the protein kinase CK2. Conjugation of these fragments with a mimic of CK2 substrate peptide resulted in bisubstrate inhibitors with increased affinity towards the kinase. Attachment of the fluorescent acceptor dye 5-TAMRA to the conjugates led to significant enhancement of intensity of long-lifetime (microsecond-scale) photoluminescence of both sulfur- and selenium-containing compounds. The developed photoluminescent probes make possible selective determination of the concentration of CK2 in cell lysates and characterization of CK2 inhibitors by means of time-gated measurement of photoluminescence.

Dihydrobenzo[4,5]imidazo[1,2-a]pyrimidine-4-ones as a new class of CK2 inhibitors

Identification of new small molecules inhibiting protein kinase CK2 is highly required for the study of this protein's functions in cell and for the further development of novel pharmaceuticals against a variety of disorders associated with CK2 activity. In this article, a virtual screening of a random small-molecule library was performed and 12 compounds were initially selected for biochemical tests toward CK2. Among them, the most active compound 1 ([Formula: see text]) belonged to dihydrobenzo[4,5]imidazo[1,2-a]pyrimidine-4-ones. The complex of this compound with CK2 was analyzed, and key ligand-enzyme interactions were determined. Then, a virtual screening of 231 dihydrobenzo[4,5]imidazo[1,2-a]pyrimidine-4-one derivatives was performed and 37 compounds were chosen for in vitro testing. It was found that 32 compounds inhibit CK2 with [Formula: see text] values from 2.5 to 7.5 [Formula: see text]. These results demonstrate that dihydrobenzo[4,5]imidazo[1,2-a]pyrimidine-4-one is a novel class of CK2 inhibitors.

Conformational dynamics of human protein kinase CK2α and its effect on function and inhibition

Protein kinase, casein kinase II (CK2), is ubiquitously expressed and highly conserved protein kinase that shows constitutive activity. It phosphorylates a diverse set of proteins and plays crucial role in several cellular processes. The catalytic subunit of this enzyme (CK2α) shows remarkable flexibility as evidenced in numerous crystal structures determined till now. Here, using analysis of multiple crystal structures and long timescale molecular dynamics simulations, we explore the conformational flexibility of CK2α. The enzyme shows considerably higher flexibility in the solution as compared to that observed in crystal structure ensemble. Multiple conformations of hinge region, located near the active site, were observed during the dynamics. We further observed that among these multiple conformations, the most populated conformational state was inadequately represented in the crystal structure ensemble. The catalytic spine, was found to be less dismantled in this state as compared to the "open" hinge/αD state crystal structures. The comparison of dynamics in unbound (Apo) state and inhibitor (CX4945) bound state exhibits inhibitor induced suppression in the overall dynamics of the enzyme. This is especially true for functionally important glycine-rich loop above the active site. Together, this work gives novel insights into the dynamics of CK2α in solution and relates it to the function. This work also explains the effect of inhibitor on the dynamics of CK2α and paves way for development of better inhibitors.

Exploring the Pivotal Role of the CK2 Hinge Region Sub-Pocket in Binding with Tricyclic Quinolone Analogues by Computational Analysis

Protein kinase CK2 has been considered as an attractive therapeutic target of cancer therapy. The tricyclic quinoline compound CX-4945 is the first representative of CK2 inhibitors used in human clinical trials. The binding of non-2,6-naphtyridine substituted compounds 27e (IC₅₀ > 500 nM) and 27h (IC₅₀ > 1000 nM) to CK2 is abolished. However, the unbinding mechanisms due to the key pharmacophore group replacement of compounds 27e and 27h are unveiled. In the present work, combined computational analysis was performed to investigate the underlying structural basis of the low-affinity of two systems. As indicated in the results, the loss of hydrogen bonds between the non-2,6-naphtyridine and the hinge region destroyed the proper recognition of the two complexes. Besides, the allosteric mechanisms between the deviated ligands and the changed regions (G-loop, C-loop and β4/β5 loop) are proposed. Furthermore, energetic analysis was evaluated by detailed energy calculation and residue-based energy decomposition. More importantly, the summary of known polar pharmacophore groups elucidates the pivotal roles of hinge region sub-pocket in the binding of CK2 inhibitors. These results provide rational clues to the fragment-based design of more potent CK2 inhibitors.

Enriching screening libraries with bioactive fragment space

By deconvoluting 238,073 bioactive molecules in the ChEMBL library into extended Murcko ring systems, we identified a set of 2245 ring systems present in at least 10 molecules. These ring systems belong to 2221 clusters by ECFP4 fingerprints with a minimum intracluster similarity of 0.8. Their overlap with ring systems in commercial libraries was further quantified. Our findings suggest that success of a small fragment library is driven by the convergence of effective coverage of bioactive ring systems (e.g., 10% coverage by 1000 fragments vs. 40% by 2million HTS compounds), high enrichment of bioactive ring systems, and low molecular complexity enhancing the probability of a match with the protein targets. Reconciling with the previous studies, bioactive ring systems are underrepresented in screening libraries. As such, we propose a library of virtual fragments with key functionalities via fragmentation of bioactive molecules. Its utility is exemplified by a prospective application on protein kinase CK2, resulting in the discovery of a series of novel inhibitors with the most potent compound having an IC50 of 0.5μM and a ligand efficiency of 0.41kcal/mol per heavy atom.