Gaining Insights into Key Structural Hotspots within the Allosteric Binding Pockets of Protein Kinases.

Protein kinases are essential regulators of cell function and represent one of the largest and most diverse protein families. They are particularly influential in signal transduction and coordinating complex processes like the cell cycle. Out of the 518 human protein kinases identified, 478 are part of a single superfamily sharing catalytic domains that are related in sequence. The dysregulation of protein kinases due to certain mutations has been associated with various diseases, including cancer. Although most of the protein kinase inhibitors identified as type I or type II primarily target the ATP-binding pockets of kinases, the structural and sequential resemblances among these pockets pose a significant challenge for selective inhibition. Therefore, targeting allosteric pockets that are beside highly conserved ATP pockets has emerged as a promising strategy to prevail current limitations, such as poor selectivity and drug resistance. In this article, we compared the binding pockets of various protein kinases for which allosteric (type III) inhibitors have already been developed. Additionally, understanding the structure and shape of existing ligands could aid in identifying key interaction sites within the allosteric pockets of kinases. This comprehensive review aims to facilitate the design of more effective and selective allosteric inhibitors.

An Innovative Approach to Address Neurodegenerative Diseases through Kinase-Targeted Therapies: Potential for Designing Covalent Inhibitors.

Owing to the dysregulation of protein kinase activity in various diseases such as cancer and autoimmune, cardiovascular, neurodegenerative, and inflammatory conditions, the protein kinase family has emerged as a crucial drug target in the 21st century. Notably, many kinases have been targeted to address cancer and neurodegenerative diseases using conventional ATP-mimicking kinase inhibitors. Likewise, irreversible covalent inhibitors have also been developed for different types of cancer. The application of covalent modification to target proteins has led to significant advancements in the treatment of cancer. However, while covalent drugs have significantly impacted medical treatment, their potential for neurodegenerative diseases remains largely unexplored. Neurodegenerative diseases present significant risks to brain function, leading to progressive deterioration in sensory, motor, and cognitive abilities. Alzheimer's disease (AD), Huntington's disease (HD), Parkinson's disease (PD), and multiple sclerosis (MS) are among the various examples of such disorders. Numerous research groups have already reported insights through reviews and research articles on FDA-approved covalent inhibitors, revealing their mechanisms and the specific covalent warheads that preferentially interact with particular amino acid residues in intricate detail. Hence, in this review, we aim to provide a concise summary of these critical topics. This summary endeavors to guide medicinal chemists in their quest to design covalent inhibitors for protein kinases, specifically targeting neurodegenerative diseases.

An Intriguing Purview on the Design of Macrocyclic Inhibitors for Unexplored Protein Kinases through Their Binding Site Comparison.

Kinases play an important role in regulating various intracellular signaling pathways that control cell proliferation, differentiation, survival, and other cellular processes, and their deregulation causes more than 400 diseases. Consequently, macrocyclization can be considered a noteworthy approach to developing new therapeutic agents for human diseases. Macrocyclization has emerged as an effective drug discovery strategy over the past decade to improve target selectivity and potency of small molecules. Small compounds with linear structures upon macrocyclization can lead to changes in their physicochemical and biological properties by firmly reducing conformational flexibility. A number of distinct protein kinases exhibit similar binding sites. Comparison of protein binding sites provides crucial insights for drug discovery and development. Binding site similarities are helpful in understanding polypharmacology, identifying potential off-targets, and repurposing known drugs. In this review, we focused on comparing the binding sites of those kinases for which macrocyclic inhibitors are available/studied so far. Furthermore, we calculated the volume of the binding site pocket for each targeted kinase and then compared it with the binding site pocket of the kinase for which only acyclic inhibitors were designed to date. Our review and analysis of several explored kinases might be useful in targeting new protein kinases for macrocyclic drug discovery.

Design of Novel IRAK4 Inhibitors Using Molecular Docking, Dynamics Simulation and 3D-QSAR Studies.

