Structural dynamics and kinase inhibitory activity of three generations of tyrosine kinase inhibitors against wild-type, L858R/T790M, and L858R/T790M/C797S forms of EGFR

Imidazole[1,5-a]pyridine derivatives as EGFR tyrosine kinase inhibitors unraveled by umbrella sampling and steered molecular dynamics simulations

Although the use of the tyrosine kinase inhibitors (TKIs) has been proved that it can save live in a cancer treatment, the currently used drugs bring in many undesirable side-effects. Therefore, the search for new drugs and an evaluation of their efficiency are intensively carried out. Recently, a series of eighteen imidazole[1,5-a]pyridine derivatives were synthetized by us, and preliminary analyses pointed out their potential to be an important platform for pharmaceutical development owing to their promising actions as anticancer agents and enzyme (kinase, HIV-protease,…) inhibitors. In the present theoretical study, we further analyzed their efficiency in using a realistic scenario of computational drug design. Our protocol has been developed to not only observe the atomistic interaction between the EGFR protein and our 18 novel compounds using both umbrella sampling and steered molecular dynamics simulations, but also determine their absolute binding free energies. Calculated properties of the 18 novel compounds were in detail compared with those of two known drugs, erlotinib and osimertinib, currently used in cancer treatment. Inspiringly the simulation results promote three imidazole[1,5-a]pyridine derivatives as promising inhibitors into a further step of clinical trials.

GraphEGFR: Multi-task and transfer learning based on molecular graph attention mechanism and fingerprints improving inhibitor bioactivity prediction for EGFR family proteins on data scarcity

The proteins within the human epidermal growth factor receptor (EGFR) family, members of the tyrosine kinase receptor family, play a pivotal role in the molecular mechanisms driving the development of various tumors. Tyrosine kinase inhibitors, key compounds in targeted therapy, encounter challenges in cancer treatment due to emerging drug resistance mutations. Consequently, machine learning has undergone significant evolution to address the challenges of cancer drug discovery related to EGFR family proteins. However, the application of deep learning in this area is hindered by inherent difficulties associated with small-scale data, particularly the risk of overfitting. Moreover, the design of a model architecture that facilitates learning through multi-task and transfer learning, coupled with appropriate molecular representation, poses substantial challenges. In this study, we introduce GraphEGFR, a deep learning regression model designed to enhance molecular representation and model architecture for predicting the bioactivity of inhibitors against both wild-type and mutant EGFR family proteins. GraphEGFR integrates a graph attention mechanism for molecular graphs with deep and convolutional neural networks for molecular fingerprints. We observed that GraphEGFR models employing multi-task and transfer learning strategies generally achieve predictive performance comparable to existing competitive methods. The integration of molecular graphs and fingerprints adeptly captures relationships between atoms and enables both global and local pattern recognition. We further validated potential multi-targeted inhibitors for wild-type and mutant HER1 kinases, exploring key amino acid residues through molecular dynamics simulations to understand molecular interactions. This predictive model offers a robust strategy that could significantly contribute to overcoming the challenges of developing deep learning models for drug discovery with limited data and exploring new frontiers in multi-targeted kinase drug discovery for EGFR family proteins.

Discovery of furopyridine-based compounds as novel inhibitors of Janus kinase 2: In silico and in vitro studies

