Functionally specific binding regions of microtubule-associated protein 2c exhibit distinct conformations and dynamics

Microtubule-associated protein 2c (MAP2c) is a 49-kDa intrinsically disordered protein regulating the dynamics of microtubules in developing neurons. MAP2c differs from its sequence homologue Tau in the pattern and kinetics of phosphorylation by cAMP-dependent protein kinase (PKA). Moreover, the mechanisms through which MAP2c interacts with its binding partners and the conformational changes and dynamics associated with these interactions remain unclear. Here, we used NMR relaxation and paramagnetic relaxation enhancement techniques to determine the dynamics and long-range interactions within MAP2c. The relaxation rates revealed large differences in flexibility of individual regions of MAP2c, with the lowest flexibility observed in the known and proposed binding sites. Quantitative conformational analyses of chemical shifts, small-angle X-ray scattering (SAXS), and paramagnetic relaxation enhancement measurements disclosed that MAP2c regions interacting with important protein partners, including Fyn tyrosine kinase, plectin, and PKA, adopt specific conformations. High populations of polyproline II and α-helices were found in Fyn- and plectin-binding sites of MAP2c, respectively. The region binding the regulatory subunit of PKA consists of two helical motifs bridged by a more extended conformation. Of note, although MAP2c and Tau did not differ substantially in their conformations in regions of high sequence identity, we found that they differ significantly in long-range interactions, dynamics, and local conformation motifs in their N-terminal domains. These results highlight that the N-terminal regions of MAP2c provide important specificity to its regulatory roles and indicate a close relationship between MAP2c's biological functions and conformational behavior.