Treatment of several autoimmune diseases and types of cancer has been an intense area of research over the past two decades. Many signaling pathways that regulate innate and/or adaptive immunity, as well as those that induce overexpression or mutation of protein kinases, have been targeted for drug discovery. One of the serine/threonine kinases, Interleukin-1 Receptor Associated Kinase 4 (IRAK4) regulates signaling through various Toll-like receptors (TLRs) and interleukin-1 receptor (IL1R). It controls diverse cellular processes including inflammation, apoptosis, and cellular differentiation. MyD88 gain-of-function mutations or overexpression of IRAK4 has been implicated in various types of malignancies such as Waldenström macroglobulinemia, B cell lymphoma, colorectal cancer, pancreatic ductal adenocarcinoma, breast cancer, etc. Moreover, over activation of IRAK4 is also associated with several autoimmune diseases. The significant role of IRAK4 makes it an interesting target for the discovery and development of potent small molecule inhibitors. A few potent IRAK4 inhibitors such as PF-06650833, RA9 and BAY1834845 have recently entered phase I/II clinical trial studies. Nevertheless, there is still a need of selective inhibitors for the treatment of cancer and various autoimmune diseases. A great need for the same intrigued us to perform molecular modeling studies on 4,6-diaminonicotinamide derivatives as IRAK4 inhibitors. We performed molecular docking and dynamics simulation of 50 ns for one of the most active compounds of the dataset. We also carried out MM-PBSA binding free energy calculation to identify the active site residues, interactions of which are contributing to the total binding energy. The final 50 ns conformation of the most active compound was selected to perform dataset alignment in a 3D-QSAR study. Generated RF-CoMFA (q2 = 0.751, ONC = 4, r2 = 0.911) model revealed reasonable statistical results. Overall results of molecular dynamics simulation, MM-PBSA binding free energy calculation and RF-CoMFA model revealed important active site residues of IRAK4 and necessary structural properties of ligand to design more potent IRAK4 inhibitors. We designed few IRAK4 inhibitors based on these results, which possessed higher activity (predicted pIC50) than the most active compounds of the dataset selected for this study. Moreover, ADMET properties of these inhibitors revealed promising results and need to be validated using experimental studies.

Design and Synthesis of Aminopyrimidinyl Pyrazole Analogs as PLK1 Inhibitors Using Hybrid 3D-QSAR and Molecular Docking.

Cancer continues to be one of the world's most severe public health issues. Polo-like kinase 1 (PLK1) is one of the most studied members of the polo-like kinase subfamily of serine/threonine protein kinases. PLK1 is a key mitotic regulator responsible for cell cycle processes, such as mitosis initiation, bipolar mitotic spindle formation, centrosome maturation, the metaphase to anaphase transition, and mitotic exit, whose overexpression is often associated with oncogenesis. Moreover, it is also involved in DNA damage response, autophagy, cytokine signaling, and apoptosis. Due to its fundamental role in cell cycle regulation, PLK1 has been linked to various types of cancer onset and progression, such as lung, colon, prostate, ovary, breast cancer, melanoma, and AML. Hence, PLK1 is recognized as a critical therapeutic target in the treatment of various proliferative diseases. PLK1 inhibitors developed in recent years have been researched and studied through clinical trials; however, most of them have failed because of their toxicity and poor therapeutic response. To design more potent and selective PLK1 inhibitors, we performed a receptor-based hybrid 3D-QSAR study of two datasets, possessing similar common scaffolds. The developed hybrid CoMFA (q2 = 0.628, r2 = 0.905) and CoMSIA (q2 = 0.580, r2 = 0.895) models showed admissible statistical results. Comprehensive, molecular docking of one of the most active compounds from the dataset and hybrid 3D-QSAR studies revealed important active site residues of PLK1 and requisite structural characteristics of ligand to design potent PLK1 inhibitors. Based on this information, we have proposed approximately 38 PLK1 inhibitors. The newly designed PLK1 inhibitors showed higher activity (predicted pIC50) than the most active compounds of all the derivatives selected for this study. We selected and synthesized two compounds, which were ultimately found to possess good IC50 values. Our design strategy provides insight into development of potent and selective PLK1 inhibitors.

Generation of Non-Nucleotide CD73 Inhibitors Using a Molecular Docking and 3D-QSAR Approach.