Janus kinase 2 (JAK2), one of the JAK isoforms participating in a JAK/STAT signaling cascade, has been considered a potential clinical target owing to its critical role in physiological processes involved in cell growth, survival, development, and differentiation of various cell types, especially immune and hematopoietic cells. Substantial studies have proven that the inhibition of this target could disrupt the JAK/STAT pathway and provide therapeutic outcomes for cancer, immune disorders, inflammation, and COVID-19. Herein, we performed docking-based virtual screening of 63 in-house furopyridine-based compounds and verified the first-round screened compounds by in vitro enzyme- and cell-based assays. By shedding light on the integration of both in silico and in vitro methods, we could elucidate two promising compounds. PD19 showed cytotoxic effects on human erythroblast cell lines (TF-1 and HEL) with IC50 values of 57.27 and 27.28 μM, respectively, while PD12 exhibited a cytotoxic effect on TF-1 with an IC50 value of 83.47 μM by suppressing JAK2/STAT5 autophosphorylation. In addition, all screened compounds were predicted to meet drug-like criteria based on Lipinski's rule of five, and none of the extreme toxicity features were found. Molecular dynamic simulations revealed that PD12 and PD19 could form stable complexes with JAK2 in an aqueous environment, and the van der Waals interactions were the main force driving the complex formation. Besides, all compounds sufficiently interacted with surrounding amino acids in all crucial regions, including glycine, catalytic, and activation loops. Altogether, PD12 and PD19 identified here could potentially be developed as novel therapeutic inhibitors disrupting the JAK/STAT pathway.

Screening of phytochemicals from Clerodendrum inerme (L.) Gaertn as potential anti-breast cancer compounds targeting EGFR: an in-silico approach

Breast cancer (BC) is the most prevalent malignancy among women around the world. The epidermal growth factor receptor (EGFR) is a tyrosine kinase receptor (RTK) of the ErbB/HER family. It is essential for triggering the cellular signaling cascades that control cell growth and survival. However, perturbations in EGFR signaling lead to cancer development and progression. Hence, EGFR is regarded as a prominent therapeutic target for breast cancer. Therefore, in the current investigation, EGFR was targeted with phytochemicals from Clerodendrum inerme (L.) Gaertn (C. inerme). A total of 121 phytochemicals identified by gas chromatography-mass spectrometry (GC-MS) analysis were screened against EGFR through molecular docking, ADMET analyses (Absorption, Distribution, Metabolism, Excretion, and Toxicity), PASS predictions, and molecular dynamics simulation, which revealed three potential hit compounds with CIDs 10586 [i.e. alpha-bisabolol (-6.4 kcal/mol)], 550281 [i.e. 2,(4,4-Trimethyl-3-hydroxymethyl-5a-(3-methyl-but-2-enyl)-cyclohexene) (-6.5 kcal/mol)], and 161271 [i.e. salvigenin (-7.4 kcal/mol)]. The FDA-approved drug gefitinib was used to compare the inhibitory effects of the phytochemicals. The top selected compounds exhibited good ADMET properties and obeyed Lipinski's rule of five (ROF). The molecular docking analysis showed that salvigenin was the best among the three compounds and formed bonds with the key residue Met 793. Furthermore, the molecular mechanics generalized born surface area (MMGBSA) calculations, molecular dynamics simulation, and normal mode analysis validated the binding affinity of the compounds and also revealed the strong stability and compactness of phytochemicals at the docked site. Additionally, DFT and DOS analyses were done to study the reactivity of the compounds and to further validate the selected phytochemicals. These results suggest that the identified phytochemicals possess high inhibitory potential against the target EGFR and can treat breast cancer. However, further in vitro and in vivo investigations are warranted towards the development of these constituents into novel anti-cancer drugs.Communicated by Ramaswamy H. Sarma.

Recruitment of hexahydroquinoline as anticancer scaffold targeting inhibition of wild and mutants EGFR (EGFRWT, EGFRT790M, and EGFRL858R)

Hexahydroquinoline (HHQ) scaffold was constructed and recruited for development of new series of anticancer agents. Thirty-two new compounds were synthesised where x-ray crystallography was performed to confirm enantiomerism. Thirteen compounds showed moderate to good activity against NCI 60 cancer cell lines, with GI % mean up to 74% for 10c. Expending erlotinib as a reference drug, target compounds were verified for their inhibiting activities against EGFRWT, EGFRT790M, and EGFRL858R where compound 10d was the best inhibitor with IC50 = 0.097, 0.280, and 0.051 µM, respectively, compared to erlotinib (IC50 = 0.082 µM, 0.342 µM, and 0.055 µM, respectively). Safety profile was validated using normal human lung (IMR-90) cells. 10c and 10d disrupted cell cycle at pre-G1 and G2/M phases in lung cancer, HOP-92, and cell line. Molecular docking study was achieved to understand the potential binding interactions and affinities in the active sites of three versions of EGFRs.

