Search

Abstract Bivalent molecules consisting of groups connected through bridging linkers often exhibit strong target binding and unique biological effects. However, developing bivalent inhibitors with the desired activity is challenging due to the dual motif architecture of these molecules and the variability that can be introduced through differing linker structures and geometries. We report a set of alternatively linked bivalent EGFR inhibitors that simultaneously occupy the ATP substrate and allosteric pockets. Crystal structures show that initial and redesigned linkers bridging a trisubstituted imidazole ATP-site inhibitor and dibenzodiazepinone allosteric-site inhibitor proved successful in spanning these sites. The re-engineered linker yielded a compound that exhibited significantly higher potency (~60 pM) against the drug-resistant EGFR L858R/T790M and L858R/T790M/C797S, which was superadditive as compared with the parent molecules. The enhanced potency is attributed to factors stemming from the linker connection to the allosteric-site group and informs strategies to engineer linkers in bivalent agent design.

Abstract Fibroblast growth factor receptor (FGFR)−2 can be inhibited by FGFR-selective or non-selective tyrosine kinase inhibitors (TKIs). Selective TKIs are approved for cholangiocarcinoma (CCA) with FGFR2 fusions; however, their application is limited by a characteristic pattern of adverse events or evocation of kinase domain mutations. A comprehensive characterization of a patient cohort treated with the non-selective TKI lenvatinib reveals promising efficacy in FGFR2-driven CCA. In a bed-to-bench approach, we investigate FGFR2 fusion proteins bearing critical tumor-relevant point mutations. These mutations confer growth advantage of tumor cells and increased resistance to selective TKIs but remain intriguingly sensitive to lenvatinib. In line with clinical observations, in-silico analyses reveal a more favorable interaction pattern of lenvatinib with FGFR2, including an increased flexibility and ligand efficacy, compared to FGFR-selective TKIs. Finally, the treatment of a patient with progressive disease and a newly developed kinase mutation during therapy with a selective inhibitor results in a striking response to lenvatinib. Our in vitro, in silico, and clinical data suggest that lenvatinib is a promising treatment option for FGFR2-driven CCA, especially when insurmountable adverse reactions of selective TKIs or acquired kinase mutations occur.

Abstract We report a series of 1,3-diphenylureido hydroxamate HDAC inhibitors evaluated against sensitive and drug-resistant P. falciparum strains. Compounds 8a–d show potent antiplasmodial activity, indicating that a phenyl spacer allows improved potency relative to cinnamyl and di-hydrocinnamyl linkers. In vitro, mechanistic studies demonstrated target activity for PfHDAC1 on a recombinant level, which agreed with cell quantification of the acetylated histone levels. Compounds 6c, 7c, and 8c, identified as the most active in phenotypic assays and PfHDAC1 enzymatic inhibition. Compound 8c stands out as a remarkable inhibitor, displaying an impressive 85% inhibition of PfHDAC1, with an IC50 value of 0.74 µM in the phenotypic screening on Pf3D7 and 0.8 µM against multidrug-resistant PfDd2 parasites. Despite its potent inhibition of PfHDAC1, 8c remains the least active on human HDAC1, displaying remarkable selectivity. In silico studies suggest that the phenyl linker has an ideal length in the series for permitting effective interactions of the hydroxamate with PfHDAC1 and that this compound series could bind as well as in HsHDAC1. Taken together, these results highlight the potential of diphenylurea hydroxamates as a privileged scaffold for the generation of potent antimalarial HDAC inhibitors with improved selectivity over human HDACs.

