Abstract
Hyperactivation caused by mutations in the Epidermal Growth Factor Receptor (EGFR) kinase domain is implicated in various diseases, including cancer. However, the structural mechanisms underlying overactivation in many EGFR mutations remain poorly understood, and exploring these mechanisms through conventional experiments or in silico simulations is often labor- and cost-intensive. Here, we establish a Molecular Dynamics (MD) protocol capable of rapidly revealing EGFR mutant modes of action using multiple short simulations. We first simulated wild-type EGFR and the well-studied oncogenic mutations L858R and T790M/L858R under different simulation conditions, to derive a protocol which could recapitulate their experimentally established behavior. We then applied this protocol to three rare EGFR mutations: S768I, S768N, and D761N. Experimental studies have suggested that S768I and D761N are oncogenic, whereas S768N is likely a neutral mutation that does not significantly alter EGFR activity. Our simulation results were consistent with these functional indications and provided the corresponding molecular bases - S768I and S768N affect the orientation and stability of the catalytically important αC-helix, while D761N introduces a new hydrogen bonding network between the αC-helix and activation loop. Collectively, the protocol presented here provides a robust and rapid framework for characterizing EGFR mutation mechanisms and is readily adaptable to novel or uncharacterized variants.