Specific phosphorylation of microtubule-associated protein 2c by extracellular signal-regulated kinase reduces interactions at its Pro-rich regions.

Microtubule-associated protein 2 (MAP2) is an important neuronal target of extracellular signal-regulated kinase 2 (ERK2) involved in Raf signaling pathways, but mechanistic details of MAP2 phosphorylation are unclear. Here, we used NMR spectroscopy to quantitatively describe the kinetics of phosphorylation of individual serines and threonines in the embryonic MAP2 variant MAP2c. We carried out real-time monitoring of phosphorylation to discover major phosphorylation sites that were not identified in previous studies relying on specific antibodies. Our comparison with the phosphorylation of MAP2c by a model cyclin-dependent kinase CDK2 and with phosphorylation of the MAP2c homolog Tau revealed differences in phosphorylation profiles that explain specificity of regulation of biological functions of MAP2c and Tau. To probe the molecular basis of the regulatory effect of ERK2, we investigated the interactions of phosphorylated and unphosphorylated MAP2c by NMR with single-residue resolution. As ERK2 phosphorylates mostly outside the regions binding microtubules, we studied the binding of proteins other than tubulin, namely regulatory subunit RIIα of cAMP-dependent PKA, adapter protein Grb2, Src homology domain 3 of tyrosine kinases Fyn and Abl, and ERK2 itself. We found ERK2 phosphorylation interfered mostly with binding to proline-rich regions of MAP2c. Furthermore, our NMR experiments in SH-SY5Y neuroblastoma cell lysates showed that the kinetics of dephosphorylation are compatible with in-cell NMR studies and that residues targeted by ERK2 and PKA are efficiently phosphorylated in the cell lysates. Taken together, our results provide a deeper characterization of MAP2c phosphorylation and its effects on interactions with other proteins.

Investigation of the structure and dynamics of gallium binding to high-affinity peptides elucidated by multi-scale simulation, quantum chemistry, NMR and ITC.

Gallium (as Ga3+) is a Group IIIa metal and its recovery from wastewaters has become increasingly important for its reuse. The use of peptides for recycling offers a low-cost and environmentally-friendly option but the structural characteristics of peptides likely to bind Ga3+ are largely unknown. Multiple computational methods, coupled with experimental verification via NMR and Isothermal Calorimetry (ITC), were used to establish that Ga3+ binds with high affinity to peptide sequences and to elucidate the structural characteristics that contributed. It was demonstrated that peptide pre-organisation is key to Ga3+ binding and that a favourable binding position is necessarily governed by the size and shape of the electrostatic environment as much as individual electrostatic interactions with peptide residues themselves. Given favourable conditions, Ga3+ retrieved plausible binding positions involving both charged and uncharged residues that greatly increases the range of bonding possibilities with other peptide sequences and offers insights for binding other metals. The addition of pH buffer substantially improved the affinity of Ga3+ and a structural role for a buffer component was demonstrated.

Choice of Force Field for Proteins Containing Structured and Intrinsically Disordered Regions.

Biomolecular force fields optimized for globular proteins fail to properly reproduce properties of intrinsically disordered proteins. In particular, parameters of the water model need to be modified to improve applicability of the force fields to both ordered and disordered proteins. Here, we compared performance of force fields recommended for intrinsically disordered proteins in molecular dynamics simulations of three proteins differing in the content of ordered and disordered regions (two proteins consisting of a well-structured domain and of a disordered region with and without a transient helical motif and one disordered protein containing a region of increased helical propensity). The obtained molecular dynamics trajectories were used to predict measurable parameters, including radii of gyration of the proteins and chemical shifts, residual dipolar couplings, paramagnetic relaxation enhancement, and NMR relaxation data of their individual residues. The predicted quantities were compared with experimental data obtained within this study or published previously. The results showed that the NMR relaxation parameters, rarely used for benchmarking, are particularly sensitive to the choice of force-field parameters, especially those defining the water model. Interestingly, the TIP3P water model, leading to an artificial structural collapse, also resulted in unrealistic relaxation properties. The TIP4P-D water model, combined with three biomolecular force-field parameters for the protein part, significantly improved reliability of the simulations. Additional analysis revealed only one particular force field capable of retaining the transient helical motif observed in NMR experiments. The benchmarking protocol used in our study, being more sensitive to imperfections than the commonly used tests, is well suited to evaluate the performance of newly developed force fields.