A current clinical challenge in cancer is multidrug resistance (MDR) mediated by ABC transporters. Breast cancer resistance protein (BCRP) or ABCG2 transporter is one of the most important ABC transporters implicated in MDR and the use of inhibitors is a promising approach to overcome the resistance in cancer. This study aimed to characterize the molecular mechanism of ABCG2 inhibitors identified by a repurposing drug strategy using antiviral, anti-inflammatory and antiparasitic agents. Lopinavir and ivermectin can be considered as pan-inhibitors of ABC transporters, since both compounds inhibited ABCG2, P-glycoprotein and MRP1. They inhibited ABCG2 activity showing IC50 values of 25.5 and 23.4 µM, respectively. These drugs were highly cytotoxic and not transported by ABCG2. Additionally, these drugs increased the 5D3 antibody binding and did not affect the mRNA and protein expression levels. Cell-based analysis of the type of inhibition suggested a non-competitive inhibition, which was further corroborated by in silico approaches of molecular docking and molecular dynamics simulations. These results showed an overlap of the lopinavir and ivermectin binding sites on ABCG2, mainly interacting with E446 residue. However, the substrate mitoxantrone occupies a different site, binding to the F436 region, closer to the L554/L555 plug. In conclusion, these results revealed the mechanistic basis of lopinavir and ivermectin interaction with ABCG2.

The Kelch-like ECH-associated protein 1 (KEAP1) – Nuclear factor erythroid 2 -related factor 2 (NRF2) pathway is the major transcriptional stress response system in cells against oxidative and electrophilic stress. NRF2 is frequently constitutively active in many cancers, rendering the cells resistant to chemo- and radiotherapy. Loss-of-function (LOF) mutations in the repressor protein KEAP1 are common in non-small cell lung cancer, particularly adenocarcinoma. While the mutations can occur throughout the gene, they are enriched in certain areas, indicating that these may have unique functional importance. In this study, we show that in the GSEA analysis of TCGA lung adenocarcinoma RNA-seq data, the KEAP1 mutations in R320 and R470 were associated with enhanced Tumor Necrosis Factor alpha (TNFα) – Nuclear Factor kappa subunit B (NFκB) signaling as well as MYC and MTORC1 pathways. To address the functional role of these hotspot mutations, affinity purification and mass spectrometry (AP-MS) analysis of wild type (wt) KEAP1 and its mutation forms, R320Q and R470C were employed to interrogate differences in the protein interactome. We identified TNF receptor associated factor 2 (TRAF2) as a putative protein interaction partner. Both mutant KEAP1 forms showed increased interaction with TRAF2 and other anti-apoptotic proteins, suggesting that apoptosis signalling could be affected by the protein interactions. A549 lung adenocarcinoma cells overexpressing mutant KEAP1 showed high TRAF2-mediated NFκB activity and increased protection against apoptosis, XIAP being one of the key proteins involved in anti-apoptotic signalling. To conclude, KEAP1 R320Q and R470C and its interaction with TRAF2 leads to activation of NFκB pathway, thereby protecting against apoptosis.

Abstract Anxiety is a troublesome symptom for many patients, especially those suffering from schizophrenia. Its regulation involves serotonin receptors, targeted e.g. by antipsychotics or psychedelics such as LSD. 5-HT2A receptors are known for an extremely long LSD residence time, enabling minute doses to exert a long-lasting effect. In this work, we explore the changes in anxiety-like processes induced by the previously reported antipsychotic, D2AAK1. In vivo studies revealed that the effect of D2AAK1 on the anxiety is mediated through serotonin 5-HT1A and 5-HT2A receptors, and that it is time-dependent (anxiogenic after 30 min, anxiolytic after 60 min) and dose-dependent. The funnel metadynamics simulations suggest complicated ligand-5HT2AR interactions, involving an allosteric site located under the third extracellular loop, which is a possible explanation of the time-dependency. The binding of D2AAK1 at the allosteric site results in a broader opening of the extracellular receptor entry, possibly altering the binding kinetics of orthosteric ligands.

Abstract SARS-CoV-2’s papain-like protease (PLpro) interaction with ligands has recently been explored with a myriad of crystal structures. We used molecular dynamics (MD) simulations to study different PLpro-ligand complexes, their ligand-induced conformational changes, and interactions. We focused on inhibitors reported with known IC50 against PLpro, namely GRL-0617, XR8-89, PLP_Snyder530, and Sander’s recently published compound 7 (CPD7), and compared these trajectories against the apostructure (Apo), with a total of around 60 µs worth simulation data. We aimed to study the conformational changes using molecular dynamics simulations for the inhibitors in the PLpro. PCA analyses and the MSM models revealed distinct conformations of PLpro in the absence/presence of ligands and proposed that BL2-loop contributes to the accessibility of these inhibitors. Further, bulkier substituents closer to Tyr268 and Gln269 could improve inhibition of SARS-CoV-2 PLpro by occupying the region between BL2-groove and BL2-loop, but we also expand on the relevance of exploring multiple PLpro sub-pockets to improve inhibition.

