Computational analysis of kinase inhibitor selectivity using structural knowledge

KinPred-RNA-kinase activity inference and cancer type classification using machine learning on RNA-seq data

Kinases as important enzymes can transfer phosphate groups from high-energy and phosphate-donating molecules to specific substrates and play essential roles in various cellular processes. Existing algorithms for kinase activity from phosphorylated proteomics data are often costly, requiring valuable samples. Moreover, methods to extract kinase activities from bulk RNA sequencing data remain undeveloped. In this study, we propose a computational framework KinPred-RNA to derive kinase activities from bulk RNA-sequencing data in cancer samples. KinPred-RNA framework, using the extreme gradient boosting (XGBoost) regression model, outperforms random forest regression, multiple linear regression, and support vector machine regression models in predicting kinase activities from cancer-related RNA sequencing data. Efficient gene signatures from the LINCS-L1000 dataset were used as inputs for KinPred-RNA. The results highlight its potential to be related to biological function. In conclusion, KinPred RNA constitutes a significant advance in cancer research by potentially facilitating the identification of cancer.

Evaluation of novel pyrazol-4-yl pyridine derivatives possessing arylsulfonamide tethers as c-Jun N-terminal kinase (JNK) inhibitors in leukemia cells

A series of 36 pyrazol-4-yl pyridine derivatives (8a-i, 9a-i, 10a-i, and 11a-i) was designed, synthesized, and evaluated for its antiproliferative activity over NCI-60 cancer cell line panel and inhibitory effect against JNK isoforms (JNK1, JNK2, and JNK3). All the synthesized compounds were tested against the NCI-60 cancer cell line panel. Compounds 11b, 11c, 11g, and 11i were selected to determine their GI50s and exerted a superior potency over the reference standard SP600125 against the tested cell lines. 11c showed a GI50 of 1.28 μM against K562 leukemic cells. Vero cells were used to assess 11c cytotoxicity compared to the tested cancer cells. The target compounds were tested against hJNK isoforms in which compound 11e exhibited the highest potency against JNK isoforms with IC50 values of 1.81, 12.7, and 10.5 nM against JNK1, JNK2, and JNK3, respectively. Kinase profiling of 11e showed higher JNK selectivity in 50 kinase panels. Compounds 11c and 11e showed cell population arrest at the G2/M phase, induced early apoptosis, and slightly inhibited beclin-1 production at higher concentrations in K562 leukemia cells relative to SP600125. NanoBRET assay of 11e showed intracellular JNK1 inhibition with an IC50 of 2.81 μM. Also, it inhibited CYP2D6 and 3A4 with different extent and its hERG activity showed little cardiac toxicity with an IC50 of 4.82 μM. hJNK3 was used as a template to generate the hJNK1 crystal structure to explore the binding mode of 11e (PDB ID: 8ENJ) with a resolution of 2.8 °A and showed a typical type I kinase inhibition against hJNK1. Binding energy scores showed that selectivity of 11e towards JNK1 could be attributed to additional hydrophobic interactions relative to JNK3.

Off-target pharmacological activity at various kinases: Potential functional and pathological side effects

In drug discovery, during the lead optimization and candidate characterization stages, novel small molecules are frequently evaluated in a battery of in vitro pharmacology assays to identify potential unintended, off-target interactions with various receptors, transporters, ion channels, and enzymes, including kinases. Furthermore, these screening panels may also provide utility at later stages of development to provide a mechanistic understanding of unexpected safety findings. Here, we present a compendium of the most likely functional and pathological outcomes associated with interaction(s) to a panel of 95 kinases based on an extensive curation of the scientific literature. This panel of kinases was designed by AbbVie based on safety-related data extracted from the literature, as well as from over 20 years of institutional knowledge generated from discovery efforts. For each kinase, the scientific literature was reviewed using online databases and the most often reported functional and pathological effects were summarized. This work should serve as a practical guide for small molecule drug discovery scientists and clinical investigators to predict and/or interpret adverse effects related to pharmacological interactions with these kinases.

AiKPro: deep learning model for kinome-wide bioactivity profiling using structure-based sequence alignments and molecular 3D conformer ensemble descriptors

