Exploration of type II binding mode: A privileged approach for kinase inhibitor focused drug discovery?

Unveiling novel type 1 inhibitors for targeting LIM kinase 2 (LIMK2) for cancer therapeutics: An integrative pharmacoinformatics approach

LIMK2 is crucial in regulating actin cytoskeleton dynamics, significantly contributing to cancer cell proliferation, invasion, and metastasis. Inhibitors like LIMKi3 effectively suppress LIMK2 kinase activity by directly affecting actin polymerization and preventing the formation of structures like filopodia and lamellipodia, which are typical of motile cancer cells. By modulating these actin dynamics, LIMKi3 inhibits cancer cell migration and invasion, reducing the potential for metastasis. Thus, this study aims to explore potential anti-cancer therapeutic LIMK2 inhibitors with properties resembling LIMKi3. Henceforth, molecular docking was utilized in this study to comprehend the ATP mimetic binding mode of LIMKi3, followed by Pharmacophore-based virtual screening to identify small molecules resembling LIMKi3. In addition, molecular dynamics simulations were performed to explore the dynamic behavior of LIMK2 and potential inhibitors. Further, network analysis and binding free energy calculations were implemented to comprehensively assess the interactions between the compounds and LIMK2. In molecular docking, LIMKi3 demonstrated an ATP mimetic hinge binding mode with hydrogen bonds at Ile408. Among the screened compounds (NCI300395, ChemDiv-8020-2508, and ChemDiv-7997-0024), three displayed "ADRH" pharmacophoric features like LIMKi3, with favorable ADMET properties, higher binding affinity, and significant hydrogen bond interactions at Ile408. LIMK2-inhibitor complexes showed lower RMSD than LIMK2-LIMKi3, indicating higher equilibrium by identified compounds. Protein-drug Complexes exhibited significant inter-domain correlation in N-lobe residues of LIMK2, including conserved β3, αC, and Hinge residues. Binding free energy analysis ranked LIMK2-NCI300395 highest, followed by LIMK2-ChemDiv-7997-0024 and LIMK2-ChemDiv-8020-2508, highlighting their potential as effective LIMK2-targeting compounds. Hence, this study emphasizes LIMKi3's significance and identifies potential candidates (NCI300395, ChemDiv-7997-0024, and ChemDiv-8020-2508) for developing cancer therapeutics targeting LIMK2. These findings open avenues for further investigations into the complex interplay between cytoskeletal dynamics and cancer progression.

Advances in reversible covalent kinase inhibitors

Reversible covalent kinase inhibitors (RCKIs) are a class of novel kinase inhibitors attracting increasing attention because they simultaneously show the selectivity of covalent kinase inhibitors yet avoid permanent protein-modification-induced adverse effects. Over the last decade, RCKIs have been reported to target different kinases, including Atypical group of kinases. Currently, three RCKIs are undergoing clinical trials. Here, advances in RCKIs are reviewed to systematically summarize the characteristics of electrophilic groups, chemical scaffolds, nucleophilic residues, and binding modes. In so doing, we integrate key insights into privileged electrophiles, the distribution of nucleophiles, and hence effective design strategies for the development of RCKIs. Finally, we provide a further perspective on future design strategies for RCKIs, including those that target proteins other than kinases.

Rational Design and Synthesis of a Novel Series of Thiosemicarbazone-Containing Quinazoline Derivatives as Potential VEGFR2 Inhibitors

Background/Objectives: Angiogenesis plays a crucial role in tumor development and is a driving force for the aggressiveness of several types of cancer. Our team developed a novel series of thiosemicarbazone-containing quinazoline derivatives, TSC1-TSC10, as potential VEGFR2 inhibitors with proven anti-angiogenic and antiproliferative potential. Methods: The TSC1-TSC10 series was synthesized and characterized by spectral data. Extensive methodology was applied both in vitro (Alamar Blue assay, Scratch assay, CAM assay, and VEGFR2 kinase assay) and in silico (docking studies, MDs, and MM-PBSA) for the confirmation of the biological potential. Results: TSC10 emerged as the most promising compound, with a favorable cytotoxic potential across the cell panel (Ea.Hy296, HaCaT, and A375) in agreement with the in vitro VEGFR2 kinase assay (IC50 = 119 nM). A comparable motility reduction in the vascular endothelial cells to that of the reference drug sorafenib was provided by TSC10, with a similar anti-angiogenic potential in the more complex in ovo model of the CAM assay. The in silico experiments confirmed the successful accommodation of the active site of the kinase domain similar to sorafenib for the entire TSC1-TSC10 series, providing valuable key insight into the complex stability driving force for the evaluated compounds. Conclusions: The in vitro evaluations of the biological potential correlated with the in silico predictions by computer-aided complex simulations provided a solid confirmation of the initial hypothesis for the TSC1-TSC10 series.

Type-II kinase inhibitors that target Parkinson's Disease-associated LRRK2

Aberrant increases in kinase activity of leucine-rich repeat kinase 2 (LRRK2) are associated with Parkinson's disease (PD). Numerous LRRK2-selective type-I kinase inhibitors have been developed and some have entered clinical trials. In this study, we present the first LRRK2-selective type-II kinase inhibitors. Targeting the inactive conformation of LRRK2 is functionally distinct from targeting the active-like conformation using type-I inhibitors. We designed these inhibitors using a combinatorial chemistry approach fusing selective LRRK2 type-I and promiscuous type-II inhibitors by iterative cycles of synthesis supported by structural biology and activity testing. Our current lead structures are selective and potent LRRK2 inhibitors. Through cellular assays, cryo-electron microscopy structural analysis, and in vitro motility assays, we show that our inhibitors stabilize the open, inactive kinase conformation. These new conformation-specific compounds will be invaluable as tools to study LRRK2's function and regulation, and expand the potential therapeutic options for PD.

AXL signaling in cancer: from molecular insights to targeted therapies

AXL, a member of the TAM receptor family, has emerged as a potential target for advanced-stage human malignancies. It is frequently overexpressed in different cancers and plays a significant role in various tumor-promoting pathways, including cancer cell proliferation, invasion, metastasis, epithelial-mesenchymal transition (EMT), angiogenesis, stemness, DNA damage response, acquired therapeutic resistance, immunosuppression, and inflammatory responses. Beyond oncology, AXL also facilitates viral infections, including SARS-CoV-2 and Zika highlighting its importance in both cancer and virology. In preclinical models, small-molecule kinase inhibitors targeting AXL have shown promising anti-tumorigenic potential. This review primarily focuses on the induction, regulation and biological functions of AXL in mediating these tumor-promoting pathways. We discuss a range of therapeutic strategies, including recently developed small-molecule tyrosine kinase inhibitors (TKIs), monoclonal antibodies, and antibody-drug conjugates (ADCs), anti-AXL-CAR, and combination therapies. These interventions are being examined in both preclinical and clinical studies, offering the potential for improved drug sensitivity and therapeutic efficacy. We further discuss the mechanisms of acquired therapeutic resistance, particularly the crosstalk between AXL and other critical receptor tyrosine kinases (RTKs) such as c-MET, EGFR, HER2/HER3, VEGFR, PDGFR, and FLT3. Finally, we highlight key research areas that require further exploration to enhance AXL-mediated therapeutic approaches for improved clinical outcomes.

Crystallographic mining driven computer-guided approach to identify the ASK1 inhibitor likely to perturb the catalytic region

The pathological levels of reactive oxygen species (ROS) and oxidative stress has been recognized as a critical driver for inflammatory disorders. Apoptosis signal-regulating kinase 1 (ASK1) has been reported to be activated by intracellular ROS and its inhibition leads to a down regulation of p38-and JNK-dependent signaling. ASK1 inhibitors are reported to have the potential to treat clinically important inflammatory pathologies including liver, pulmonary and renal disorders. In view of its biological and pathological significance, inhibition of ASK1 with small molecules has been pursued as an attractive strategy to combat human diseases such as non-alcoholic steatohepatitis (NASH). Despite several ASK1 inhibitors being developed, the failure in Phase 3 clinical trials of most advanced candidate selonsertib's, underscores to discover therapeutic agents with diverse chemical moiety. Here, by using structural pharmacophore and enumeration strategy on mining co-crystals of ASK1, different scaffolds were generated to enhance the chemical diversity keeping the critical molecular interaction in the catalytic site intact. A total of 15,772 compounds were generated from diverse chemical scaffolds and were evaluated using a virtual screening pipeline. Based on docking and MM-GBSA scores, a lead candidate, S3C-1-D424 was identified from top hits. A comparative molecular dynamics simulations (MD) of APO, Selonsertib and shortlisted potential candidates combined with pharmacokinetics profiling and thermodynamic analysis, demonstrating their suitability as potential ASK1 inhibitors to explore further for establishment towards hit-to-lead campaign.Communicated by Ramaswamy H. Sarma.

First-in-class ultralong-target-residence-time p38α inhibitors as a mitosis-targeted therapy for colorectal cancer