Abstract Intrahepatic cholangiocarcinoma (iCCA) has emerged as a promising candidate for precision medicine, especially in the case of activating FGFR2 gene fusions. In addition to fusions, a considerable fraction of iCCA patients reveals FGFR2 mutations, which might lead to uncontrolled activation of the FGFR2 pathway but are mostly of unknown functional significance. A current challenge for molecular tumor boards (MTB) is to predict the functional consequences of such FGFR2 alterations to guide potential treatment decisions. We report two iCCA patients with extracellular and juxtamembrane FGFR2 mutations. After in silico investigation of the alterations and identification of activated FGFR2 downstream targets in tumor specimens by immunohistochemistry and transcriptome analysis, the MTB recommended treatment with an FGFR-inhibiting tyrosine kinase inhibitor. Both patients developed a rapidly detectable and prolonged partial response to treatment. These two cases suggest an approach to characterize further detected FGFR2 mutations in iCCA to enable patients´ selection for a successful application of the FGFR -inhibiting drugs.

Inhibition of ABC transporters is a promising approach to overcome multidrug resistance in cancer. Herein, we report the characterization of a potent ABCG2 inhibitor, namely, chromone 4a (C4a). Molecular docking and in vitro assays using ABCG2 and P-glycoprotein (P-gp) expressing membrane vesicles of insect cells revealed that C4a interacts with both transporters, while showing selectivity toward ABCG2 using cell-based transport assays. C4a inhibited the ABCG2-mediated efflux of different substrates and molecular dynamic simulations demonstrated that C4a binds in the Ko143-binding pocket. Liposomes and extracellular vesicles (EVs) of Giardia intestinalis and human blood were used to successfully bypass the poor water solubility and delivery of C4a as assessed by inhibition of the ABCG2 function. Human blood EVs also promoted delivery of the well-known P-gp inhibitor, elacridar. Here, for the first time, we demonstrated the potential use of plasma circulating EVs for drug delivery of hydrophobic drugs targeting membrane proteins.

Trypanosoma cruzi, the etiological agent of Chagas disease, relies on finely coordinated epigenetic regulation during the transition between hosts. Herein we targeted the silent information regulator 2 (Sir2) enzyme, a NAD+-dependent class III histone deacetylase, to interfere with the parasites’ cell cycle. A combination of molecular modelling with on-target experimental validation was used to discover new inhibitors from commercially available compound libraries. We selected six inhibitors from the virtual screening, which were validated on the recombinant Sir2 enzyme. The most potent inhibitor (CDMS-01, IC50 = 40 μM) was chosen as a potential lead compound.

Dry age-related macular degeneration (AMD) is a currently untreatable vision threatening disease. Impaired proteasomal clearance and autophagy in the retinal pigment epithelium (RPE) and subsequent photoreceptor damage are connected with dry AMD, but detailed pathophysiology is still unclear. In this paper, we discover inhibition of cytosolic protease, prolyl oligopeptidase (PREP), as a potential pathway to treat dry AMD. We showed that PREP inhibitor exposure induced autophagy in the RPE cells, shown by increased LC3-II levels and decreased p62 levels. PREP inhibitor treatment increased total levels of autophagic vacuoles in the RPE cells. Global proteomics was used to examine the phenotype of a commonly used cell model displaying AMD characteristics, oxidative stress and altered protein metabolism, in vitro. These RPE cells displayed induced protein aggregation and clear alterations in macromolecule metabolism, confirming the relevance of the cell model. Differences in intracellular target engagement of PREP inhibitors were observed with cellular thermal shift assay (CETSA). These differences were explained by intracellular drug exposure (the unbound cellular partition coefficient, Kpuu). Importantly, our data is in line with previous observations regarding the discrepancy between PREP’s cleaving activity and outcomes in autophagy. This highlights the need to further explore PREP’s role in autophagy so that more effective compounds can be designed to battle diseases in which autophagy induction is needed. The present work is the first report investigating the PREP pathway in the RPE and we predict that the PREP inhibitors can be further optimized for treatment of dry AMD.