The discovery of selective and potent kinase inhibitors is crucial for the treatment of various diseases, but the process is challenging due to the high structural similarity among kinases. Efficient kinome-wide bioactivity profiling is essential for understanding kinase function and identifying selective inhibitors. In this study, we propose AiKPro, a deep learning model that combines structure-validated multiple sequence alignments and molecular 3D conformer ensemble descriptors to predict kinase-ligand binding affinities. Our deep learning model uses an attention-based mechanism to capture complex patterns in the interactions between the kinase and the ligand. To assess the performance of AiKPro, we evaluated the impact of descriptors, the predictability for untrained kinases and compounds, and kinase activity profiling based on odd ratios. Our model, AiKPro, shows good Pearson's correlation coefficients of 0.88 and 0.87 for the test set and for the untrained sets of compounds, respectively, which also shows the robustness of the model. AiKPro shows good kinase-activity profiles across the kinome, potentially facilitating the discovery of novel interactions and selective inhibitors. Our approach holds potential implications for the discovery of novel, selective kinase inhibitors and guiding rational drug design.

Selectivity and Ranking of Tight-Binding JAK-STAT Inhibitors Using Markovian Milestoning with Voronoi Tessellations

Janus kinases (JAK), a group of proteins in the nonreceptor tyrosine kinase (NRTKs) family, play a crucial role in growth, survival, and angiogenesis. They are activated by cytokines through the Janus kinase-signal transducer and activator of a transcription (JAK-STAT) signaling pathway. JAK-STAT signaling pathways have significant roles in the regulation of cell division, apoptosis, and immunity. Identification of the V617F mutation in the Janus homology 2 (JH2) domain of JAK2 leading to myeloproliferative disorders has stimulated great interest in the drug discovery community to develop JAK2-specific inhibitors. However, such inhibitors should be selective toward JAK2 over other JAKs and display an extended residence time. Recently, novel JAK2/STAT5 axis inhibitors (N-(1H-pyrazol-3-yl)pyrimidin-2-amino derivatives) have displayed extended residence times (hours or longer) on target and adequate selectivity excluding JAK3. To facilitate a deeper understanding of the kinase-inhibitor interactions and advance the development of such inhibitors, we utilize a multiscale Markovian milestoning with Voronoi tessellations (MMVT) approach within the Simulation-Enabled Estimation of Kinetic Rates v.2 (SEEKR2) program to rank order these inhibitors based on their kinetic properties and further explain the selectivity of JAK2 inhibitors over JAK3. Our approach investigates the kinetic and thermodynamic properties of JAK-inhibitor complexes in a user-friendly, fast, efficient, and accurate manner compared to other brute force and hybrid-enhanced sampling approaches.

Computational Approaches to Enzyme Inhibition by Marine Natural Products in the Search for New Drugs

The exploration of biologically relevant chemical space for the discovery of small bioactive molecules present in marine organisms has led not only to important advances in certain therapeutic areas, but also to a better understanding of many life processes. The still largely untapped reservoir of countless metabolites that play biological roles in marine invertebrates and microorganisms opens new avenues and poses new challenges for research. Computational technologies provide the means to (i) organize chemical and biological information in easily searchable and hyperlinked databases and knowledgebases; (ii) carry out cheminformatic analyses on natural products; (iii) mine microbial genomes for known and cryptic biosynthetic pathways; (iv) explore global networks that connect active compounds to their targets (often including enzymes); (v) solve structures of ligands, targets, and their respective complexes using X-ray crystallography and NMR techniques, thus enabling virtual screening and structure-based drug design; and (vi) build molecular models to simulate ligand binding and understand mechanisms of action in atomic detail. Marine natural products are viewed today not only as potential drugs, but also as an invaluable source of chemical inspiration for the development of novel chemotypes to be used in chemical biology and medicinal chemistry research.

Target-specific compound selectivity for multi-target drug discovery and repurposing