Colorectal cancer (CRC) constitutes the second leading cause of cancer-related death worldwide and advanced CRCs are resistant to targeted therapies, chemotherapies and immunotherapies. p38α (Mapk14) has been suggested as a therapeutic target in CRC; however, available p38α inhibitors only allow for insufficient target inhibition. Here we describe a unique class of p38α inhibitors with ultralong target residence times (designated ULTR-p38i) that robustly inhibit p38α downstream signaling and induce distinct biological phenotypes. ULTR-p38i monotherapy triggers an uncontrolled mitotic entry by activating Cdc25 and simultaneously blocking Wee1. Consequently, CRC cells undergo mitotic catastrophe, resulting in apoptosis or senescence. ULTR-p38i exhibit high selectivity, good pharmaco-kinetic properties and no measurable toxicity with strong therapeutic effects in patient-derived CRC organoids and syngeneic CRC mouse models. Conceptually, our study suggests ultralong-target-residence-time kinase inhibitors as an alternative to covalent inhibitors, which, because of the lack of cysteine residues, cannot be generated for many kinase cancer targets.

Harnessing allosteric inhibition: prioritizing LIMK2 inhibitors for targeted cancer therapy through pharmacophore-based virtual screening and essential molecular dynamics

The therapeutic potential of small molecule kinase inhibitors in cancer treatment is well recognized. However, achieving selectivity remains a formidable challenge, primarily due to the structural similarity of ATP binding pockets among kinases. Allosteric inhibition, which involves targeting binding pockets beyond the ATP-binding site, provides a promising alternative to overcome this challenge. In this study, a meticulous approach was implemented to prioritize type 3 inhibitors for LIMK2, employing a range of techniques including Molecular Dynamics (MD) simulations, e-pharmacophore-guided High Throughput Virtual Screening (HTVS), MM/GBSA and ADMETox analyses, Density Functional Theory (DFT) calculations, and MM/PBSA investigations. The e-pharmacophore model identifies a hypothesis featuring five essential pharmacophoric elements (RRRAH). Through virtual screening of the ZINC compound database, we identified only five compounds that align with all four pharmacophoric features: ZINC1044382792, ZINC1433610865, ZINC1044109145, ZINC952869440, and ZINC490621334. These compounds not only exhibit higher binding affinity but also demonstrate favorable ADME/Tox profiles. Molecular dynamics simulations underscore the stability of hydrogen bond interactions with critical cryptic LIMK2 pocket residues, Asp469 and Arg474, only for two compounds: ZINC143361086 and ZINC1044382792. These compounds also exhibit superior occupancy interactions, as indicated by HOMO-LUMO analysis. Additionally, binding free energy calculations highlight the significant affinities of these two compounds when complexed with LIMK2: -83.491 ± 1.230 kJ/mol and -90.122 ± 1.248 kJ/mol for ZINC1044382792 and ZINC1433610862, respectively. Hence, this comprehensive investigation identifies ZINC1433610862 and ZINC1044382792 as prospective hits, representing promising leads for targeting LIMK2 in cancer therapeutics.Communicated by Ramaswamy H. Sarma.

Improving the Reliability of Language Model-Predicted Structures as Docking Targets through Geometric Graph Learning

Applying artificial intelligence techniques to flexibly model the binding between the ligand and protein has attracted extensive interest in recent years, but their applicability remains improved. In this study, we have developed CarsiDock-Flex, a novel two-step flexible docking paradigm that generates binding poses directly from predicted structures. CarsiDock-Flex consists of an equivariant deep learning-based model termed CarsiInduce to refine ESMFold-predicted protein pockets with the induction of specific ligands and our existing CarsiDock algorithm to redock the ligand into the induced binding pockets. Extensive evaluations demonstrate the effectiveness of CarsiInduce, which can successfully guide the transition of ESMFold-predicted pockets into their holo-like conformations for numerous cases, thus leading to the superior docking accuracy of CarsiDock-Flex even on unseen sequences. Overall, our approach offers a novel design for flexible modeling of protein-ligand binding poses, paving the way for a deeper understanding of protein-ligand interactions that account for protein flexibility.

Discovery of a potent anticancer agent against pancreatic ductal adenocarcinoma targeting FAK with DFG-out state and JAK/Aurora kinases

Pancreatic ductal adenocarcinoma (PDAC) is a clinically challenging cancer because of the difficulty in diagnosis and its resistance to chemotherapy. Focal adhesion kinase (FAK) is found overexpressed in PDAC, and targeting FAK has been proved to impede the progress of PDAC. However, most of FAK inhibitors were reported to bind with FAK in a DFG-in conformation, leading to a limited anti-tumor effect in clinical studies. Herein, to develop FAK inhibitors targeting the inactive DFG-out conformation, a series of large aromatic rings were selected to improve the interaction with Phe565 of the DFG motif. Compound 26 was designed to effectively inhibit FAK and the proliferation of PANC-1 cells with IC50 of 50.94 nM and 0.15 μM, respectively. Besides, compound 26 was proved to strongly suppress the proliferation, colony formation, migration, and invasion in FAK-overexpressing PDAC cells. This inhibitor was confirmed to induce the apoptosis and G2/M arrest in PANC-1 cells through the suppression of FAK/PI3K/Akt signal pathway. Meanwhile, compound 26 was found to simultaneously inhibit FAK with DFG-out conformation and JAK3/Aurora B (IC50 of 9.99 nM and 0.49 nM, respectively). In vivo, compound 26 effectively inhibited the tumorigenesis and metastasis of PDAC with desirable biosafety. Overall, these results suggested that compound 26 was a promising candidate for the treatment of PDAC.

New thiazolidin-4-ones as anti-cervical cancer agents targeting EGFR: design, synthesis, and computational studies

AIM

A new series of 3,4-dihydronaphthalen-1(2 h)-ylidene)hydrazineylidene)-5-substituted thiazolidin-4-one derivatives were designed and synthesized.

RESULTS & METHODOLOGY

The new compounds were screened for in vitro antitumor activity against Hela cancer cell line. The compounds 7b, 7 h, and 7i produced more potent cytotoxicity than doxorubicin with IC50 values of 1.83 ± 0.1, 2.54 ± 0.14, 2.75 ± 0.15, and 3.63 ± 0.2 μM, respectively. They also showed a promising safety profile against WI-38 normal cells. In addition, compound 7b produced a promising multi-kinase inhibition against EGFR (WT) while being very selective toward the mutant forms (L858R and T790M) with IC50 values of 0.099 ± 0.006, 0.064 ± 0.006, and 0.026 ± 0.007 μM, respectively, in comparison to gefitinib and osimertinib. A study of the cell cycle in Hela cells showed that 7b arrests cell cycle in the pre-G1 phase and causes early and late apoptosis. Eventually, the molecular docking results showed that 7b had good-binding interactions with EGFRWT, EGFRL858R, and EGFRT790M.

CONCLUSION

Compound 7b was predicted to have promising oral absorption, good drug-likeness, and low toxicity risks in humans. Moreover, MD simulations confirmed the stable complexes of 7b with EGFRWT, EGFRL858R, and EGFRT790M (with RMSD 0.12-0.35 nm, RMSF 0.2-0.55 nm, SASA 140-150, and Rg 1.80-2.00 nm).

Exploring Extended Warheads toward Developing Cysteine-Targeted Covalent Kinase Inhibitors

In designing covalent kinase inhibitors (CKIs), the inclusion of electrophiles as attacking warheads demands careful choreography, ensuring not only their presence on the scaffold moiety but also their precise interaction with nucleophiles in the binding sites. Given the limited number of known electrophiles, exploring adjacent chemical space to broaden the palette of available electrophiles capable of covalent inhibition is desirable. Here, we systematically analyze the characteristics of warheads and the corresponding adjacent fragments for use in CKI design. We first collect all the released cysteine-targeted CKIs from multiple databases and create one CKI data set containing 16,961 kinase-inhibitor data points from 12,381 unique CKIs covering 146 kinases with accessible cysteines in their binding pockets. Then, we analyze this data set, focusing on the extended warheads (i.e., warheads + adjacent fragments)─including 30 common warheads and 1344 unique adjacent fragments. In so doing, we provide structural insights and delineate chemical properties and patterns in these extended warheads. Notably, we highlight the popular patterns observed within reversible CKIs for the popular warheads cyanoacrylamide and aldehyde. This study provides medicinal chemists with novel insights into extended warheads and a comprehensive source of adjacent fragments, thus guiding the design, synthesis, and optimization of CKIs.

Tilting the Scales toward EGFR Mutant Selectivity: Expanding the Scope of Bivalent "Type V" Kinase Inhibitors

Binding multiple sites within proteins with bivalent compounds is a strategy for developing uniquely active agents. A new class of dual-site inhibitors has emerged targeting the epidermal growth factor receptor (EGFR) anchored to both the orthosteric (ATP) and allosteric sites. Despite proof-of-concept successes, enabling selectivity against oncogenic activating mutations has not been achieved and classifying these inhibitors among kinase inhibitors remains underexplored. This study investigates the structure-activity relationships, binding modes, and biological activity of ATP-allosteric bivalent inhibitors (AABIs). We find that AABIs selectively inhibit drug-resistant EGFR mutants (L858R/T790M and L858R/T790M/C797S) by anchoring a methyl isoindolinone moiety along the αC-helix channel of the allosteric site. In contrast, related Type I1/2 inhibitors target wild-type EGFR but are less effective against resistant mutants. This shift in selectivity demonstrates that mutant-selective AABIs classify as "Type V" bivalent inhibitors.

FragGen: towards 3D geometry reliable fragment-based molecular generation

