Active EGFR-TKD was modelled on the basis of PDB entries 1M17 (which also provided the initial erlotinib conformation) [10 (link)] and 2ITX [11 (link)], and the L834R mutant was modelled on the basis of PDB entry 2ITV [11 (link)]. Inactive EGFR-TKD was modelled based on PDB entries 2GS7 [12 (link)] and 1XKK [13 (link)]. Protein and ligand conformations were prepared using the Protein Preparation Wizard and LigPrep protocols from Schrödinger Software. All docking simulations used the OPLS (Optimized Potentials for Liquid Simulations) force field [30 ], and used Schrödinger's IFD (Induced Fit Docking) package [31 (link)]. Ligand was first docked to rigid protein using Glide XP [32 (link)]. For the resulting top 20 complex conformations, the protein side chains within 5.0 Å of the ligand in that pose were subjected to conformational search and minimized using Prime [33 (link)] and the ligand was redocked to the 20 new receptor conformations.