Protein kinases are crucial components of the cell-signalling machinery that orchestrate and convey messages to their downstream targets. Most often, kinases are activated upon a phosphorylation to their activation loop, which will shift the kinase into the active conformation. The Dual specificity mitogen-activated protein kinase kinase 4 (MKK4) exists in a unique conformation in its inactive unphosphorylated state, where its activation segment appears in a stable α-helical conformation. However, the precise role of this unique conformational state of MKK4 is unknown. Here, by all-atom molecular dynamics simulations (MD simulations), we show that this inactive state is unstable as monomer even when unphosphorylated and that the phosphorylation of the activation segment further destabilizes the autoinhibited α-helix. The specific phosphorylation pattern of the activation segment has also a unique influence on MKK4 dynamics. Furthermore, we observed that this specific inactive state is stable as a dimer, which becomes destabilized upon phosphorylation. Finally, we noticed that the most frequent MKK4 mutation observed in cancer, R134W, which role has not been disclosed to date, contributes to the dimer stability. Based on these data we postulate that MKK4 occurs as a dimer in its inactive autoinhibited state, providing an additional layer for its activity regulation.

Abstract Treatment options for COVID‐19, caused by SARS‐CoV‐2, remain limited. Understanding viral pathogenesis at the molecular level is critical to develop effective therapy. Some recent studies have explored SARS‐CoV‐2–host interactomes and provided great resources for understanding viral replication. However, host proteins that functionally associate with SARS‐CoV‐2 are localized in the corresponding subnetwork within the comprehensive human interactome. Therefore, constructing a downstream network including all potential viral receptors, host cell proteases, and cofactors is necessary and should be used as an additional criterion for the validation of critical host machineries used for viral processing. This study applied both affinity purification mass spectrometry (AP‐MS) and the complementary proximity‐based labeling MS method (BioID‐MS) on 29 viral ORFs and 18 host proteins with potential roles in viral replication to map the interactions relevant to viral processing. The analysis yields a list of 693 hub proteins sharing interactions with both viral baits and host baits and revealed their biological significance for SARS‐CoV‐2. Those hub proteins then served as a rational resource for drug repurposing via a virtual screening approach. The overall process resulted in the suggested repurposing of 59 compounds for 15 protein targets. Furthermore, antiviral effects of some candidate drugs were observed in vitro validation using image‐based drug screen with infectious SARS‐CoV‐2. In addition, our results suggest that the antiviral activity of methotrexate could be associated with its inhibitory effect on specific protein–protein interactions.

To better understand the functionality of organic anion transporting polypeptides (OATPs) and to design new ligands, reliable structural data of each OATP is needed. In this work, we used a combination of homology model with molecular dynamics simulations to generate a comprehensive structural dataset, that encompasses a diverse set of OATPs but also their relevant conformations. Our OATP models share a conserved transmembrane helix folding harbouring a druggable binding pocket in the shape of an inner pore. Our simulations suggest that the conserved salt bridges at the extracellular region between residues on TM1 and TM7 might influence the entrance of substrates. Interactions between residues on TM1 and TM4 within OATP1 family shown their importance in transport of substrates. Additionally, in transmembrane (TM) 1/2, a known conserved element, interact with two identified motifs in the TM7 and TM11. Our simulations suggest that TM1/2-TM7 interaction influence the inner pocket accessibility, while TM1/2-TM11 salt bridges control the substrate binding stability.

