Computational Analysis of Activation of Dimerized Epidermal Growth Factor Receptor Kinase Using the String Method and Markov State Model

Epidermal growth factor receptor (EGFR) activation is accompanied by dimerization. During the activation of the intracellular kinase domain, two EGFR kinases form an asymmetric dimer, and one side of the dimer (receiver) is activated. Using the string method and Markov state model (MSM), we performed a computational analysis of the structural changes in the activation of the EGFR dimer in this study. The string method reveals the minimum free-energy pathway (MFEP) from the inactive to active structure. The MSM was constructed from numerous trajectories of molecular dynamics simulations around the MFEP, which revealed the free-energy map of structural changes. In the activation of the receiver kinase, the unfolding of the activation loop (A-loop) is followed by the rearrangement of the C-helix, as observed in other kinases. However, unlike other kinases, the free-energy map of EGFR at the asymmetric dimer showed that the active state yielded the highest stability and revealed how interactions at the dimer interface induced receiver activation. As the H-helix of the activator approaches the C-helix of the receiver during activation, the A-loop unfolds. Subsequently, L782 of the receiver enters the pocket between the G- and H-helices of the activator, leading to a rearrangement of the hydrophobic residues around L782 of the receiver, which constitutes a structural rearrangement of the C-helix of the receiver from an outward to an inner position. The MSM analysis revealed long-time scale trajectories via kinetic Monte Carlo.