Epidermal growth factor receptor (EGFR) has been recognized as one of the attractive targets for anticancer drug development. Herein, a set of anilino-1,4-naphthoquinone derivatives (3-18) was synthesized and investigated for their anticancer and EGFR inhibitory potentials. Among all tested compounds, three derivatives (3, 8, and 10) were selected for studying EGFR inhibitory activity (in vitro and in silico) due to their most potent cytotoxic activities against six tested cancer cell lines (i.e., HuCCA-1, HepG2, A549, MOLT-3, MDA-MB-231, and T47D; IC₅₀ values = 1.75-27.91 μM), high selectivity index (>20), and good predicted drug-like properties. The experimental results showed that these three promising compounds are potent EGFR inhibitors with nanomolar IC₅₀ values (3.96-18.64 nM). Interestingly, the most potent compound 3 bearing 4-methyl substituent on the phenyl ring displayed 4-fold higher potency than the known EGFR inhibitor, erlotinib. Molecular docking, molecular dynamics simulation, and MM/GBSA-based free energy calculation revealed that van der Waals force played a major role in the accommodations of compound 3 within the ATP-binding pocket of EGFR. Additionally, the 4-CH₃ moiety of the compound was noted to be a key chemical feature contributing to the highly potent EGFR inhibitory activity via its formations of alkyl interactions with A743, K745, M766, and L788 residues as well as additional interactions with M766 and T790.

Epidermal growth factor receptor (EGFR), overexpressed in many types of cancer, has been proved as a high potential target for targeted cancer therapy due to its role in regulating proliferation and survival of cancer cells. In the present study, a series of designed vinyl sulfone derivatives was screened against EGFR tyrosine kinase (EGFR-TK) using in silico and in vitro studies. The molecular docking results suggested that, among 78 vinyl sulfones, there were eight compounds that could interact well with the EGFR-TK at the ATP-binding site. Afterwards, these screened compounds were tested for the inhibitory activity towards EGFR-TK using ADP-Glo™ kinase assay, and we found that only VF16 compound exhibited promising inhibitory activity against EGFR-TK with the IC₅₀ value of 7.85 ± 0.88 nM. In addition, VF16 showed a high cytotoxicity with IC₅₀ values of 33.52 ± 2.57, 54.63 ± 0.09, and 30.38 ± 1.37 µM against the A431, A549, and H1975 cancer cell lines, respectively. From 500-ns MD simulation, the structural stability of VF16 in complex with EGFR-TK was quite stable, suggesting that this compound could be a novel small molecule inhibitor targeting EGFR-TK.

Mutations in the tyrosine kinase domain of epidermal growth factor receptor (EGFR), including L858R/T790M double and L858R/T790M/C797S triple mutations, are major causes of acquired resistance towards EGFR targeted drugs. In this work, a combination of comprehensive molecular modeling and in vitro kinase inhibition assay was used to unravel the mutational effects of EGFR on the susceptibility of three generations of EGFR tyrosine kinase inhibitors (erlotinib, gefitinib, afatinib, dacomitinib, and osimertinib) in comparison with the wild-type EGFR. The binding affinity of all studied inhibitors towards the double and triple EGFR mutations was in good agreement with the experimental data, ranked in the order of osimertinib > afatinib > dacomitinib > erlotinib > gefitinib. Three hot-spot residues at the hinge region (M790, M793, and C797) were involved in the binding of osimertinib and afatinib, enhancing their inhibitory activity towards mutated EGFRs. Both double and triple EGFR mutations causing erlotinib and gefitinib resistance are mainly caused by the low number of H-bond occupations, the low number of surrounding atoms, and the high number of water molecules accessible to the enzyme active site. According to principal component analysis, the molecular complexation of osimertinib against the two mutated EGFRs was in a closed conformation, whereas that against wild-type EGFR was in an open conformation, resulting in drug resistance. This work paves the way for further design of the novel EGFR inhibitors to overcome drug resistance mechanisms.

Sorafenib (SOR) is an oral multikinase inhibitor that effectively hampers the growth and spread of cancer cells by targeting angiogenesis and proliferation. However, SOR tablets (Nexavar) have limited oral bioavailability, ranging from 38% to 49%, due to their low water solubility. To address this issue, cyclodextrins (CDs), widely used to enhance the solubility and stability of lipophilic drugs by encapsulating them within their molecular structure, were considered in this study. We focused on β-cyclodextrin (βCD) and its derivatives, including hydroxypropyl-β-cyclodextrin (HPβCD), dimethyl-β-cyclodextrin (DMβCD), sulfobutylether-β-cyclodextrin (SBEβCD), and compared them with γ-cyclodextrin (γCD) for generating inclusion complexes with SOR. The 200 ns molecular dynamics simulations revealed that SOR could form inclusion complexes with all CDs in two possible orientations: pyridine group insertion (P-form) and chlorobenzotrifluoride group insertion (C-form), primarily driven by van der Waals interactions. Among the four βCD derivatives studied, SOR exhibited the highest number of atom contacts with SBEβCD and demonstrated the lowest solvent accessibility within the hydrophobic cavity of SBEβCD. These findings correlated with the highest binding affinity of SOR/SBEβCD complex determined by SIE, MM/GBSA, and MM/PBSA methods. Experimental results further supported our computational predictions, in which SBEβCD exhibited a stability constant of 940 M⁻¹ at 25 °C, surpassing βCD's stability constant of 210 M⁻¹. Taken together, our results suggest that the modified CDs, particularly SBEβCD, hold promising potential as an efficient molecular encapsulating agent for SOR, offering improved solubility and stability for this lipophilic drug.

Targeting L858R/T790M and L858R/T790M/C797S mutant EGFR is a critical challenge in developing EGFR tyrosine kinase inhibitors to overcome drug resistance in non-small cell lung cancer (NSCLC). The discovery of next-generation EGFR tyrosine kinase inhibitors (TKIs) is therefore necessary. To this end, a series of furopyridine derivatives were evaluated for their EGFR-based inhibition and antiproliferative activities using computational and biological approaches. We found that several compounds derived from virtual screening based on a molecular docking and solvated interaction energy (SIE) method showed the potential to suppress wild-type and mutant EGFR. The most promising PD13 displayed strong inhibitory activity against wild-type (IC₅₀ of 11.64 ± 1.30 nM), L858R/T790M (IC₅₀ of 10.51 ± 0.71 nM), which are more significant than known drugs. In addition, PD13 revealed a potent cytotoxic effect on A549 and H1975 cell lines with IC₅₀ values of 18.09 ± 1.57 and 33.87 ± 0.86 µM, respectively. The 500-ns MD simulations indicated that PD13 formed a hydrogen bond with Met793 at the hinge region, thus creating excellent EGFR inhibitory activity. Moreover, the binding of PD13 in the hinge region of EGFR was the major determining factor in stabilizing the interactions via hydrogen bonds and van der Waals (vdW). Altogether, PD13 is a promising novel EGFR inhibitor that could be further clinically developed as fourth-generation EGFR-TKIs.