Radiotherapy and chemotherapy are conventional cancer treatments. Around 60% of all patients who are diagnosed with cancer receive radio- or chemotherapy in combination with surgery during their disease. Only a few patients respond to the blockage of immune checkpoints alone, or in combination therapy, because their tumours might not be immunogenic. Under these circumstances, an increasing level of extracellular adenosine via the activation of ecto-5'-nucleotidase (CD73) and consequent adenosine receptor signalling is a typical mechanism that tumours use to evade immune surveillance. CD73 is responsible for the conversion of adenosine monophosphate to adenosine. CD73 is overexpressed in various tumour types. Hence, targetting CD73's signalling is important for the reversal of adenosine-facilitated immune suppression. In this study, we selected a potent series of the non-nucleotide small molecule inhibitors of CD73. Molecular docking studies were performed in order to examine the binding mode of the inhibitors inside the active site of CD73 and 3D-QSAR was used to study the structure-activity relationship. The obtained CoMFA (q2 = 0.844, ONC = 5, r2 = 0.947) and CoMSIA (q2 = 0.804, ONC = 4, r2 = 0.954) models showed reasonable statistical values. The 3D-QSAR contour map analysis revealed useful structural characteristics that were needed to modify non-nucleotide small molecule inhibitors. We used the structural information from the overall docking and 3D-QSAR results to design new, potent CD73 non-nucleotide inhibitors. The newly designed CD73 inhibitors exhibited higher activity (predicted pIC50) than the most active compound of all of the derivatives that were selected for this study. Further experimental studies are needed in order to validate the new CD73 inhibitors.

Molecular Modelling Studies on Pyrazole Derivatives for the Design of Potent Rearranged during Transfection Kinase Inhibitors.

RET (rearranged during transfection) kinase, one of the receptor tyrosine kinases, plays a crucial role in the development of the human nervous system. It is also involved in various cell signaling networks responsible for the normal cell division, growth, migration, and survival. Previously reported clinical studies revealed that deregulation or aberrant activation of RET signaling can cause several types of human cancer. For example, medullary thyroid carcinoma (MTC) and multiple endocrine neoplasia (MEN2A, MEN2B) occur due to sporadic mutation or germline RET mutation. A number of RET kinase inhibitors have been approved by the FDA for the treatment of cancer, such as cabozantinib, vandetanib, lenvatinib, and sorafenib. However, each of these drugs is a multikinase inhibitor. Hence, RET is an important therapeutic target for cancer drug design. In this work, we have performed various molecular modelling studies, such as molecular docking and dynamics simulation for the most active compound of the pyrazole series as RET kinase inhibitors. Furthermore, molecular mechanics Poisson-Boltzmann surface area (MM/PBSA) free energy calculation and 3-dimensional quantitative structure-activity relationship (3D-QSAR) were performed using g_mmpbsa and SYBYL-X 2.1 package. The results of this study revealed the crucial binding site residues at the active site of RET kinase and contour map analysis showed important structural characteristics for the design of new highly active inhibitors. Therefore, we have designed ten RET kinase inhibitors, which showed higher inhibitory activity than the most active compound of the series. The results of our study provide insights to design more potent and selective RET kinase inhibitors.

Molecular modeling studies of pyrrolo[2,3-d]pyrimidin-4-amine derivatives as JAK1 inhibitors based on 3D-QSAR, molecular docking, molecular dynamics (MD) and MM-PBSA calculations.

Rheumatoid Arthritis (RA) is an autoimmune disease caused by overproduction of pro-inflammatory cytokines. Janus Kinases (JAKs) mediate cytokines signaling through the Janus Kinase (JAK)/signal transducer and activator of transcription (STAT) signaling pathways. Clinical studies have shown that Janus kinase 1 (JAK1) mediated signaling plays a key role in synovial response in rheumatoid arthritis. Hence, the inhibition JAK1 is considered as an important therapeutic route for treatment of rheumatoid arthritis. In this study, we have performed three-dimensional quantitative structure-activity relationship (3 D-QSAR), molecular docking, molecular dynamics (MD) and free energy calculations on a series of pyrrolo[2,3-d]pyrimidin-4-amine JAK1 inhibitors. Molecular docking studies of the compounds 03, 13, 36 and 49 with JAK1 were performed to study the binding interactions. The binding conformations of the compounds from docking studies were selected based on binding energy and H-bond interactions and were used as initial structure for MD simulations. Using 3 D-QSAR techniques, a ligand-based comparative molecular field analysis (CoMFA) model (q2 = 0.5, r2 = 0.96) and a receptor-based CoMFA model (q2 = 0.78, r2 = 0.98) were developed. Analysis of the MD results of the most active compound (compound 49) with JAK1 showed the formation of H-bond interactions with residues Glu957, Leu959 and Gly887 and water-mediated H-bond interaction with Gly887 and His885. Based on the contour map analyses of the receptor-based CoMFA, a design strategy was proposed and was used for designing new JAK1 inhibitors. Four of the designed compounds (D57, D58, D98 and D99) showed predicted activity values (pIC50> 8.8) greater than the most active compound for JAK1. MM-PBSA based free energy calculations indicated that the designed compounds were able to form stable binding with JAK1 primarily through electrostatic interactions and van der Waal interactions. Collectively, the outcome of this study can be used to further the progress of JAK1 inhibition for the treatment of rheumatoid arthritis. Communicated by Ramaswamy H. Sarma.