3D structure-based molecular generation is a successful application of generative AI in drug discovery. Most earlier models follow an atom-wise paradigm, generating molecules with good docking scores but poor molecular properties (like synthesizability and drugability). In contrast, fragment-wise generation offers a promising alternative by assembling chemically viable fragments. However, the co-design of plausible chemical and geometrical structures is still challenging, as evidenced by existing models. To address this, we introduce the Deep Geometry Handling protocol, which decomposes the entire geometry into multiple sets of geometric variables, looking beyond model architecture design. Drawing from a newly defined six-category taxonomy, we propose FragGen, a novel hybrid strategy as the first geometry-reliable, fragment-wise molecular generation method. FragGen significantly enhances both the geometric quality and synthesizability of the generated molecules, overcoming major limitations of previous models. Moreover, FragGen has been successfully applied in real-world scenarios, notably in designing type II kinase inhibitors at the ∼nM level, establishing it as the first validated 3D fragment-based drug design algorithm. We believe that this concept-algorithm-application cycle will not only inspire researchers working on other geometry-centric tasks to move beyond architecture designs but also provide a solid example of how generative AI can be customized for drug design.

The evolution of BRAF-targeted therapies in melanoma: overcoming hurdles and unleashing novel strategies

Melanoma, a highly aggressive form of skin cancer, poses a significant global health burden, with 331,647 new cases and 58,645 deaths reported in 2022. The development of melanoma is influenced by various factors, including sunlight exposure and BRAFV600 mutations that activate the MAPK/ERK pathway. The introduction of BRAF and MEK inhibitors has revolutionized the treatment landscape for melanoma patients. However, innate and acquired therapeutic resistance remains a significant challenge. This review provides a comprehensive overview of the current state of BRAF-targeted therapies in melanoma, highlighting the efficacy and limitations of FDA-approved combinations of BRAF and MEK inhibitors such as vemurafenib, dabrafenib, trametinib, and cobimetinib. The review also explores the off-target effects of BRAF inhibitors on endothelial cells, emphasizing the need for more selective therapies to minimize vascular complications and metastatic potential. The article also discusses potential druggable targets, including ERK5, CD73, ALDH1A1, PLA1A, and DMKN, which are promising in addressing diagnostic hurdles and guiding personalized therapeutic decisions. Recent studies on regorafenib, ERK5 signaling, and CD73 inhibition are highlighted as novel strategies to overcome resistance and improve treatment outcomes. The review also delves into the role of advanced therapeutic tools, such as mRNA vaccines and CRISPR-Cas9, in revolutionizing personalized oncology by targeting specific genetic mutations and enhancing immune responses against melanoma. The ongoing synergy between advancing research, targeted interventions, strategic treatment combinations, and cost-effectiveness evaluations offers a promising pathway to elevate patient outcomes in the persistent battle against melanoma significantly.

Design, synthesis and biological evaluation of thieno[3,2-c]pyrazol-urea derivatives as potent glycogen synthase kinase 3β inhibitors based on the DFG-out conformation

Glycogen synthase kinase 3β (GSK-3β) is a potential therapeutic target for the treatment of a variety of human diseases. Here, we report the design and synthesis of a series of thieno[3,2-c]pyrazol-urea derivatives and evaluation of their GSK-3β inhibitory activity. Among these analogues, the compound without substitution on terminal phenyl ring (3a) was found to be the most potent GSK-3β inhibitor with an IC50 of 74.4 nM, while substitution on the terminal phenyl (3b-3p) led to decreased potency, independent of the position, size, or electronic properties of the substituents. Kinase selectivity assay revealed that 3a showed good selectivity over a panel of kinases, but was less selective over CDK1, CDK2 and CDK5. Additionally, the pharmacological properties of the synthesized compounds were investigated computationally by the SwissADME and the results showed that most of the compounds have good ADME profiles.

AXL as immune regulator and therapeutic target in Acute Myeloid Leukemia: from current progress to novel strategies

Until recently, treatment options for patients diagnosed with Acute Myeloid Leukemia (AML) were limited and predominantly relied on various combinations, dosages, or schedules of traditional chemotherapeutic agents. Patients with advanced age, relapsed/refractory disease or comorbidities were often left without effective treatment options. Novel advances in the understanding of leukemogenesis at the molecular and genetic levels, alongside recent progress in drug development, have resulted in the emergence of novel therapeutic agents and strategies for AML patients. Among these innovations, the receptor tyrosine kinase AXL has been established as a promising therapeutic target for AML. AXL is a key regulator of several cellular functions, including epithelial-to-mesenchymal transition in tumor cells, immune regulation, apoptosis, angiogenesis and the development of chemoresistance. Clinical studies of AXL inhibitors, as single agents and in combination therapy, have demonstrated promising efficacy in treating AML. Additionally, novel AXL-targeted therapies, such as AXL-specific antibodies or antibody fragments, present potential solutions to overcome the limitations associated with traditional small-molecule AXL inhibitors or multikinase inhibitors. This review provides a comprehensive overview of the structure and biological functions of AXL under normal physiological conditions, including its role in immune regulation. We also summarize AXL's involvement in cancer, with a specific emphasis on its role in the pathogenesis of AML, its contribution to immune evasion and drug resistance. Moreover, we discuss the AXL inhibitors currently undergoing (pre)clinical evaluation for the treatment of AML.

Recent developments in receptor tyrosine kinase inhibitors: A promising mainstay in targeted cancer therapy

During the past two decades, significant advances have been made in the discovery and development of targeted inhibitors aimed at improving the survival rates of cancer patients. Among the multitude of potential therapeutic targets identified thus far, Receptor Tyrosine Kinases (RTKs) are of particular importance. Dysregulation of RTKs has been implicated in numerous human diseases, particularly cancer, where aberrant signaling pathways contribute to disease progression. RTKs have a profound impact on intra and intercellular communication, and they also facilitate post-translational modifications, notably phosphorylation, which intricately regulates a multitude of cellular processes. Prolonged phosphorylation or the disruption of kinase regulation may lead to significant alterations in cell signaling. The emergence of small molecule kinase inhibitors has revolutionized cancer therapy by offering a targeted and strategic approach that surpasses the efficacy of traditional chemotherapeutic drugs. Over the last two decades, a plethora of targeted inhibitors have been identified or engineered and have undergone clinical evaluation to enhance the survival rates of cancer patients. In this review, we have compared the expression of different RTKs, including Met, KDR/VEGFR2, EGFR, BRAF, BCR, and ALK across different cancer types in TCGA samples. Additionally, we have summarized the recent development of small molecule inhibitors and their potential in treating various malignancies. Lastly, we have discussed the mechanisms of acquired therapeutic resistance with a focus on kinase inhibitors in EGFR mutant and ALK-rearranged non-small cell lung cancer and BCR-ABL positive chronic myeloid leukemia.

Discovery of Pyrazolopyrazines as Selective, Potent, and Mutant-Active MET Inhibitors with Intracranial Efficacy

Mesenchymal-epithelial transition factor (MET) is a receptor tyrosine kinase that serves a critical function in numerous developmental, morphogenic, and proliferative signaling pathways. If dysregulated, MET has been shown to be involved in the development and survival of several cancers, including non-small cell lung cancer (NSCLC), renal cancer, and other epithelial tumors. Currently, the clinical efficacy of FDA approved MET inhibitors is limited by on-target acquired resistance, dose-limiting toxicities, and less than optimal efficacy against brain metastasis. Therefore, there is still an unmet medical need for the development of MET inhibitors to address these issues. Herein we report the application of structure-based design for the discovery and development of a novel class of brain-penetrant MET inhibitors with enhanced activity against clinically relevant mutations and improved selectivity. Compound 13 with a MET D1228N cell line IC50 value of 23 nM showed good efficacy in an intracranial tumor model and increased the median overall survival of the animals to 100% when dosed orally at 100 mg/kg daily for 21 days.

Novel therapeutic agents in clinical trials: emerging approaches in cancer therapy

Novel therapeutic agents in clinical trials offer a paradigm shift in the approach to battling this prevalent and destructive disease, and the area of cancer therapy is on the precipice of a trans formative revolution. Despite the importance of tried-and-true cancer treatments like surgery, radiation, and chemotherapy, the disease continues to evolve and adapt, making new, more potent methods necessary. The field of cancer therapy is currently witnessing the emergence of a wide range of innovative approaches. Immunotherapy, including checkpoint inhibitors, CAR-T cell treatment, and cancer vaccines, utilizes the host's immune system to selectively target and eradicate malignant cells while minimizing harm to normal tissue. The development of targeted medicines like kinase inhibitors and monoclonal antibodies has allowed for more targeted and less harmful approaches to treating cancer. With the help of genomics and molecular profiling, "precision medicine" customizes therapies to each patient's unique genetic makeup to maximize therapeutic efficacy while minimizing unwanted side effects. Epigenetic therapies, metabolic interventions, radio-pharmaceuticals, and an increasing emphasis on combination therapy with synergistic effects further broaden the therapeutic landscape. Multiple-stage clinical trials are essential for determining the safety and efficacy of these novel drugs, allowing patients to gain access to novel treatments while also furthering scientific understanding. The future of cancer therapy is rife with promise, as the integration of artificial intelligence and big data has the potential to revolutionize early detection and prevention. Collaboration among researchers, and healthcare providers, and the active involvement of patients remain the bedrock of the ongoing battle against cancer. In conclusion, the dynamic and evolving landscape of cancer therapy provides hope for improved treatment outcomes, emphasizing a patient-centered, data-driven, and ethically grounded approach as we collectively strive towards a cancer-free world.

