Structural features underlying selective inhibition of protein kinase CK2 by ATP site-directed tetrabromo-2-benzotriazole

Synthesis and evaluation of chemical linchpins for highly selective CK2α targeting

Casein kinase-2 (CK2) are serine/threonine kinases with dual co-factor (ATP and GTP) specificity, that are involved in the regulation of a wide variety of cellular functions. Small molecules targeting CK2 have been described in the literature targeting different binding pockets of the kinase with a focus on type I inhibitors such as the recently published chemical probe SGC-CK2-1. In this study, we investigated whether known allosteric inhibitors binding to a pocket adjacent to helix αD could be combined with ATP mimetic moieties defining a novel class of ATP competitive compounds with a unique binding mode. Linking both binding sites requires a chemical linking moiety that would introduce a 90-degree angle between the ATP mimetic ring system and the αD targeting moiety, which was realized using a sulfonamide. The synthesized inhibitors were highly selective for CK2 with binding constants in the nM range and low micromolar activity. While these inhibitors need to be further improved, the present work provides a structure-based design strategy for highly selective CK2 inhibitors.

Molecular Interaction Fields Describing Halogen Bond Formable Areas on Protein Surfaces

Molecular interaction fields (MIFs) are three-dimensional interaction maps that describe the intermolecular interactions expected to be formed around target molecules. In this paper, a method for the fast computation of MIFs using the approximation functions of quantum mechanics-level MIFs of small model molecules is proposed. MIF functions of N-methylacetamide with chlorobenzene, bromobenzene, and iodobenzene probes were precisely approximated and used to calculate the MIFs on protein surfaces. This method appropriately reproduced halogen-bond-formable areas around the ligand-binding sites of proteins, where halogen bond formation was suggested in a previous study.

Serine/Threonine Protein Kinases as Attractive Targets for Anti-Cancer Drugs-An Innovative Approach to Ligand Tuning Using Combined Quantum Chemical Calculations, Molecular Docking, Molecular Dynamic Simulations, and Network-like Similarity Graphs

Serine/threonine protein kinases (CK2, PIM-1, RIO1) are constitutively active, highly conserved, pleiotropic, and multifunctional kinases, which control several signaling pathways and regulate many cellular functions, such as cell activity, survival, proliferation, and apoptosis. Over the past decades, they have gained increasing attention as potential therapeutic targets, ranging from various cancers and neurological, inflammation, and autoimmune disorders to viral diseases, including COVID-19. Despite the accumulation of a vast amount of experimental data, there is still no "recipe" that would facilitate the search for new effective kinase inhibitors. The aim of our study was to develop an effective screening method that would be useful for this purpose. A combination of Density Functional Theory calculations and molecular docking, supplemented with newly developed quantitative methods for the comparison of the binding modes, provided deep insight into the set of desirable properties responsible for their inhibition. The mathematical metrics helped assess the distance between the binding modes, while heatmaps revealed the locations in the ligand that should be modified according to binding site requirements. The Structure-Binding Affinity Index and Structural-Binding Affinity Landscape proposed in this paper helped to measure the extent to which binding affinity is gained or lost in response to a relatively small change in the ligand's structure. The combination of the physico-chemical profile with the aforementioned factors enabled the identification of both "dead" and "promising" search directions. Tests carried out on experimental data have validated and demonstrated the high efficiency of the proposed innovative approach. Our method for quantifying differences between the ligands and their binding capabilities holds promise for guiding future research on new anti-cancer agents.

Recent Advances in the Discovery of CK2 Inhibitors

CK2 is a vital enzyme that phosphorylates a large number of substrates and thereby controls many processes in the body. Its upregulation was reported in many cancer types. Inhibitors of CK2 might have anticancer activity, and two compounds are currently under clinical trials. However, both compounds are ATP-competitive inhibitors that may have off-target side effects. The development of allosteric and dual inhibitors can overcome this drawback. These inhibitors showed higher selectivity and specificity for the CK2 enzyme compared to the ATP-competitive inhibitors. The present review summarizes the efforts exerted in the last five years in the design of 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.

Structural Insights on the Role of Halogen Bonding in Protein MEK Kinase-Inhibitor Complexes

Kinases are enzymes that play a critical role in governing essential biological processes. Due to their pivotal involvement in cancer cell signaling, they have become key targets in the development of anti-cancer drugs. Among these drugs, those containing the 2,4-dihalophenyl moiety demonstrated significant potential. Here we show how this moiety, particularly the 2-fluoro-4-iodophenyl one, is crucial for the structural stability of the formed drug-enzyme complexes. Crystallographic analysis of reported kinase-inhibitor complex structures highlights the role of the halogen bonding that this moiety forms with specific residues of the kinase binding site. This interaction is not limited to FDA-approved MEK inhibitors, but it is also relevant for other kinase inhibitors, indicating its broad relevance in the design of this class of drugs.

Effect of histidine protonation state on ligand binding at the ATP-binding site of human protein kinase CK2

Histidine residues contribute to numerous molecular interactions, owing to their structure with the ionizable aromatic side chain with pKa close to the physiological pH. Herein, we studied how the two histidine residues, His115 and His160 of the catalytic subunit of human protein kinase CK2, affect the binding of the halogenated heterocyclic ligands at the ATP-binding site. Thermodynamic studies on the interaction between five variants of hCK2α (WT protein and four histidine mutants) and three ionizable bromo-benzotriazoles and their conditionally non-ionizable benzimidazole counterparts were performed with nanoDSF, MST, and ITC. The results allowed us to identify the contribution of interactions involving the particular histidine residues to ligand binding. We showed that despite the well-documented hydrogen bonding/salt bridge formation dragging the anionic ligands towards Lys68, the protonated His160 also contributes to the binding of such ligands by long-range electrostatic interactions. Simultaneously, His 115 indirectly affects ligand binding, placing the hinge region in open/closed conformations.

