Development of triple mutant T790M/C797S allosteric EGFR inhibitors: a computational approach

Computational investigation of Moringa oleifera phytochemicals targeting EGFR: molecular docking, molecular dynamics simulation and density functional theory studies

Epidermal growth factor receptor (EGFR) is a prominent target for anticancer therapy due to its role in activating several cell signaling cascades. Clinically approved EGFR inhibitors are reported to show treatment resistance and toxicity, this study, therefore, investigates Moringa oleifera phytochemicals to find potent and safe anti-EGFR compounds. For that, phytochemicals were screened based on drug-likeness and molecular docking analysis followed by molecular dynamics simulation, density functional theory analysis and ADMET analysis to identify the effective inhibitors of EGFR tyrosine kinase (EGFR-TK) domain. Known EGFR-TK inhibitors (1-4 generations) were used as control. Among 146 phytochemicals, 136 compounds showed drug-likeness, of which Delta 7-Avenasterol was the most potential EGFR-TK inhibitor with a binding energy of -9.2 kcal/mol followed by 24-Methylenecholesterol (-9.1 kcal/mol), Campesterol (-9.0 kcal/mol) and Ellagic acid (-9.0 kcal/mol). In comparison, the highest binding affinity from control drugs was displayed by Rociletinib (-9.0 kcal/mol). The molecular dynamics simulation (100 ns) exhibited the structural stability of native EGFR-TK and protein-inhibitor complexes. Further, MM/PBSA computed the binding free energies of protein complex with Delta 7-Avenasterol, 24-Methylenecholesterol, Campesterol and Ellagic acid as -154.559 ± 18.591 kJ/mol, -139.176 ± 19.236 kJ/mol, -136.212 ± 17.598 kJ/mol and -139.513 ± 23.832 kJ/mol, respectively. Non-polar interactions were the major contributors to these energies. The density functional theory analysis also established the stability of these inhibitor compounds. ADMET analysis depicted acceptable outcomes for all top phytochemicals without displaying any toxicity. In conclusion, this report has identified promising EGFR-TK inhibitors to treat several cancers that can be further investigated through laboratory and clinical tests.

Target-Based Small Molecule Drug Discovery for Colorectal Cancer: A Review of Molecular Pathways and In Silico Studies

Colorectal cancer is one of the most prevalent cancer types. Although there have been breakthroughs in its treatments, a better understanding of the molecular mechanisms and genetic involvement in colorectal cancer will have a substantial role in producing novel and targeted treatments with better safety profiles. In this review, the main molecular pathways and driver genes that are responsible for initiating and propagating the cascade of signaling molecules reaching carcinoma and the aggressive metastatic stages of colorectal cancer were presented. Protein kinases involved in colorectal cancer, as much as other cancers, have seen much focus and committed efforts due to their crucial role in subsidizing, inhibiting, or changing the disease course. Moreover, notable improvements in colorectal cancer treatments with in silico studies and the enhanced selectivity on specific macromolecular targets were discussed. Besides, the selective multi-target agents have been made easier by employing in silico methods in molecular de novo synthesis or target identification and drug repurposing.

Free energy perturbation guided Synthesis with Biological Evaluation of Substituted Quinoline derivatives as small molecule L858R/T790M/C797S mutant EGFR inhibitors targeting resistance in Non-Small Cell Lung Cancer (NSCLC)

Two different schemes of novel substituted quinoline derivatives were designed and synthesized via simple reaction steps and conditions. A comparative molecular docking study was carried out on two different types of EGFR enzymes which include wild-type (PDB: 4I23) and T790M mutated (PDB: 2JIV) respectively. Compounds were also validated upon T790M/C797S mutated (PDB ID: 5D41) EGFR enzyme at the allosteric binding site. Free energy perturbations were carried out to determine the absolute binding free energy of a protein-ligand complex in the form of ΔG_binding, which in turn provided 4ab and 5ad as the most potential contenders through the structural enhancement in the determined initial scaffolds. Anticancer activity of the synthesized derivatives was examined against HCC827, H1975 (L858R/T790M), A549, and HT-29 cell lines by standard MTT assay. Compound 4ad (6-chloro-2-(isoindolin-2-yl)-4-methylquinoline) has shown excellent inhibitory activities against mutant EGFR kinase with IC₅₀ value 0.91 µM. The potency of compounds 4ab, 4ad and 5adwas compared throughan insilicoADMET study.

Computational and Synthetic approach with Biological Evaluation of Substituted Quinoline derivatives as small molecule L858R/T790M/C797S triple mutant EGFR inhibitors targeting resistance in Non-Small Cell Lung Cancer (NSCLC)

New substituted quinoline derivatives were designed and synthesized via a five-step modified Suzuki coupling reaction. A comparative molecular docking study was carried out on two different types of EGFR enzymes which include wild-type (PDB: 4I23) and T790M mutated (PDB: 2JIV) respectively. Compounds were also validated upon T790M/C797S mutated (PDB ID: 5D41) EGFR enzyme at the allosteric binding site. All docking studies confirmed high potency and flexibility towards wild type as well as a mutated enzyme. Anticancer activity of the synthesized derivatives was examined against HCC827, H1975 (L858R/T790M/C797S and L858R/T790M), A549, and HT-29 cell lines by standard MTT assay. Most of the quinoline derivatives revealed a significant cytotoxic effect. The IC50 values of 4-(4-methylquinolin-2-yl)phenyl 4-(chloromethyl)benzoate (5j) were found to be 0.0042 µM, 0.02 µM, 1.91 µM, 3.82 µM and 3.67 µM while IC50 values of osimertinib were 0.0040 µM, 0.02 µM, ND, 0.99 µM and 1.22 µM, respectively. Compound 5j has shownexcellent inhibitory activities against EGFR kinases triple mutant with IC 50 value 1.91 µM. It was observed that, compared to H1975, A549 and A431 cell lines, synthesized compounds significantly inhibited proliferation of the HCC827 cell line. These data suggested that synthesized compounds showed promising selective anticancer activity against tumor cells harboring EGFR Del E746-A750. The potency of compound 5j was compared through molecular dynamic simulations andan insilicoADMET study. QSAR models were generated and the best model was correctly compared with respect to predicted and observed activity of compounds. The built model will assist to design, refine and construct novel substituted quinoline derivatives as potent EGFR inhibitors in near future.

Toward structure-based drug design against the epidermal growth factor receptor (EGFR)

•

Structural variations in EGFR should not be ignored in structure-based drug design.

•

Main variations involve inward and outward folding of C-helix in the kinase N-lobe.

•

Origins of variations are mutations and drug R-groups but not the drug core.

•

Comparative modeling, fitting and clustering are imperative steps in EGFR drug design.

•

Alternatively, volume and shape of binding site can be used to filter ligands against structures.

Most of the available crystal structures of epidermal growth factor receptor (EGFR) kinase domain, bound to drug inhibitors, originated from ligand-based drug design studies. Here, we used variations in 110 crystal structures to assemble eight distinct families highlighting the C-helix orientation in the N-lobe of the EGFR kinase domain. The families shared similar mutational profiles and similarity in the ligand R-groups (chemical composition, geometry, and charge) facing the C-helix, mutation sites, and DFG domain. For structure-based drug design, we recommend a systematic decision-making process for choice of template, guided by appropriate pairwise fitting and clustering before the molecular docking step. Alternatively, the binding site shape/volume can be used to filter and select the compound libraries.