A comprehensive exploration of the druggable conformational space of protein kinases using AI-predicted structures

Protein kinase function and interactions with drugs are controlled in part by the movement of the DFG and ɑC-Helix motifs that are related to the catalytic activity of the kinase. Small molecule ligands elicit therapeutic effects with distinct selectivity profiles and residence times that often depend on the active or inactive kinase conformation(s) they bind. Modern AI-based structural modeling methods have the potential to expand upon the limited availability of experimentally determined kinase structures in inactive states. Here, we first explored the conformational space of kinases in the PDB and models generated by AlphaFold2 (AF2) and ESMFold, two prominent AI-based protein structure prediction methods. Our investigation of AF2's ability to explore the conformational diversity of the kinome at various multiple sequence alignment (MSA) depths showed a bias within the predicted structures of kinases in DFG-in conformations, particularly those controlled by the DFG motif, based on their overabundance in the PDB. We demonstrate that predicting kinase structures using AF2 at lower MSA depths explored these alternative conformations more extensively, including identifying previously unobserved conformations for 398 kinases. Ligand enrichment analyses for 23 kinases showed that, on average, docked models distinguished between active molecules and decoys better than random (average AUC (avgAUC) of 64.58), but select models perform well (e.g., avgAUCs for PTK2 and JAK2 were 79.28 and 80.16, respectively). Further analysis explained the ligand enrichment discrepancy between low- and high-performing kinase models as binding site occlusions that would preclude docking. The overall results of our analyses suggested that, although AF2 explored previously uncharted regions of the kinase conformational space and select models exhibited enrichment scores suitable for rational drug discovery, rigorous refinement of AF2 models is likely still necessary for drug discovery campaigns.

Overcoming Secondary Mutations of Type II Kinase Inhibitors

Type II kinase inhibitors bind in the "DFG-out" kinase conformation and are generally considered to be more potent and selective than type I inhibitors, which target a DFG-in conformation. Nine type II inhibitors are currently clinically approved, with more undergoing clinical development. Resistance-conferring secondary mutations emerged with the first series of type II inhibitors, most commonly at residues within the kinase activation loop and at the "gatekeeper" position. Recently, new inhibitors have been developed to overcome such mutations; however, mutations activating other pathways (and/or other targets) have subsequently emerged on occasion. Here, we systematically summarize the secondary mutations that confer resistance to type II inhibitors, the structural basis for resistance, newer inhibitors designed to overcome resistance, as well as the challenges and opportunities for the development of new inhibitors to overcome secondary kinase domain mutations.

Small molecule anticancer drugs approved during 2021-2022: Synthesis and clinical applications

Drugs have structural homology across similar biological targets. Small molecule drugs have the efficacy to target specific molecular targets within the cancer cells with enhanced cell membrane permeability, oral administration, selectivity, and specific affinity. The objective of this review is to highlight the clinical importance and synthetic routes of new small molecule oncology drugs approved by the FDA during the period 2021-2022. These marketed drugs are listed based on the month and year of approval in chronological order. We believed that an in-depth insight into the synthetic approaches for the construction of these chemical entities would enhance the ability to develop new drugs more efficiently.

Vorolanib, sunitinib, and axitinib: A comparative study of vascular endothelial growth factor receptor inhibitors and their anti-angiogenic effects

PURPOSE

Pathological angiogenesis and vascular instability are observed in diabetic retinopathy (DR), diabetic macular edema (DME), and wet age-related macular degeneration (wAMD). Many receptor tyrosine kinases (RTKs) including vascular endothelial growth factor receptors (VEGFRs) contribute to angiogenesis, whereas the RTK TIE2 is important for vascular stability. Pan-VEGFR tyrosine kinase inhibitors (TKIs) such as vorolanib, sunitinib, and axitinib are of therapeutic interest over current antibody treatments that target only one or two ligands. This study compared the anti-angiogenic potential of these TKIs.

METHODS

A kinase HotSpot™ assay was conducted to identify TKIs inhibiting RTKs associated with angiogenesis and vascular stability. Half-maximal inhibitory concentration (IC50) for VEGFRs and TIE2 was determined for each TKI. In vitro angiogenesis inhibition was investigated using a human umbilical vein endothelial cell sprouting assay, and in vivo angiogenesis was studied using the chorioallantoic membrane assay. Melanin binding was assessed using a melanin-binding assay. Computer modeling was conducted to understand the TIE2-axitinib complex as well as interactions between vorolanib and VEGFRs.

RESULTS

Vorolanib, sunitinib, and axitinib inhibited RTKs of interest in angiogenesis and exhibited pan-VEGFR inhibition. HotSpot™ assay and TIE2 IC50 values showed that only axitinib potently inhibited TIE2 (up to 89%). All three TKIs effectively inhibited angiogenesis in vitro. In vivo, TKIs were more effective at inhibiting VEGF-induced angiogenesis than the anti-VEGF antibody bevacizumab. Of the three TKIs, only sunitinib bound melanin. TKIs differ in their classification and binding to VEGFRs, which is important because type II inhibitors have greater selectivity than type I TKIs.

CONCLUSIONS

Vorolanib, sunitinib, and axitinib exhibited pan-VEGFR inhibition and inhibited RTKs associated with pathological angiogenesis. Of the three TKIs, only axitinib potently inhibited TIE2 which is an undesired trait as TIE2 is essential for vascular stability. The findings support the use of vorolanib for therapeutic inhibition of angiogenesis observed in DR, DME, and wAMD.

New apoptotic anti-triple-negative breast cancer theobromine derivative inhibiting EGFRWT and EGFRT790M: in silico and in vitro evaluation

A new theobromine-derived EGFR inhibitor (2-(3,7-Dimethyl-2,6-dioxo-2,3,6,7-tetrahydro-1H-purin-1-yl)-N-(2,6-dimethylphenyl)acetamide) has been developed that has the essential structural characteristics to interact with EGFR's pocket. The designed compound is 2,6-di ortho methylphenyl)acetamide derivative of the well-known alkaloid, theobromine, (T-1-DOMPA). Firstly, deep DFT studies have been conducted to study the optimized chemical structure, molecular orbital and chemical reactivity analysis of T-1-DOMPA. Then, T-1-DOMPA's anticancer potentialities were estimated first through a structure-based computational approach. Utilizing molecular docking, molecular dynamics, MD, simulations over 100 ns, MM-PBSA and PLIP studies, T-1-DOMPA bonded to and inhibited the EGFR protein effectively. Subsequently, the ADMET profiles of T-1-DOMPA were computed before preparation, and its drug-likeness was anticipated. Therefore, T-1-DOMPA was prepared for the purposes of scrutinizing both the design and the results obtained in silico. The in vitro potential of T-1-DOMPA against triple-negative breast cancer cell lines, MDA- MB-231, was very promising with an IC50 value of1.8 µM, comparable to the reference drug (0.9 µM), and a much higher selectivity index of 2.6. Interestingly, T-1-DOMPA inhibited three other cancer cell lines (CaCO-2, HepG-2, and A549) with IC50 values of 1.98, 2.53, and 2.39 µM exhibiting selectivity index values of 2,4, 1.9, and 2, respectively. Additionally, T-1-DOMPA prevented effectively the MDA-MB-231cell line's healing and migration abilities. Also, T-1-DOMPA's abilities to induce apoptosis were confirmed by acridine orange/ethidium bromide (AO/EB) staining assay. Finally, T-1-DOMPA caused an up-regulation of the gene expression of the apoptotic gene, Caspase-3, in the treated MDA-MB-231cell.

Discovery of novel pyrazole based Urea/Thiourea derivatives as multiple targeting VEGFR-2, EGFRWT, EGFRT790M tyrosine kinases and COX-2 Inhibitors, with anti-cancer and anti-inflammatory activities

A novel series of pyrazole derivatives with urea/thiourea scaffolds 16a-l as hybrid sorafenib/erlotinib/celecoxib analogs was designed, synthesized and tested for its VEGFR-2, EGFRWT, EGFRT790M tyrosine kinases and COX-2, pro-inflammatory cytokines TNF-α and IL-6 inhibitory activities. All the tested compounds showed excellent COX-2 selectivity index in range of 18.04-47.87 compared to celecoxib (S.I. = 26.17) and TNF-α and IL-6 inhibitory activities (IC50 = 5.0-7.50, 6.23-8.93 respectively, compared to celecoxib IC50 = 8.40 and 8.50, respectively). Screening was carried out against 60 human cancer cell lines by National Cancer Institute (NCI), compounds 16a, 16c, 16d and 16 g were the most potent inhibitors with GI% ranges of 100 %, 99.63-87.02 %, 98.98-43.10 % and 98.68-23.62 % respectively, and with GI50 values of 1.76-15.50 µM, 1.60-5.38 µM, 1.68-7.39 µM and 1.81-11.0 µM respectively, in addition, of showing good safety profile against normal cell line (F180). Moreover, compounds 16a, 16c, 16d and 16 g had cell cycle arrest at G2/M phase with induced necrotic percentage compared to sorafenib of 2.06 %, 2.47 %, 1.57 %, 0.88 % and 1.83 % respectively. Amusingly, compounds 16a, 16c, 16d and 16 g inhibited VEGFR-2 with IC50 of 25 nM, 52 nM, 324 nM and 110 nM respectively, compared to sorafenib (IC50 = 85 nM), and had excellent EGFRWT and EGFRT790M kinase inhibitory activities (IC50 = 94 nM, 128 nM, 160 nM, 297 nM), (10 nM, 25 nM, 36 nM and 48 nM) respectively, compared to both erlotinib and osimertinib (IC50 = 114 nM, 56 nM) and (70 nM, 37 nM) respectively and showed (EGFRwt/EGFRT790M S.I.) of (range: 4.44-9.40) compared to erlotinib (2.03) and osmertinib (1.89).