Strategies of Targeting CK2 in Drug Discovery: Challenges, Opportunities, and Emerging Prospects

CK2 (casein kinase 2) is a serine/threonine protein kinase that is ubiquitous in eukaryotic cells and plays important roles in a variety of cellular functions, including cell growth, apoptosis, circadian rhythms, DNA damage repair, transcription, and translation. CK2 is involved in cancer pathogenesis and the occurrence of many diseases. Therefore, targeting CK2 is a promising therapeutic strategy. Although many CK2-specific small-molecule inhibitors have been developed, only CX-4945 has progressed to clinical trials. In recent years, novel CK2 inhibitors have gradually become a research hotspot, which is expected to overcome the limitations of traditional inhibitors. Herein, we summarize the structure, biological functions, and disease relevance of CK2 and emphatically analyze the structure-activity relationship (SAR) and binding modes of small-molecule CK2 inhibitors. We also discuss the latest progress of novel strategies, providing insights into new drugs targeting CK2 for clinical practice.

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.

Competition between electrostatic interactions and halogen bonding in the protein-ligand system: structural and thermodynamic studies of 5,6-dibromobenzotriazole-hCK2α complexes

CK2 is a member of the CMGC group of eukaryotic protein kinases and a cancer drug target. It can be efficiently inhibited by halogenated benzotriazoles and benzimidazoles. Depending on the scaffold, substitution pattern, and pH, these compounds are either neutral or anionic. Their binding poses are dictated by a hydrophobic effect (desolvation) and a tug of war between a salt bridge/hydrogen bond (to K68) and halogen bonding (to E114 and V116 backbone oxygens). Here, we test the idea that binding poses might be controllable by pH for ligands with near-neutral pK_a, using the conditionally anionic 5,6-DBBt and constitutively anionic TBBt as our models. We characterize the binding by low-volume Differential Scanning Fluorimetry (nanoDSF), Isothermal Calorimetry (ITC), Hydrogen/Deuterium eXchange (HDX), and X-ray crystallography (MX). The data indicate that the ligand pose away from the hinge dominates for the entire tested pH range (5.5-8.5). The insensitivity of the binding mode to pH is attributed to the perturbation of ligand pK_a upon binding that keeps it anionic in the ligand binding pocket at all tested pH values. However, a minor population of the ligand, detectable only by HDX, shifts towards the hinge in acidic conditions. Our findings demonstrate that electrostatic (ionic) interactions predominate over halogen bonding.

Casein kinase CK2 structure and activities in plants