Microtubule-associated protein 2c (MAP2c) is a 49-kDa
intrinsically disordered protein regulating the dynamics of microtubules in developing neurons. MAP2c differs from its sequence homologue Tau in the pattern and kinetics of phosphorylation by cAMP-dependent protein kinase (PKA). Moreover, the mechanisms through which MAP2c interacts with its binding partners and the conformational changes and dynamics associated with these interactions remain unclear. Here, we used NMR relaxation and paramagnetic relaxation enhancement techniques to determine the dynamics and long-range interactions within MAP2c. The relaxation rates revealed large differences in flexibility of individual regions of MAP2c, with the lowest flexibility observed in the known and proposed binding sites. Quantitative conformational analyses of chemical shifts, small-angle X-ray scattering (SAXS), and paramagnetic relaxation enhancement measurements disclosed that MAP2c regions interacting with important protein partners, including Fyn tyrosine kinase, plectin, and PKA, adopt specific conformations. High populations of polyproline II and ␣-helices were found in Fyn-and plectin-binding sites of MAP2c, respectively. The region binding the regulatory subunit of PKA consists of two helical motifs bridged by a more extended conformation. Of note, although MAP2c and Tau did not differ substantially in their conformations in regions of high sequence identity, we found that they differ significantly in long-range interactions, dynamics, and local conformation motifs in their N-terminal domains. These results highlight that the N-terminal regions of MAP2c provide important specificity to its regulatory roles and indicate a close relationship between MAP2c's biological functions and conformational behavior.
Cytoskeletal microtubule-associated proteins (MAPs) 6 bind, stabilize, and regulate dynamics of microtubules, in a phosphorylation-dependent manner. MAP2 and Tau are neuronal MAPs, MAP2 being expressed mainly in dendrites, whereas Tau is found in axons (1). Both MAP2 and Tau are expressed as different spliced variants. Tau isoforms expressed in human brain differ in the number of near-N-terminal inserts as well as in the number of repeats in the microtubule-binding domain (MTBD), whereas MAP2 isoforms differ in the length of the N-terminal projection domain (2,3). Expression of both MAP2 and Tau isoforms is regulated during development. The highmolecular weight isoforms, MAP2a and MAP2b, contain 1830 amino acids, and the two low-molecular weight variants, MAP2c and MAP2d, consist of 467 and 498 amino acids, respectively. MAP2c is the smallest functional isoform, expressed mainly during embryonal brain development. After birth, its expression is restricted to regions exhibiting post-natal plasticity, such as the olfactory bulb (4), suggesting a role in neuronal development.
Tau and MAP2c share a high sequence homology in the C-terminal part, containing MTBD, but differ in the N-terminal part, comprising the acidic and the proline-rich subdomains ( Fig. 1) (5). The sequence diversity is further magnified by different phosphorylation patterns (6 -13).
Tau, but not MAP2c, forms paired helical filaments observed in brains of patients suffering from Alzheimer's disease. The major difference in the aggregation properties seems to be encoded in a few amino acids of the aggregation seed motifs in the otherwise highly homologous C-terminal domains of Tau and MAP2c (14). In addition, interactions between the MTBD and the N-terminal domain of Tau are present in paired helical filaments (15).
Identification of structural motifs responsible for specific functions of MAP2c and Tau is complicated by the fact that MAP2c and Tau belong to the class of intrinsically disordered proteins (IDPs) lacking a unique structure and exist in multiple, quickly interconverting conformations (16 -20). NMR has proved to be a reliable tool to study their dynamics and structural features (11,13,(21)(22)(23)(24)(25). Despite their highly flexible nature, Tau (23) and MAP2c (25) adopt transient secondary structures, indicating that the structures of Tau and MAP2c are far from a random coil. Also, long-range intramolecular contacts have been observed in Tau (15), supporting the "paperclip" model of the Tau molecule.
Direct studies of structure-function relations rely on the availability of atomistic structural models. Quantitative computational analysis of experimental data reflecting the local conformations and long-range contacts provided an atomicresolution ensemble model describing molecular properties of Tau (26). In this paper, we built a similar model for MAP2c and investigated molecular motions of MAP2c using NMR relaxation. The results allowed us to correlate structural features and dynamics of MAP2c with its known and proposed binding sites. In particular, we addressed questions related to interactions of MAP2c with three important binding partners. The first of them was the Src-homology domain 3 (SH3) domain of the Fyn tyrosine kinase, recognizing different motifs of MAP2c and Tau (27,28). It was not clear whether the conformation optimal for the SH3 binding is highly populated already in the free form of intrinsically disordered MAP2c or if it is formed only in the bound state. The second examined interaction partner was plectin. Its noncanonical SH3 domain (29) was reported to bind MAP2c, but the recognized motif was unknown (30). The third inspected interaction region was the site recognized by the RII regulatory subunit of cAMP-dependent protein kinase (PKA). This site was expected to adopt ␣-helical conformation, but populations of ␣-helix in MAP2c fragments covering the corresponding sequence were low (31).

Overall shape from small-angle X-ray scattering
Sampling of the volume space by MAP2c was probed by small-angle X-ray scattering (SAXS). The scattering curves were acquired for various concentrations of unphosphorylated MAP2c and of MAP2c phosphorylated by PKA. A concentration dependence of the scattering intensity at low angles (corresponding to the momentum transfer q Ͻ 0.5 nm Ϫ1 ) was observed at concentrations of unphosphorylated MAP2c exceeding 0.1 mM (Fig. 2a). The concentration effects were greatly reduced when MAP2c was phosphorylated by PKA (Fig.  2b). This finding indicates that intermolecular interactions existing at submillimolar concentrations of MAP2c are suppressed by the PKA phosphorylation. The scattering data were further analyzed quantitatively using the ASTEROIDS approach (see below).

Paramagnetic relaxation enhancement
Long-range interactions of MAP2c were monitored using paramagnetic relaxation enhancement (PRE) of NMR signals (32). The native MAP2c contains a single cysteine at position 348, located in MTBD, providing a convenient site for attachment of the paramagnetic nitroxide radical label (1-oxyl-2,2,5,5-tetramethyl-⌬3-pyrroline-3-methyl)-methanethiosulfonate (MTSL). Six mutants, with Cys 348 replaced by serine and with a single cysteine located in the N-terminal acidic domain (E52C and Q110C), the proline-rich domain (T153C, E203C, Regions associated with functions of MAP2c, htau40, and both proteins are shown as cyan, pink, and gray boxes, respectively, and labeled as follows. a-d, f, and l, proposed neurosteroid-binding site (57); e, RII site; g, proline-rich region coded by exon 7; h, proline-rich region coded by exon 14 and phosphorylated by PKA; i, short proline-rich region P1 of MAP2c; j, proline-rich region P2; k, Class I Fyn-binding site; m, MTBR1; n, MTBR3; o, MTBR4; p, region RЈ; q, region homologous to the muscarinic receptor binding site of Tau; r, C-terminal ␣-helix; s, near-N-terminal insert I1; t, near-N-terminal insert I2; u, region including ␣-helix 114 LEDEAAGHVT 123 ; v, long proline-rich region P1 of Tau; w, Class II Fyn-binding site; x, MTBR2. Regions of Tau with high homology with MAP2c are drawn closer to the scheme of MAP2c. Red and yellow circles, residues phosphorylated by PKA with high and medium rate, respectively (11,13). SH3 recognition motifs are displayed in white boxes. Colored segments of the middle bar refer to clusters of spectral density values plotted in Fig. 7

Conformation and dynamics of MAP2c
and K282C), and in the C-terminal domain (L402C), were designed and labeled with MTSL.
To achieve high resolution, the peak intensities were monitored in 3D nonuniformly sampled HNCO spectra. Considering the lower sensitivity of the HNCO experiment compared with the 2D 1 H-15 N heteronuclear single quantum coherence spectra, we tested whether the lower tendency of MAP2c to aggregate allows us to work at higher concentrations than used in the PRE studies of Tau (23,24). PRE observed for 70 M unphosphorylated, 70 M phosphorylated, and 0.3 mM phosphorylated WT MAP2c were comparable. In contrast, a strong concentration dependence of PRE was observed for 0.3 mM unphosphorylated WT MAP2c, in agreement with the SAXS data indicating possible intermolecular interactions at high concentrations of unphosphorylated MAP2c (Fig. 2c). Therefore, only data obtained at the 70 M concentration were used for the quantitative conformational analysis (see below) in the case of unphosphorylated MAP2c. At this concentration, the sensitivity of the NMR experiment is still acceptable, and the concentration dependence of the SAXS curves is already weak.
The obtained PRE profiles (Fig. 3d) revealed the following interactions. Within the N-terminal domain, intramolecular interactions were observed between the vicinity of Glu 52 (in the negatively charged stretch of glutamates 52 EEEE 55 ) and residues in the vicinity of Trp 14 and Val 95 . In the C-terminal domain, weak intramolecular interactions were noticed between microtubule-binding repeats and proline-rich regions. Finally, long-range contacts were observed between the C-terminal domain and the N-terminal domain (namely in the vicinity of Trp 14 , Glu 52 , and Val 95 ), with a strong contribution of intermolecular (concentration-dependent) interactions. The long-range interactions are in agreement with previous observations of bent conformations and antiparallel dimers of MAP2c molecules (33).

Conformational analysis
To convert the obtained experimental data into a structural model, conformational analysis was performed using the ASTEROIDS algorithm (34). Assigned chemical shifts (13,25), SAXS, and PRE data served as an experimental input, reflecting local conformation, long-range contacts, and overall shape of the molecule, respectively. The potential presence of dimers (or higher oligomers) has a negligible effect on the description of local conformational behavior but may influence analysis of long-range interactions. To minimize the effect of intermolecular interactions, PRE and SAXS data measured at low concentrations were used in the case of unphosphorylated MAP2c. From a large set of 84,000 starting structures generated by the flexible-meccano program (see "Experimental procedures"), a set of 600 structures was selected to reproduce the experimental data. Despite the large conformational space of the 467residue polypeptide chain of MAP2c, good agreement between the experimental data and values calculated from the selected ensemble was achieved, as documented in Fig. 3.
To validate the ensemble of the selected MAP2c conformers, (i) values of 15 N secondary chemical shifts, ⌬␦(N), were removed from the selection procedure and back-calculated from the selected ensembles; and (ii) values of 1 D(NH) residual dipolar couplings (RDCs), measured as described under "Experimental procedures" but not used in the ASTEROID analysis, were calculated from the ensembles selected based on chemical shift, PRE, and SAXS data of MAP2c phosphorylated by PKA. In both cases, the independent data were well reproduced by the selected ensembles (Table 1 and Fig. 3 (a and b)). Similarly to the data published for htau40, 2 for the comparison of 259 experimental RDCs with the selected ensembles was a factor of 2 lower than 2 for the comparison with the pool of statistical coil structures. The difference between 2 values for

Conformation and dynamics of MAP2c
the selected and statistical coil ensembles further increased when the compared values were limited to a subset of residues with a higher difference in local sampling between the selected ensemble and the statistical coil pool (measured by the parameter ⌬ Rama , defined under "Experimental procedures"). The differences in 2 were even greater for ⌬␦(N), measured with a high precision in this study. The validation confirmed that the conformational description provided by ASTEROIDS is meaningful and predictive, as reported for other proteins (26). Local conformational propensities of MAP2c derived from the distribution of the and torsion angles, expressed as populations of three typical regions of the Ramachandran diagram in Fig. 4, were in general agreement with previously reported (13, 25) secondary structure propensities calculated using the SSP program (35). Similarly to SSP, ASTEROIDS derives the local conformations mostly from chemical shifts. However, the ASTEROIDS analysis of chemical shifts, based on the SPARTA program (36), is more detailed. Of particular interest is the population of the polyproline II conformation, that can be estimated on the basis of the currently available experimental data (37). The propensities to form particular secondary structures were often identified for continuous sequences of residues, although the ensemble selection treats each amino acid independently. To better visualize the continuous segments of secondary structures, we also plotted populations of short stretches of amino acids in the same conformation (orange and red bars in Fig. 4). Note that the important information is not the absolute value of population of the continuous stretch (which reflects composition of the starting pool of structures with the statistical distribution of torsion angles), but the relative increase of the population in the selected ensemble.
Highly populated continuous stretches of four and seven ␣-helical residues (corresponding to one and two turns of an ␣-helix, respectively) were identified in the C-terminal region (box r in Fig. 1). In addition, two regions were identified where a complete ␣-helical turn was observed in more than 5% structures of the selected ensembles, and at least five structures

Conformation and dynamics of MAP2c
(0.8%) with two turns were selected. These regions correspond to the binding site of the regulatory RII subunit of PKA (RII site, Asp 83 -Glu 113 , box e in Fig. 1) (5,31,38) and the region 206 SFSLNSSISSARR 218 (box h in Fig. 1). Most of the extended regions (i.e. exhibiting a negative SSP score (35)) were found to highly populate the polyproline II region of the Ramachandran plot. This applies not only to the proline-rich sequences, but also to a major portion of MTBD and residues between Leu 402 and Gln 421 . The highest population of continuous stretches of residues in the polyproline II conformation was observed for a short proline-rich segment encoded by exon 7 (box g in Fig. 1), both in the selected ensembles and in the statistical coil pool. In addition, two regions were identified where populations of four-and seven-residue stretches of polyproline II conformation exceed 30 and 5%, respectively, in the selected ensemble only, whereas the corresponding populations do not exceed 10 and 1% in the statistical coil pool, respectively. These regions are the sequence 120 PAALPLAA 127 and a binding site for the SH3 domain of the Fyn kinase (box k in Fig. 1) (27). The differences between secondary structures of unphosphorylated and phosphorylated MAP2c were very small.
We compared results of our conformational analysis with secondary structure predictions provided by several publicly available software tools. The predictors correctly identified the position of the C-terminal helix (box r in Fig. 1) and predicted the extended conformation in MTBD, but provided inconsistent results for the first 300 residues of MAP2c (data not shown). This finding is not surprising, considering that the predictors were developed for well-structured proteins, and it documents that experimental data are currently essential for conformational analysis of IDPs.
Long-range intramolecular interactions in the ensemble of structures selected by ASTEROIDS are presented in Fig. 5 as a map of relative distances between individual amino acids. The map shows that the average distances in the N-terminal region (including the vicinity of Trp 14 , the vicinity of Glu 52 , and the region Asp 83 -Glu 113 ) and between the N-terminal region and MTBD are shorter in the selected ensembles than in the statistical coil structures, reflecting only distribution of torsion angles typical for individual amino acids. As mentioned above, the analyzed SAXS and PRE data were obtained at conditions where formation of oligomeric species is greatly reduced but cannot be completely excluded. Therefore, we do not interpret the ASTEROIDS distance maps quantitatively under "Discussion," but we take into account also possible dimerization.

NMR relaxation analysis
Backbone dynamics of MAP2c were investigated by 15 N spin relaxation. To achieve sufficient resolution, relaxation rates were measured using 13 C, 15 N-labeled MAP2c and nonuniformly sampled 3D relaxation experiments based on the HNCO correlation (39). The R 1 , R 2 relaxation rates and steady-state heteronuclear Overhauser enhancement were measured with 600-and 950-MHz spectrometers. Because the R 2 rate is Population a α-helix P P P P P P P P P P PP P PP P P P P P P P P P PPPPPP P P PPP PP P P P

Conformation and dynamics of MAP2c
potentially influenced by chemical and/or conformational exchange and by the presence of 13 C (39), the R 2 experiment was complemented by a less sensitive, but more accurate, measurement of the transverse cross-correlated relaxation rate ⌫ x (a result of interference between 15 N chemical shift anisotropy and 1 H-15 N dipole-dipole interaction (40)) at 950 MHz. The relaxation rates obtained for unphosphorylated and PKA-phosphorylated MAP2c are presented in Fig. 6. Because MAP2c is an IDP, the measured relaxation rates represent average values for ensembles of all conformational and oligomeric states present at the given temperature. With the exception of the region Asp 83 -Glu 113 , phosphorylation by PKA has little effect on the measured rates. This indicates that the obtained values reflect mostly internal dynamics and not intermolecular interactions. The more ordered (more rigid) regions can be distinguished by higher values of the R 2 and ⌫ x relaxation rates. Such regions include the vicinity of Trp 14 , the vicinity of Glu 52 , residues Asp 83 -Glu 113 , serines in the proline-rich sequence 134 PPSPPPSP 141 , residues Arg 187 -Gly 197 and Glu 204 -Arg 226 , N-terminal regions of microtubule binding repeats and a homologous sequence following the microtubule-binding repeat 4 (MTBR4), the vicinity of the major phosphorylation site Ser 435 , and the C-terminal helix 456 EDVTAALAK 464 . The obtained relaxation data are in a good agreement with outputs of some disorder predictors (bottom panel in Fig. 6).
The measured relaxation rates were further analyzed to obtain values of the spectral density function J(), which is a Fourier transform of the time-correlation function (TCF) directly describing stochastic motions of the molecule (41)(42)(43)(44)(45). TCF can be rather complex in the case of disordered proteins, but it can be approximated as a series of weighted mono-exponential functions characterized by correlation times. Various relaxation rates are different linear combinations of values of the spectral density function J(), where the values of are given by the Larmor frequencies of the observed nuclei ( 1 H and 15 N). In this work, we calculated J(0) from R 1 , ⌫ x , and NOE measured at 950 MHz (46); J() for 61 and 515 MHz from R 1 and NOE measured at 600 MHz; and J() for 96 and 815 MHz from R 1 and NOE measured at 950 MHz (47,48). Note that the obtained set of spectral density values is not affected by the aforementioned inaccuracy of the measured R 2 rates.

Interaction with plectin SH3 domain
Regions of MAP2c responsible for interactions with several proteins have been characterized in the past. To extend the knowledge of binding sites, we examined interaction of MAP2c with a region of plectin containing a noncanonical SH3 domain. Interactions of disordered MAP2c with a well-folded fragment of plectin, consisting of spectrin repeats SR4 and SR5 and of the SH3 domain (29), were monitored as a decrease of MAP2c peak height in HNCO spectra after the addition of the SR4-SR5-SH3 plectin fragment (Fig. 9). Peak heights were reduced in several regions, most notably in MTBR3 and the region 208 SLNSSIS-SARR 218 . Canonical SH3 domains are known to interact with PXXP motifs and positively charged region (RXXK) motifs and via tertiary contacts involving hydrophobic residues (52). Comparison of the peak height decrease with distribution of prolines, hydrophobic residues, charged amino acids, and secondary structure motifs (Fig. 9) shows that positively charged residues, but not the PXXP motifs, are often found in the inter-

Conformation and dynamics of MAP2c
action sites. Interestingly, the 208 SLNSSISSA 216 sequence contains no prolines or positively charged amino acids but tends to form an amphiphilic helical structure.

Discussion
Intrinsic conformational behavior of MAP2c is an important factor in a complex network of intermolecular interactions con-  Fig. 1). The limits corresponding to mono-exponential TCF for 600-and 950-MHz spectrometers are plotted as dotted and solid gray curves, respectively.