Natural Compounds for Bone Remodeling: A Computational and Experimental Approach Targeting Bone Metabolism-Related Proteins

Osteoporosis, characterized by reduced bone density and increased fracture risk, affects over 200 million people worldwide, predominantly older adults and postmenopausal women. The disruption of the balance between bone-forming osteoblasts and bone-resorbing osteoclasts underlies osteoporosis pathophysiology. Standard treatment includes lifestyle modifications, calcium and vitamin D supplementation and specific drugs that either inhibit osteoclasts or stimulate osteoblasts. However, these treatments have limitations, including side effects and compliance issues. Natural products have emerged as potential osteoporosis therapeutics, but their mechanisms of action remain poorly understood. In this study, we investigate the efficacy of natural compounds in modulating molecular targets relevant to osteoporosis, focusing on the Mitogen-Activated Protein Kinase (MAPK) pathway and the gut microbiome's influence on bone homeostasis. Using an in silico and in vitro methodology, we have identified quercetin as a promising candidate in modulating MAPK activity, offering a potential therapeutic perspective for osteoporosis treatment.

Predicting the functional state of protein kinases using interpretable graph neural networks from sequence and structural data

Protein kinases are central to cellular activities and are actively pursued as drug targets for several conditions including cancer and autoimmune diseases. Despite the availability of a large structural database for kinases, methodologies to elucidate the structure-function relationship of these proteins (without manual intervention) are lacking. Such techniques are essential in structural biology and to accelerate drug discovery efforts. Here, we implement an interpretable graph neural network (GNN) framework for classifying the functionally active and inactive states of a large set of protein kinases by only using their tertiary structure and amino acid sequence. We show that the GNN models can classify kinase structures with high accuracy (>97%). We implement the Gradient-weighted Class Activation Mapping for graphs (Graph Grad-CAM) to automatically identify structurally important residues and residue-residue contacts of the kinases without any a priori input. We show that the motifs identified through the Graph Grad-CAM methodology are functionally critical, consistent with the existing kinase literature. Notably, the highly conserved DFG and HRD motifs of the well-known hydrophobic spine are identified by the interpretable framework in addition to some of the lesser known motifs. Further, using Grad-CAM maps as the vector embedding of the protein structures, we identify the subtle differences in the crystal structures among different sub-classes of kinases in the Protein Data Bank (PDB). Frameworks such as the one implemented here, for high-throughput identification of protein structure-function relationships are essential in designing targeted small molecules therapies as well as in engineering new proteins for novel applications.

Structural and molecular insights from dual inhibitors of EGFR and VEGFR2 as a strategy to improve the efficacy of cancer therapy

Epidermal growth factor receptor (EGFR) and vascular endothelial growth factor 2 (VEGFR2) are known as valid targets for cancer therapy. Overexpression of EGFR induces uncontrolled cell proliferation and VEGF expression triggering angiogenesis via VEGFR2 signaling. On the other hand, VEGF expression independent of EGFR signaling is already known as one of the mechanisms of resistance to anti-EGFR therapy. Therefore, drugs that act as dual inhibitors of EGFR and VEGFR2 can be a solution to the problem of drug resistance and increase the effectiveness of therapy. In this review, we summarize the relationship between EGFR and VEGFR2 signal transduction in promoting cancer growth and how their kinase domain structures can affect the selectivity of an inhibitor as the basis for designing dual inhibitors. In addition, several recent studies on the development of dual EGFR and VEGFR2 inhibitors involving docking simulations were highlighted in this paper to provide some references such as pharmacophore features of inhibitors and key residues for further research, especially in computer-aided drug design.

Discovery of Potent Multikinase Type-II Inhibitors Targeting CDK5 in the DFG-out Inactive State with Promising Potential against Glioblastoma

Kinases have proven valuable targets in successful cancer drug discovery projects, but not yet for malignant brain tumors where type-II inhibition of cyclin-dependent kinase 5 (CDK5) stabilizing the DFG-out inactive state has potential for design of selective and clinically efficient drug candidates. In the absence of crystallographic evidence for a CDK5 DFG-out inactive state protein-ligand complex, for the first time, a model was designed using metadynamics/molecular dynamics simulations. Glide docking of the ZINC15 biogenic database identified [pyrimidin-2-yl]amino-furo[3,2-b]-furyl-urea/amide hit chemical scaffolds. For four selected analogues (4, 27, 36, and 42), potent effects on glioblastoma cell viability in U87-MG, T98G, and U251-MG cell lines and patient-derived cultures were generally observed (IC50s ∼ 10-40 μM at 72 h). Selectivity profiling against 11 homologous kinases revealed multikinase inhibition (CDK2, CDK5, CDK9, and GSK-3α/β), most potent for GSK-3α in the nanomolar range (IC50s ∼ 0.23-0.98 μM). These compounds may therefore have diverse anticancer mechanisms of action and are of considerable interest for lead optimization.

Type II & III inhibitors of tropomyosin receptor kinase (Trk): a 2020-2022 patent update

INTRODUCTION

The Trk family proteins are membrane-bound kinases predominantly expressed in neuronal tissues. Activated by neurotrophins, they regulate critical cellular processes through downstream signaling pathways. Dysregulation of Trk signaling can drive a range of diseases, making the design and study of Trk inhibitors a vital area of research. This review explores recent advances in the development of type II and III Trk inhibitors, with implications for various therapeutic applications.

AREAS COVERED

Patents covering type II and III inhibitors targeting the Trk family are discussed as a complement of the previous review, Type I inhibitors of tropomyosin receptor kinase (Trk): a 2020-2022 patent update. Relevant patents were identified using the Web of Science database, Google, and Google Patents.

EXPERT OPINION

While type II and III Trk inhibitor development has advanced more gradually compared to their type I counterparts, they hold significant promise in overcoming resistance mutations and achieving enhanced subtype selectivity - a critical factor in reducing adverse effects associated with pan-Trk inhibition. Recent interdisciplinary endeavors have marked substantial progress in the design of subtype selective Trk inhibitors, with impressive success heralded by the type III inhibitors. Notably, the emergence of mutant-selective Trk inhibitors introduces an intriguing dimension to the field, offering precise treatment possibilities.

Novel Pyrimidine-5-Carbonitriles as potential apoptotic and antiproliferative agents by dual inhibition of EGFRWT/T790M and PI3k enzymes; Design, Synthesis, biological Evaluation, and docking studies

A new series of 6-(4-fluorophenyl)-2-(methylthio) pyrimidine-5-carbonitrile derivatives were designed and synthesized as EGFR/PI3K dual inhibitors, and potential antiproliferative agents. The new 22 compounds were screened by DTP-NCI against all NCI60 cell lines. Almost all compounds showed cytotoxic activity. Compound 7c showed a promising antitumour activity on CNS cancer (SNB-75), and ovarian cancer (OVAR-4) with IC50 < 0.01, and 0.64 µM, respectively. Fortunately, 7c exhibited a better safety profile on normal cells (WI-38) than doxorubicin by 2.2-fold. Compound 7c displayed selective inhibitory activity on EGFRt790m over EGFRWT with IC50 = 0.08, and 0.13 µM, respectively, wherefore it might overcome EGFR-TKIs resistance. In addition to its remarkable inhibitory activity on all PI3K isoforms, specifically PI3K-δ with IC50 = 0.64 µM Compared with LY294002 IC50 = 7.6 µM. Compound 7c arrested the cell cycle of SNB-75 & OVAR-4 at the G0-G1 phase coupled with apoptosis induction. The western blotting analysis approved decreasing the expression level of p-AKT coupled with an increase in Casp3, Casp9, and BAX proteins in the SNB-75 & OVAR-4 after being treated with 7c which may support the suggested mechanism of action of 7c as EGFR/PI3K dual inhibitor. Physicochemical parameters were forecasted using SwissADME online tool. MD showed the interaction of 7c with the crucial amino acids of the active domain of both EGFR/PI3K which may explain its potent inhibitory activities. In vivo study disclosed a significant decrease in tumor weight and the number of nodules in the group of mice treated with 7c compared with the control group.

Development of 1,5-diarylpyrazoles as EGFR/JNK-2 dual inhibitors: design, synthesis, moleecular docking, and bioactivity evaluation