Transcription and translation are fundamental cellular processes that govern the protein production of cells. These processes are generally up regulated in cancer cells, to maintain the enhanced metabolism and proliferative state of these cells. As such cancerous cells can be susceptible to transcription and translation inhibitors. There are numerous druggable proteins involved in transcription and translation which make lucrative targets for cancer drug development. In addition to proteins, recent years have shown that the “undruggable” transcription factors and RNA molecules can also be targeted to hamper the transcription or translation in cancer. In this review, we summarize the properties and function of the transcription and translation inhibitors that have been tested and developed, focusing on the advances of the last 5 years. To complement this, we also discuss some of the recent advances in targeting oncogenes tightly controlling transcription including transcription factors and KRAS. In addition to natural and synthetic compounds, we review DNA and RNA based approaches to develop cancer drugs. Finally, we conclude with the outlook to the future of the development of transcription and translation inhibitors.

In this study, we investigated the delivery of synthetic neurosteroids into MCF-7 human breast adenocarcinoma cells via Organic Anionic Transporting Polypeptides (OATPs) (pH 7.4 and 5.5) to identify the structural components required for OATP-mediated cellular uptake and to get insight into brain drug delivery. Then, we identified structure-uptake relationships using in-house developed OATP1A2 homology model to predict binding sites and modes for the ligands. These binding modes were studied by molecular dynamics simulations to rationalize the experimental results. Our results show that carboxylic acid needs to be at least at 3 carbon-carbon bonds distance from amide bond at the C-3 position of the androstane skeleton and have an amino group to avoid efflux transport. Replacement of hydroxyl group at C-3 with any of the 3, 4, and 5-carbon chained terminal carboxylic groups improved the affinity. We attribute this to polar interactions between carboxylic acid and side-chains of Lys33 and Arg556. The additional amine group showed interactions with Glu172 and Glu200. Based on transporter capacities and efficacies, it could be speculated that the functionalization of acetyl group at the C-17 position of the steroidal skeleton might be explored further to enable OAT1A2-mediated delivery of neurosteroids into the cells and also across the blood-brain barrier.

This work aimed to construct 3D-QSAR CoMFA and CoMSIA models for a series of 31 FAAH inhibitors, containing the 1,3,4-oxadiazol-2-one moiety. The obtained models were characterized by good statistical parameters: CoMFA Q2 = 0.61, R2 = 0.98; CoMSIA Q2 = 0.64, R2 = 0.93. The CoMFA model field contributions were 54.1% and 45.9% for steric and electrostatic fields, respectively. In the CoMSIA model, electrostatic, steric, hydrogen bond donor, and hydrogen acceptor properties were equal to 34.6%, 23.9%, 23.4%, and 18.0%, respectively. These models were validated by applying the leave-one-out technique, the seven-element test set (CoMFA r2test-set = 0.91; CoMSIA r2test-set = 0.91), a progressive scrambling test, and external validation criteria developed by Golbraikh and Tropsha (CoMFA r20 = 0.98, k = 0.95; CoMSIA r20 = 0.98, k = 0.89). As the statistical significance of the obtained model was confirmed, the results of the CoMFA and CoMSIA field calculation were mapped onto the enzyme binding site. It gave us the opportunity to discuss the structure–activity relationship based on the ligand–enzyme interactions. In particular, examination of the electrostatic properties of the established CoMFA model revealed fields that correspond to the regions where electropositive substituents are not desired, e.g., in the neighborhood of the 1,3,4-oxadiazol-2-one moiety. This highlights the importance of heterocycle, a highly electronegative moiety in this area of each ligand. Examination of hydrogen bond donor and acceptor properties contour maps revealed several spots where the implementation of another hydrogen-bond-donating moiety will positively impact molecules’ binding affinity, e.g., in the neighborhood of the 1,3,4-oxadiazol-2-one ring. On the other hand, there is a large isopleth that refers to the favorable H-bond properties close to the terminal phenoxy group of a ligand, which means that, generally speaking, H-bond acceptors are desired in this area.