Casein kinase CK2 is a highly conserved serine/threonine protein kinase and exists in all eukaryotes. It has been demonstrated to be widely involved in the biological processes of plants. The CK2 holoenzyme is a heterotetramer consisting of two catalytic subunits (α and/or α') and two regulatory subunits (β). CK2 in plants is generally encoded by multiple genes, with monomeric and oligomeric forms present in the tissue. Various subunit genes of CK2 have been cloned and characterized from Arabidopsis thaliana, tobacco, maize, wheat, tomato, and other plants. This paper reviews the structural features of CK2, provides a clear classification of its physiological functions and mechanisms of action, and elaborates on the regulation of CK2 activity to provide a knowledge base for subsequent studies of CK2 in plants.

The Immune Regulatory Role of Protein Kinase CK2 and Its Implications for Treatment of Cancer

Protein Kinase CK2, a constitutively active serine/threonine kinase, fulfills its functions via phosphorylating hundreds of proteins in nearly all cells. It regulates a variety of cellular signaling pathways and contributes to cell survival, proliferation and inflammation. CK2 has been implicated in the pathogenesis of hematologic and solid cancers. Recent data have documented that CK2 has unique functions in both innate and adaptive immune cells. In this article, we review aspects of CK2 biology, functions of the major innate and adaptive immune cells, and how CK2 regulates the function of immune cells. Finally, we provide perspectives on how CK2 effects in immune cells, particularly T-cells, may impact the treatment of cancers via targeting CK2.

Synthesis of Novel Halogenated Heterocycles Based on o-Phenylenediamine and Their Interactions with the Catalytic Subunit of Protein Kinase CK2

Protein kinase CK2 is a highly pleiotropic protein kinase capable of phosphorylating hundreds of protein substrates. It is involved in numerous cellular functions, including cell viability, apoptosis, cell proliferation and survival, angiogenesis, or ER-stress response. As CK2 activity is found perturbed in many pathological states, including cancers, it becomes an attractive target for the pharma. A large number of low-mass ATP-competitive inhibitors have already been developed, the majority of them halogenated. We tested the binding of six series of halogenated heterocyclic ligands derived from the commercially available 4,5-dihalo-benzene-1,2-diamines. These ligand series were selected to enable the separation of the scaffold effect from the hydrophobic interactions attributed directly to the presence of halogen atoms. In silico molecular docking was initially applied to test the capability of each ligand for binding at the ATP-binding site of CK2. HPLC-derived ligand hydrophobicity data are compared with the binding affinity assessed by low-volume differential scanning fluorimetry (nanoDSF). We identified three promising ligand scaffolds, two of which have not yet been described as CK2 inhibitors but may lead to potent CK2 kinase inhibitors. The inhibitory activity against CK2α and toxicity against four reference cell lines have been determined for eight compounds identified as the most promising in nanoDSF assay.

Protein kinase CK2: a potential therapeutic target for diverse human diseases

CK2 is a constitutively active Ser/Thr protein kinase, which phosphorylates hundreds of substrates, controls several signaling pathways, and is implicated in a plethora of human diseases. Its best documented role is in cancer, where it regulates practically all malignant hallmarks. Other well-known functions of CK2 are in human infections; in particular, several viruses exploit host cell CK2 for their life cycle. Very recently, also SARS-CoV-2, the virus responsible for the COVID-19 pandemic, has been found to enhance CK2 activity and to induce the phosphorylation of several CK2 substrates (either viral and host proteins). CK2 is also considered an emerging target for neurological diseases, inflammation and autoimmune disorders, diverse ophthalmic pathologies, diabetes, and obesity. In addition, CK2 activity has been associated with cardiovascular diseases, as cardiac ischemia-reperfusion injury, atherosclerosis, and cardiac hypertrophy. The hypothesis of considering CK2 inhibition for cystic fibrosis therapies has been also entertained for many years. Moreover, psychiatric disorders and syndromes due to CK2 mutations have been recently identified. On these bases, CK2 is emerging as an increasingly attractive target in various fields of human medicine, with the advantage that several very specific and effective inhibitors are already available. Here, we review the literature on CK2 implication in different human pathologies and evaluate its potential as a pharmacological target in the light of the most recent findings.

Halogen Atoms in the Protein-Ligand System. Structural and Thermodynamic Studies of the Binding of Bromobenzotriazoles by the Catalytic Subunit of Human Protein Kinase CK2

Binding of a family of brominated benzotriazoles to the catalytic subunit of human protein kinase CK2 (hCK2α) was used as a model system to assess the contribution of halogen bonding to protein–ligand interaction. CK2 is a constitutively active pleiotropic serine/threonine protein kinase that belongs to the CMGC group of eukaryotic protein kinases (EPKs). Due to the addiction of some cancer cells, CK2 is an attractive and well-characterized drug target. Halogenated benzotriazoles act as ATP-competitive inhibitors with unexpectedly good selectivity for CK2 over other EPKs. We have characterized the interaction of bromobenzotriazoles with hCK2α by X-ray crystallography, low-volume differential scanning fluorimetry, and isothermal titration calorimetry. Properties of free ligands in solution were additionally characterized by volumetric and RT-HPLC measurements. Thermodynamic data indicate that the affinity increases with bromo substitution, with greater contributions from 5- and 6-substituents than 4- and 7-substituents. Except for 4,7-disubstituted compounds, the bromobenzotriazoles adopt a canonical pose with the triazole close to lysine 68, which precludes halogen bonding. More highly substituted benzotriazoles adopt many additional noncanonical poses, presumably driven by a large hydrophobic contribution to binding. Some noncanonical ligand orientations allow the formation of halogen bonds with the hinge region. Consistent with a predominantly hydrophobic interaction, the isobaric heat capacity decreases upon ligand binding, the more so the higher the substitution.

Structural basis for the design of bisubstrate inhibitors of protein kinase CK2 provided by complex structures with the substrate-competitive inhibitor heparin

The Ser/Thr kinase CK2, a member of the superfamily of eukaryotic protein kinases, has an acidophilic substrate profile with the substrate recognition sequence S/T-D/E-X-D/E, and it is inhibited by polyanionic substances like heparin. The latter, a highly sulphated glucosamino glycan composed mainly of repeating 2-O-sulpho-α-l-idopyranuronic acid/N,O6-disulpho-α-d-glucosamine disaccharide units, is the longest known substrate-competitive CK2 inhibitor. The structural basis of CK2's preference for anionic substrates and substrate-competitive inhibitors is only vaguely known which limits the value of the substrate-binding region for the structure-based development of CK2 bisubstrate inhibitors. Here, a tetragonal and a monoclinic co-crystal structure of CK2α, the catalytic subunit of CK2, with a decameric heparin fragment are described. In the tetragonal structure, the heparin molecule binds to the polybasic stretch at the beginning of CK2α's helix αC, whereas in the monoclinic structure it occupies the central substrate-recognition region around the P+1 loop. Together, the structures rationalize the inhibitory efficacy of heparin fragments as a function of chain length. The monoclinic CK2α/heparin structure, in which the heparin fragment is particularly well defined, is the first CK2 structure with an anionic inhibitor of considerable size at the central part of the substrate-recognition site. The bound heparin fragment is so close to the binding site of ATP-competitive inhibitors that it can guide the design of linkers and pave the way to efficient CK2 bisubstrate inhibitors in the future.

Protein kinase CK2 inhibition as a pharmacological strategy

CK2 is a constitutively active Ser/Thr protein kinase which phosphorylates hundreds of substrates. Since they are primarily related to survival and proliferation pathways, the best-known pathological roles of CK2 are in cancer, where its targeting is currently being considered as a possible therapy. However, CK2 activity has been found instrumental in many other human pathologies, and its inhibition will expectably be extended to different purposes in the near future. Here, after a description of CK2 features and implications in diseases, we analyze the different inhibitors and strategies available to target CK2, and update the results so far obtained by their in vivo application.

Synthesis, biological properties and structural study of new halogenated azolo[4,5-b]pyridines as inhibitors of CK2 kinase

The new halogenated 1H-triazolo[4,5-b]pyridines and 1H-imidazo[4,5-b]pyridines were synthesised as analogues of known CK2 inhibitors: 4,5,6,7-tetrabromo-1H-benzotriazole (TBBt) and 4,5,6,7-tetrabromo-1H-benzimidazole (TBBi). Their influence on the activity of recombinant human CK2α, CK2α' and PIM1 kinases was determined. The most active inhibitors were di- and trihalogenated 1H-triazolo[4,5-b]pyridines (4a, 5a and 10a) with IC₅₀ values 2.56, 3.82 and 3.26 μM respectively for CK2α. Furthermore, effect on viability of cancer cell lines MCF-7 (human breast adenocarcinoma) and CCRF-CEM (T lymphoblast leukemia) of all final compounds was evaluated. Finally, three crystal structures of complexes of CK2α^1-335 with inhibitors 4a, 5a and 10a were obtained. In addition, new protocol was used to obtain high-resolution crystal structures of CK2α'^Cys336Ser in complex with four inhibitors (4a, 5a, 5b, 10a).

Protein kinase CK2 and ion channels (Review)

Protein kinase CK2 appears as a tetramer or higher molecular weight oligomer composed of catalytic CK2α, CK2α' subunits and non-catalytic regulatory CK2β subunits or as individual subunits. It is implicated in a variety of different regulatory processes, such as Akt signalling, splicing and DNA repair within eukaryotic cells. The present review evaluates the influence of CK2 on ion channels in the plasma membrane. CK2 phosphorylates platform proteins such as calmodulin and ankyrin G, which bind to channel proteins for a physiological transport to and positioning into the membrane. In addition, CK2 directly phosphorylates a variety of channel proteins directly to regulate opening and closing of the channels. Thus, modulation of CK2 activities by specific inhibitors, by siRNA technology or by CRISPR/Cas technology has an influence on intracellular ion concentrations and thereby on cellular signalling. The physiological regulation of the intracellular ion concentration is important for cell survival and correct intracellular signalling. Disturbance of this regulation results in a variety of different diseases including epilepsy, heart failure, cystic fibrosis and diabetes. Therefore, these effects should be considered when using CK2 inhibition as a treatment option for cancer.

Emodin reduces Breast Cancer Lung Metastasis by suppressing Macrophage-induced Breast Cancer Cell Epithelial-mesenchymal transition and Cancer Stem Cell formation

Our previous studies demonstrated that the natural compound emodin blocks the tumor-promoting feedforward interactions between cancer cells and macrophages, and thus ameliorates the immunosuppressive state of the tumor microenvironment. Since tumor-associated macrophages (TAMs) also affect epithelial mesenchymal-transition (EMT) and cancer stem cell (CSC) formation, here we aimed to test if emodin as a neoadjuvant therapy halts breast cancer metastasis by attenuating TAM-induced EMT and CSC formation of breast cancer cells.

Methods: Bioinformatical analysis was performed to examine the correlation between macrophage abundance and EMT/CSC markers in human breast tumors. Cell culture and co-culture studies were performed to test if emodin suppresses TGF-β1 or macrophage-induced EMT and CSC formation of breast cancer cells, and if it inhibits breast cancer cell migration and invasion. Using mouse models, we tested if short-term administration of emodin before surgical removal of breast tumors halts breast cancer post-surgery metastatic recurrence in the lungs. The effects of emodin on TGF-β1 signaling pathways in breast cancer cells were examined by western blots and immunofluorescent imaging.

Results: Macrophage abundance positively correlates with EMT and CSC markers in human breast tumors. Emodin suppressed TGF-β1 production in breast cancer cells and macrophages and attenuated TGF-β1 or macrophage-induced EMT and CSC formation of breast cancer cells. Short-term administration of emodin before surgery halted breast cancer post-surgery metastatic recurrence in the lungs by reducing tumor-promoting macrophages and suppressing EMT and CSC formation in the primary tumors. Mechanistic studies revealed that emodin inhibited both canonical and noncanonical TGF-β1 signaling pathways in breast cancer cells and suppressed transcription factors key to EMT and CSC.

Conclusion: Natural compound emodin suppresses EMT and CSC formation of breast cancer cells by blocking TGF-β1-mediated crosstalk between TAMs and breast cancer cells. Our study provides evidence suggesting that emodin harbors the potential for clinical development as a new effective and safe agent to halt metastatic recurrence of breast cancer.

Crystal structure of Arabidopsis thaliana casein kinase 2 α1

Casein kinase 2 (CK2) is a ubiquitous pleiotropic enzyme that is highly conserved across eukaryotic kingdoms. CK2 is singular amongst kinases as it is highly rigid and constitutively active. Arabidopsis thaliana is widely used as a model system in molecular plant research; the biological functions of A. thaliana CK2 are well studied in vivo and many of its substrates have been identified. Here, crystal structures of the α subunit of A. thaliana CK2 in three crystal forms and of its complex with the nonhydrolyzable ATP analog AMppNHp are presented. While the C-lobe of the enzyme is highly rigid, structural plasticity is observed for the N-lobe. Small but significant displacements within the active cleft are necessary in order to avoid steric clashes with the AMppNHp molecule. Binding of AMppNHp is influenced by a rigid-body motion of the N-lobe that was not previously recognized in maize CK2.

Multitarget Anticancer Agents Based on Histone Deacetylase and Protein Kinase CK2 inhibitors

The design of multitarget drugs (MTDs) has become an innovative approach for the search of effective treatments in complex diseases such as cancer. In this work, we communicate our efforts in the design of multi-targeting histone deacetylase (HDAC) and protein kinase CK2 inhibitors as a novel therapeutic strategy against cancer. Using tetrabromobenzotriazole (TBB) and 2-dimethylamino-4,5,6,7-tetrabromo-benzimidazole (DMAT) as scaffolds for CK2 inhibition, and a hydroxamate to coordinate the zinc atom present in the active site of HDAC (zinc binding group, ZBG), new multitarget inhibitors have been designed and synthesized. According to the in vitro assays, N-Hydroxy-6-(4,5,6,7-tetrabromo-2-(dimethylamino)-1H-benzo[d]imidazol-1-yl)hexanamide (11b) is the most interesting compound, with IC₅₀ values of 0.66; 1.46 and 3.67 µM. for HDAC6; HDAC1 and CK2; respectively. Cellular assays on different cancer cell lines rendered promising results for N-Hydroxy-8-(4,5,6,7-tetrabromo-2-(dimethylamino)-1H-benzo[d]imidazol-1-yl)octanamide (11d). This inhibitor presented the highest cytotoxic activity, proapoptotic capability, and the best mitochondria-targeting and multidrug-circumventing properties, thus being the most promising drug candidate for further in vivo studies.

The Marine Natural Product Manzamine A Inhibits Cervical Cancer by Targeting the SIX1 Protein

Natural products remain an important source of drug leads covering unique chemical space and providing significant therapeutic value for the control of cancer and infectious diseases resistant to current drugs. Here, we determined the antiproliferative activity of a natural product manzamine A (1) from an Indo-Pacific sponge following various in vitro cellular assays targeting cervical cancer (C33A, HeLa, SiHa, and CaSki). Our data demonstrated the antiproliferative effects of 1 at relatively low and non-cytotoxic concentrations (up to 4 μM). Mechanistic investigations confirmed that 1 blocked cell cycle progression in SiHa and CaSki cells at G1/S phase and regulated cell cycle-related genes, including restoration of p21 and p53 expression. In apoptotic assays, HeLa cells showed the highest sensitivity to 1 as compared to other cell types (C33A, SiHa, and CaSki). Interestingly, 1 decreased the levels of the oncoprotein SIX1, which is associated with oncogenesis in cervical cancer. To further investigate the structure-activity relationship among manzamine A (1) class with potential antiproliferative activity, molecular networking facilitated the efficient identification, dereplication, and assignment of structures from the manzamine class and revealed the significant potential in the design of optimized molecules for the treatment of cervical cancer. These data suggest that this sponge-derived natural product class warrants further attention regarding the design and development of novel manzamine analogues, which may be efficacious for preventive and therapeutic treatment of cancer. Additionally, this study reveals the significance of protecting fragile marine ecosystems from climate change-induced loss of species diversity.

Thermodynamic contribution of iodine atom to the binding of heterogeneously polyhalogenated benzotriazoles by the catalytic subunit of human protein kinase CK2

A series of novel benzotriazole derivatives containing iodine atom(s) were synthesized. The binding of these compounds to the catalytic subunit of human protein kinase CK2 was evaluated using differential scanning fluorimetry. The obtained thermodynamic data were compared with those determined previously for the brominated and chlorinated benzotriazole analogues to get a deeper insight into the thermodynamic contribution of iodine substitution to the free energy of ligand binding. We have shown that iodine atom(s) attached to the benzene ring of benzotriazole enhance(s) its binding by the target protein. This effect is the strongest when two iodine atoms are attached at positions peripheral to the triazole ring, which according to the structures deposited in protein data bank may be indicative for the formation of the halogen bond between iodine and carbonyl groups of residues located in the hinge region of the protein. Finally, quantitative structure-activity relationship analysis pointed the solute hydrophobicity as the main factor contributing to the binding affinity.

CK2 inhibition, lipophilicity and anticancer activity of new N1versus N2-substituted tetrabromobenzotriazole regioisomers

Both the type and position of polar group substitutions in polybrominated benzotriazoles dramatically change their lipophilicity, kinase inhibition and anticancer activity.

Making endo-cyclizations favorable again: a conceptually new synthetic approach to benzotriazoles via azide group directed lithiation/cyclization of 2-azidoaryl bromides

Although benzotriazoles are important and ubiquitous, currently there is only one conceptual approach to their synthesis: bridging the two ortho-amino groups with an electrophilic nitrogen atom. Herein, we disclose a new practical alternative - the endo-cyclization of 2-azidoaryl lithiums obtained in situ from 2-azido-aryl bromides. The scope of the reaction is illustrated using twenty-four examples with a variety of alkyl, alkoxy, perfluoroalkyl, and halogen substituents. We found that the directing effect of the azide group allows selective metal-halogen exchange in aryl azides containing several bromine atoms. Furthermore, (2-bromophenyl)diazomethane undergoes similar cyclization to give an indazole. Thus, cyclizations of aryl lithiums containing an ortho-X = Y = Z group emerge as a new general approach for the synthesis of aromatic heterocycles. DFT computations suggested that the observed endo-selectivity applies to the anionic cyclizations of other functionalities that undergo "1,1-additions" (i.e., azides, diazo compounds, and isonitriles). In contrast, cyclizations with the heteroatomic functionalities that follow the "1,2-addition" pattern (cyanates, thiocyanates, isocyanates, isothiocyanates, and nitriles) prefer the exo-cyclization path. Hence, such reactions expand the current understanding of stereoelectronic factors in anionic cyclizations.

Rational drug-design approach supported with thermodynamic studies - a peptide leader for the efficient bi-substrate inhibitor of protein kinase CK2

Numerous inhibitors of protein kinases act on the basis of competition, targeting the ATP binding site. In this work, we present a procedure of rational design of a bi-substrate inhibitor, complemented with biophysical assays. The inhibitors of this type are commonly engineered by combining ligands carrying an ATP-like part with a peptide or peptide-mimicking fragment that determines specificity. Approach presented in this paper led to generation of a specific system for independent screening for efficient ligands and peptides, by means of thermodynamic measurements, that assessed the ability of the identified ligand and peptide to combine into a bi-substrate inhibitor. The catalytic subunit of human protein kinase CK2 was used as the model target. Peptide sequence was optimized using peptide libraries [KGDE]-[DE]-[ST]-[DE]_3-4-NH_2, originated from the consensus CK2 sequence. We identified KESEEE-NH₂ peptide as the most promising one, whose binding affinity is substantially higher than that of the reference RRRDDDSDDD peptide. We assessed its potency to form an efficient bi-substrate inhibitor using tetrabromobenzotriazole (TBBt) as the model ATP-competitive inhibitor. The formation of ternary complex was monitored using Differential Scanning Fluorimetry (DSF), Microscale Thermophoresis (MST) and Isothermal Titration Calorimetry (ITC).

Diacritic Binding of an Indenoindole Inhibitor by CK2α Paralogs Explored by a Reliable Path to Atomic Resolution CK2α' Structures

CK2α and CK2α' are the two isoforms of the catalytic subunit of human protein kinase CK2, an important target for cancer therapy. They have similar, albeit not identical functional and structural properties, and were occasionally reported to be inhibited with distinct efficacies by certain ATP-competitive ligands. Here, we present THN27, an indeno[1,2-b]indole derivative, as a further inhibitor with basal isoform selectivity. The selectivity disappears when measured using CK2α/CK2α' complexes with CK2β, the regulatory CK2 subunit. Co-crystal structures of THN27 with CK2α and CK2α' reveal that subtle differences in the conformational variability of the interdomain hinge region are correlated with the observed effect. In the case of CK2α', a crystallographically problematic protein so far, this comparative structural analysis required the development of an experimental strategy that finally enables atomic resolution structure determinations with ab initio phasing of potentially any ATP-competitive CK2 inhibitor and possibly many non-ATP-competitive ligands as well bound to CK2α'.

From a MMP2/CK2 multitarget approach to the identification of potent and selective MMP13 inhibitors

In this article, we describe our efforts in the search of MMP2/CK2 dual targeting inhibitors. We have followed a rational drug design approach based on our experience in the selective inhibition of these two enzymes. We have successfully obtained highly active MMP2 (10, IC50 = 70 nM; 11, IC50 = 100 nM) and CK2 (16a, IC50 = 500 nM) inhibitors. However, structural fine tuning of these small molecules to simultaneously target both enzymes turned out to be an unattainable goal. Unexpectedly, we were lucky to identify new and selective MMP13 inhibitors (10, IC50 = 3.7 nM and 11, IC50 = 5.6 nM) with a novel TBB-derived scaffold. These compounds constitute an interesting starting point for further optimization.

Biological properties and structural study of new aminoalkyl derivatives of benzimidazole and benzotriazole, dual inhibitors of CK2 and PIM1 kinases

The new aminoalkyl-substituted derivatives of known CK2 inhibitors 4,5,6,7-tetrabromo-1H-benzimidazole (TBBi) and 4,5,6,7-tetrabromo-1H-benzotriazole (TBBt) were synthesized, and their influence on the activity of recombinant human CK2 α, CK2 holoenzyme and PIM1 kinases was evaluated. All derivatives inhibited the activity of studied kinases and the most efficient were aminopropyl-derivatives 8b and 14b. These compounds also exerted inhibition of cancer cell lines - CCRF-CEM (acute lymphoblastoid leukemia), MCF-7 (human breast cancer), and PC-3 (prostate cancer) proliferation and their EC₅₀ is comparable with the value for clinically studied CK2 inhibitor CX-4945. Preliminary structure activity relationship analysis indicated that the spacer length affected antitumor potency, and two to three methylene units were more favorable. The complex of CK2 α^1-335/8b was crystallized, both under high-salt conditions and under low-salt conditions giving crystals which diffracted X-rays to about 2.4 Å resolution, what enabled the determination of the corresponding 3D-structures.

Changes in potassium channel modulation may underlie afterhyperpolarization plasticity in oxytocin neurons during late pregnancy

Oxytocin (OT) neurons exhibit larger afterhyperpolarizations (AHPs) following spike trains during late pregnancy and lactation, times when these neurons fire in bursts and release more OT associated with labor and lactation. Calcium-dependent AHPs mediated by SK channels show this plasticity, and are reduced when the channel complex is phosphorylated by casein kinase 2 (CK2), and increased when dephosphorylated by protein phosphatase (PP)2A, by altering Ca²⁺ sensitivity. We compared AHP currents in supraoptic OT neurons after CK2 inhibition with 4,5,6,7-tetrabromobenzotriazole (TBB), or PP1-PP2A inhibition with okadaic acid (OA), to determine the roles of these enzymes in AHP plasticity, focusing on the peak current at 100 ms representing the SK-mediated, medium AHP (I_mAHP). In slices from virgin and two groups of pregnant rats [embryonic days (E18-19, or E20-21], I_mAHPs were evoked with 3-, 10-, and 17-spike trains (20 Hz). With 3-spike trains, TBB increased the I_mAHP to the greatest extent in virgin compared with both groups of pregnant animals. A difference between virgins and E20-21 rats was also evident with a 10-spike train but the increases in I_mAHPs were similar among groups with 17-spike trains. In contrast, OA, while consistently reducing the I_mAHP in all cases, showed no differential effects among groups. In Western blots, CK2α, CK2β, PP2A-A, PP2A-B, and PP2A-C were found in supraoptic lysates, and expression of CK2α and CK2β was reduced in E20-21 rats. Coimmunoprecipitation revealed that calmodulin, CK2α, and PP2A-C were associated with SK3 protein. The results suggest that a downregulation of SK3-associated CK2α during late pregnancy may increase the sensitivity of the SK calmodulin (Ca²⁺) sensor for I_mAHP, contributing to the enhanced I_mAHP. NEW & NOTEWORTHY The article demonstrates for the first time that enhancement in spike afterhyperpolarizations in oxytocin neurons during pregnancy may be related to a downregulation in the small-conductance Ca²⁺-activated potassium channels (SK)/calmodulin binding protein casein kinase 2, which phosphorylates the SK channel complex and reduces its Ca²⁺ sensitivity.

Unexpected Binding Mode of a Potent Indeno[1,2-b]indole-Type Inhibitor of Protein Kinase CK2 Revealed by Complex Structures with the Catalytic Subunit CK2α and Its Paralog CK2α'

Protein kinase CK2, a member of the eukaryotic protein kinase superfamily, is associated with cancer and other human pathologies and thus an attractive drug target. The indeno[1,2-b]indole scaffold is a novel lead structure to develop ATP-competitive CK2 inhibitors. Some indeno[1,2-b]indole-based CK2 inhibitors additionally obstruct ABCG2, an ABC half transporter overexpressed in breast cancer and co-responsible for drug efflux and resistance. Comprehensive derivatization studies revealed substitutions of the indeno[1,2-b]indole framework that boost either the CK2 or the ABCG2 selectivity or even support the dual inhibition potential. The best indeno[1,2-b]indole-based CK2 inhibitor described yet (IC₅₀ = 25 nM) is 5-isopropyl-4-(3-methylbut-2-enyl-oxy)-5,6,7,8-tetrahydroindeno[1,2-b]indole-9,10-dione (4p). Herein, we demonstrate the membrane permeability of 4p and describe co-crystal structures of 4p with CK2α and CK2α′, the paralogs of human CK2 catalytic subunit. As expected, 4p occupies the narrow, hydrophobic ATP site of CK2α/CK2α′, but surprisingly with a unique orientation: its hydrophobic substituents point towards the solvent while its two oxo groups are hydrogen-bonded to a hidden water molecule. An equivalent water molecule was found in many CK2α structures, but never as a critical mediator of ligand binding. This unexpected binding mode is independent of the interdomain hinge/helix αD region conformation and of the salt content in the crystallization medium.

Ultrafast dynamics induced by the interaction of molecules with electromagnetic fields: Several quantum, semiclassical, and classical approaches

Several strategies for simulating the ultrafast dynamics of molecules induced by interactions with electromagnetic fields are presented. After a brief overview of the theory of molecule-field interaction, we present several representative examples of quantum, semiclassical, and classical approaches to describe the ultrafast molecular dynamics, including the multiconfiguration time-dependent Hartree method, Bohmian dynamics, local control theory, semiclassical thawed Gaussian approximation, phase averaging, dephasing representation, molecular mechanics with proton transfer, and multipolar force fields. In addition to the general overview, some focus is given to the description of nuclear quantum effects and to the direct dynamics, in which the ab initio energies and forces acting on the nuclei are evaluated on the fly. Several practical applications, performed within the framework of the Swiss National Center of Competence in Research "Molecular Ultrafast Science and Technology," are presented: These include Bohmian dynamics description of the collision of H with H2, local control theory applied to the photoinduced ultrafast intramolecular proton transfer, semiclassical evaluation of vibrationally resolved electronic absorption, emission, photoelectron, and time-resolved stimulated emission spectra, infrared spectroscopy of H-bonding systems, and multipolar force fields applications in the condensed phase.

Looking Back, Looking Forward at Halogen Bonding in Drug Discovery

Halogen bonding has emerged at the forefront of advances in improving ligand: receptor interactions. In particular the newfound ability of this extant non-covalent-bonding phenomena has revolutionized computational approaches to drug discovery while simultaneously reenergizing synthetic approaches to the field. Here we survey, via examples of classical applications involving halogen atoms in pharmaceutical compounds and their biological hosts, the unique advantages that halogen atoms offer as both Lewis acids and Lewis bases.

Ensemble-based virtual screening: identification of a potential allosteric inhibitor of Bcr-Abl

Ensemble-based virtual screening using different conformations of a target protein is gaining popularity, as it can leverage information from target flexibility for effective lead identification. In this paper, molecular dynamics simulation followed by RMSD-based clustering was employed to generate and choose distinct conformations of Bcr-Abl. Three representative structures from the most-populated clusters along with the crystal structure conformation (PDBID: 3K5V) were used to perform docking-based virtual screening of 14,400 compounds (in the Maybridge database) in order to identify potential allosteric site binders. Seven compounds found as hits in at least three of the four virtual screenings had higher Glide docking scores than the co-crystallized allosteric inhibitor GNF-2. Detailed computational analyses of the seven hits identified SEW02675 (ΔG _bind = -164.92 kJ/mol with the wild-type (wt) Bcr-Abl and -167.37 kJ/mol with the T334I Bcr-Abl mutant) as a better allosteric site binder with both the wt and the mutant Bcr-Abl protein than the reference allosteric inhibitor GNF-2 (ΔG _bind = -103.12 with wt and -142.96 kJ/mol with T334I). Moreover, the presence of SEW02675 in the allosteric site enhanced the binding of imatinib (ΔG _bind = -367.58 with wt and -294.56 kJ/mol with T334I) to the ATP sites of the wt and the mutant Bcr-Abl. However, when GNF-2 was present in the allosteric site, the binding of imatinib (ΔG _bind = -351.76 with wt and -273.94 kJ/mol with T334I) to the ATP site was weaker. The in silico findings suggest that SEW02675 could be used in combination with imatinib to treat chronic myeloid leukemia, and that it could help to overcome resistance due to T334I Bcr-Abl mutation. Graphical abstract Virtual screening strategy to identify allosteric inhbitors of Bcr-Abl for the treatment of Chronic myeloid leukemia.

ITC-derived binding affinity may be biased due to titrant (nano)-aggregation. Binding of halogenated benzotriazoles to the catalytic domain of human protein kinase CK2

The binding of four bromobenzotriazoles to the catalytic subunit of human protein kinase CK2 was assessed by two complementary methods: Microscale Thermophoresis (MST) and Isothermal Titration Calorimetry (ITC). New algorithm proposed for the global analysis of MST pseudo-titration data enabled reliable determination of binding affinities for two distinct sites, a relatively strong one with the K_d of the order of 100 nM and a substantially weaker one (K_d > 1 μM). The affinities for the strong binding site determined for the same protein-ligand systems using ITC were in most cases approximately 10-fold underestimated. The discrepancy was assigned directly to the kinetics of ligand nano-aggregates decay occurring upon injection of the concentrated ligand solution to the protein sample. The binding affinities determined in the reverse ITC experiment, in which ligands were titrated with a concentrated protein solution, agreed with the MST-derived data. Our analysis suggests that some ITC-derived K_d values, routinely reported together with PDB structures of protein-ligand complexes, may be biased due to the uncontrolled ligand (nano)-aggregation, which may occur even substantially below the solubility limit.

Impact of Quadrupolar Electrostatics on Atoms Adjacent to the Sigma-Hole in Condensed-Phase Simulations

Halogenation is one of the cases for which advanced molecular simulation methods are mandatory for quantitative and predictive studies. The present work provides a systematic investigation of the importance of higher-order multipoles on specific sites of halobenzenes, other than the halogen, for static and dynamic properties in condensed-phase simulations. For that purpose, solute-solvent interactions using point charge (PC), multipole (MTP), and hybrid point charge/multipole (HYB) electrostatic models are analyzed in regions of halogen bonding and extended to regions of π orbitals of phenyl carbons. Using molecular dynamics simulations and quantum chemical methods, it is found that the sigma-hole does not only affect the halogen and the carbon bound to it but its effect extends to the carbons adjacent to the CX bond. This effect increases with the magnitude of the positive potential of the sigma-hole. With the MTP and HYB3 models, all hydration free energies of the PhX compounds are reproduced within 0.1 kcal/mol. Analysis of pair distribution functions and hydration free energies of halogenated benzenes provides a microscopic explanation why "point charge"-based representations with off-site charges fail in reproducing thermodynamic properties of the sigma-hole. Application of the hybrid models to study protein-ligand binding demonstrates both their accuracy and computational efficiency.

The Halogen Bond

The halogen bond occurs when there is evidence of a net attractive interaction between an electrophilic region associated with a halogen atom in a molecular entity and a nucleophilic region in another, or the same, molecular entity. In this fairly extensive review, after a brief history of the interaction, we will provide the reader with a snapshot of where the research on the halogen bond is now, and, perhaps, where it is going. The specific advantages brought up by a design based on the use of the halogen bond will be demonstrated in quite different fields spanning from material sciences to biomolecular recognition and drug design.

Design and synthesis of novel dual-target agents for HDAC1 and CK2 inhibition

Drug entities able to address multiple targets can be more effective than those directed to just one biological target.

Benzotriazole: An overview on its versatile biological behavior

Discovered in late 1960, azoles are heterocyclic compounds class which constitute the largest group of available antifungal drugs. Particularly, the imidazole ring is the chemical component that confers activity to azoles. Triazoles are obtained by a slight modification of this ring and similar or improved activities as well as less adverse effects are reported for triazole derivatives. Consequently, it is not surprising that benzimidazole/benzotriazole derivatives have been found to be biologically active. Since benzimidazole has been widely investigated, this review is focused on defining the place of benzotriazole derivatives in biomedical research, highlighting their versatile biological properties, the mode of action and Structure Activity Relationship (SAR) studies for a variety of antimicrobial, antiparasitic, and even antitumor, choleretic, cholesterol-lowering agents.

Thermodynamics parameters for binding of halogenated benzotriazole inhibitors of human protein kinase CK2α

The interaction of human CK2α (hCK2α) with nine halogenated benzotriazoles, TBBt and its analogues representing all possible patterns of halogenation on the benzene ring of benzotriazole, was studied by biophysical methods. Thermal stability of protein-ligand complexes, monitored by calorimetric (DSC) and optical (DSF) methods, showed that the increase in the mid-point temperature for unfolding of protein-ligand complexes (i.e. potency of ligand binding to hCK2α) follow the inhibitory activities determined by biochemical assays. The dissociation constant for the ATP-hCK2α complex was estimated with the aid of microscale thermophoresis (MST) as 4.3±1.8 μM, and MST-derived dissociation constants determined for halogenated benzotriazoles, when converted according to known ATP concentrations, perfectly reconstruct IC50 values determined by the biochemical assays. Ligand-dependent quenching of tyrosine fluorescence, together with molecular modeling and DSC-derived heats of unfolding, support the hypothesis that halogenated benzotriazoles bind in at least two alternative orientations, and those that are efficient hCK2α inhibitors bind in the orientation which TBBt adopts in its complex with maize CK2α. DSC-derived apparent heat for ligand binding (ΔΔHbind) is driven by intermolecular electrostatic interactions between Lys68 and the triazole ring of the ligand, as indicated by a good correlation between ΔΔHbind and ligand pKa. Overall results, additionally supported by molecular modeling, confirm that a balance of hydrophobic and electrostatic interactions contribute predominantly (~40 kJ/mol), relative to possible intermolecular halogen/hydrogen bonding (less than 10 kJ/mol), in binding of halogenated benzotriazoles to the ATP-binding site of hCK2α. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases.

Identification of protein kinase CK2 inhibitors using solvent dipole ordering virtual screening

Novel protein kinase CK2 inhibitors were identified using the solvent dipole ordering virtual screening method. A total of 26 compounds categorized in 15 distinct scaffold classes inhibited greater than 50% of enzyme activity at 50 μM, and eight exhibited IC50 values less than 10 μM. Most of the identified compounds are lead-like and dissimilar to known inhibitors. The crystal structures of two of the CK2 complexes revealed the high accuracy of the predicted binding modes.