The eradication of multifactorial diseases, such as cancer, requires the design of drug candidates that attack multiple targets that contribute to the progression and proliferation of such diseases. Here, 1,5-diarylpyrazole derivatives bearing vanillin or sulfanilamide are developed as potential dual inhibitors of epidermal growth factor receptor (EGFR)/c-Jun N-terminal kinase 2 (JNK-2) for possible anticancer activity. These derivatives inhibited the growths of DLD-1, HeLa, K-562, SUIT-2 and HepG2 cancer cell lines, with minimum concentration required to inhibit half of the cellular growth (IC50) values of 2.7-63 μM. The tests confirmed that 5b and 5d were potent JNK-2 inhibitors, with IC50 of 2.0 and 0.9 μM, respectively, whereas 6 h selectively inhibited EGFR protein kinase (EGFR-PK) (IC50 = 1.7 μM). Notably, 6c inhibited both kinases, with IC50 values of 2.7 and 3.0 μM against EGFR-PK and JNK-2, respectively, offering a reference for designing mutual inhibitors of EGFR/JNK-2. The docking studies revealed the ability of the pyrazole ring to bind to the hinge region of the ATP binding site, thereby supporting the experimental inhibitory results. Furthermore, the developed compounds could induce apoptosis and induce cell cycle arrest at different cell phases.

Expanding FTMap for Fragment-Based Identification of Pharmacophore Regions in Ligand Binding Sites

The knowledge of ligand binding hot spots and of the important interactions within such hot spots is crucial for the design of lead compounds in the early stages of structure-based drug discovery. The computational solvent mapping server FTMap can reliably identify binding hot spots as consensus clusters, free energy minima that bind a variety of organic probe molecules. However, in its current implementation, FTMap provides limited information on regions within the hot spots that tend to interact with specific pharmacophoric features of potential ligands. E-FTMap is a new server that expands on the original FTMap protocol. E-FTMap uses 119 organic probes, rather than the 16 in the original FTMap, to exhaustively map binding sites, and identifies pharmacophore features as atomic consensus sites where similar chemical groups bind. We validate E-FTMap against a set of 109 experimentally derived structures of fragment-lead pairs, finding that highly ranked pharmacophore features overlap with the corresponding atoms in both fragments and lead compounds. Additionally, comparisons of mapping results to ensembles of bound ligands reveal that pharmacophores generated with E-FTMap tend to sample highly conserved protein-ligand interactions. E-FTMap is available as a web server at https://eftmap.bu.edu.

Computer aided drug discovery (CADD) of a thieno[2,3-d]pyrimidine derivative as a new EGFR inhibitor targeting the ribose pocket

Depending on the pharmacophoric characteristics of EGFR inhibitors, a new thieno[2,3-d]pyrimidine derivative has been developed. Firstly, the potential inhibitory effect of the designed compound against EGFR has been proven by docking experiments that showed correct binding modes and excellent binding energies of -98.44 and -88.00 kcal/mol, against EGFR wild-type and mutant type, respectively. Furthermore, MD simulations studies confirmed the precise energetic, conformational, and dynamic alterations that occurred after binding to EGFR. The correct binding was also confirmed by essential dynamics studies. To further investigate the general drug-like properties of the developed candidate, in silico ADME and toxicity studies have also been carried out. The thieno[2,3-d]pyrimidine derivative was synthesized following the earlier promising findings. Fascinatingly, the synthesized compound (4) showed promising inhibitory effects against EGFRWT and EGFRT790M with IC50 values of 25.8 and 182.3 nM, respectively. Also, it exhibited anticancer potentialities against A549 and MCF-7cell lines with IC50 values of 13.06 and 20.13 µM, respectively. Interestingly, these strong activities were combined with selectivity indices of 2.8 and 1.8 against the two cancer cell lines, respectively. Further investigations indicated the ability of compound 4 to arrest the cancer cells' growth at the G2/M phase and to increase early and late apoptosis percentages from 2.52% and 2.80 to 17.99% and 16.72%, respectively. Additionally, it was observed that compound 4 markedly increased the levels of caspase-3 and caspase-9 by 4 and 3-fold compared to the control cells. Moreover, it up-regulated the level of BAX by 3-fold and down-regulated the level of Bcl-2 by 3-fold affording a BAX/Bcl-2 ratio of 9.Communicated by Ramaswamy H. Sarma.

Receptor-based pharmacophore modeling, molecular docking, synthesis and biological evaluation of novel VEGFR-2, FGFR-1, and BRAF multi-kinase inhibitors

A receptor-based pharmacophore model describing the binding features required for the multi-kinase inhibition of the target kinases (VEGFR-2, FGFR-1, and BRAF) were constructed and validated. It showed a good overall quality in discriminating between the active and the inactive in a compiled test set compounds with F1 score of 0.502 and Mathew's correlation coefficient of 0.513. It described the ligand binding to the hinge region Cys or Ala, the glutamate residue of the Glu-Lys αC helix conserved pair, the DFG motif Asp at the activation loop, and the allosteric back pocket next to the ATP binding site. Moreover, excluded volumes were used to define the steric extent of the binding sites. The application of the developed pharmacophore model in virtual screening of an in-house scaffold dataset resulted in the identification of a benzimidazole-based scaffold as a promising hit within the dataset. Compounds 8a-u were designed through structural optimization of the hit benzimidazole-based scaffold through (un)substituted aryl substitution on 2 and 5 positions of the benzimidazole ring. Molecular docking simulations and ADME properties predictions confirmed the promising characteristics of the designed compounds in terms of binding affinity and pharmacokinetic properties, respectively. The designed compounds 8a-u were synthesized, and they demonstrated moderate to potent VEGFR-2 inhibitory activity at 10 µM. Compound 8u exhibited a potent inhibitory activity against the target kinases (VEGFR-2, FGFR-1, and BRAF) with IC50 values of 0.93, 3.74, 0.25 µM, respectively. The benzimidazole derivatives 8a-u were all selected by the NCI (USA) to conduct their anti-proliferation screening. Compounds 8a and 8d resulted in a potent mean growth inhibition % (GI%) of 97.73% and 92.51%, respectively. Whereas compounds 8h, 8j, 8k, 8o, 8q, 8r, and 8u showed a mean GI% > 100% (lethal effect). The most potent compounds on the NCI panel of 60 different cancer cell lines were progressed further to NCI five-dose testing. The benzimidazole derivatives 8a, 8d, 8h, 8j, 8k, 8o, 8q, 8r and 8u exhibited potent anticancer activity on the tested cell lines reaching sub-micromolar range. Moreover, 8u was found to induce cell cycle arrest of MCF-7 cell line at the G2/M phase and accumulating cells at the sub-G1 phase as a result of cell apoptosis.

Discovery of A Novel Series of Quinazoline-Thiazole Hybrids as Potential Antiproliferative and Anti-Angiogenic Agents

Considering the pivotal role of angiogenesis in solid tumor progression, we developed a novel series of quinazoline-thiazole hybrids (SA01-SA07) as antiproliferative and anti-angiogenic agents. Four out of the seven compounds displayed superior antiproliferative activity (IC50 =1.83-4.24 µM) on HepG2 cells compared to sorafenib (IC50 = 6.28 µM). The affinity towards the VEGFR2 kinase domain was assessed through in silico prediction by molecular docking, molecular dynamics studies, and MM-PBSA. The series displayed a high degree of similarity to sorafenib regarding the binding pose within the active site of VEGFR2, with a different orientation of the 4-substituted-thiazole moieties in the allosteric pocket. Molecular dynamics and MM-PBSA evaluations identified SA05 as the hybrid forming the most stable complex with VEGFR2 compared to sorafenib. The impact of the compounds on vascular cell proliferation was assessed on EA.hy926 cells. Six compounds (SA01-SA05, SA07) displayed superior anti-proliferative activity (IC50 = 0.79-5.85 µM) compared to sorafenib (IC50 = 6.62 µM). The toxicity was evaluated on BJ cells. Further studies of the anti-angiogenic effect of the most promising compounds, SA04 and SA05, through the assessment of impact on EA.hy296 motility using a wound healing assay and in ovo potential in a CAM assay compared to sorafenib, led to the confirmation of the anti-angiogenic potential.

Deciphering the molecular choreography of Janus kinase 2 inhibition via Gaussian accelerated molecular dynamics simulations: a dynamic odyssey

The Janus kinases (JAK) are crucial targets in drug development for several diseases. However, accounting for the impact of possible structural rearrangements on the binding of different kinase inhibitors is complicated by the extensive conformational variability of their catalytic kinase domain (KD). The dynamic KD contains mainly four prominent mobile structural motifs: the phosphate-binding loop (P-loop), the αC-helix within the N-lobe, the Asp-Phe-Gly (DFG) motif, and the activation loop (A-loop) within the C-lobe. These distinct structural orientations imply a complex signal transmission path for regulating the A-loop's flexibility and conformational preference for optimal JAK function. Nevertheless, the precise dynamical features of the JAK induced by different types of inhibitors still remain elusive. We performed comparative, microsecond-long, Gaussian accelerated molecular dynamics simulations in triplicate of three phosphorylated JAK2 systems: the KD alone, type-I ATP-competitive inhibitor (CI) bound KD in the catalytically active DFG-in conformation, and the type-II inhibitor (AI) bound KD in the catalytically inactive DFG-out conformation. Our results indicate significant conformational variations observed in the A-loop and αC helix motions upon inhibitor binding. Our studies also reveal that the DFG-out inactive conformation is characterized by the closed A-loop rearrangement, open catalytic cleft of N and C-lobe, the outward movement of the αC helix, and open P-loop states. Moreover, the outward positioning of the αC helix impacts the hallmark salt bridge formation between Lys882 and Glu898 in an inactive conformation. Finally, we compared their ligand binding poses and free energy by the MM/PBSA approach. The free energy calculations suggested that the AI's binding affinity is higher than CI against JAK2 due to an increased favorable contribution from the total non-polar interactions and the involvement of the αC helix. Overall, our study provides the structural and energetic insights crucial for developing more promising type I/II JAK2 inhibitors for treating JAK-related diseases.