Conformation and dynamics of MAP2c
trolling microtubule dynamics essential for neuron development. Therefore, we examined how the observed regions of distinctive conformations and dynamics correlate with known and proposed functional motifs of MAP2c. We also compared the properties of MAP2c with its close homologue Tau (clone htau40, splicing variant 2N4R), exhibiting different expression and phosphorylation patterns and cellular localization and playing an important role in neurodegenerative processes associated with Alzheimer's disease. The determined populations of conformations and dynamics did not differ substantially between MAP2c and Tau in regions of high sequence homology (boxes m-r in Fig. 1), including MTBD. Therefore, we can conclude that the structural basis of interactions with microtubules and their competitors (such as 14-3-3 proteins) is similar for MAP2c and Tau. Results providing new insight into the binding of MAP2c to its important interaction partners were obtained for regions of lower sequence similarity with Tau and are discussed below.

Site of interaction with canonical Fyn SH3 domain
Fyn is an important kinase containing an SH3 domain and phosphorylating single tyrosine residues in different regions of the compared proteins: Tyr 18 in Tau (12) and Tyr 67 in MAP2c (53). In general, SH3 domains bind to proteins with the sequential motif PXXP, preferably with positively charged amino acids Relative peak height P P P P P P P P P P PP P PP P P P P P P P P P PPPPPP P P PPP PP P P P P P P 20   Secondary structure motifs determined by the ASTEROIDS analysis, hydropathicity, and charge distribution are shown above the plot. The symbols used to describe the secondary structure motifs are explained in the legend to Fig. 4. The hydrophaticity index according to Kyte and Doolittle (85) is shown as darkness of the top bar above the plot (white and black correspond to values of Ϫ4.5 and ϩ4.5, respectively). The charge distribution is represented by the color in the bottom bar above the plot, corresponding to a relative electrostatic potential approximated by ⌺ j CQ i /(d 0 ϩ d 1 ͉n i Ϫ n j ͉), where Q i and n i are charge and sequential number of the ith residue, C is a constant including the electric permitivity, and d k are distance constants. The ratio d 1 /d 0 was set to 2.0, and the colors were chosen so that red and blue correspond to the highest negative and positive potential, respectively, which makes the color code independent of C/d 0 (13). Letters P above the plot indicate positions of prolines.

