Discovery of Novel EGFR Inhibitor Targeting Wild-Type and Mutant Forms of EGFR: In Silico and In Vitro Study

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 (IC50 of 11.64 ± 1.30 nM), L858R/T790M (IC50 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 IC50 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.


Introduction
EGFR is a transmembrane protein tyrosine kinase that plays an essential role in cellular signaling for cell proliferation, invasion, metastasis, apoptosis, and angiogenesis [1][2][3][4][5][6]. Structurally, EGFR is composed of three domains-a transmembrane domain, an extracellular domain, and an intracellular tyrosine kinase (TK) domain-which are constructed of four important conserved regions: (i) the activation loop (A-loop), (ii) the glycine-rich loop (G-loop), (iii) the catalytic loop (C-loop), and (iv) the hinge region [7,8]. Deregulation of EGFR is caused by increased EGFR activity such as overexpression and mutations of EGFR [9][10][11]. Increased EGFR activity is associated with numerous malignant tumors, including esophageal cancers, glioblastoma, anal cancers, epithelial cancers of the head  , Table 1), a furo [2,3-c]pyridine compound, is a product of a structural modification of vitamin B6, which plays an extremely important biological role in living organisms, since it is involved in a wide range of biochemical reactions necessary for cellular metabolism and function. The latter explains why vitamin B6 and its derivatives are considered biologically privileged molecules with which to develop new therapeutics.  , Table 1), a furo [2,3-c]pyridine compound, is a product of a structural modification of vitamin B 6 , which plays an extremely important biological role in living organisms, since it is involved in a wide range of biochemical reactions necessary for cellular metabolism and function. The latter explains why vitamin B 6 and its derivatives are considered biologically privileged molecules with which to develop new therapeutics.
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3-c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI 50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI 50 -23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC 50 0.93 µM as a CDK2 Molecules 2023, 28, 3014 3 of 20 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type,

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type,

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type,

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and anti-cancer activities against NSCLC cell lines in comparison with normal cell lines were Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and anti-cancer activities against NSCLC cell lines in comparison with normal cell lines were

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and anti-cancer activities against NSCLC cell lines in comparison with normal cell lines were

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro [2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and anti-cancer activities against NSCLC cell lines in comparison with normal cell lines were

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro[2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and anti-cancer activities against NSCLC cell lines in comparison with normal cell lines were

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro[2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and anti-cancer activities against NSCLC cell lines in comparison with normal cell lines were

PD8
Furthermore, PDs display a wide range of biological activities, including anticancer [32,33], antiviral [34], anti-inflammatory [35], and HIV protease inhibitory activities [36], and also inhibit α-glucosidase [37]. The products of the furan cyclization of pyridoxal under the action of hetarylamines 2-alkylamino-3-hetarylamino-4-hydroxymethylfuro[2,3c]pyridines have an immunostimulating effect due to their selective activation of Toll-like TLR8 receptors [38]. Even though anticancer activity of PD has been reported, an understanding of the binding modes of action of PD and its derivatives against NSCLC remains elusive. Hung et al. [33] reported the antiproliferative activity of synthesized thienopyridine derivatives using the panel of NCI-60 cell lines with GI50 values in the low nanomolarity range, especially the melanoma cell line MDA-MD-435 (GI50-23 nM). Rahman et al. also evaluated the CDK2 inhibitory and antiproliferative activity of pyrazolopyridine, furopyridine, and pyridine derivatives. Furopyridine derivative (14) showed IC50 0.93 µM as a CDK2 inhibitor as well as significant inhibition on different human cancer cell lines (HCT-116, MCF-7, HepG2, and A549).
The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and anti-cancer activities against NSCLC cell lines in comparison with normal cell lines were The main aim of this work was the synthesis and evaluation of novel classes of anti-EGFR agents in order to understand the mechanism of anticancer activity as well as their binding modes to the receptor. Herein, we have reported the synthesis of a number of novel furopyridine derivatives with pharmacologically relevant substituents and also used molecular docking and molecular dynamics (MD) simulations with free energy calculations to virtually screen 14 in-house PD compounds towards wild-type, L858R/T790M, and L858R/T790M/C797S EGFR. Moreover, in vitro EGFR inhibitory and anti-cancer activities against NSCLC cell lines in comparison with normal cell lines were studied. Subsequently, 500-ns MD simulations were carried out to investigate the binding behavior of the most promising PD against wild-type and mutant forms of EGFR. The most promising PD from this work could be further developed as a novel anticancer drug by targeting both wild-type and mutant forms of EGFR.

Virtual Screening
To assess the potency of EGFR inhibitor(s) against wild-type, L858R/T790M, and L858R/T790M/C797S EGFR, molecular docking was carried out on 14 PD compounds.
Three drugs-erlotinib, afatinib, and Osimertinib-were selected as a positive control for wild-type, L858R/T790M, and L858R/T790M/C797S EGFR, respectively. The GOLD docking results were plotted in Figure 2. We found that nine PD compounds, including PD3, PD4, PD6, PD8, and PD10-PD13, showed a higher fitness score than erlotinib. Moreover, those PDs were found to inhibit L858R/T790M double mutant EGFR. Interestingly, a series of four PDs (PD4, PD8, PD13, and PD14) exhibited a strong affinity with L858R/T790M/C797S, which is better than osimertinib. The first-round screening suggested that several PD derivatives could be effective in inhibiting both wild-type and mutant forms of EGFR.

Virtual Screening
To assess the potency of EGFR inhibitor(s) against wild-type, L858R/T L858R/T790M/C797S EGFR, molecular docking was carried out on 14 PD co Three drugs -erlotinib, afatinib, and Osimertinib-were selected as a positive wild-type, L858R/T790M, and L858R/T790M/C797S EGFR, respectively. The GO ing results were plotted in Figure 2. We found that nine PD compounds, inclu PD4, PD6, PD8, and PD10-PD13, showed a higher fitness score than erlotinib. those PDs were found to inhibit L858R/T790M double mutant EGFR. Interest ries of four PDs (PD4, PD8, PD13, and PD14) exhibited a strong aff L858R/T790M/C797S, which is better than osimertinib. The first-round scre gested that several PD derivatives could be effective in inhibiting both wild-typ tant forms of EGFR. To verify the first-round screening, the screened PDs/EGFRs complexes de molecular docking were subsequently studied on 100 snapshots extracted from  To verify the first-round screening, the screened PDs/EGFRs complexes derived from molecular docking were subsequently studied on 100 snapshots extracted from the last 10 ns of 100-ns molecular dynamics (MD) simulations. The predicted binding affinity based on the solvated interaction energy (SIE) method of screened PD compounds and three positive controls (erlotinib, afatinib, and osimertinib) is shown in Figure 2. PD6 (∆G bind of −9.96 kcal/mol) and PD11 (−11.30 kcal/mol) displayed greater binding affinity than erlotinib (−8.50 kcal/mol) in the wild-type EGFR system, suggesting that both PD6 and PD11 were specific for the wild-type EGFR. Moreover, the three PDs (PD3, PD4, and PD12) exhibited the most efficient binding to both wild-type and L858R/T790M EGFR, whereas that against L858R/T790M/C797S EGFR showed poor binding affinities. Surprisingly, the most promising, PD13, could be able to suppress both wild-type and mutant EGFR. The ∆G bind calculations of PD13 based on the SIE method gave ∆G bind of wild-type, L858R/T790M, and L858R/T790M/C797S EGFR as follows: −11.81, −9.39 and −9.7 kcal/mol. Hence, the screened compounds with estimated binding free energies based on the SIE method that were higher than those of the positive controls from each strain of EGFR were synthesized and chosen for in vitro evaluation.

Chemistry
All the target compounds were synthesized as outlined in Scheme 1.
To verify the first-round screening, the screened PDs/EGFRs complexes derived from molecular docking were subsequently studied on 100 snapshots extracted from the last 10 ns of 100-ns molecular dynamics (MD) simulations. The predicted binding affinity based on the solvated interaction energy (SIE) method of screened PD compounds and three positive controls (erlotinib, afatinib, and osimertinib) is shown in Figure 2. PD6 (ΔGbind of −9.96 kcal/mol) and PD11 (−11.30 kcal/mol) displayed greater binding affinity than erlotinib (−8.50 kcal/mol) in the wild-type EGFR system, suggesting that both PD6 and PD11 were specific for the wild-type EGFR. Moreover, the three PDs (PD3, PD4, and PD12) exhibited the most efficient binding to both wild-type and L858R/T790M EGFR, whereas that against L858R/T790M/C797S EGFR showed poor binding affinities. Surprisingly, the most promising, PD13, could be able to suppress both wild-type and mutant EGFR. The ΔGbind calculations of PD13 based on the SIE method gave ΔGbind of wild-type, L858R/T790M, and L858R/T790M/C797S EGFR as follows: −11.81, −9.39 and −9.7 kcal/mol. Hence, the screened compounds with estimated binding free energies based on the SIE method that were higher than those of the positive controls from each strain of EGFR were synthesized and chosen for in vitro evaluation.

Chemistry
All the target compounds were synthesized as outlined in Scheme 1. The cyclization of pyridoxal hydrochloride 1 in aqueous acetonitrile at 50-60 • C used potassium carbonate as a catalyst soft base, and proceeded fairly quickly in all cases-within 2-3 h-and led to the synthesis of 4-hydroxymethyl-7-methylfuro [2,3-c]pyridines (PD1, PD2). The yields of compounds PD1, PD2 were 55% and 79%, respectively.
The R 1 and R 2 substituents were chosen considering that the biological activity of the compounds, as well as their pharmacodynamics and pharmacokinetics, are significantly affected by both electronic factors (electron donation or electron acceptor of substituents) and their physical properties (lipophilicity, hydrophilicity, polarizability, etc.) The structure of the obtained furo [2,3-c]pyridines was confirmed by 1 H NMR spectroscopy. 1 H NMR data for compounds PD1, PD2 can be illustrated by the example of 4-hydroxymethylfuro[2,3-c]pyridines PD1. The 1 H NMR of this ketone showed the proton signals of H-3 and H-5 of the bicyclic system (singlets at 7.79 and 8.24 ppm, respectively), the benzoyl group (two-proton, somewhat distorted doublet at 8.10 ppm (H-2 , H-6 , J 8.7 Hz) as well as multiplate at 7.08-7.14 (2H, H-3 , H-5 )), the hydroxymethyl group (two-proton and single-proton singlets at 4.79 and 5.24 ppm, respectively), and the methyl group (singlet at 2.76 ppm) (see Supplementary Files, 1 H-NMR spectra) (SI).
In the spectra of amino derivatives PD3-PD14, there was no signal of the CH 2 OH group, while in the hydrochlorides of these compounds (PD3-PD5, PD9, PD11-PD14), a broadened singlet of the NH + group appeared in the region of 5.01-5.71 ppm. The aromatic proton signals of these compounds (PD3-PD14) correspond to the proposed structures and differ from the signals of 4-hydroxymethylfuro[2,3-c]pyridine (PD1) by 01-03 ppm (see Supplementary Files, 1 H-NMR spectra). The 13 C-NMR of all these compounds was in agreement with their structure (see Supplementary Materials, 13 C-NMR).

Kinase Inhibitory Activities of EGFR
Potent compounds obtained from in silico screening were selected to elucidate their inhibitory EGFR activity against different types of EGFR (wild-type, L858R/T790M and L858R/T790M/C797S). Note that erlotinib, afatinib, and osimertinib, which are FDAapproved drugs targeting EGFR, were used as the positive control. The selected compounds and drugs were first screened to observe the sign of EGFR using a single 1 µM dose concentration. The PD compounds that present a percentage value of EGFRs inhibition <50 were selected to calculate the IC 50 values. The relative inhibition of EGFR and IC 50 curves were analyzed and are plotted in Figure 3, and the IC 50 values of focused PD and drugs are summarized in Table 2. It was found that PD13 was active in wild-type and L858R/T790M forms of EGFR, whereas PD4 was only against wild-type EGFR.
8, x FOR PEER REVIEW 7 of 20 Figure 3. Three forms of EGFR (wild-type, L858R/T790M, and L858R/T790M/C797S) inhibitory screening of PD derivatives at 1 µM concentration (top) and IC50 curves of compounds that present more than 50% PD4 and PD13 in comparison with reference drugs (bottom).

Antiproliferative Activities on NSCLC and Normal Cell Lines
The antiproliferative activity of two focused PD compounds and drugs were evaluated against NSCLC cell lines, namely the wild-type NSCLC cell line A549 as well as L858R/T790M mutant NSCLC cell line H1975. Erlotinib, afatinib, and osimertinib were used as positive controls. The results are summarized in Table 3. We found that both PD4 and PD13 exhibited strong anti-cancer activity by inhibiting the proliferation of A549 rather than H1975 cell lines. In addition, PD13 exhibited stronger antiproliferative activities than afatinib and osimertinib against the A549 cell line. Moreover, PD13 showed cell growth inhibition in H1975 cell cultures at a similar level to afatinib. Furthermore, PD4, PD13, and positive controls were then investigated for their cytotoxicity to the normal cell line Vero (monkey kidney epithelial cells). We found that the IC 50 value of PD13 was >100 µM in Vero cells, whereas for PD4 IC 50 was <100 µM. These results suggested that PD13 showed low cytotoxicity against the normal cell. Therefore, PD13 was chosen in order to study its binding affinity by MD simulations.

System Stability
The stability of PD13 and the reference drugs erlotinib, afatinib, and osimertinib was determined using all-atom root mean square displacement (RMSD) to observe the system stability. The RMSD values of each system were obtained from three independent simulations. As shown in Figure 4, the RMSD values of the erlotinib system dramatically increased in the first 100 ns and were maintained at the fluctuation of~2.0-3.0 Å until the end of the simulation for all independent runs. The wild-type/PD13 complexes reached three runs after 300 ns. The RMSD values of afatinib complexed with L858R/T790M EGFR fluctuated 1.5-3.0 Å and were maintained at the fluctuation of~2.0-3.0 Å until the end of the simulation. In the case of PD13/double mutant, EGFR showed similar patterns to the three independent simulations. The RMSD value of osimertinib continuously increased in the first 100 ns and was stably maintained at a fluctuation of 3.0-4.0 Å. L858R/T790M/C797S complexed with PD13 appeared to be quite stable after 400 ns for all simulated systems, while the RMSD value of PD13 and known drugs complexed with EGFRs of all systems tended to be stable after 200 ns. Altogether, the MD trajectories from 400 to 500 ns of each simulation were extracted for further analysis in terms of: (i) key binding residues; (ii) protein-ligand hydrogen bonding; and (iii) binding affinity. Molecules 2023, 28, x FOR PEER REVIEW 9 of 20

Key Binding Residues
To explore the binding affinity of PD13 and reference drugs, the per-residue free energy decomposition (∆G ) based on the MM/PBSA method was calculated on 100 snapshots extracted from the last 100 ns of all simulations. Note that, among the 695-1018 residues of EGFR, only the energy stabilization of ≤−1.0 kcal/mol for the 700-900 residues of all systems was plotted in Figure 5. The compounds were shaded based on their highest and lowest energies, with white representing the highest energies and blue representing the lowest energies.
The wild-type/erlotinib showed important residue for binding ten residues including L718, V726, A743, K745, T790, L792, M793, G796, C797, and L844. Interestingly, the PD13/wild-type complex displayed a similar binding pattern to that of reported EGFR inhibitors and was stabilized by L718 (−3.74 kcal/mol, dark red) at the hinge region, which is in good agreement with a previous study that identified L718, A743, L792, M793, G796, and L844 as key residues involved in the main binding of EGFR inhibitors [43][44][45]. In addition, the residue M793 of L858R/T790M EGFR exhibited the highest stabilization for PD13 (−3.43 kcal/mol, dark red) compared to the afatinib (−3.21 kcal/mol, red) system. This may be one of the reasons that explains why PD13 could inhibit L858R/T790M EGFR better than afatinib. The binding of PD13 with L858R/T790M/C797S EGFR displayed a higher number of amino acid residues than Osimertinib; L718, V726, A743, K745, M766, L788, M790, L792, M793, G796, S797, and L844 are all involved in ligand stabilization. Among all the residues, the moderate L858R/T790M/C797S EGFR inhibitory potency of PD13 may be due to steric clashing of M790 with the phenyl ring of the inhibitor.

Key Binding Residues
To explore the binding affinity of PD13 and reference drugs, the per-residue free energy decomposition (∆G residue bind ) based on the MM/PBSA method was calculated on 100 snapshots extracted from the last 100 ns of all simulations. Note that, among the 695-1018 residues of EGFR, only the energy stabilization of ≤−1.0 kcal/mol for the 700-900 residues of all systems was plotted in Figure 5. The compounds were shaded based on their highest and lowest energies, with white representing the highest energies and blue representing the lowest energies.
The wild-type/erlotinib showed important residue for binding ten residues including L718, V726, A743, K745, T790, L792, M793, G796, C797, and L844. Interestingly, the PD13/wild-type complex displayed a similar binding pattern to that of reported EGFR inhibitors and was stabilized by L718 (−3.74 kcal/mol, dark red) at the hinge region, which is in good agreement with a previous study that identified L718, A743, L792, M793, G796, and L844 as key residues involved in the main binding of EGFR inhibitors [43][44][45]. In addition, the residue M793 of L858R/T790M EGFR exhibited the highest stabilization for PD13 (−3.43 kcal/mol, dark red) compared to the afatinib (−3.21 kcal/mol, red) system. This may be one of the reasons that explains why PD13 could inhibit L858R/T790M EGFR better than afatinib. The binding of PD13 with L858R/T790M/C797S EGFR displayed a higher number of amino acid residues than Osimertinib; L718, V726, A743, K745, M766, L788, M790, L792, M793, G796, S797, and L844 are all involved in ligand stabilization. Among all the residues, the moderate L858R/T790M/C797S EGFR inhibitory potency of PD13 may be due to steric clashing of M790 with the phenyl ring of the inhibitor.

Protein-Ligand Hydrogen Bonding
The hydrogen bond is one of the main factors for determining the binding strength of the protein-ligand complex. We further calculated the hydrogen bond occupations between EGFRs with PD13 and reference drugs using the two following criteria: (i) distance between hydrogen donor (HD) and acceptor (HA) ≤ 3.5 Å; and (ii) the angle between HD and HA being over 120° (HD−H⋯HA ≥ 120°). The percentages of hydrogen bond occupations for the studied ligands were observed and are compared in Figure 6.

Protein-Ligand Hydrogen Bonding
The hydrogen bond is one of the main factors for determining the binding strength of the protein-ligand complex. We further calculated the hydrogen bond occupations between EGFRs with PD13 and reference drugs using the two following criteria: (i) distance between hydrogen donor (HD) and acceptor (HA) ≤ 3.5 Å; and (ii) the angle between HD and HA being over 120 • (HD−H· · · HA ≥ 120 • ). The percentages of hydrogen bond occupations for the studied ligands were observed and are compared in Figure 6. The hydrogen bond formations of wild-type/PD13 indicated strong stabilization with the M793 (78-98%) and L718 (24%), and are similar to that found in wild-type/erlotinib complex. In the case of L858R/T790M double EGFR, this compound formed a hydrogen bond with M793 (67-95%) and L178 (42%), in good agreement with the per-residue decomposition analysis mentioned above. PD13 complexed with triple mutant EGFR showed similar H-bond formations to those of other forms of EGFR (except S797). We speculate that S797 residue that is mutated EGFR could lead to moderate inhibitory activity of L858R/T790M/C797S EGFR due to more favorable interactions with S797 via a hydrogen bond. Therefore, the inhibitory activity of PD13/EGFRs was induced by the two key residues, which are M793 and S797, at the hinge region. The main EGFR residues (L718, M793, and S797) evidenced crucial binding via hydrogen bonds, which is similar to other reported EGFR tyrosine kinases such as tupichinol E [46], quinazoline [47], and 9-heterocyclyl substituted 9H-purine derivatives [48].

Binding Affinity
The binding efficiency of PD13 and reference drugs was calculated using the MM/PBSA approach on 100 snapshots extracted from the last 100 ns of three independent simulations. The values of the averaged ∆Gbind over the three independent simulations are listed in Table 4 in comparison with the half-maximal inhibitory concentration (IC50). The hydrogen bond formations of wild-type/PD13 indicated strong stabilization with the M793 (78-98%) and L718 (24%), and are similar to that found in wild-type/erlotinib complex. In the case of L858R/T790M double EGFR, this compound formed a hydrogen bond with M793 (67-95%) and L178 (42%), in good agreement with the per-residue decomposition analysis mentioned above. PD13 complexed with triple mutant EGFR showed similar H-bond formations to those of other forms of EGFR (except S797). We speculate that S797 residue that is mutated EGFR could lead to moderate inhibitory activity of L858R/T790M/C797S EGFR due to more favorable interactions with S797 via a hydrogen bond. Therefore, the inhibitory activity of PD13/EGFRs was induced by the two key residues, which are M793 and S797, at the hinge region. The main EGFR residues (L718, M793, and S797) evidenced crucial binding via hydrogen bonds, which is similar to other reported EGFR tyrosine kinases such as tupichinol E [46], quinazoline [47], and 9-heterocyclyl substituted 9H-purine derivatives [48].

Binding Affinity
The binding efficiency of PD13 and reference drugs was calculated using the MM/PBSA approach on 100 snapshots extracted from the last 100 ns of three independent simulations. The values of the averaged ∆G bind over the three independent simulations are listed in Table 4 in comparison with the half-maximal inhibitory concentration (IC 50 ). In the gas term, we found that vdW interactions were three-to fourfold higher than electrostatic interactions in all studied ligands. Therefore, vdW interactions play an important role in the recognition of PD13 and reference drugs. These interaction findings are strongly consistent with previous studies on several small-molecule therapeutic kinase inhibitors such as olmutinib, lapatinib, and icotinib [49][50][51]. From the ∆G bind calculation, the predicted binding values produced similar phenomena to that of experimental ∆G bind values converted from the IC 50 values, indicating that the MM/PBSA approach produced reliable ∆G bind results. It was found that PD13 complexed with the wild-type and L858R/T790M system had a better binding affinity than erlotinib or afatinib, respectively. While PD13 complexed with L858R/T790M/C797S displayed lower binding affinity than osimertinib, the results agree well with in vitro L858R/T790M/C797S EGFR inhibitory activities, showing that this compound could inhibit wild-type and L858R/T790M better than L858R/T790M/C797S EGFR.

General Information
NMR 1 H (300 MHz) and 13 C (125 MHz) spectra of newly synthesized compounds were recorded on a spectrometer Bruker AC-300 (300 MHz) and BRUKER AVANCE NEO 500 MHz FT spectrometers in DMSO-d6. Chemical shifts of nuclei 1 H were measured relative to the residual signals of deuterosolvent (δ = 2.50 ppm). Coupling constants (J) are reported in Hz. Melting points were determined by using Fisher-Johns Melting Point Apparatus (Fisher Scientific) and are uncorrected. Elemental analysis was performed using the classical method of microanalysis. The reaction and purity of the obtained compounds were monitored by TLC (plates with Al 2 O 3 III activity grade, eluent CHCl 3 , development of TLC plates by exposition to iodine vapors in "iodine chamber").

Chemistry
The reaction conditions and yields of compounds PD3-PD14 are presented in Table 5. Compound PD2 was synthesized as described previously [39].

System Preparation
The three-dimensional structure of wild-type (PDB ID: 1M17) [52], L858R/T790M (PDB ID: 4I22) [53] and L858R/T790M/C797S EGFR (PDB ID: 6LUD) [52] were obtained from the Protein Data Bank. These crystal structures complexed with known drugs and subsequently known drugs were extracted. A series of 14 furo[2,3-c]pyridine derivatives was built using the Gaussview 09 program. All the ligands were optimized using the Gaussian 09 program (HF/6-31G*) and, subsequently, the protonation states of all ionizable amino acids were characterized using ChemAxon at pH 7.0. Molecular docking studies were performed using the GOLD program. The setting of the docking was studied using the following parameter as 10 Å with a sphere and 100 docking poses. The docking runs in each complex were sorted and selected based on GOLD fitness score. Then, the docking results were visualized for interaction using the UCSF Chimera package. The docked PDs/EGFRs with the highest GOLD fitness scores were chosen as the initial structures for performing MD simulations.

Molecular Dynamics Simulation and Free Energy Calculation Based on Solvated Interaction Energy (SIE) Method
The initial structures of PDs/EGFRs derived from molecular docking were carried out using pmemd CUDA in the AMBER 16 program [54] under periodic boundary conditions with the isobaric isothermal (NPT) ensemble in triplicate. The general AMBER force fields GAFF [55] and FF14SB [56] were applied for the ligand and protein, respectively. The partial charges of ligand were generated as follows: (i) the 3D structure of PD derivatives were optimized using the HF/6-31G* method [57] in the Gaussian09 program; (ii) the electrostatic potential (ESP) charge and restrained ESP (RESP) charge of the ligand were obtained by the AMBER16 program. All studied complexes were solvated by the TIP3P water model [58]. All missing hydrogen atoms of the protein and ligand were added by the LEaP module and subsequently minimized. The sodium ions were randomly neutralized in the simulated systems by Cl − counter ions. The added hydrogen atoms and water molecules were then minimized using 2500 steps of the steepest descents (SD), followed by 2500 steps of the conjugated gradient (CG), before starting the MD simulations. The SHAKE algorithm was used to fix all covalent bonds involving hydrogen atoms [59]. The systems were heated up to 310 K for 100 ps. Finally, the whole system was performed under the NPT ensemble (310 K, 1 atm) until reaching 100 ns. The root mean square deviation (RMSD) of each system was calculated using all-atom, and then the MD trajectories from the last 10-ns simulation were extracted for analysis in terms of total binding free energy (∆G bind ) based on the solvated interaction energy (SIE) method. The potent PD derivatives with the lowest binding free energy were selected to determine EGFR inhibitory as well as anticancer activity. After finding the more potent PD derivative obtained from in silico screening and in vitro evaluation, the 500-ns MD simulation was performed. In addition, the per-residue decomposition free energy (∆G bind,residue ) calculations were carried out using the MMPBSA.py module in AMBER16 [60]. The CPPTRJ module was used to calculate intermolecular hydrogen bonding.

EGFR Inhibitory ACTIVITY
The IC 50 values of the potent PD derivatives obtained from the SIE method were measured using the ADP-Glo TM kinase assay kit as previously reported [19]. Firstly, 8 µL of buffer (40 mM Tris-HCl pH 7.5, 20 mM MgCl 2 , and 0.1 mg/mL bovine serum albumin) were added to a 384-well plate (Promega, solid white). Secondly, 5 µL of 1.25 ng/µL EGFRs enzyme (wild-type, L858R/T790M, and L858R/T790M/C797S EGFR) and 2 µL of 1 µM of compounds/known drugs were added separately to individual cell cultures, followed by 10 µL of a mixture of 5 µM ATP and 2.5 µM poly(glu-tyr), and were incubated for 1 h at room temperature. Thirdly, 5 µL of the ADP-Glo reagent was added and incubated for 40 min.
Finally, 10 µL of kinase detection reagent was added and incubated at room temperature for 30 min and was detected by measuring the luminescence using a microplate reader (Infinite M200 microplate reader, Tecan, Männedorf, Switzerland). A triplicate assay was performed on all potent PD derivatives and known drugs. The percentage relative inhibition (%) of potent compounds was measured compared to the control with no inhibitor as shown in Equation (1):

Cell Viability Assay
The cytotoxicity activities of potent PD derivatives obtained from the SIE method against A549 and H1975, as well as Vero (which are adherent cells), were assessed using MTT (3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2H-tetrazolium bromide) assays. The A549 and Vero cells were grown in complete DMEM medium whereas H1975 cells were grown in complete RPMI medium supplemented with 10% (v/v) FBS, 100 U/mL penicillin and 100 µg/mL streptomycin. All cells were routinely cultured at 37 • C in a 5% (v/v) CO 2 , 95% (v/v) air humidified incubator. A549 (5 × 10 3 cells per well), H1975 (5 × 10 3 cells per well), and Vero (2 × 10 3 cells per well) cells were seeded into 96-well plates and incubated overnight. After cell attachment, cells were replenished using fresh medium with 0.1% FBS containing various concentrations of potent PD derivatives and known drugs and then incubated at 37 • C for 72 h. After that, fresh medium containing MTT solution (5 mg/mL) was added and incubated at 37 • C. After 3 h, the formed formazan crystal was dissolved with 100 µL of dimethyl sulfoxide (DMSO). Finally, the absorbance of formazan solution was measured at a wavelength 570 nm using a microplate reader (Infinite M200 microplate reader, Tecan, Männedorf, Switzerland). The independent treatments were repeated in triplicate. The IC 50 determination of compounds was achieved using GraphPad Prism version 9.0.

Conclusions
In this work, a combination of in silico and in vitro studies were performed to screen for novel EGFR inhibitors that could inhibit both wild-type and mutant forms of EGFR. We found that potent compounds from virtual screening exhibited promising EGFR inhibitory activity. In particular, PD13 showed strong inhibition of wild-type and L858R/T790M double mutant EGFR, while L858R/T790M/C797S showed moderate inhibition. In addition, the most potent PD13 revealed high cytotoxicity against A549 and H1975. In 500-ns MD simulations, the model of PD13 and known drugs/EGFRs was stable. The binding affinity of PD13 complexed with wild-type and mutant forms of EGFR was significantly higher than that of known drugs (except for the L858R/T790M/C797S system), which was consistent with the hot-spot residues and EGFR inhibitory activity. The ∆G residue bind calculations revealed that M793 and S797 are the main residues in the hinge region for the binding of PD13 with three EGFR strains via H-bond formations. In addition, the binding of all studied PD compounds is driven mainly by vdW interaction. Our results showed that PD13 could potentially act as a promising non-mutant and mutant EGFR drug for cancer therapy.