Tetrahydrobenzothiophene derivatives ameliorate Mia PaCa-2 cell progression and induces apoptosis via inhibiting EGFR2 tyrosine kinase signal

A series of new thiophene analogues with acarbonitrile-basedmoiety were designed and synthesized via structural optimization. The conjugates were assessed for their in-vitro cytotoxic activity against a human pancreatic cancer cell line (Mia PaCa-2) and among them compound 5b showed IC50 value of 13.37 ± 2.37 μM. The compounds 5b (20 µM & 25 µM) and 7c (30 & 35 µM) also showed reduced clonogenicity, enhanced ROS and decreased mitochondrial membrane potential in Mia PaCa-2 cells. Treatment with these compounds also increased apoptotic population as evident with the double staining assay. Among the evaluated series, compounds 5b, 5g, 7c, and 9a attained a greater inhibitory potency than first generation's reversible EGFR inhibitor, Gefitinib. EGFR2 enzyme inhibitory studies revealed that 5b efficiently and arbitrarily suppressed the development of EGFR2 dependent cells and inhibited the enzymatic activity with an IC50 value of 0.68 µM; interestingly, the most effective molecule 5b with N-methyl piperazine substitution, has 1.29-fold greater potency than well-known EGFR inhibitor Gefitinib and increased Gefitinib's anti-growth impact with 2.04 folds greater against Mia PaCa-2. The in-vitro studies were validated with in-silico docking studies wherein compounds 5b and 7c exhibited binding energies of -8.2 and -7.4 Kcal/mol respectively. The present study reveals that tetrahydrobenzothiophene based analogues could be a promising lead for the evolution of potent chemo preventives over pancreatic cancer.

N-(3-Methoxyphenyl)-6-(7-(1-methyl-1H-pyrazol-4-yl)imidazo[1,2-a]pyridin-3-yl)pyridin-2-amine is an inhibitor of the FLT3-ITD and BCR-ABL pathways, and potently inhibits FLT3-ITD/D835Y and FLT3-ITD/F691L secondary mutants

Activating mutations within FLT3 make up 30 % of all newly diagnosed acute myeloid leukemia (AML) cases, with the most common mutation being an internal tandem duplication (FLT3-ITD) in the juxtamembrane region (25 %). Currently, two generations of FLT3 kinase inhibitors have been developed, with three inhibitors clinically approved. However, treatment of FLT3-ITD mutated AML is limited due to the emergence of secondary clinical resistance, caused by multiple mechanism including on-target FLT3 secondary mutations - FLT3-ITD/D835Y and FLT3-ITD/F691L being the most common, as well as the off-target activation of alternative pathways including the BCR-ABL pathway. Through the screening of imidazo[1,2-a]pyridine derivatives, N-(3-methoxyphenyl)-6-(7-(1-methyl-1H-pyrazol-4-yl)imidazo[1,2-a]pyridin-3-yl)pyridin-2-amine (compound 1) was identified as an inhibitor of both the FLT3-ITD and BCR-ABL pathways. Compound 1 potently inhibits clinically related leukemia cell lines driven by FLT3-ITD, FLT3-ITD/D835Y, FLT3-ITD/F691L, or BCR-ABL. Studies indicate that it mediates proapoptotic effects on cells by inhibiting FLT3 and BCR-ABL pathways, and other possible targets. Compound 1 is more potent against FLT3-ITD than BCR-ABL, and it may have other possible targets; however, compound 1 is first step for further optimization for the development of a balanced FLT3-ITD/BCR-ABL dual inhibitor for the treatment of relapsed FLT3-ITD mutated AML with multiple secondary clinical resistant subtypes such as FLT3-ITD/D835Y, FLT3-ITD/F691L, and cells co-expressing FLT3-ITD and BCR-ABL.

Potential inhibitors of RPS6KB2 and NRF2 in head and neck squamous cell carcinoma

Among the major altered pathways in head and neck squamous cell carcinoma, AKT/mTORC1/S6K and NRF2/KEAP1 pathway are quite significant. The overexpression and overstimulation of proteins from both these pathways makes them the promising candidates in cancer therapeutics. Inhibiting mTOR has been in research from past several decades but the tumour heterogeneity, and upregulation of several compensatory feed-back mechanisms, encourages to explore other downstream targets for inhibiting the pathway. One such downstream effectors of mTOR is S6K2. It is reported to be overexpressed in cancers such as head and neck cancer, breast cancer and prostate cancer. In case of NRF2/KEAP1 pathway, nuclear factor erythroid 2-related factor 2 (NFE2L2 or NRF2) is overexpressed in ∼90% of head and neck squamous cell carcinoma (HNSCC) cases. It associates with poor survival rate and therapeutic resistance in HNSCC treatment. NRF2 pathway is the primary antioxidant pathway in the cell which also serves pro-tumorigenic functions, such as repression of apoptosis, cell proliferation support and chemoresistance. The aim of this work was to explore S6K2 and NRF2 and identify novel and potential inhibitors against them for treating head and neck squamous cell carcinoma. Since the crystal structure of S6K2 was not available at the time of this study, we modelled its structure using homology modelling and performed high throughput screening, molecular dynamics simulations, free energy calculations and protein-ligand interaction studies to identify the inhibitors. We identified natural compounds Crocin and Gypenoside XVII against S6K2 and Chebulinic acid and Sennoside A against NRF2. This study provides a significant in-depth understanding of the two studied pathways and therefore can be used in the development of potential therapeutics against HNSCC.Communicated by Ramaswamy H. Sarma.

An imidazo[1,2-a]pyridine-pyridine derivative potently inhibits FLT3-ITD and FLT3-ITD secondary mutants, including gilteritinib-resistant FLT3-ITD/F691L

FLT3 activating mutations are detected in approximately 30 % of newly diagnosed acute myeloid leukemia (AML) cases, most commonly consisting of internal tandem duplication (ITD) mutations in the juxtamembrane region. Recently, several FLT3 inhibitors have demonstrated clinical activity and three are currently approved - midostaurin, quizartinib, and gilteritinib. Midostaurin is a first-generation FLT3 inhibitor with minimal activity as monotherapy. Midostaurin lacks selectivity and is only approved by the USFDA for use in combination with other chemotherapy agents. The second-generation inhibitors quizartinib and gilteritinib display improved specificity and selectivity, and have been approved for use as monotherapy. However, their clinical efficacies are limited in part due to the emergence of drug-resistant FLT3 secondary mutations in the tyrosine kinase domain at positions D835 and F691. Therefore, in order to overcome drug resistance and further improve outcomes, new compounds targeting FLT3-ITD with secondary mutants are urgently needed. In this study, through the structural modification of a reported compound Ling-5e, we identified compound 24 as a FLT3 inhibitor that is equally potent against FLT3-ITD and the clinically relevant mutants FLT3-ITD/D835Y, and FLT3-ITD/F691L. Its inhibitory effects were demonstrated in both cell viability assays and western blots analyses. When tested against cell lines lacking activating mutations in FLT3, no non-specific cytotoxicity effect was observed. Interestingly, molecular docking results showed that compound 24 may adopt different binding conformations with FLT3-F691L compared to FLT3, which may explain its retained activity against FLT3-ITD/F691L. In summary, compound 24 has inhibition potency on FLT3 comparable to gilteritinib, but a more balanced inhibition on FLT3 secondary mutations, especially FLT3-ITD/F691L which is gilteritinib resistant. Compound 24 may serve as a promising lead for the drug development of either primary or relapsed AML with FLT3 secondary mutations.

Novel insights into the role of Discoidin domain receptor 2 (DDR2) in cancer progression: a new avenue of therapeutic intervention

Discoidin domain receptors (DDRs), including DDR1 and DDR2, are a unique class of receptor tyrosine kinases (RTKs) activated by collagens at the cell-matrix boundary interface. The peculiar mode of activation makes DDRs as key cellular sensors of microenvironmental changes, with a critical role in all physiological and pathological processes governed by collagen remodeling. DDRs are widely expressed in fetal and adult tissues, and experimental and clinical evidence has shown that their expression is deregulated in cancer. Strong findings supporting the role of collagens in tumor progression and metastasis have led to renewed interest in DDRs.  However, despite an increasing number of studies, DDR biology remains poorly understood, particularly the less studied DDR2, whose involvement in cancer progression mechanisms is undoubted. Thus, the understanding of a wider range of DDR2 functions and related molecular mechanisms is expected. To date, several lines of evidence support DDR2 as a promising target in cancer therapy. Its involvement in key functions in the tumor microenvironment makes DDR2 inhibition particularly attractive to achieve simultaneous targeting of tumor and stromal cells, and tumor regression, which is beneficial for improving the response to different types of anti-cancer therapies, including chemo- and immunotherapy. This review summarizes current research on DDR2, focusing on its role in cancer progression through its involvement in tumor and stromal cell functions, and discusses findings that support the rationale for future development of direct clinical strategies targeting DDR2.

Novel pyrazolo[3,4-d]pyrimidine derivatives: design, synthesis, anticancer evaluation, VEGFR-2 inhibition, and antiangiogenic activity