Most drug molecules modulate multiple target proteins, leading either to therapeutic effects or unwanted side effects. Such target promiscuity partly contributes to high attrition rates and leads to wasted costs and time in the current drug discovery process, and makes the assessment of compound selectivity an important factor in drug development and repurposing efforts. Traditionally, selectivity of a compound is characterized in terms of its target activity profile (wide or narrow), which can be quantified using various statistical and information theoretic metrics. Even though the existing selectivity metrics are widely used for characterizing the overall selectivity of a compound, they fall short in quantifying how selective the compound is against a particular target protein (e.g., disease target of interest). We therefore extended the concept of compound selectivity towards target-specific selectivity, defined as the potency of a compound to bind to the particular protein in comparison to the other potential targets. We decompose the target-specific selectivity into two components: 1) the compound’s potency against the target of interest (absolute potency), and 2) the compound’s potency against the other targets (relative potency). The maximally selective compound-target pairs are then identified as a solution of a bi-objective optimization problem that simultaneously optimizes these two potency metrics. In computational experiments carried out using large-scale kinase inhibitor dataset, which represents a wide range of polypharmacological activities, we show how the optimization-based selectivity scoring offers a systematic approach to finding both potent and selective compounds against given kinase targets. Compared to the existing selectivity metrics, we show how the target-specific selectivity provides additional insights into the target selectivity and promiscuity of multi-targeting kinase inhibitors. Even though the selectivity score is shown to be relatively robust against both missing bioactivity values and the dataset size, we further developed a permutation-based procedure to calculate empirical p-values to assess the statistical significance of the observed selectivity of a compound-target pair in the given bioactivity dataset. We present several case studies that show how the target-specific selectivity can distinguish between highly selective and broadly-active kinase inhibitors, hence facilitating the discovery or repurposing of multi-targeting drugs.

An overview of kinase downregulators and recent advances in discovery approaches

Since the clinical approval of imatinib, the discovery of protein kinase downregulators entered a prosperous age. However, challenges still exist in the discovery of kinase downregulator drugs, such as the high failure rate during development, side effects, and drug-resistance problems. With the progress made through multidisciplinary efforts, an increasing number of new approaches have been applied to solve the above problems during the discovery process of kinase downregulators. In terms of in vitro and in vivo drug evaluation, progress was also made in cellular and animal model platforms for better and more clinically relevant drug assessment. Here, we review the advances in drug design strategies, drug property evaluation technologies, and efficacy evaluation models and technologies. Finally, we discuss the challenges and perspectives in the development of kinase downregulator drugs.

Active Site Sequence Representations of Human Kinases Outperform Full Sequence Representations for Affinity Prediction and Inhibitor Generation: 3D Effects in a 1D Model

Recent advances in deep learning have enabled the development of large-scale multimodal models for virtual screening and de novo molecular design. The human kinome with its abundant sequence and inhibitor data presents an attractive opportunity to develop proteochemometric models that exploit the size and internal diversity of this family of targets. Here, we challenge a standard practice in sequence-based affinity prediction models: instead of leveraging the full primary structure of proteins, each target is represented by a sequence of 29 discontiguous residues defining the ATP binding site. In kinase-ligand binding affinity prediction, our results show that the reduced active site sequence representation is not only computationally more efficient but consistently yields significantly higher performance than the full primary structure. This trend persists across different models, data sets, and performance metrics and holds true when predicting pIC₅₀ for both unseen ligands and kinases. Our interpretability analysis reveals a potential explanation for the superiority of the active site models: whereas only mild statistical effects about the extraction of three-dimensional (3D) interaction sites take place in the full sequence models, the active site models are equipped with an implicit but strong inductive bias about the 3D structure stemming from the discontiguity of the active sites. Moreover, in direct comparisons, our models perform similarly or better than previous state-of-the-art approaches in affinity prediction. We then investigate a de novo molecular design task and find that the active site provides benefits in the computational efficiency, but otherwise, both kinase representations yield similar optimized affinities (for both SMILES- and SELFIES-based molecular generators). Our work challenges the assumption that the full primary structure is indispensable for modeling human kinases.

ProfKin: A comprehensive web server for structure-based kinase profiling

Protein kinases are central mediators of signal-transduction cascades and attractive drug targets for therapeutic intervention. Since kinases are structurally and mechanistically related to each other, kinase inhibitor selectivity is often investigated by kinase profiling and considered as an important index for drug discovery. We here describe a versatile web server termed ProfKin for structure-based kinase profiling, which is based on a kinase-ligand focused database (KinLigDB). It provides all ready-to-use 3D structure coordinates of 4219 kinase-ligand complex structures covering 297 human kinases and the associated information, particularly including binding site type, binding ligand type, interaction fingerprints, downstream molecules and related human diseases. The web server works via predicting possible binding modes for the query molecule, prioritizing the binding modes guided by an interaction fingerprint analysis method, and giving a list of ranked kinases by a comprehensive index. Users can freely select entire or part of the KinLigDB database, e.g. via subfamily and binding site type, to customize the profiling contents. The superimpositions of the predicted binding poses of the query molecule with reference binding modes can be visually inspected on the website. The additional classification attributes and phylogenetic tree are also given for each top-ranked kinase.