Rational approach toward COVID-19 main protease inhibitors via molecular docking, molecular dynamics simulation and free energy calculation.

In the rapidly evolving coronavirus disease (COVID-19) pandemic, repurposing existing drugs and evaluating commercially available inhibitors against druggable targets of the virus could be an effective strategy to accelerate the drug discovery process. The 3C-Like proteinase (3CLpro) of the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has been identified as an important drug target due to its role in viral replication. The lack of a potent 3CLpro inhibitor and the availability of the X-ray crystal structure of 3CLpro (PDB-ID 6LU7) motivated us to perform computational studies to identify commercially available potential inhibitors. A combination of modeling studies was performed to identify potential 3CLpro inhibitors from the protease inhibitor database MEROPS ( https://www.ebi.ac.uk/merops/index.shtml ). Binding energy evaluation identified key residues for inhibitor design. We found 15 potential 3CLpro inhibitors with higher binding affinity than that of an α-ketoamide inhibitor determined via X-ray structure. Among them, saquinavir and three other investigational drugs aclarubicin, TMC-310911, and faldaprevir could be suggested as potential 3CLpro inhibitors. We recommend further experimental investigation of these compounds.

Computational study of paroxetine-like inhibitors reveals new molecular insight to inhibit GRK2 with selectivity over ROCK1.

The G-protein coupled receptor kinase 2 (GRK2) regulates the desensitization of beta-adrenergic receptors (ß-AR), and its overexpression has been implicated in heart failure. Hence, the inhibition of GRK2 is considered to be an important drug target for the treatment of heart failure. Due to the high sequence similarity of GRK2 with the A, G, and C family (AGC family) of kinases, the inhibition of GRK2 also leads to the inhibition of AGC kinases such as Rho-associated coiled-coil kinase 1 (ROCK1). Therefore, unraveling the mechanisms to selectively inhibit GRK2 poses an important challenge. We have performed molecular docking, three dimensional quantitative structure activity relationship (3D-QSAR), molecular dynamics (MD) simulation, and free energy calculations techniques on a series of 53 paroxetine-like compounds to understand the structural properties desirable for enhancing the inhibitory activity for GRK2 with selectivity over ROCK1. The formation of stable hydrogen bond interactions with the residues Phe202 and Lys220 of GRK2 seems to be important for selective inhibition of GRK2. Electropositive substituents at the piperidine ring and electronegative substituents near the amide linker between the benzene ring and pyrazole ring showed a higher inhibitory preference for GRK2 over ROCK1. This study may be used in designing more potent and selective GRK2 inhibitors for therapeutic intervention of heart failure.

3D-QSAR-aided design of potent c-Met inhibitors using molecular dynamics simulation and binding free energy calculation.

Mesenchymal-epithelial transition factor (c-Met) is a member of receptor tyrosine kinase. It involves in various cellular signaling pathways which includes proliferation, motility, migration, and invasion. Over-expression of c-Met has been reported in various cancers. Hence, it is an ideal therapeutic target for cancer. The main objective of the study is to identify crucial residues involved in the inhibition of c-Met kinase and to design a series of potent imidazo [4,5-b] pyrazine derivatives as c-Met inhibitors. Docking was used to identify important active site residues involved in the inhibition of c-Met kinase which was further validated by 100 ns of molecular dynamics simulation and free energy calculation using molecular mechanics generalized born surface area. Furthermore, binding energy decomposition identified that residues Tyr1230, Met1211, Asp1222, Tyr1159, Met1160, Val1092, Ala1108, and Leu1157 contributed favorably to the binding stability of compound 32. Receptor-guided Comparative Molecular Field Analysis (CoMFA) (q2 = 0.751, NOC = 6, r2 = 0.933) and Comparative Molecular Similarity Indices Analysis (COMSIA) (q2 = 0.744, NOC = 6, r2 = 0.950) models were generated based on the docked conformation of the most active compound 32. The robustness of these models was tested using various validation techniques and found to be predictive. The results of CoMFA and CoMSIA contour maps exposed the regions favorable to enhance the activity. Based on this information, 27 novel c-Met inhibitors were designed. These designed compounds exhibited potent activity than the most active compound of the existing dataset. Communicated by Ramaswamy H. Sarma.