Abstract
The epidermal growth factor receptor (EGFR) is a tyrosine kinase receptor that plays a fundamental role in regulating cellular proliferation, survival, and differentiation. EGFR's kinase domain catalyses the autophosphorylation that drives downstream signalling and is a known hotspot for cancer-associated mutations. Such substitutions often enhance receptor activation, driving oncogenesis and therapeutic resistance, with major implications for prognosis and treatment choice; consequently, a detailed understanding of this domain's functional dynamics, and how mutations reshape them, is essential. Here, we apply dynamical-nonequilibrium molecular dynamics (D-NEMD) simulations to probe how changes at the ATP-binding site are transmitted throughout the kinase domain, ultimately shaping the behaviour of distal regions in the wild-type protein and in two of the most clinically observed non-small cell lung cancer variants, L858R and L858R+T790M. D-NEMD reveals previously uncharacterised communication networks linking the ATP pocket to key functional motifs, including the P-loop, αC-helix, hinge region, activation loop, and αF-helix, showing the intricate allosteric connectivity within this domain. Notably, numerous cancer-related mutation sites, despite being distant from the ATP site, map onto or lie close to these communication pathways, indicating that D-NEMD simulations can identify functionally relevant allosteric positions and potential mutation sites in EGFR. Our simulations further reveal that the two variants rewire EGFR's internal signal propagation network in distinct ways: L858R diminishes the transmission of structural changes to distal regions of the C-lobe, whereas introducing T790M in addition to L858R partially restores wild-type-like behaviour through compensatory rerouting of early-stage dynamical pathways.