Conformation and dynamics of MAP2c
in its vicinity. Tau and MAP2c contain 7 and 13 PXXP motifs, respectively, most of them in the proline-rich region P2 (box j in Fig. 1) with a high sequence homology. P2 of Tau contains two Fyn-binding sites, the classical class II motif 216 PTTPTR 221 (box w in Fig. 1) and the classical Class I motif 230 RTPPKSP 236 (box k).
Only the latter site is present in MAP2c ( 288 RTPPKSP 294 ). Binding assays performed with synthetic biotinylated peptides showed that the Fyn SH3 domain binds preferentially to the Class II site of Tau (28) and Class I site of MAP2c (27). Our analysis (Fig. 4) revealed that the selected ensembles contained 6% MAP2c structures with all residues of the 288 RTPPKSP 294 motif in the polyproline II conformation, optimal for binding to the Fyn SH3 domain. This number is 25 times higher than the corresponding population in the statistical coil pool. Only two other sites outside the proline-rich regions P1 and P2 have populations of seven-residue polyproline II stretches higher than 5% (one of them being 134 PPSPPPSP 141 with the discussed conformation highly populated already in the statistical coil pool). This shows that (i) the 288 RTPPKSP 294 sequence preferentially adopts the conformation optimal for binding already in the free form of MAP2c, and (ii) the 288 RTPPKSP 294 site differs from the other PXXP motifs not only in charge distribution, but also by well-defined conformation. The relaxation data (Fig. 6) also show somewhat higher ordering for the 288 RTPPKSP 294 site.