We report our investigations into the underlying differences between 1,2,3-dithiazole and their ultra-rare counterpart, 1,2,3-thiaselenazole. This rare 1,2,3-thiaselenazole chemotype was afforded by sulfur extrusion and selenium insertion into the preconstructed 1,2,3-dithiazoles. We built a library of matched paired compounds to compare and contrast the two ring systems. This led to the development of both narrow and broad-spectrum antimicrobial compounds with sub-micro molar potency, limited to no toxicity and a further understanding of the transition state electronics through molecular simulations. We also identified the potent 4,5,6-trichlorocyclopenta[d][1,2,3]thiaselenazole 11a, for use against Candida albicans, Cryptococcus neoformans var. grubii, Staphylococcus aureus and Acinetobacter baumannii, all of which have limited clinical treatment options. The 1,2,3-thiaselenazole represents a new class of potential compounds for the treatment of a host of multi-resistant hospital derived infections.

A mutated KRAS protein is frequently observed in human cancers. Traditionally, the oncogenic properties of KRAS missense mutants at position 12 (G12X) have been considered as equal. Here, by assessing the probabilities of occurrence of all KRAS G12X mutations and KRAS dynamics we show that this assumption does not hold true. Instead, our findings revealed an outstanding mutational bias. We conducted a thorough mutational analysis of KRAS G12X mutations and assessed to what extent the observed mutation frequencies follow a random distribution. Unique tissue-specific frequencies are displayed with specific mutations, especially with G12R, which cannot be explained by random probabilities. To clarify the underlying causes for the nonrandom probabilities, we conducted extensive atomistic molecular dynamics simulations (170 μs) to study the differences of G12X mutations on a molecular level. The simulations revealed an allosteric hydrophobic signaling network in KRAS, and that protein dynamics is altered among the G12X mutants and as such differs from the wild-type and is mutation-specific. The shift in long-timescale conformational dynamics was confirmed with Markov state modeling. A G12X mutation was found to modify KRAS dynamics in an allosteric way, which is especially manifested in the switch regions that are responsible for the effector protein binding. The findings provide a basis to understand better the oncogenic properties of KRAS G12X mutants and the consequences of the observed nonrandom frequencies of specific G12X mutations.

The human family of type II transmembrane serine proteases includes 17 members. The defining features of these proteases are an N-terminal transmembrane domain and a C-terminal serine protease of the chymotrypsin (S1) fold, separated from each other by a variable stem region. Recently accumulated evidence suggests a critical role for these proteases in development of cancer and metastatic capacity. Both the cancer relevance and the accessibility of the extracellularly oriented catalytic domain for therapeutic and imaging agents have fueled drug discovery interest in the type II class of transmembrane serine proteases. Typically, the initial hit discovery processes aim to identify molecules with verifiable activity at the drug target and with sufficient drug-like characters. We present here protocols for structure-based virtual screening of candidate ligands for transmembrane serine protease hepsin. The methods describe use of the 3D structure of the catalytic site of hepsin for molecular docking with ZINC, which is a molecular database of > 30 million purchasable compounds. Small candidate subsets were experimentally tested with demonstrable hits, which provided meaningful cues of the ligand structures for further lead development.

In order to search for novel antipsychotics acting through the D2 receptor, it is necessary to know the structure−activity relationships for dopamine D2 receptor antagonists. In this context, we constructed the universal three-dimensional quantitative structure−activity relationship (3D- QSAR) model for competitive dopamine D2 receptor antagonists. We took 176 compounds from chemically different groups characterized by the half maximal inhibitory concentration (IC50)from the CHEMBL database and docked them to the X-ray structure of the human D2 receptor in the inactive state. Selected docking poses were applied for Comparative Molecular Field Analysis (CoMFA) alignment. The obtained CoMFA model is characterized by a cross-validated coefficient Q2 of 0.76 with an optimal component of 5, R2 of 0.92, and an F value of 338.9. The steric and electrostatic field contributions are 67.4% and 32.6%, respectively. The statistics obtained prove that the CoMFA model is significant. Next, the IC50 of the 16 compounds from the test set was predicted with R2 of 0.95. Finally, a progressive scrambling test was carried out for additional validation. The CoMFA fields were mapped onto the dopamine D2 receptor binding site, which enabled a discussion of the structure−activity relationship based on ligand−receptor interactions. In particular, it was found that one of the desired steric interactions covers the area of a putative common allosteric pocket suggested for some other G protein-coupled receptors (GPCRs), which would suggest that some of the known dopamine receptor antagonists are bitopic in their essence. The CoMFA model can be applied to predict the potential activity of novel dopamine D2 receptor antagonists.

