The stability and dynamics of cytoskeleton in brain nerve cells are regulated by microtubule associated proteins (MAPs), tau and MAP2. Both proteins are intrinsically disordered and involved in multiple molecular interactions important for normal physiology and pathology of chronic neurodegenerative diseases. Nuclear magnetic resonance and cryo-electron microscopy recently revealed propensities of MAPs to form transient local structures and long-range contacts in the free state, and conformations adopted in complexes with microtubules and filamentous actin, as well as in pathological aggregates. In this paper, we compare the longest, 441-residue brain isoform of tau (tau40), and a 467-residue isoform of MAP2, known as MAP2c. For both molecules, we present transient structural motifs revealed by conformational analysis of experimental data obtained for free soluble forms of the proteins. We show that many of the short sequence motifs that exhibit transient structural features are linked to functional properties, manifested by specific interactions. The transient structural motifs can be therefore classified as molecular recognition elements of tau40 and MAP2c. Their interactions are further regulated by post-translational modifications, in particular phosphorylation. The structure-function analysis also explains differences between biological activities of tau40 and MAP2c.

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.