A novel series of 12 pyrazolo[3,4-d]pyrimidine derivatives were created and evaluated in vitro for their antiproliferative activity against the NCI 60 human tumor cell line panel. Compounds 12a-d displayed significant antitumor activity against MDA-MB-468 and T-47D (breast cancer cell lines), especially compound 12b, which exhibited the highest anticancer activity against MDA-MB-468 and T-47D cell lines with IC50 values of 3.343 ± 0.13 and 4.792 ± 0.21 μM, respectively compared to staurosporine with IC50 values of 6.358 ± 0.24 and 4.849 ± 0.22 μM. The most potent cytotoxic derivatives 12a-d were studied for their VEGFR-2 inhibitory activity to explore the mechanism of action of these substances. Compound 12b had potent activity against VEGFR-2 with an IC50 value of 0.063 ± 0.003 μM, compared to sunitinib with IC50 = 0.035 ± 0.012 μM. Moreover, there was an excellent reduction in HUVEC migratory potential that resulted in a significant disruption of wound healing patterns by 23% after 72 h of treatment with compound 12b. Cell cycle and apoptosis investigations showed that compound 12b could stop the cell cycle at the S phase and significantly increase total apoptosis in the MDA-MB-468 cell line by 18.98-fold compared to the control. Moreover, compound 12b increased the caspase-3 level in the MDA-MB-468 cell line by 7.32-fold as compared to the control.

Computer-Assisted Drug Discovery of a Novel Theobromine Derivative as an EGFR Protein-Targeted Apoptosis Inducer

The overexpression of the Epidermal Growth Factor Receptor (EGFR) marks it as a pivotal target in cancer treatment, with the aim of reducing its proliferation and inducing apoptosis. This study aimed at the CADD of a new apoptotic EGFR inhibitor. The natural alkaloid, theobromine, was used as a starting point to obtain a new semisynthetic (di-ortho-chloro acetamide) derivative (T-1-DOCA). Firstly, T-1-DOCA's total electron density, energy gap, reactivity indices, and electrostatic surface potential were determined by DFT calculations, Then, molecular docking studies were carried out to predict the potential of T-1-DOCA against wild and mutant EGFR proteins. T-1-DOCA's correct binding was further confirmed by molecular dynamics (MD) over 100 ns, MM-GPSA, and PLIP experiments. In vitro, T-1-DOCA showed noticeable efficacy compared to erlotinib by suppressing EGFRWT and EGFRT790M with IC50 values of 56.94 and 269.01 nM, respectively. T-1-DOCA inhibited also the proliferation of H1975 and HCT-116 malignant cell lines, exhibiting IC50 values of 14.12 and 23.39 µM, with selectivity indices of 6.8 and 4.1, respectively, indicating its anticancer potential and general safety. The apoptotic effects of T-1-DOCA were indicated by flow cytometric analysis and were further confirmed through its potential to increase the levels of BAX, Casp3, and Casp9, and decrease Bcl-2 levels. In conclusion, T-1-DOCA, a new apoptotic EGFR inhibitor, was designed and evaluated both computationally and experimentally. The results suggest that T-1-DOCA is a promising candidate for further development as an anti-cancer drug.

How Ligands Interact with the Kinase Hinge

ATP-competitive kinase inhibitors form hydrogen bond interactions with the kinase hinge region at the adenine binding site. Thus, it is crucial to explore hinge-ligand recognition as part of a rational drug design strategy. Here, harnessing known ligand-bound kinase structures and experimental assay resources, we first created a kinase structure-assay database (KSAD) containing 2705 nM ligand-bound kinase complexes. Then, using KSAD, we systematically investigate hinge-ligand binding patterns using interaction fingerprints, thereby delineating 15 different hydrogen-bond interaction modes. We believe these results will be valuable for de novo drug design and/or scaffold hopping of kinase-targeted drugs.

A novel in silico scaffold-hopping method for drug repositioning in rare and intractable diseases

In the field of rare and intractable diseases, new drug development is difficult and drug repositioning (DR) is a key method to improve this situation. In this study, we present a new method for finding DR candidates utilizing virtual screening, which integrates amino acid interaction mapping into scaffold-hopping (AI-AAM). At first, we used a spleen associated tyrosine kinase inhibitor as a reference to evaluate the technique, and succeeded in scaffold-hopping maintaining the pharmacological activity. Then we applied this method to five drugs and obtained 144 compounds with diverse structures. Among these, 31 compounds were known to target the same proteins as their reference compounds and 113 compounds were known to target different proteins. We found that AI-AAM dominantly selected functionally similar compounds; thus, these selected compounds may represent improved alternatives to their reference compounds. Moreover, the latter compounds were presumed to bind to the targets of their references as well. This new "compound-target" information provided DR candidates that could be utilized for future drug development.

Development of Kinase-Centric Drugs: A Computational Perspective

Kinases are prominent drug targets in the pharmaceutical and research community due to their involvement in signal transduction, physiological responses, and upon dysregulation, in diseases such as cancer, neurological and autoimmune disorders. Several FDA-approved small-molecule drugs have been developed to combat human diseases since Gleevec was approved for the treatment of chronic myelogenous leukemia. Kinases were considered "undruggable" in the beginning. Several FDA-approved small-molecule drugs have become available in recent years. Most of these drugs target ATP-binding sites, but a few target allosteric sites. Among kinases that belong to the same family, the catalytic domain shows high structural and sequence conservation. Inhibitors of ATP-binding sites can cause off-target binding. Because members of the same family have similar sequences and structural patterns, often complex relationships between kinases and inhibitors are observed. To design and develop drugs with desired selectivity, it is essential to understand the target selectivity for kinase inhibitors. To create new inhibitors with the desired selectivity, several experimental methods have been designed to profile the kinase selectivity of small molecules. Experimental approaches are often expensive, laborious, time-consuming, and limited by the available kinases. Researchers have used computational methodologies to address these limitations in the design and development of effective therapeutics. Many computational methods have been developed over the last few decades, either to complement experimental findings or to forecast kinase inhibitor activity and selectivity. The purpose of this review is to provide insight into recent advances in theoretical/computational approaches for the design of new kinase inhibitors with the desired selectivity and optimization of existing inhibitors.

Exploring the Druggable Conformational Space of Protein Kinases Using AI-Generated Structures

Protein kinase function and interactions with drugs are controlled in part by the movement of the DFG and ɑC-Helix motifs, which enable kinases to adopt various conformational states. Small molecule ligands elicit therapeutic effects with distinct selectivity profiles and residence times that often depend on the kinase conformation(s) they bind. However, the limited availability of experimentally determined structural data for kinases in inactive states restricts drug discovery efforts for this major protein family. Modern AI-based structural modeling methods hold potential for exploring the previously experimentally uncharted druggable conformational space for kinases. Here, we first evaluated the currently explored conformational space of kinases in the PDB and models generated by AlphaFold2 (AF2) (1) and ESMFold (2), two prominent AI-based structure prediction methods. We then investigated AF2's ability to predict kinase structures in different conformations at various multiple sequence alignment (MSA) depths, based on this parameter's ability to explore conformational diversity. Our results showed a bias within the PDB and predicted structural models generated by AF2 and ESMFold toward structures of kinases in the active state over alternative conformations, particularly those conformations controlled by the DFG motif. Finally, we demonstrate that predicting kinase structures using AF2 at lower MSA depths allows the exploration of the space of these alternative conformations, including identifying previously unobserved conformations for 398 kinases. The results of our analysis of structural modeling by AF2 create a new avenue for the pursuit of new therapeutic agents against a notoriously difficult-to-target family of proteins.

Significance Statement

Greater abundance of kinase structural data in inactive conformations, currently lacking in structural databases, would improve our understanding of how protein kinases function and expand drug discovery and development for this family of therapeutic targets. Modern approaches utilizing artificial intelligence and machine learning have potential for efficiently capturing novel protein conformations. We provide evidence for a bias within AlphaFold2 and ESMFold to predict structures of kinases in their active states, similar to their overrepresentation in the PDB. We show that lowering the AlphaFold2 algorithm's multiple sequence alignment depth can help explore kinase conformational space more broadly. It can also enable the prediction of hundreds of kinase structures in novel conformations, many of whose models are likely viable for drug discovery.

Next-Generation JAK2 Inhibitors for the Treatment of Myeloproliferative Neoplasms: Lessons from Structure-Based Drug Discovery Approaches

Selective inhibitors of Janus kinase (JAK) 2 have been in demand since the discovery of the JAK2 V617F mutation present in patients with myeloproliferative neoplasms (MPN); however, the structural basis of V617F oncogenicity has only recently been elucidated. New structural studies reveal a role for other JAK2 domains, beyond the kinase domain, that contribute to pathogenic signaling. Here we evaluate the structure-based approaches that led to recently-approved type I JAK2 inhibitors (fedratinib and pacritinib), as well as type II (BBT594 and CHZ868) and pseudokinase inhibitors under development (JNJ7706621). With full-length JAK homodimeric structures now available, superior selective and mutation-specific JAK2 inhibitors are foreseeable.

SIGNIFICANCE

The JAK inhibitors currently used for the treatment of MPNs are effective for symptom management but not for disease eradication, primarily because they are not strongly selective for the mutant clone. The rise of computational and structure-based drug discovery approaches together with the knowledge of full-length JAK dimer complexes provides a unique opportunity to develop better targeted therapies for a range of conditions driven by pathologic JAK2 signaling.