Site of interaction with the noncanonical plectin SH3 domain
Interactions with plectin where proposed to regulate the microtubule-stabilizing activity of MAP2c. Binding of MAP2c to the plectin SH3 domain and longer plectin fragments was observed experimentally (30). X-ray crystallography revealed that the plectin SH3 domain adopts the standard fold, but it lacks some aromatic side chains involved in canonical interactions with PXXP motifs, and its PXXP-binding site is partially occluded (29). It was therefore proposed that the plectin SH3 domain does not recognize the classical PXXP motifs. Here, we tested this hypothesis for intrinsically disordered MAP2c by solution NMR spectroscopy. Decrease in NMR peak height, indicative for MAP2c interactions with the well-folded SR4-SR5-SH3 plectin fragment, was not observed for the PXXP motifs, but in other regions (Fig. 9). Some of the observed decrease can be explained by electrostatic interactions of positively charged MAP2c residues in MTBD (e.g. of 298 KQLR 301 in MTBR1 and 352 KNIR 355 in MTBR3, resembling the alternative SH3-binding motifs RXXK (52)). However, significant binding was observed also in the region 208 SLNSSISSARR 218 , preceding a positively charged sequence, 217 RRTTRSEPIRRAGK 230 . We used the conformational analysis to examine whether this region can be distinguished by its structure. The analysis revealed (Fig. 4) that the sequence 206 SFSLNSSISSARR 218 is one of three regions with the highest population of ␣-helical stretches (together with the C-terminal helix and with the RII site). Based on this finding, we propose that the SR4-SR5-SH3 plectin fragment recognizes also an amphiphilic helix consisting of residues 206 SFSLNSSISSARR 218 and present in free MAP2c.

