Average geometrical values (Å) over the last 50 ns of 100-ns-long MD simulations.
The selective α1-adrenergic receptor antagonist doxazosin is used for the treatment of hypertension. More recently, an experimental report demonstrated that this compound exhibits antiproliferative activity in breast cancer cell lines with similar inhibitory activity to gefitinib, a selective inhibitor of EGFR in the active state (EGFRAC). This experimental study provided evidence that doxazosin can be employed as an anticancer compound, however, the structural basis for its inhibitory properties is poorly understood at the atomic level. To gain insight about this molecule, molecular dynamics (MD) simulation with the molecular mechanics generalized Born surface area (MMGBSA) approach was employed to explore the structural and energetic features that guide the inhibitory properties of doxazosin and gefitinib in overexpressing EGFR/HER2 cell lines. Our result suggest that doxazosin exerts its inhibitory properties in breast cancer cell lines by targeting EGFR/HER2 but mainly HER2 in the inactive state (HER2IN), whereas gefitinib by targeting mainly EGFRAC, in line with previous literature. Decomposition of the binding affinity into individual contributions of HER2IN-doxazosin and EGFRAC-gefitinib systems detected hot spot residues but also showed polar interactions of Met801/Met793 with the quinazoline ring of both compounds. Principal component (PC) analysis revealed that the molecular recognition of the HER2IN-doxazosin system was linked to conformational changes but EGFRAC-gefitinib was not.
- MD simulations
Human epidermal growth factor receptor 1 (EGFR) and 2 (HER2) form part of a family of human epidermal grown factor receptors (EGFRs), and whose phosphorylation impacts cell proliferation, differentiation, and migration . The cytoplasmic tyrosine kinase domain (TKD) is considered one of the most studied receptors for developing new anticancer drugs . Activation of EGFR starts the molecular recognition of endogenous growth factors at the extracellular domain that, at the time, promotes the formation of homo- and heterodimers among the different members of EGFR, with HER2 the preferred member of EGFRs to form heterodimers [3, 4, 5]. The transition from a monomeric to dimeric state in EGFR is coupled to a conformational change in the TKD from an inactive to active state [6, 7, 8], whereas that for HER2 transitions from inactive, intermediate, and inactive states [9, 10, 11]. Generally, the signaling activity regulated by EGFR/HER2 is under control, however, mutations in TKD give place to constitutive activation of these receptors, which results in the development of different types of cancer, such as lung  and breast cancer . In addition, overexpression of EGFR/HER2 also happens with radiotherapy and chemotherapy resistance [14, 15, 16].
Based on the ability of tyrosine kinase inhibitors (TKIs) to inhibit EGFR, they can be divided in two types, those targeting the active state, such as Iressa, and those targeting the inactive state, such as erlotinib and lapatinib [17, 18, 19, 20]. Lapatinib showed dual activity on EGFR/HER2 [21, 22, 23, 24, 25, 26]. Despite the benefits of using these TKIs, the employment of them has been linked to severe side effects and drug resistance [27, 28, 29, 30]. Therefore, it is necessary to identify new compounds, either through drug design or drug repurposing, that target EGFR and/or HER2 receptors and are effective for cancer therapy. In this context, the combination of docking and molecular dynamics (MD) simulations has been widely exploited to generate new information about the binding properties between natural or synthetic TKIs and EGFR/HER2 [10, 11, 19, 20, 31, 32, 33, 34, 35, 36]. In a previous study, Hui et al. explored the inhibitory properties of doxazosin, an α-1 antagonist used for the treatment of hypertension, in two human breast cancer cell lines: BCC MDA-MB-231 and MCF-7 cells . MDA-MB-231 and MCF-7 cells are estrogen receptor (ER) positive and ER negative, respectively , and both cell lines also expressed EGFR and HER2; however, MDA-MB-231 expressed both receptors in higher concentrations than MCF-7 . Although EGFR and HER2 are important regulators for normal cellular processes, their dysregulation has been associated to protein overexpression that leads to the development of different types of cancer [1, 5]. They demonstrated that doxazosin induces apoptosis in breast cancer cell lines similar to Iressa (Gefitinib), reducing phosphorylated EGFR by a mechanism that does not involve the α1-adrenergic receptor, however, the structural and energetic basis for its inhibitory properties is poorly understood. In addition, Sharkawi et al. identified similar experimental antiproliferative activity of doxazosin in an MCF-7 cell line through the inhibition of EGFR . Thus, more robust structural and energetic studies are required to provide structural insight into the affinity of doxazosin for EGFR/HER2 compared with gefitinib. Structural data, docking, and molecular dynamics (MD) simulations combined with the MMGBSA approach were used to elucidate the molecular mechanism through which doxazosin and gefitinib inhibit EGFR/HER2.
2.1 Structural modeling
The free forms of EGFR in the inactive (EGFRIN) and active (EGFRAC) states were taken from the crystallographic structures of EGFRIN (PDB entry 1XKK) and EGFRAC (PDB entry 1 M17) conformations. The free forms of HER2 in the inactive (HER2IN) and active (HER2AC) states were taken from previous MD simulation studies; HER2IN  and HER2AC  conformations. Amino acid residues missing in the electron density map of EGFR structures were built with MODELER Version 9.14 .
2.2 Docking studies
Docking calculations were carried out using AutoDock 4.2 and AutoDock Tools 1.5.6 software . The ligand structures were built and optimized with the Gaussian package . The initial geometries of ligands were optimized at the AM1 level. Hydrogen atoms were added to ligands and receptors, and Kollman and Gasteiger partial charges were assigned for ligand and proteins, respectively. The affinity grid maps were constructed on the receptor using a grid size of 70 × 70 × 70 Å and 0.370 Å of spacing. Due to the stochastic nature of the Lamarckian algorithm, 20 runs were performed for each compound, and 30 conformations of the ligand (binding poses) were observed between ligand and protein. The best binding poses were selected using the criteria of having the lowest energetic conformations at the receptor binding site.
2.3 Molecular dynamics simulations
The protein-ligand results obtained by docking were checked through MD simulation studies. MD simulations were carried out using the AMBER16 package , in conjunction with the ff14SB force field . The systems simulated were put into a space-filling dodecahedric box of 12 Å, solvated with TIP3P water model , and neutralized with sodium and chloride ions (0.10 M) to create a physiological concentration. The parameterizations of the ligands were performed assigning AM1-BCC atomic charges and matching the atoms with the general Amber force field (GAFF) . Once the systems were constructed, they were minimized using steepest descent with position restraint of the ligands, followed by steepest descent without position restraint and conjugate gradients. The minimized systems were then submitted to 100 ns-long MD simulations using an NPT ensemble with the velocity rescaling arrangement to simulate a constant temperature at 310 K. A constant temperature and pressure (1 atm) were maintained using the weak-coupling algorithm , with coupling constants τT and τP of 1.0 and 0.2 ps, respectively. The electrostatic term was described by the PME method , and a 10 Å cut-off was selected for the van der Waals interactions. The time step for the MD simulations was set to 2.0 fs. The SHAKE algorithm  was employed to reset bonds to their right lengths after an unconstrained update. The conformations obtained from MD simulations at intervals of 20 picoseconds (ps) were analyzed using the cpptraj tool in Amber16. Plots of variation of root mean squared deviation (RMSD) and radius of gyration (RG) were generated to evaluate convergence. Clustering analysis using a cutoff of 2.5 Å was performed to identify the most populated conformation in the simulation. Principal components (PC) analysis along the most essential eigenvectors was carried out to evaluate total flexibility. A map of interactions was generated using Maestro Version 10.1, 2015–1 .
2.4 Affinity prediction and per-residue decomposition
The binding free energy (Δ
3. Results and discussion
3.1 Convergence and equilibrium
The stability of the evaluated systems was observed by measuring two geometrical parameters. The root mean squared deviation (RMSD) and the radius of gyration (RG) were determined to identify the time at which the systems reached convergence (Table 1). RMSD analysis showed that free and bound EGFR/HER2 systems reached stability between 20 to 50 ns with RMSD values, which oscillated between 1.40 and 4.20 Å. RG examination revealed that free and bound EGFR/HER2 systems exhibited stability from 20 to 50 ns with values oscillating between 18.8 and 20.2 Å. Based on this result, further analysis was carried out discarding the first 50 ns.
|EGFRAC||2.1 ± 0.20||19.0 ± 0.12|
|EGFRAC-doxazosin||1.7 ± 0.10||19.1 ± 0.10|
|EGFRAC-gefitinib||2.2 ± 0.20||19.2 ± 0.10|
|EGFRIN||1.8 ± 0.20||18.8 ± 0.10|
|EGFRIN-doxazosin||2.2 ± 0.20||19.0 ± 0.10|
|EGFRIN-gefitinib||2.7 ± 0.21||19.0 ± 0.10|
|HER2AC||3.6 ± 0.17||20.0 ± 0.10|
|HER2AC-doxazosin||1.7 ± 0.20||20.0 ± 0.14|
|HER2AC-gefitinib||1.4 ± 0.20||19.9 ± 0.10|
|HER2IN||3.9 ± 0.40||20.0 ± 0.01|
|HER2IN-doxazosin||3.4 ± 0.40||20.2 ± 0.10|
|HER2IN-gefitinib||4.2 ± 0.22||19.6 ± 0.13|
3.2 Structural analysis of complexes between doxazosin and EGFRAC/EGFRIN
To explore the structural differences between doxazosin and gefitinib on EGFR/HER2, the most populated receptor-ligand conformations were retrieved over the equilibrated simulation time (the last 50 ns) through clustering analysis. Analysis of the complex between doxazosin and EGFRAC showed that the ligand was stabilized through van der Waals interactions with Phe723, Val726, Leu792, Met793, Pro794, Phe795, Cys797, Leu799, and Leu844, and polar interactions with Gly719, Lys745, Gly796, Asp800, Asp837, Arg841, Asn842, and Asp855 (Figure 1A). In contrast, Val718 formed both van der Waals interactions and a hydrogen bond with the quinazoline ring of doxazosin. In the complex with EGFRIN, dozaxosin was bound through van der Waals interactions with Phe723, Val718, Val726, Ala743, Leu792, Phe795, Tyr801, and Leu844. The polar contact took place with Lys745, Gln791, Gly796, Glu804, and His805. Met793 formed both van der Waals interactions and one hydrogen bond with the benzodioxin moiety of doxazosin. Pro794 established both van der Waals contacts and one hydrogen bond with the quinazoline ring of doxazosin, whereas Glu804 formed hydrogen bonds with the quinazoline ring of doxazosin (Figure 1B). Stabilization of doxazosin did not establish interactions with Thr790 and Met766, two residues whose mutations have been linked to EGFR drug resistance [57, 58]. In addition, the characteristic interactions between Met793 and the quinazoline moiety were not observed, which has previously been observed for other TKIs of EGFR .
3.3 Structural analysis of complexes between doxazosin and HER2AC/HER2IN
Dozaxosin in complex with HER2AC was bound through van der Waals interactions with Leu726, Val734, Ala751, Phe864, Leu852, and Leu800 and polar interactions with Gly727, Ser783, Thr798, Gln799, Arg849, Thr862, Asp863, and Lys921 (Figure 1C). For the complex between doxazosin and HER2IN, the ligand was stabilized by Val754, Leu755, Met774, Leu785, Leu796, Pro802, Cys805, Leu852, and Phe864 and polar contacts with Lys753, Arg756, Gln799, Thr862, and Asp863. Tyr803 and Met801 formed van der Waals and hydrogen bonds with polar groups of the quinazoline ring, whereas Ser783 and Thr798 formed polar contacts with the linker between piperazine and the benzodioxin moiety (Figure 1D). Structural comparison of the complexes of doxazosin with HER2AC/HER2IN showed that doxazosin was better coordinated on HER2IN than HER2AC through more well-adjusted types of van der Waals and hydrogen bonds. In addition, the characteristic hinge hydrogen bond between Met801 and the polar atoms of the quinazoline moiety of several TKIs [31, 59] was present only in for the complex with doxazosin and HER2IN.
3.4 Structural analysis of complexes between gefitinib and EGFRAC/EGFRIN
Analysis of complexes between gefitinib and EGFRAC illustrated that the ligand was bound through van der Waals interactions with Leu718, Val726, Ala743, Met766, Leu788, Ile789, Leu792, Pro794, Phe795, Cys797, Leu844, and Phe856. The polar interactions were through contacts with Gly719, Ser720, Gly721, Lys745, Glu762, Thr790, Gln791, Gly796, Arg803, Thr854, and Asp855. Met793 formed van der Waals interactions and one polar interaction with the quinazoline ring of gefitinib, whereas Asp800 formed a hydrogen bond with one of the substituents at the quinazoline ring (Figure 2A). In complex with EGFRIN, gefitinib was stabilized by van der Waals contacts with Phe723, Val718, Val726, Ala743, Cys775, Leu792, and Cys797. The polar contacts were through Gly719, Ser720, Gly721, Gly724, Thr790, Gln791, Arg841, Asn842, Thr854, and Asp 855. Met793 established both van der Waals interactions and one hydrogen bond with the quinazoline ring, whereas Lys745 made a hydrogen bond with one of the substituents of the quinazoline ring (Figure 2B). Comparison of the map of interactions of both complexes showed that gefitinib was better coordinated on EGFRIN than EGFRAC. In both complexes, gefitinib established interactions with Thr790, a residue whose mutation is linked to EGFR drug resistance [57, 58]. In addition, in both complexes, the characteristic interactions between Met793 and the quinazoline moiety of ligand were observed, which has been reported elsewhere .
3.5 Structural analysis of complexes between gefitinib and HER2AC/HER2IN
Gefitinib in complex with HER2AC was bound through van der Waals interactions by Leu726, Val734, Leu800, Cys805, Leu807, Val851, and Leu852 (Figure 2C). Polar interactions were stabilized by Gly727, Ser728, Gly729, Asp808, Asp845, Arg849, Asn850, Thr862, and Gly865 residues, whereas Asp863 formed a hydrogen bond with one of the substituents of the quinazoline ring (Figure 2C). Gefitinib formed a complex with HER2IN, coordinated by Leu726, Val734, Val731, Cys805, Leu807, Leu866, Ala867, Leu869, Tyr923, Ala920, and Pro922 residues through van der Waals interactions. Polar interactions took place by Gly727, Ser728, Lys753, Arg811, Arg849, Asp863, and Lys921 residues, whereas Asp808 formed a hydrogen bond with one of the quinazoline ring substituents (Figure 2D). Structural comparison of both systems depicted that gefitinib was better stabilized on HER2IN than HER2AC. In addition, in both complexes, the characteristic polar interaction between Met801 and polar atoms of the quinazoline moiety of ligand was not observed [31, 59], as observed for the complex between doxasozin and HER2IN (Figure 1C).
3.6 Binding free energy
Determination of the Δ
|EGFRAC-doxazosin||−43.52 (5.13)||22.33 (4.91)||−2.69 (0.80)||−5.46 (0.50)||−29.33 (6.14)|
|EGFRIN-doxazosin||−38.78 (3.90)||3.25 (0.91)||9.68 (2.80)||−4.62 (0.37)||−30.46 (4.20)|
|HER2AC-doxazosin||−36.77 (6.0)||40.15 (12.0)||−21.06 (11.0)||−4.36 (0.55)||−22.05 (4.60)|
|HER2IN-doxazosin||−48.93 (4.0)||23.99 (9.0)||−7.39 (1.50)||−5.99 (0.30)||−38.32 (4.0)|
|EGFRAC-gefitinib||−55.01 (0.16)||−11.11 (0.65)||−25.40 (0.60)||−7.16 (0.01)||−47.88 (0.17)|
|EGFRIN-gefitinib||−44.10 (0.15)||27.00 (0.69)||−16.48 (0.65)||−5.74 (0.01)||−39.32 (0.23)|
|HER2AC-gefitinib||−34.76 (4.0)||38.46 (16.0)||−22.98 (5.0)||−4.90 (0.70)||−24.18 (5.0)|
|HER2IN-gefitinib||−38.78 (5.7)||−31.92 (13.0)||49.70 (1.50)||−4.90 (0.50)||−25.91 (6.0)|
3.7 Decomposition of the per-residue free energy
This analysis identified the residues that contributed the most to the Δ
Table 4 shows that Leu726, Gly727, Ser728, Val734, Asp863, and Phe864 contributed the most to the Δ
3.8 Principal component analysis
We evaluated the differences in mobility for the free and bound EGFR/HER2 systems via PC analysis. Evaluation of the quantification of the diagonalized covariance matrix based on covariance showed the following values: HER2AC, 15.0 nm2; HER2AC-doxazosin, 10.26 nm2; HER2AC-gefitinib, 8.83 nm2; HER2IN, 28.9 nm2; HER2IN-doxazosin, 18.21 nm2; HER2IN-gefitinib, 20.41 nm2; EGFRAC, 10.79 nm2; EGFRACAC-doxazosin, 9.46 nm2; EGFRAC-gefitinib, 7.69 nm2; EGFRIN, 8.17 nm2; EGFRIN-doxazosin, 11.55 nm2; and EGFRAC-gefitinib, 10.0 nm2. This result indicates that the molecular recognition of doxazosin or gefitinib on HER2AC, HER2IN, and EGFRAC decreased the number of conformational states compared to that of free HER2AC, HER2IN, and EGFRAC states. However, this conformational reduction was more significant for the free and bound HER2IN system. In contrast, a small increase in the conformational mobility was experienced upon the coupling of gefitinib by EGFRIN. This indicated that doxazosin and gefitinib binding to HER2IN was linked with reduced heterogeneity, which suggests that this molecular recognition was associated with an unfavorable entropy contribution that could contribute to decrease the favorable Δ
In this chapter, we explored the structural and energetic features that guide the similar inhibitory properties of doxazosin with gefitinib in overexpressing EGFR/HER2 cell lines combining docking and MD simulation with the MMGBSA approach. Based on these studies, we identified that doxazosin was able to target the active and inactive states of EGFR and HER2, however, its inhibitory activity against breast cancer cell lines was mainly by targeting HER2IN. Similarly, although gefitinib was able to target the inactive and inactive states of EGFR and HER2, its activity mainly targeted EGFRAC, in line with previous reports. Per-residue free energy analysis identified the key residues stabilizing HER2IN-DOX and EGFRAC-GEF systems, showing that in the stabilization of both systems, Met793 and Met801 were involved for EGFR and HER2, respectively. These residues stabilized HER2IN-DOX and EGFRAC-GEF systems through the formation of hydrogen bonds with the quinazoline ring, as reported for other TKIs. This study provides structural and energetic information that can be used to design new inhibitors for HER2IN or EGFRAC using doxazosin or gefitinib, respectively, as a pharmacophoric model.
This is an unpublished original article. The work was supported by grants from CONACYT (CB-A1-S-21278) and SIP/IPN (20201015).
M.A. Olayioye, R.M. Neve, H.A. Lane, N.E. Hynes, The ErbB signaling network: receptor heterodimerization in development and cancer, EMBO J. 19 (2000) 3159-3167.
G. Lurje, H.-J. Lenz, EGFR signaling and drug discovery, Oncology 77 (6) (2009) 400-410.
E. Tzahar, H. Waterman, X. Chen, G. Levkowitz, D. Karunagaran, S. Lavi, B.J. Ratzkin, Y. Yarden. A hierarchical network of interreceptor interactions determines signal transduction by Neu differentiation factor/neuregulin and epidermal growth factor. Mol. Cell. Biol. 16 (1996), pp. 5276-5287
X. Qian, C.M. LeVea, J.K., Freeman, W. C. Dougall, M.I. Greene. Heterodimerization of epidermal growth factor receptor and wild-type or kinase-deficient Neu: a mechanism of interreceptor kinase activation and transphosphorylation. Proc. Natl. Acad. Sci. U. S. A., 91 (1994), pp. 1500-1504
D. Graus-Porta, R.R. Beerli, J. M. Daly, N.E. Hynes. ErbB-2, the preferred heterodimerization partner of all ErbB receptors, is a mediator of lateral signaling. EMBO J., 16 (1997), pp. 1647-1655
D.J. Riese 2nd, R.M. Gallo, J. Settleman, Mutational activation of ErbB family receptor tyrosine kinases: insights into mechanisms of signal transduction and tumorigenesis, BioEssays 29 (2007) 558-565.
X. Zhang, J. Gureasko, K. Shen, P.A. Cole, J. Kuriyan, An allosteric mechanism for activation of the kinase domain of epidermal growth factor receptor, Cell 125 (2006) 1137-1149.
S.R. Hubbard, W.T. Miller, Receptor tyrosine kinases: mechanisms of activation and signaling, Curr. Opin. Cell Biol. 19 (2007) 117-123.
K. Aertgeerts, R. Skene, J. Yano, B.-C. Sang, H. Zou, G. Snell, A. Jennings, K. Iwamoto, N. Habuka, A. Hirokawa, T. Ishikawa, T. Tanaka, H. Miki, Y. Ohta, S. 15 Sogabe. Structural analysis of the mechanism of inhibition and allosteric activation of the kinase domain of HER2 protein. J. Biol. Chem., 286 (2011), pp. 18756-18765
Bello, M., Guadarrama-García, C. & Rodriguez-Fonseca, R.A. Dissecting the Molecular Recognition of Dual Lapatinib Derivatives for EGFR/HER2. J Comput Aided Mol Des 2020 Mar; 34(3):293-303.
Bello, Martiniano, et al. "Structural and energetic basis for the molecular recognition of dual synthetic vs. natural inhibitors of EGFR/HER2." International journal of biological macromolecules 111 (2018): 569-586.
S.V. Sharma, D.W. Bell, J. Settleman, D.A. Haber, Epidermal growth factor receptor mutations in lung cancer, Nat. Rev. Cancer 7 (2007) 169-181.
J.R.C. Sainsbury, J.R. Farndon, G.K. Needham, A.J. Malcolm, A.L. Harris, Epidermalgrowth-factor receptor status as predictor of early recurrence of and death from breast-cancer, Lancet 1 (1987) 1398-1402
F. Ciardiello, G. Tortora. A novel approach in the treatment of cancer: targeting the epidermal growth factor receptor. Clin Cancer Res., 7(2001), pp. 2958-2970
B.N. Rexer, R. Ghosh, A. Narasanna. Human breast cancer cells harboring a gatekeeper T798M mutation in HER2 overexpress EGFR ligands and are sensitive to dual inhibition of EGFR and HER2. Clin Cancer Res., 19(2013), pp. 5390-401
I.M. Gonzaga, S.C. Soares-Lima, P.T. de Santos. Alterations in epidermal growth factor receptors 1 and 2 in esophageal squamous cell carcinomas. BMC Cancer., 12(2012), pp. 56
E.R. Wood, A.T. Truesdale, O.B. McDonald, D. Yuan, A. Hassell, S.H. Dickerson, B. Ellis, C. Pennisi, E. Horne, K. Lackey, K.J. Alligood, D.W. Rusnak, T.M. Gilmer, L. Shewchuk. A unique structure for epidermal growth factor receptor bound to GW572016 (Lapatinib): relationships among protein conformation, inhibitor off-rate, and receptor activity in tumor cells. Cancer Res., 64 (2004), pp. 6652-6659
M.A. Seeliger, P. Ranjitkar, C. Kasap, Y. Shan, D.E. Shaw, N.P. Shah, J. Kuriyan, D.J. Maly. Equally potent inhibition of c-Src and Abl by compounds that recognize inactive kinase conformations. Cancer Res., 69 (2009), pp. 2384-2392
M. Bello. Binding mechanism of kinase inhibitors to EGFR and T790M, L858R and L858R/T790M mutants through structural and energetic analysis. International journal of biological macromolecules, 118 2018, pp. 1948-1962
L. Saldaña-Rivera, M. Bello, D. Méndez-Luna. Structural insight into the binding mechanism of ATP to EGFR and L858R, and T790M and L858R/T790 mutants. Journal of Biomolecular Structure and Dynamics, 1 (2019), pp. 1-14
D.W. Rusnak, K. Lackey, K. Affleck. The effects of the novel, reversible epidermal growth factor receptor/ErbB-2 tyrosine kinase inhibitor, GW2016, on the growth of human normal and tumor-derived cell lines in vitro and in vivo. Mol. Cancer Ther., 1 (2001), pp. 85-9
G.E. Konecny, M.D. Pegram, N. Venkatesan, R. Finn, G. Yang. Activity of the dual kinase inhibitor lapatinib (GW572016) against HER-2- overexpressing and trastuzumab-treated breast cancer cells. Cancer Res., 66 (2006), pp. 1630-1639
S.R. Johnston, A. Leary. Lapatinib: a novel EGFR/HER2 tyrosine kinase inhibitor for cancer. Drugs Today Barc, 42 (2006), pp. 441-453
Y. Zhou, S. Li, Hu YP, J. Wang, J. Hauser. Blockade of EGFR and ErbB2 by the novel dual EGFR and ErbB2 tyrosine kinase inhibitor GW572016 sensitizes human colon carcinoma GEO cells to apoptosis. Cancer Res., 66 (2006), pp. 404-411
P.J. Medina, S. Goodin. Lapatinib: a dual inhibitor of human epidermal growth factor receptor tyrosine kinases. Clin. Ther., 30 (2008), pp. 1426-1447
S.R. Johnston, A. Leary, Lapatinib: a novel EGFR/HER2 tyrosine kinase inhibitor for cancer, Drugs Today (Barc) 42 (2006) 441-453.
B. Forsythe, K. Faulkner, Overview of the tolerability of gefitinib (IRESSA) monotherapy: clinical experience in non-small-cell lung cancer, Drug Saf. 27 (2004) 1081-1092.
D.M. Jackman, V.A. Miller, L.A. Cioffredi, B.Y. Yeap, P.A. Jänne, G.J. Riely, M.G. Ruiz, G. Giaccone, L.V. Sequist, B.E. Johnson, Clin. Cancer Res. 15 (2009) 5267.
W. Pao, V.A. Miller, K.A. Politi, G.J. Riely, R. Somwar, M.F. Zakowski, M.G. Kris, H. Varmus, Acquired resistance of lung adenocarcinomas to gefitinib or erlotinib is associated with a second mutation in the EGFR kinase domain, PLoS Med. 2 (2005), e73.
C.H. Yun, K.E. Mengwasser, A.V. Toms, M.S. Woo, H. Greulich, K.K. Wong, M. Meyerson, M.J. Eck, The T790M mutation in EGFR kinase causes drug resistance by increasing the affinity for ATP, PNAS 105 (2008) 2070-2075.
M. Ahmed, M. Sadek, K. A. Abouzid, F. Wang. In silico design: extended molecular dynamic simulations of a new series of dually acting inhibitors against EGFR and HER2. Journal of Molecular Graphics and Modelling, 44 (2013), pp. 220-231
M.M. Sadek, R.A. Serrya, A.-H.N. Kafafy, M. Ahmed, F. Wang, K.A.M. Abouzid, Discovery of new HER2/EGFR dual kinase inhibitors based on the anilinoquinazoline scaffold as potential anticancer agents, J. Enzyme Inhib. Med. Chem. 1 (2013) 1-8.
S.C. Yang, S.S. Chang, H.Y. Chen, C.Y.C. Chen, Identification of potent EGFR inhibitors from TCM database@ Taiwan, PLoS Comput. Biol. 7 (10) (2011), e1002189.
Y. Liu, N.S. Gray, Rational design of inhibitors that bind to inactive kinase conformations, Nat. Chem. Biol. 2 (2006) 358-364.
C.N. Cavasotto, R.A. Abagyan, Protein flexibility in ligand docking and virtual screening to protein kinases, J. Mol. Biol. 337 (2004) 209-226.
B. Chandrika, J. Subramanian, S.D. Sharma, Managing protein flexibility in docking and its applications, Drug Discov. Today 14 (2009) 394-400.
Hui, Hongxiang, Manory A. Fernando, and Anthony P. Heaney. "The α1-adrenergic receptor antagonist doxazosin inhibits EGFR and NF-κB signalling to induce breast cancer cell apoptosis." European journal of cancer44.1 (2008): 160-166.
Ford, Christopher HJ, et al. "Reassessment of estrogen receptor expression in human breast cancer cell lines." Anticancer research31.2 (2011): 521-527.
Stanley, Aryan, et al. "Synergistic effects of various Her inhibitors in combination with IGF-1R, C-MET and Src targeting agents in breast cancer cell lines." Scientific Reports7.1 (2017): 1-15.
El Sharkawi, Fathia Zaky, Hany Abdelaziz El Shemy, and Hussein Moustafa Khaled. "Possible anticancer activity of rosuvastatine, doxazosin, repaglinide and oxcarbazepin." Asian Pacific Journal of Cancer Prevention15.1 (2014): 199-203.
A. Sali, T.L. Blundell, Comparative protein modelling by satisfaction of spatial restraints, J. Mol. Biol. 234 (3) (1993) 779-815.
G.M. Morris, R. Huey, W. Lindstrom, M.F. Sanner, R.K. Belew, D.S. Goodsell, A.J. Olson, Computational, Chemistry 16 (2009) 2785-2791.
M.J. Frisch, G.W. Trucks, H.B. Schlegel, G.E. Scuseria, M.A. Robb, J.R. Cheeseman, J.A. Montgomery, Jr., T. Vreven, K.N. Kudin, J.C. Burant, J.M. Millam, S.S. Iyengar, J. Tomasi, V. Barone, B. Mennucci, M. Cossi, G. Scalmani, N. Rega, G.A. Petersson, H. Nakatsuji, M. Hada, M. Ehara, K. Toyota, R. Fukuda, J. Hasegawa, M. Ishida, T. Nakajima, Y. Honda, O. Kitao, H. Nakai, M. Klene, X. Li, J.E. Knox, H.P. Hratchian, J.B. Cross, V. Bakken, C. Adamo, J. Jaramillo, R. Gomperts, R.E. Stratmann, O. Yazyev, A.J. Austin, R. Cammi, C. Pomelli, J.W. Ochterski, P.Y. Ayala, K. Morokuma, G.A. Voth, P. Salvador, J.J. Dannenberg, V.G. Zakrzewski, S. Dapprich, A.D. Daniels, M.C. Strain, O. Farkas, D.K. Malick, A.D. Rabuck, K. Raghavachari, J.B. Foresman, J.V. Ortiz, Q. Cui, A.G. Baboul, S. Clifford, J. Cioslowski, B.B. Stefanov, G. Liu, A. Liashenko, P. Piskorz, I. Komaromi, R.L. Martin, D.J. Fox, T. Keith, M.A. Al-Laham, C.Y. Peng, A. Nanayakkara, M. Challacombe, P.M.W. Gill, B. Johnson, W. Chen, M.W. Wong, C. Gonzalez, and J.A. Pople, Gaussian 03, Revision B.04, Gaussian, Inc., Wallingford, CT, 2004
D.A. Case, T.E. Cheatham, T. Darden, H. Gohlke, R. Luo, K.M. Merz Jr., R.J. Woods, The Amber biomolecular simulation programs, J. Comput. Chem. 26 (2005) 1668-1688
Y. Duan, C. Wu, S. Chowdhury, M.C. Lee, G. Xiong, W. Zhang, et al., Apoint-charge force field for molecular mechanics simulations of proteinsbased on condensedphase quantum mechanical calculations, J. Comput. Chem. 24 (16) (2003) 1999-2012
W.L. Jorgensen, J. Chandrasekhar, J.D. Madura, R.W. Impey, M.L. Klein, Comparison of simple potential functions for simulating liquid water, J. Chem. Phys. 79 (1983) 926-935.
J. Wang, R.M. Wolf, J.W. Caldwell, P.A. Kollman, D.A. Case, Development andtesting of a general amber force field, J. Comput. Chem. 25 (9) (2004) 1157-1174
H.J.C. Berendsen, J.P.M. Postma, W.F. van Gunsteren, A. DiNola, J.R. Haak, Moleculardynamics with coupling to an external bath, J. Chem. Phys. 81 (8) (1984) 3684-3690.
T. Darden, D. York, L. Pedersen, Particle mesh Ewald-an N. log(N) method forewald sums in large systems, J. Chem. Phys. 98 (1993) 10089-10092.
W.F. Van Gunsteren, H.J.C. Berendsen, Algorithms for macromoleculardynamics and constraint dynamics, Mol. Phys. 34 (1977) 1311-1327.
Maestro, version 10.1. Schrödinger, LLC; New York, NY, USA: 2015-1.
B.R. Miller, T.D. McGee, J.M. Swails, N. Homeyer, H. Gohlke, A.E. Roitberg. MMPBSA.py: an efficient program for end-state free energy calculations. J.Chem. Theory Comput., 8 (2012), pp. 3314-3321
H. Gohlke, C. Kiel, D. A. J. Case. Insights into protein-protein binding by binding free energy calculation and free energy decomposition for the Ras-Raf and Ras-RalGDS complexes. Mol. Biol., 330 (2003), pp. 891-913
P.A. Kollman, I. Massova, C. Reyes, B. Kuhn, S. Huo, L. Chong. Calculating structures and free energies of complex molecules: combining molecular mechanics and continuum models. Acc. Chem. Res., 33 (2000), pp. 889-897
J. M. Wang, T. J. Hou, X. J. Xu. Recent advances in free energy calculations with a combination of molecular mechanics and continuum model. Drug Des., 2 (2006), pp. 287-306
A. Onufriev, V. Bashford, D. A. Case. Exploring protein native states and large-scale conformational changes with a modified generalized born model. Proteins., 55 (2004), pp. 383-394
C.H. Yun, K.E. Mengwasser, A.V. Toms, M.S. Woo, H. Greulich, K.K. Wong, M. Meyerson, M.J. Eck. The T790M mutation in EGFR kinase causes drug resistance by increasing the affinity for ATP. PNAS, 105 (2008), pp. 2070-2075
Z. Ruan, S. Katiyar, N. Kannan. Computational and Experimental Characterization of Patient Derived Mutations Reveal an Unusual Mode of Regulatory Spine Assembly and Drug Sensitivity in EGFR Kinase. Biochemistry, 56(2017), pp. 22-32
Xing, L., Klug-Mcleod, J., Rai, B., & Lunney, E. A. (2015). Kinase hinge binding scaffolds and their hydrogen bond patterns. Bioorganic & medicinal chemistry, 23(19), 6520-6527.
J.H. Park, Y. Liu, M.A. Lemmon, R. Radhakrishnan, Erlotinib binds both inactive and active conformations of the EGFR tyrosine kinase domain, Biochem. J. 448 (3) (2012) 417-423