Unveiling Potential Targeted Therapeutic Opportunities for Co-Overexpressed Targeting Protein for Xklp2 and Aurora-A Kinase in Lung Adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the most prevalent and leading causes of cancer deaths globally, with limited diagnostic and clinically significant therapeutic targets. Identifying the genes and processes involved in developing and progressing LUAD is crucial for developing effective targeted therapeutics and improving patient outcomes. Therefore, the study aimed to explore the RNA sequencing data of LUAD from The Cancer Genome Atlas (TCGA) and gene expression profile datasets involving GSE10072, GSE31210, and GSE32863 from the Gene Expression Omnibus (GEO) databases. The differential gene expression and the downstream analysis determined clinically significant biomarkers using a network-based approach. These therapeutic targets predominantly enriched the dysregulation of mitotic cell cycle regulation and revealed the co-overexpression of Aurora-A Kinase (AURKA) and Targeting Protein for Xklp2 (TPX2) with high survival risk in LUAD patients. The hydrophobic residues of the AURKA-TPX2 interaction were considered as the target site to block the autophosphorylation of AURKA during the mitotic cell cycle. The tyrosine kinase inhibitor (TKI) dacomitinib demonstrated the strong binding potential to hinder TPX2, shielding the AURKA destabilization. This in silico study lays the foundation for repurposing targeted therapeutic options to impede the Protein-Protein Interactions (PPIs) in LUAD progression and aid in future translational investigations.

Gefitinib derivatives and drug-resistance: A perspective from molecular dynamics simulations

Epidermal-growth factor receptor (EGFR) is a transmembrane tyrosine kinase (TK) with a significant role in cell survival. EGFR is upregulated in various cancer cells and known as a druggable target. Gefitinib is a first-line TK inhibitor used against metastatic non-small cell lung cancer (NSCLC). Despite initial clinical response, a conserved therapeutic effect could not be achieved due to the occurrence of resistance mechanisms. Point mutations in EGFR genes are one of the major causes of rendered tumor sensitivity. To aid in the development of more efficient TKIs, chemical structures of prevailing drugs and their target binding patterns are very important. The aim of the present study was to propose synthetically-accessible gefitinib congeners with enhanced binding fitness to clinically frequent EGFR mutants. Docking simulations of intended molecules identified 1-(4-(3-chloro-4-fluorophenylamino)-7-methoxyquinazolin-6-yl)-3-(oxazolidin-2-ylmethyl) thiourea (23) as a top-binder structure inside G719S, T790 M, L858R and T790 M/L858R-EGFR active sites. Superior docked complexes were subjected to the entire 400 ns molecular dynamics (MD) simulations. Analysis of data revealed the stability of mutant enzymes upon binding to molecule 23. All mutant complexes with the exception of a T790 M/L858R-EGFR, were majorly stabilized through cooperative hydrophobic contacts. Pairwise analysis of H-bonds proved Met793 as the conserved residue with stable H-bond participations as hydrogen bond donor (Frequency 63-96%). Amino acid decomposition analysis confirmed the probable role of Met793 in complex stabilization. Estimated binding free energies indicated the proper accommodation of molecule 23 inside target active sites. Pairwise energy decompositions of stable binding modes revealed the energetic contribution of key residues. Although wet lab experiments are required to unravel the mechanistic details of mEGFR inhibition, MD results provide structural basis for those events that are difficult to address experimentally. The outputs of the current study may assist to design small molecules with high potency to mEGFRs.

SUMOylation of AnxA6 facilitates EGFR-PKCα complex formation to suppress epithelial cancer growth

BACKGROUND

The Annexin A6 (AnxA6) protein is known to inhibit the epidermal growth factor receptor (EGFR)-extracellular signal regulated kinase (ERK)1/2 signaling upon EGF stimulation. While the biochemical mechanism of AnxA6 inactivating phosphorylation of EGFR and ERK1/2 is not completely explored in cancer cells.

