Activating mutations in the epidermal growth factor receptor (EGFR) tyrosine kinase are frequently found in many cancers. It has been suggested that changes in the equilibrium between its active and inactive conformations are linked to its oncogenic potential. Here, we quantify the effects of some of the most common single (L858R and T790M) and double (T790M-L858R) oncogenic mutations on the conformational free-energy landscape of the EGFR kinase domain by using massive molecular dynamics simulations together with parallel tempering, metadynamics, and one of the best force-fields available. Whereas the wild-type EGFR catalytic domain monomer is mostly found in an inactive conformation, our results show a clear shift toward the active conformation for all of the mutants. The L858R mutation stabilizes the active conformation at the expense of the inactive conformation and rigidifies the αC-helix. The T790M gatekeeper mutant favors activation by stabilizing a hydrophobic cluster. Finally, T790M with L858R shows a significant positive epistasis effect. This combination not only stabilizes the active conformation, but in nontrivial ways changes the free-energy landscape lowering the transition barriers.conformational changes | resistance-causing mutations

Metadynamics is an algorithm for accelerating rare events and reconstructing the associated free energy surface. It works by biasing the evolution of the system by a history‐dependent potential that is adaptively constructed in the space of a suitably chosen set of collective variables. Since its first appearance, the method has been successfully applied in several domains of science. Its widespread adoption is not only due to its efficiency, flexibility, and availability but also to its continuous evolution and its combination with complementary enhanced sampling algorithms. Here, we focus on the progress made in the development of more general and powerful collective variables and on the very recent and exciting evolutions of the method. © 2012 John Wiley & Sons, Ltd.This article is categorized under: Molecular and Statistical Mechanics > Molecular Dynamics and Monte-Carlo Methods

c-Src and c-Abl are two closely related protein kinases that constitute important anticancer targets. Despite their high sequence identity, they show different sensitivities to the anticancer drug imatinib, which binds specifically to a particular inactive conformation in which the Asp of the conserved DFG motif points outward (DFG-out). We have analyzed the DFG conformational transition of the two kinases using massive molecular dynamics simulations, free energy calculations, and isothermal titration calorimetry. On the basis of the reconstruction of the free energy surfaces for the DFG-in to DFG-out conformational changes of c-Src and c-Abl, we propose that the different flexibility of the two kinases results in a different stability of the DFG-out conformation and might be the main determinant of imatinib selectivity.