Protein kinase CK2 is an emerging target for therapeutic intervention in human diseases, particularly in cancer. Inhibitors of this enzyme are currently in clinical trials, indicating the druggability of human CK2. By virtual screening of the ZINC database, we found that the natural compound bikaverin can fit well in the ATP binding site of the target enzyme CK2. By further in vitro evaluation using CK2 holoenzyme, bikaverin turned to be a potent inhibitor with an IC50 value of 1.24 µM. In this work, the cell permeability of bikaverin was determined using a Caco-2 cell permeability assay as a prerequisite for cellular evaluation and the compound turned out to be cell permeable with a Papp- value of 4.46 × 10−6 cm/s. Bikaverin was tested for its effect on cell viability using a MTT assay and cell proliferation using an EdU assay in different cancer cell lines (MCF7, A427 and A431 cells). Cell viability and cell proliferation were reduced dramatically after treatment with 10 µM bikaverin for 24 h. Additionally the IncuCyte® live-cell imaging system was applied for monitoring the cytotoxicity of bikaverin in the three tested cancer cell lines. Finally, molecular dynamic studies were performed to clarify the ligand binding mode of bikaverin at the ATP binding site of CK2 and to identify the amino acids involved.

In 1982, Kuntz et al. published an article with the title “A Geometric Approach to Macromolecule-Ligand Interactions”, where they described a method “to explore geometrically feasible alignment of ligands and receptors of known structure”. Since then, small molecule docking has been employed as a fast way to estimate the binding pose of a given compound within a specific target protein and also to predict binding affinity. Remarkably, the first docking method suggested by Kuntz and colleagues aimed to predict binding poses but very little was specified about binding affinity. This raises the question as to whether docking is the right tool to estimate binding affinity. The short answer is no, and this has been concluded in several comprehensive analyses. However, in this opinion paper we discuss several critical aspects that need to be reconsidered before a reliable binding affinity prediction through docking is realistic. These are not the only issues that need to be considered, but they are perhaps the most critical ones. We also consider that in spite of the huge efforts to enhance scoring functions, the accuracy of binding affinity predictions is perhaps only as good as it was 10–20 years ago. There are several underlying reasons for this poor performance and these are analyzed. In particular, we focus on the role of the solvent (water), the poor description of H-bonding and the lack of the systems’ true dynamics. We hope to provide readers with potential insights and tools to overcome the challenging issues related to binding affinity prediction via docking.

Protein kinase CK2, initially designated as casein kinase 2, is an ubiquitously expressed serine/threonine kinase. This enzyme, implicated in many cellular processes, is highly expressed and active in many tumor cells. A large number of compounds has been developed as inhibitors comprising different backbones. Beside others, structures with an indeno[1,2-b]indole scaffold turned out to be potent new leads. With the aim of developing new inhibitors of human protein kinase CK2, we report here on the generation of common feature pharmacophore model to further explain the binding requirements for human CK2 inhibitors. Nine common chemical features of indeno[1,2-b]indole-type CK2 inhibitors were determined using MOE software (Chemical Computing Group, Montreal, Canada). This pharmacophore model was used for database mining with the aim to identify novel scaffolds for developing new potent and selective CK2 inhibitors. Using this strategy several structures were selected by searching inside the ZINC compound database. One of the selected compounds was bikaverin (6,11-dihydroxy-3,8-dimethoxy-1-methylbenzo[b]xanthene-7,10,12-trione), a natural compound which is produced by several kinds of fungi. This compound was tested on human recombinant CK2 and turned out to be an active inhibitor with an IC50 value of 1.24 µM.