METHODS

Cells were transiently co-transfected with pFlag-AnxA6, pHA-UBC9 and pHis-SUMO1 plasmids to enrich the SUMOylated AnxA6 by immunoprecipitation, and the modification level of AnxA6 by SUMO1 was detected by Western blot against SUMO1 antibody. The SUMOylation level of AnxA6 was compared in response to chemical SUMOylation inhibitor treatment. AnxA6 SUMOylation sites were further identified by LC-MS/MS and amino acid site mutation validation. AnxA6 gene was silenced through AnxA6 targeting shRNA-containing pLKO.1 lentiviral transfection in HeLa cells, while AnxA6 gene was over-expressed within the Lenti-Vector carrying AnxA6 or mutant AnxA6^K299R plasmid in A431 cells using lentiviral infections. Moreover, the mutant plasmid pGFP-EGFR^T790M/L858R was constructed to test AnxA6 regulation on EGFR mutation-induced signal transduction. Moreover, cell proliferation, migration, and gefitinib chemotherapy sensitivity were evaluated in HeLa and A431 cells under AnxA6 konckdown or AnxA6 overexpression by CCK8, colony form and wound healing assays. And tumorigenicity in vivo was measured in epithelial cancer cells-xenografted nude mouse model.

RESULTS

AnxA6 was obviously modified by SUMO1 conjugation within Lys (K) residues, and the K299 was one key SUMOylation site of AnxA6 in epithelial cancer cells. Compared to the wild type AnxA6, AnxA6 knockdown and its SUMO site mutant AnxA6^K299R showed less suppression of dephosphorylation of EGFR-ERK1/2 under EGF stimulation. The SUMOylated AnxA6 was prone to bind EGFR in response to EGF inducement, which facilitated EGFR-PKCα complex formation to decrease the EGF-induced phosphorylation of EGFR-ERK1/2 and cyclin D1 expression. Similarly, AnxA6 SUMOylation inhibited dephosphorylation of the mutant EGFR, thereby impeding EGFR mutation-involved signal transduction. Moreover, AnxA6 knockdown or the K299 mutant AnxA6^K299R conferred AnxA6 inability to suppress tumor progression, resulting in drug resistance to gefitinib in epithelial cancer cells. And in epithelial cancer cells-xenografted nude mouse model, both the weight and size of tumors derived from AnxA6 knockdown or AnxA6^K299R mutation-expressing cells were much greater than that of AnxA6-expressing cells.

CONCLUSIONS

Besides EGFR gene mutation, protein SUMOylation modification of EGFR-binding protein AnxA6 also functions pivotal roles in mediating epithelial cancer cell growth and gefitinib drug effect. Video Abstract.

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 IC₅₀ 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 IC₅₀ 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.

A Closer Look at EGFR Inhibitor Resistance in Non-Small Cell Lung Cancer through the Lens of Precision Medicine

The development of EGFR small-molecule inhibitors has provided significant benefit for the affected patient population. Unfortunately, current inhibitors are no curative therapy, and their development has been driven by on-target mutations that interfere with binding and thus inhibitory activity. Genomic studies have revealed that, in addition to these on-target mutations, there are also multiple off-target mechanisms of EGFR inhibitor resistance and novel therapeutics that can overcome these challenges are sought. Resistance to competitive 1st-generation and covalent 2nd- and 3rd-generation EGFR inhibitors is overall more complex than initially thought, and novel 4th-generation allosteric inhibitors are expected to suffer from a similar fate. Additional nongenetic mechanisms of resistance are significant and can include up to 50% of the escape pathways. These potential targets have gained recent interest and are usually not part of cancer panels that look for alterations in resistant patient specimen. We discuss the duality between genetic and nongenetic EGFR inhibitor drug resistance and summarize current team medicine approaches, wherein clinical developments, hand in hand with drug development research, drive potential opportunities for combination therapy.

