High throughput virtual screening and validation of Plant-Based EGFR L858R kinase inhibitors against Non-Small cell lung Cancer: An integrated approach Utilizing GC–MS, network Pharmacology, Docking, and molecular dynamics

Abstract

Lung cancer ranks as the 2nd most common cancer globally. It’s the most prevalent cancer in men and the 2nd most common in women. The prominent events in EGFR-mutated non-small-cell lung cancer (NSCLC) include the emergence of the L858R mutation within EGFR exon 21. Despite the promising efficacy of EGFR inhibitors in managing lung cancer, the development of acquired resistance poses a significant hurdle. In the current investigation, we focused on the screening of two phytochemicals, namely Dehydrocostus lactone and Mokkolactone, derived from the Saussurea lappa plant, as potential inhibitors targeting EGFR L858R mutant lung cancer. The chloroform and ethanol extract of the plant demonstrated anti-proliferative activity through the Resazurin chemosensitivity assay, exhibiting an IC50 value of 37.90 ± 0.29 µg/ml with selectivity index 2.4. Through a GC–MS study, we identified 11 phytochemicals for further insilico analysis. These compounds underwent ADMET assessment followed by drug likeliness analysis before being subjected to molecular docking against EGFR L858R, identified through protein–protein interaction network analysis. All phytochemicals exhibited binding energy scores ranging from −6.9 to −8.1 kcal/mol. Dehydrocostus lactone and Mokkolactone were specifically identified for their binding profile. Findings from 100 ns molecular dynamics simulations demonstrated their enhanced stability compared to the reference ligand DJK. This was evident in the root mean square deviation (RMSD) values, ranging from 0.23 ± 0.01 nm to 0.30 ± 0.05 nm, the radius of gyration values, from 1.71 ± 0.01 nm to 1.72 ± 0.01 nm, and the solvent accessible surface area values, from 155.39 ± 2.40 nm² to 159.32 ± 2.14 nm². Additionally, favourable characteristics were observed in terms of hydrogen bonding, principal component analysis, and free energy landscape analysis. Examination of their electronic structure via density functional theory revealed efficient properties, with the highest occupied molecular orbital-least unoccupied molecular orbital energy gap values ranging from −3.984 eV to −6.547 eV. Further, in vivo analysis is required to gain a more comprehensive understanding and efficacy of these identified phytochemicals against lung cancer.