The epidermal growth factor receptor (EGFR) is a key protein in cellular signaling, and its kinase domain (EGFR kinase) is an intensely pursued target of small-molecule drugs. Although both catalytically active and inactive conformations of EGFR kinase have been resolved crystallographically, experimental characterization of the transitions between these conformations remains difficult. Using unbiased, all-atom molecular dynamics simulations, we observed EGFR kinase spontaneously transition from the active to the so-called "Src-like inactive" conformation by way of two sets of intermediate conformations: One corresponds to a previously identified locally disordered state and the other to previously undescribed "extended" conformations, marked by the opening of the ATP-binding site between the two lobes of the kinase domain. We also simulated the protonation-dependent transition of EGFR kinase to another ["Asp-Phe-Gly-out" ("DFG-out")] inactive conformation and observed similar intermediate conformations. A key element observed in the simulated transitions is local unfolding, or "cracking," which supports a prediction of energy landscape theory. We used hydrogen-deuterium (H/D) exchange measurements to corroborate our simulations and found that the simulated intermediate conformations correlate better with the H/D exchange data than existing active or inactive EGFR kinase crystal structures. The intermediate conformations revealed by our simulations of the transition process differ significantly from the existing crystal structures and may provide unique possibilities for structure-based drug discovery.