High-throughput kinome-RNAi screen identifies protein kinase R activator (PACT) as a novel genetic modifier of CUG foci integrity in myotonic dystrophy type 1 (DM1)

Myotonic Dystrophy Type 1 (DM1) is the most common form of adult muscular dystrophy (~1:8000). In DM1, expansion of CTG trinucleotide repeats in the 3' untranslated region of the dystrophia myotonica protein kinase (DMPK) gene results in DMPK mRNA hairpin structures which aggregate as insoluble ribonuclear foci and sequester several RNA-binding proteins. The resulting sequestration and misregulation of important splicing factors, such as muscleblind-like 1 (MBNL1), causes the aberrant expression of fetal transcripts for several genes that contribute to the disease phenotype. Previous work has shown that antisense oligonucleotide-mediated disaggregation of the intranuclear foci has the potential to reverse downstream anomalies. To explore whether the nuclear foci are, to some extent, controlled by cell signalling pathways, we have performed a screen using a small interfering RNA (siRNA) library targeting 518 protein kinases to look at kinomic modulation of foci integrity. RNA foci were visualized by in situ hybridization of a fluorescent-tagged (CAG)10 probe directed towards the expanded DMPK mRNA and the cross-sectional area and number of foci per nuclei were recorded. From our screen, we have identified PACT (protein kinase R (PKR) activator) as a novel modulator of foci integrity and have shown that PACT knockdown can both increase MBNL1 protein levels; however, these changes are not suffcient for significant correction of downstream spliceopathies.

Homology modeling of TMPRSS2 yields candidate drugs that may inhibit entry of SARS-CoV-2 into human cells

The most rapid path to discovering treatment options for the novel coronavirus SARS-CoV-2 is to find existing medications that are active against the virus. We have focused on identifying repurposing candidates for the transmembrane serine protease family member II (TMPRSS2), which is critical for entry of coronaviruses into cells. Using known 3D structures of close homologs, we created seven homology models. We also identified a set of serine protease inhibitor drugs, generated several conformations of each, and docked them into our models. We used three known chemical (non-drug) inhibitors and one validated inhibitor of TMPRSS2 in MERS as benchmark compounds and found six compounds with predicted high binding affinity in the range of the known inhibitors. We also showed that a previously published weak inhibitor, Camostat, had a significantly lower binding score than our six compounds. All six compounds are anticoagulants with significant and potentially dangerous clinical effects and side effects. Nonetheless, if these compounds significantly inhibit SARS-CoV-2 infection, they could represent a potentially useful clinical tool.

Pocket similarity identifies selective estrogen receptor modulators as microtubule modulators at the taxane site

Taxanes are a family of natural products with a broad spectrum of anticancer activity. This activity is mediated by interaction with the taxane site of beta-tubulin, leading to microtubule stabilization and cell death. Although widely used in the treatment of breast cancer and other malignancies, existing taxane-based therapies including paclitaxel and the second-generation docetaxel are currently limited by severe adverse effects and dose-limiting toxicity. To discover taxane site modulators, we employ a computational binding site similarity screen of > 14,000 drug-like pockets from PDB, revealing an unexpected similarity between the estrogen receptor and the beta-tubulin taxane binding pocket. Evaluation of nine selective estrogen receptor modulators (SERMs) via cellular and biochemical assays confirms taxane site interaction, microtubule stabilization, and cell proliferation inhibition. Our study demonstrates that SERMs can modulate microtubule assembly and raises the possibility of an estrogen receptor-independent mechanism for inhibiting cell proliferation.

Taxanes are natural products which bind beta-tubulin, stabilize microtubules and have a broad spectrum of anticancer activity. Here authors employ a computational binding site similarity screen and cell-based assays to reveal a SERM cross-reactivity between the estrogen receptor and the beta-tubulin taxane binding pocket.