It is now clear that proteins are flexible entities that in solution switch between conformations to achieve their function. Hydrogen/Deuterium Exchange Mass Spectrometry (HX/MS) is an invaluable tool to understand dynamic changes in proteins modulated by cofactor binding, post-transductional modifications, or protein-protein interactions. ERK2MAPK, a protein involved in highly conserved signal transduction pathways of paramount importance for normal cellular function, has been extensively studied by HX/MS. Experiments of the ERK2MAPK in the inactive and active states (in the presence or absence of bound ATP) have provided valuable information on the plasticity of the MAPK domain. However, interpretation of the HX/MS data is difficult, and changes are mostly explained in relation to available X-ray structures, precluding a complete atomic picture of protein dynamics. In the present work, we have used all atom Molecular Dynamics simulations (MD) to provide a theoretical framework for the interpretation of HX/MS data. Our results show that detailed analysis of protein-solvent interaction along the MD simulations allows (i) prediction of the number of protons exchanged for each peptide in the HX/MS experiments, (ii) rationalization of the experimentally observed changes in exchange rates in different protein conditions at the residue level, and (iii) that at least for ERK2MAPK, most of the functionally observed differences in protein dynamics are related to what can be considered the native state conformational ensemble. In summary, the combination of HX/MS experiments with all atom MD simulations emerges as a powerful approach to study protein native state dynamics with atomic resolution.

Mitogen-activated protein kinases (MAPKs) are serine-threonine kinases that participate in signal transduction pathways. p38 MAPKs have four isoforms (p38α, p38β, p38γ, and p38δ) which are involved in multiple cellular functions such as proliferation, differentiation, survival, and migration. MAPK kinases phosphorylate p38s in the dual-phosphorylation motif, Thr-Gly-Tyr, located in their activation loop, which induces a conformational change that increases ATP binding affinity and catalytic activity. Several works have proposed that MAPK dynamics is a key factor in determining their function. However, we still do not understand the dynamical changes that lead to MAPK activation. In this work we have used molecular dynamics techniques to study the dynamical changes associated with p38γ activation, the only fully active MAPK crystallized so far. We performed MD simulations of p38γ in three different states, fully active with ATP, active without ATP, and inactive. We found that the dynamical fluctuations of the docking sites, important for protein-protein interactions, are regulated allosterically by changes in the active site. Interestingly, in the phosphorylated and ATP-bound states the whole protein dynamics lead to concerted motions of whole protein domains in contrast to the inactive state. The binding/unbinding of ATP participates in the reorientation of the two domains and in the regulation of protein plasticity. Our study shows that beyond the conformational changes associated with MAPK activation their correlated dynamics are highly regulated by phosphorylation and ATP binding. This means that MAPK plasticity may have a role in their catalytic activity, specificity, and protein-protein interactions and, therefore, in the outcome of the signaling network.