Structural basis for the selectivity of 3rd generation EGFR inhibitors: a molecular dynamics study

Activating mutations in the EGFR kinase domain are known to be a common cause of Non-Small Cell Lung Cancer (NSCLC) and are thus targeted for treatment. First generation Tyrosine Kinase Inhibitors (TKIs) were used to treat NSCLC patients with the known activating mutations L858R and exon 19 deletion but were resisted by a second mutation T790M in the active site of the kinase domain. Second generation members of TKIs have an electrophilic moiety that can form a covalent bond with Cys797 and are effective against T790M EGFR but are toxic because they inhibit WT EGFR as well. Third generation TKIs, like Osimertinib, can bind to and irreversibly inhibit T790M mutants selectively, while sparing the wild-type enzyme. Thus, they possess a better safety profile and a wider therapeutic window. However, the reason behind their selectivity is still not well understood. In this study, computational MD simulations were carried out on Osimertinib in complex with both WT and L858R/T790M Double Mutant (DM) EGFR to provide an insight into the selectivity of Osimertinib and its molecular interactions within the active site. A high-resolution trajectory analysis showed that the key selectivity residues are Val726, met793, and Cys797. Interaction of Osimertinib with these residues is improved due to the T790M mutation which optimizes the ligand orientation for binding, as evident from the RMSD and the distances monitored. These results can provide guidance for the development of more selective 3rd generation EGFR TKIs.Communicated by Ramaswamy H. Sarma