Site of interaction with the PKA regulatory domain
The sequence between Asp 83 and Glu 113 of MAP2c (box e in Fig. 1) aligns with a motif binding both the RI-and RII-regula-tory subunit of PKA (54,55), and the binding to the RII subunit is well documented (5,31,38). Unlike the canonical SH3 domains, dimers of the regulatory subunits of PKA form X-type helical bundles, recognizing relatively long amphiphilic ␣helices (at least five turns) of well-folded protein kinase A-anchoring proteins (55). Therefore, the obvious question was whether such a helix is also significantly populated among conformers of free MAP2c. Increased population of ␣-helical conformation was observed in the selected ensembles for residues Thr 80 -Ser 91 and Lys 107 -Gln 110 . In the middle part of the motif, the helical population in the selected ensembles did not exceed the statistical coil values, whereas the polyproline II conformation was more populated. A particularly low population of ␣-helix was observed for Gln 96 and Val 98 . It shows that the RII site does not form a highly populated, uniform secondary structure in free MAP2c, as observed for the SH3-binding sites. It should be emphasized that the lack of a uniform secondary structure does not imply higher flexibility of the RII site. On the contrary, relaxation data revealed that the whole RII site is significantly more ordered than other parts of MAP2c. Contacts observed in PRE data suggest that hydrophobic interactions with aromatic residues in the N-terminal region may be involved. Taken together, our results indicate that the RII site of MAP2c represents a relatively complex local structure, qualitatively different from the SH3-binding sites. Further investigation will be needed to find out whether a regular helical structure is formed in complex with the regulatory subunit of PKA. In any case, it is interesting that the ␣-helical conformation is more populated in terminal regions of the motif, whereas the polyproline II conformation prevails in the interior part, which is supposed to form the amphiphilic helix critical for the interaction.