Overall survival in advanced epidermal growth factor receptor mutated non-small cell lung cancer using different tyrosine kinase inhibitors in The Netherlands: a retrospective, nationwide registry study

Background

Clinical guidelines advise osimertinib as preferred first line treatment for advanced epidermal growth factor receptor (EGFR) mutated non-small cell lung cancer (NSCLC) with deletions in exon 19 (del19) or exon 21 L858R mutation. However, for first-line osimertinib the real world overall survival (OS) in mutation subgroups remains unknown. Therefore, the aim of this study was to evaluate the real-world OS of those patients treated with different generations of EGFR-tyrosine kinase inhibitors (TKI), and to identify predictors of survival.

Methods

Using real-world data from the Netherlands Cancer Registry (NCR) we assessed patients diagnosed with stage IV NSCLC with del19 or L858R mutation between January 1, 2015, and December 31, 2020, primarily treated with then regularly available TKIs (including osimertinib).

Findings

Between January 1, 2015, and December 31, 2020, 57,592 patients were included in the NCR. Within this cohort we identified 1109 patients, 654 (59%) with del19 and 455 (41%) with L858R mutations, respectively; 230 (21%) patients were diagnosed with baseline brain metastases (BM). Patients were treated with gefitinib (19%, 213/1109), erlotinib (42%, 470/1109), afatinib (15%, 161/1109) or osimertinib (24%, 265/1109). Median OS was superior for del19 versus L858R (28.4 months (95% CI 25.6-30.6) versus 17.7 months (95% CI 16.1-19.5), p < 0.001. In multivariable analysis, no difference in survival was observed between various TKIs in both groups. Only in the subgroup of patients with del19 and baseline BM, a benefit was observed for treatment with osimertinib.

Interpretation

In this nationwide real-world cohort, survival of Dutch patients with advanced NSCLC and an EGFR del19 mutation was superior versus those harboring an L858R mutation. Osimertinib performed only better as first-line treatment in patients with del19 and BM.

Funding

None.

Quinoxalinones as A Novel Inhibitor Scaffold for EGFR (L858R/T790M/C797S) Tyrosine Kinase: Molecular Docking, Biological Evaluations, and Computational Insights

Combating acquired drug resistance of EGFR tyrosine kinase (TK) is a great challenge and an urgent necessity in the management of non-small cell lung cancers. The advanced EGFR (L858R/T790M/C797S) triple mutation has been recently reported, and there have been no specific drugs approved for this strain. Therefore, our research aimed to search for effective agents that could impede the function of EGFR (L858R/T790M/C797S) TK by the integration of in silico and in vitro approaches. Our in-house quinoxalinone-containing compounds were screened through molecular docking and their biological activity was then verified by enzyme- and cell-based assay. We found that the four quinoxalinone-containing compounds including CPD4, CPD15, CPD16, and CPD21 were promising to be novel EGFR (L858R/T790M/C797S) TK inhibitors. The IC₅₀ values measured by the enzyme-based assay were 3.04 ± 1.24 nM; 6.50 ± 3.02 nM,10.50 ± 1.10 nM; and 3.81 ± 1.80 nM, respectively, which are at a similar level to a reference drug; osimertinib (8.93 ± 3.01 nM). Besides that, they displayed cytotoxic effects on a lung cancer cell line (H1975) with IC₅₀ values in the range of 3.47 to 79.43 μM. In this proposed study, we found that all screened compounds could interact with M793 at the hinge regions and two mutated residues including M790 and S797; which may be the main reason supporting the inhibitory activity in vitro. The structural dynamics revealed that the screened compounds have sufficient non-native contacts with surrounding amino acids and could be well-buried in the binding site's cleft. In addition, all predicted physicochemical parameters were favorable to be drug-like based on Lipinski's rule of five, and no extreme violation of toxicity features was found. Altogether, this study proposes a novel EGFR (L858R/T790M/C797S) TK inhibitor scaffold and provides a detailed understanding of compounds' recognition and susceptibility at the molecular level.