Global conformation
The obtained experimental data not only describe local conformational motifs, discussed above, but also provide information about global structural features of MAP2c. The overall shape of the MAP2c molecule is mostly given by electrostatic interactions between acidic N-terminal and positively charged C-terminal domains. The structural effect of the intramolecular electrostatic interactions is formation of bent "paper clip" conformations (56). At higher concentration, formation of antiparallel dimers is expected based on the charge distribution. Both types of structures have been observed earlier by Wille et al. (33). The antiparallel nature of the dimers is important for the interpretation of the PRE data. Whereas PRE observed between the N-terminal and C-terminal regions can be due to intramolecular and/or intermolecular interactions, the observed PRE inside the N-terminal region can be reliably attributed to the intramolecular contacts, discussed above. Regardless of the possible contribution of intermolecular contacts to the PRE data, it is clear that the actual interacting regions of MAP2c and htau40 differ. The contacts of MTSLlabeled E52C with the vicinity of Trp 14 and the RII site are in agreement with the proposed model of the hydrophobic neurosteroid binding site and its interactions with the RII site (57). Most of the other differences can be explained by different charge distributions in htau40 and MAP2c (13). In the cellular environment, the intramolecular contacts contribute to the

Conformation and dynamics of MAP2c
delicate equilibrium of interactions related to microtubule dynamics (58), aggregation of Tau (23), and interactions with other partners. The biological relevance of the intermolecular interactions is less clear. Quantitative data (53, 59 -63) show that prenatal cytosolic concentrations of MAP2c are comparable with the concentrations of Tau in adult neurons (5-10 M) (63), but formation of dimers can be expected in regions with locally increased MAP2c concentration (64). Interestingly, intermolecular interactions in MAP2c are reduced by PKA phosphorylation. Notably, PKA phosphorylation sites differ between MAP2c and Tau. For example, the major PKA phosphorylation site in MAP2c, Ser 435 , is localized close to the C terminus in a slightly positively charged region, and its phosphorylation makes interactions with the acidic N-terminal region in the antiparallel dimer less favorable.
In summary, high-resolution NMR experiments allowed us to study structural features and motions of a vast majority of MAP2c residues. Despite the high complexity of the conformational ensembles of MAP2c, the obtained data could be converted into parameters describing populations of conformations, long-range contacts, and average dynamics of individual amino acids using calibrated analytical procedures (37) to derive representative ensemble descriptions of the conformational space sampled by AP2c. Results of the analysis and NMR relaxation data revealed lower flexibility and increased populations of specific conformations in regions of MAP2c interacting with other proteins. They indicate that the local conformational propensities may play an important role in the protein partner recognition. As tested here for MAP2c, these conformational propensities cannot be reliably predicted for IDPs by currently available bioinformatics tools. Therefore, the presented type of structural ensemble analysis based on a set of experimental (NMR) data is essential for discovering yet unexplored binding motifs.

Spin labeling
After purification, the cysteines were reduced by incubating MAP2c with 0.5 mM TCEP. For spin labeling, TCEP was removed by size-exclusion chromatography using the HiTrap desalting column (GE Healthcare), and the protein was incubated in 20 mM potassium phosphate, 150 mM potassium chloride, pH 7.4. A 10-fold molar excess of MTSL (Toronto Research Chemicals) dissolved in DMSO was added, and MTSL was allowed to bind overnight at 4°C. The protein was then dialyzed into NMR buffer (50 mM MOPS, 150 mM sodium chloride, pH 6.9) before measurement. Paramagnetic MTSL label was converted into diamagnetic MTS by adding a 3-fold molar excess of ascorbic acid to the sample and incubating for 15 min at room temperature. 3D HNCO spectra were measured for the paramagnetic and diamagnetic MAP2c mutants. The PRE profiles were obtained by calculating the ratio of peak intensities of the paramagnetic and diamagnetic samples.

NMR spectroscopy
NMR experiments were acquired using 600-, 850-, and 950-MHz Bruker Avance III spectrometers equipped with TCI cryogenic probe heads with z axis gradients and a 700-MHz Bruker Avance III spectrometer equipped with a 1 H/ 13 C/ 15 N TXO cryogenic probe head with z axis gradients. All experiments were performed at 27°C with the temperature calibrated according to the chemical shift differences of pure methanol peaks. The indirect dimensions in 3D and 5D experiments were acquired in a nonuniformly sampled manner. On-grid Poisson disk sampling with a Gaussian probability distribution (66) was applied.
The 3D (CACO)NCACO spectrum and the 5D CACONCACO spectrum (25) were acquired as described earlier (13) to confirm assignment of mutants used in the PRE study. The 3D HNCO (67) spectra were acquired with spectral widths set to 18,939 (direct) ϫ 2000 ( 15 N) ϫ 2000 ( 13 CЈ) Hz and maximal evolution times of 120 ms for 15 N and 80 ms for 13 CЈ indirectly detected dimensions. The overall number of 2048 complex points was acquired in the acquisition dimension, and 2000 hypercomplex points were randomly distributed over the indirectly detected dimensions. Interactions with the plectin SH3 domain were observed by comparing nonuniformly sampled HNCO spectra of samples containing 70 M MAP2c with and without 70 M SR4-SR5-SH3 plectin fragment in 50 mM MOPS, pH 6.9, 150 mM NaCl, and 0.7 mM TCEP. The spectra were recorded at 850 MHz and 10°C (due to the limited thermal stability of the SR4-SR5-SH3 plectin fragment). Peaks with low intensity in the spectrum of free MAP2c were excluded from the analysis. 1 D(NH) RDCs in PKA-phosphorylated MAP2c were measured at 600 MHz using a nonuniformly sampled 3D in-phase anti-phase HNCO experiment (68) with spectral widths of 2000 ( 15 N) ϫ 2000 ( 13 CЈ) Hz in the indirect dimensions, with maximal evolution times of 180 ms for 15 N and 120 ms for 13 CЈ, and with 600 hypercomplex points randomly distributed over the indirectly detected dimensions. The 1 D(NH) values were calculated as the difference between the splitting in spectra of phosphorylated MAP2c in the isotropic buffer and in a stretched 5% polyacrylamide gel. The relaxation rates of 0.9 mM unphosphorylated and phosphorylated MAP2c were measured using 3D HNCO-based relaxation experiments as described recently (39). The spectral density values J( N ) and J(0.86 H ) (i.e. J(61 MHz) and J(515 MHz), respectively, on the 600-MHz spectrometer and J(96 MHz) and J(815 MHz), respectively, on the 950 MHz spectrometer) were calculated from R 1 and steady-state NOE using reduced spectral density mapping (47,48). The J(0) values were calculated from R 1 , ⌫ x , and steady-state NOE was measured at 950 MHz using the LNX protocol described previously (46). The graphical spectral density analysis was performed as described earlier (46,49,51).

SAXS measurement
SAXS data were collected on the beamline BM29 BioSAXS ESFR in Grenoble, France. The highest concentrations of MAP2c in the dilution series were 12.0 and 13.6 mg/ml for unphosphorylated and phosphorylated sample, respectively. The data were recorded using the pixel 1M PILATUS detector at a sample-detector distance of 2.867 m, and a wavelength () of 0.099 nm, covering the range of momentum transfer 0.025 nm Ϫ1 Ͻ q Ͻ 5 nm Ϫ1 (q ϭ 4sin()/, where 2 is the scattering angle). 50 l of the MAP2c protein solution was loaded into a flow-through quartz capillary cell at 20.12°C. The total exposure time was 10 s/sample. Several successive frames (10 for MAP2c and 20 for the buffer) were recorded for both the sample and the buffer. Each frame was inspected to check for the possible presence of protein damage before calculating average intensities. Solvent contributions (buffer backgrounds collected before and after the protein sample) were averaged and subtracted from the associated protein sample. No radiation damage was observed during the data collection. The data were processed using standard procedures with PRIMUS (69).

Conformational analysis
The flexible-meccano program (70,71) was used to generate a large ensemble of conformers of the protein. For each conformer in the ensemble, the side-chain and hydrogen atoms were added to the model using SCCOMP (72). The randomcoil torsion angle database was used by flexible-meccano to create an ensemble of 10,000 conformers. Moreover, a database of torsion angles filtered to be in agreement with experimental chemical shifts of 15 N, 13 CЈ, 13 C ␣ , and 13 C ␤ was prepared to enhance convergence of ASTEROIDS using flexible-meccano/ ASTEROIDS selection and used by flexible-meccano to create another ensemble of 74,000 conformers. Both ensembles were joined together to define a starting pool of 84,000 conformers, and experimental observables were calculated for each conformer. Chemical shifts were predicted using SPARTA (36), SAXS using CRYSOL (73)(74)(75), and PRE as described previously (76). Local alignment windows of 15 amino acids in length were used to calculate averaged RDCs (77). These were then combined with a generic baseline to account for long-range effects (34,76). The starting pool was used for selection of a sub-ensemble representing the disorder protein by ASTEROIDS (34). The optimal number of structures in the selection was estimated from cross-validation procedures to be ϳ600. To achieve convergence, 200,000 iterations of the genetic algorithm were required. The ASTEROIDS selection was repeated five times with the same set of data to confirm reproducibility of torsion angle distribution in the selected ensembles and of the back-calculated chemical shift, PRE, and SAXS data. Differences between statistical coil (COIL) sampling and selected ensembles (AST) were characterized as described previously (26) with the metric ⌬ Rama defined for each residue as follows, where p COIL and p AST are statistical coil and selected populations, respectively, and k covers four regions of the Ramachan-dran diagram (␣-helix, ␤-sheet, polyproline II, and the remaining region including the left-handed ␣-helical conformation).