PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of plosbiolPLoS BiologySubmit to PLoSGet E-mail AlertsContact UsPublic Library of Science (PLoS)View this Article
 
PLoS Biol. Feb 2009; 7(2): e1000034.
Published online Feb 17, 2009. doi:  10.1371/journal.pbio.1000034
PMCID: PMC2642882
Structural Polymorphism of 441-Residue Tau at Single Residue Resolution
Marco D Mukrasch,#1 Stefan Bibow,#1 Jegannath Korukottu,1 Sadasivam Jeganathan,2 Jacek Biernat,2 Christian Griesinger,1 Eckhard Mandelkow,2 and Markus Zweckstetter1,3*
1 Department for Nuclear Magnetic Resonance (NMR)-Based Structural Biology, Max Planck Institute for Biophysical Chemistry, Göttingen, Germany
2 Max Planck Unit for Structural Molecular Biology, Hamburg, Germany
3 Deutsche Forschungsgemeinschaft (DFG) Research Center for the Molecular Physiology of the Brain (CMPB), Göttingen, Germany
Gregory A Petsko, Academic Editor
Brandeis University, United States of America
#Contributed equally.
* To whom correspondence should be addressed. E-mail: mzwecks/at/gwdg.de
Received May 21, 2008; Accepted January 7, 2009.
Alzheimer disease is characterized by abnormal protein deposits in the brain, such as extracellular amyloid plaques and intracellular neurofibrillary tangles. The tangles are made of a protein called tau comprising 441 residues in its longest isoform. Tau belongs to the class of natively unfolded proteins, binds to and stabilizes microtubules, and partially folds into an ordered β-structure during aggregation to Alzheimer paired helical filaments (PHFs). Here we show that it is possible to overcome the size limitations that have traditionally hampered detailed nuclear magnetic resonance (NMR) spectroscopy studies of such large nonglobular proteins. This is achieved using optimal NMR pulse sequences and matching of chemical shifts from smaller segments in a divide and conquer strategy. The methodology reveals that 441-residue tau is highly dynamic in solution with a distinct domain character and an intricate network of transient long-range contacts important for pathogenic aggregation. Moreover, the single-residue view provided by the NMR analysis reveals unique insights into the interaction of tau with microtubules. Our results establish that NMR spectroscopy can provide detailed insight into the structural polymorphism of very large nonglobular proteins.
The Tau protein, which plays a central role in the progression of Alzheimer disease, is normally expressed in nerve axons, where it stabilizes microtubules (MTs), supports the outgrowth of axons, and modulates the transport of vesicles and organelles along MTs. In Alzheimer disease, Tau becomes excessively phosphorylated, loses its ability to bind to MTs, and aggregates into intracellular abnormal protein deposits. Many efforts have been made over the years to understand Tau structure as a way to understand Tau function and its mechanisms of action, but these efforts have primarily used traditional biochemistry and molecular biology approaches and therefore have addressed structure and function at a relatively primitive level. Here, we show that it is possible to characterize the structure and dynamics of 441-residue Tau at single residue resolution using nuclear magnetic resonance (NMR) spectroscopy. NMR spectroscopy demonstrates that 441-residue Tau is highly dynamic in solution with a distinct domain character and an intricate network of transient long-range contacts important for pathogenic aggregation. Moreover, the single-residue view provided by the NMR analysis reveals unique insights into the interaction of Tau with MTs.
Tau protein was originally discovered as a neuronal microtubule-associated protein (MAP) that stabilizes microtubules (MTs) and supports the outgrowth of axons [1,2]. The protein can modulate the transport of vesicles and organelles along MTs, serves as an anchor for enzymes, and regulates the dynamics of MTs [3,4]. In Alzheimer disease (AD), tau becomes excessively phosphorylated, looses its ability to bind to MTs, and aggregates into neurofibrillary tangles that consist of paired helical filaments (PHFs) of tau. Mutations in the tau gene cause tau aggregation and frontotemporal dementia with parkinsonism linked to Chromosome 17 [5,6].
The human central nervous system contains six isoforms of tau, generated from a single gene by alternative splicing and ranging between 352 to 441 amino acid residues [7,8]. The isoforms differ by two inserts near the N-terminal end and the presence of either four or three imperfect repeat sequences in the C-terminal half of the protein (Figure 1A). The repeat domain represents the core of the MT interaction [9] and also forms the core of PHFs [10]. For PHF aggregation two hexapeptides at the beginning of the second and third repeats (275VQIINK280 and 306VQIVYK311) are crucial because they are able to initiate the aggregation process [11].
Figure 1
Figure 1
Primary Sequence of 441-Residue Tau
It was recognized early on that tau has an unusual character as a protein, because it was resistant to heat and acid treatment without loosing its function and had a very low content of secondary structure [12]. These properties can be traced back to the high fraction of basic and hydrophilic amino acid residues (Figure 1B), which resist the compact folding typical of most proteins. In fact, a number of biophysical studies revealed that tau is a prototypical “natively unfolded” protein [13]. In recent years, this type of protein emerged as a major fraction in the human proteome (termed “natively unfolded” or “intrinsically unstructured proteins” [IUPs] [14]). Apart from tau, most “fibrous” MAPs have the signature of the natively unfolded state, whereas other MT-binding proteins show conventional folding (e.g., motor proteins).
Since disordered proteins tend to be highly flexible and have variable conformations, they have not been amenable for structure analysis by crystallography. Thus nuclear magnetic resonance (NMR) spectroscopy is the only method that allows a description of their conformations and dynamics with high resolution [15]. The lack of an ordered structure, however, causes dramatic signal overlap. Therefore, we and others have previously performed NMR studies on fragments of tau or studied full-length tau but only on the basis of partial assignment that was scattered throughout the sequence [1625]. In particular by studying tau fragments that contain only the repeat domain (K19 or K18) or the repeat domain and the flanking regions (K32), we and others showed that the hexapeptides in repeats R2 and R3 populate β-structure and bind to MTs and polyanions [1618], that short stretches in the repeat domain assume highly populated turn conformations [19], that the repeat domain of tau folds into an α-helical conformation upon binding to lipid surfaces [20], and that PHFs formed in vitro by the three-repeat-domain (K19) of tau consist of three major β-strands [21]. Moreover, using a partial assignment of full-length tau (less than 40%) Lippens and coworkers investigated the binding of tau to MTs [22], the phosphorylation pattern of tau as induced by cAMP dependent kinase [23], tau aggregated into PHFs [24], and the impact of binding of heparin to tau [25]. In addition, small angle x-ray scattering (SAXS) and Förster resonance energy transfer (FRET) was used to obtain insight into the structure of the tau protein [26,27].
Despite the wealth of information from previous studies on the conformational properties of tau, they were always limited because they were either not of high resolution (SAXS, circular dichcroism, electron microscopy), were restricted to fragments of tau (liquid-state and solid-state NMR, x-ray crystallography of a complexed tau peptide [28]), or were limited to a subset of residues (NMR, FRET). In contrast, we show here that it is possible to obtain the complete backbone assignment of 441-residue tau (the longest tau isoform found in the human central nervous system, htau40; Figure 1B) and thus to overcome the size limitation that in the past has limited detailed NMR studies of unfolded proteins to fewer than 200 amino acids (Figure 2A) [16,29]. The complete backbone assignment of htau40 allowed us to probe the structure and dynamics of the full-length soluble protein, including the 198 residues of the N-terminal half and the 47 residues of the C-terminal domain, and determine at single residue-resolution the residues involved in the interaction between tau and MTs. Most importantly, the data provide unique insights into long-range interactions between remote regions of tau that can be studied only in the context of the full-length protein.
Figure 2
Figure 2
441-Residue Tau at Single-Residue Resolution
Backbone Resonance Assignment of 441-Residue Tau
For 441-residue htau40, we observed a narrow, highly congested cluster of amide proton signals in a 1-D NMR spectrum (Figure 2B). Correlation with the directly attached 15N-amides in a two-dimensional 1H-15N heteronuclear single quantum coherence (HSQC) spectrum, only partially resolved the degeneracy (Figure 2B): the number of overlapping signals was still a factor of 3.5 higher than in the largest currently assigned globular protein, the 731-residue malate synthase G (Figure 2C). The large number of proline residues (43 out of 441 residues) and the strongly repetitive primary sequence in the repeat domain of htau40 also complicated the analysis of sequential connectivity. Within a range of +/−0.2 ppm we observed 20.1 residues on average (Figure S1).
To obtain the sequence-specific assignment of the backbone resonances of htau40, we recorded 3-D (HA)CANNH [30] and HNN [31] experiments. For nitrogen and Cα nuclei, very high resolution was obtained at the second highest currently commercially available magnetic field (21.14 T) (Figure 2E), and more than 98% of non-proline backbone resonances for the full length htau40 protein were assigned. Thus, htau40 exceeds more than 2-fold the largest currently assigned disordered protein. Only Gly272, Gly303, Gly334, and Gly366, which are found at the C-terminal end of each repeat region and are surrounded by two glycines in the sequence motif PGGG, as well as Gly304 and Gly335 could not be assigned unambiguously owing to severe signal overlap. In addition, the resonances originating from Met1 and Ala2 were not observed in 1H-15N HSQC spectra. In case of proline, more than 83% of Cα frequencies were assigned. The assignment was corroborated by producing three overlapping fragments (Figure 2D): (i) a 185-residue fragment comprising the N-terminal half up to the repeat region, but excluding the two inserts that are affected by alternative splicing (I1, I2, encoded by exons 2 and 3); (ii) a 198-residue fragment containing the repeat region and the two proline-rich flanking regions; (iii) a 168-residue fragment covering most of the C-terminal half except for the second repeat R2. The two 29-residue inserts in the N-terminal half were only present in htau40. Superposition of HSQCs of the three fragments with that of htau40 showed that many resonances observed in the three fragments were found at identical positions as in htau40 (Figure 2D), in agreement with its high flexibility. This dataset forms the basis for probing intramolecular interactions and studying the interactions between tau and its binding partners.
Secondary Structure Propensity in 441-Residue Tau
NMR spectroscopy provides a variety of probes that are highly sensitive to backbone dihedral angles in both globular and disordered proteins [32]. We used experimental Cα chemical shifts and 3J(HNHα) couplings, from which random coil values were substracted to reveal the presence of helical or β-structure. For all residues of htau40, Cα secondary chemical shifts were below 1.5 ppm (positive or negative) (Figure 3A), indicating that rigid secondary structure elements are not present in htau40. However, several continuous stretches (containing 7–11 residues) with negative Cα secondary chemical shifts were observed in the repeat region (Figure 3A and and3C),3C), indicative of a propensity to adopt β-structure. According to a quantitative analysis, the β-structure-like conformations are populated 12%, 22%, 25%, 19%, and 12% of the time for residues 256VKSKIG262 (in R1), 274KVQIINKKLDL284 (in R2), 305SVQIVYKPVDL315 (in R3), 336QVEVKSEKLD345 (in R4), and 351QSKIGSL357 (in R4), respectively. In addition, stretches of negative Cα secondary chemical shifts consisting of more than five amino acids were found for 86GKQAAAQ92 (17% in I2), 161GQKGQA166 (17% in P1), and 224KKVAVVR230 (18% in P2). Thus, the highest β-structure content was found for residues 274KVQIINKKLDL284 and 305SVQIVYKPVDL315 in repeats R2 and R3, comprising the two aggregation-prone hexapeptides 275VQIINK280 and 306VQIVYK311. Formation of β-structure in the homologues region of R1 (243LQTAVMPDL253) is hindered by the presence of three proline residues (Figure 1B).
Figure 3
Figure 3
Residual Secondary Structure of Tau in Solution
Continuous stretches of positive secondary chemical shifts report on α-helical propensity and were observed for 114LEDEAAGHVT123 (between insert 2 [I2] and the proline-rich region P2) and 428LADEVSASLA437 in immediate proximity to the C terminus (Figure 3A, A,3D,3D, and and3E).3E). Quantitative analysis revealed 18% of α-helical population for 114LEDEAAGHVT123 and 25% for 428LADEVSASLA437. Mapping of these two residue stretches onto a helical wheel reveals two amphiphatic helices with more hydrophobic residues on one side of the helical cylinder and an excess of negative charges on the opposite side (Figure 3D and and33E).
The 3JNH-αH-coupling of a residue depends on its [var phi] backbone torsion angle. Positive Δ3JNH-αH (Jexp−Jrandom coil) values indicate a tendency towards extended and β-structure, low or even negative values indicate turns or helical propensity. In htau40, positive Δ3JNH-αH values dominate along the entire sequence indicating its overall extended chain conformation (Figure 3B). The largest 3JNH-αH values were detected for residues 305SVQIVYKPVDL315, in agreement with the highest β-structure propensity (25%) as estimated from Cα secondary chemical shifts. 274KVQIINKKLDL284 also showed increased positive Δ3JNH-αH values, but the effect was less pronounced. In contrast, negative Δ3JNH-αH values or values close to zero were found for 114LEDEAAGHVT123 and 428LADEVSASLA437 (Figure 3B, B,3D,3D, and and3E),3E), supporting the preferential population of α-helical conformations by these residues. In addition, small 3JNH-αH-couplings were detected for several non-proline residues in the two proline-rich regions P1 and P2 (Figure 3B and and3F):3F): Thr175, Ala178, Thr181, Ser184, Thr217, Thr220, Glu222, Lys224, Lys234, and Ser235. Most of these residues also showed negative or very small positive Cα secondary chemical shifts (Figure 3F), suggesting that 175TPPAPKTPPS184, 216PTPPTREP223, and 232PPKSPSSA239 transiently assume polyproline II helical conformations (Figure 3G). Two of these motifs (216PTPPTREP223 and 232PPKSPSSA239) are separated by a stretch of positive Δ3JNH-αH values (Figure 3F), consistent with the β-structure propensity of 224KKVAVVR230 that was suggested by Cα chemical shifts.
Residual dipolar couplings (RDCs) [33] report on time and ensemble-averaged conformations [34] and can be used to understand both the structure and dynamics of disordered proteins [35]. By weakly aligning htau40 in Pf1 bacteriophage, we could determine 262 residual one-bond H-N dipolar couplings (Figure 4). For other residues, peak overlap prohibited a quantitative analysis. In addition, most residues in the proline-rich region P2 (residues 198–244) as well as residues 171–183 in the proline-rich region P1 showed very strong alignment prohibiting a quantitative analysis of their dipolar couplings. Large residual one-bond H-N dipolar couplings were observed in the repeat domain and the proline-rich regions P1 and P2. In the repeat region, the largest values were found for the hexapeptide 306VQIVYK311 in the beginning of R3 (Figure 4A). Large positive H-N RDCs are associated with locally more extended conformations [35], in qualitative agreement with β-structure propensity of 306VQIVYK311.
Figure 4
Figure 4
RDCs Observed in Weakly Aligned Tau
Negative H-N RDCs were observed for 43 residues (Figure 4A). In particular, residues 430–438 showed negative RDC values, with the most negative value (−12.4 Hz) found for S433 (Figure 4A and and4C).4C). The sign inverted RDCs indicate that the H-N internuclear vectors of 428LADEVSASLA437 are parallel to the long axis of this segment [36], in agreement with the presence of a helical conformation. Isolated residues with sign-inverted RDCs are characteristic for the presence of turns in disordered proteins [19]. In htau40, most of the negative H-N RDCs belong to residues in the N-terminal part (residues 1–200), suggesting a higher number of turns in this region. Previously, we showed by a combination of H-N RDCs and molecular dynamics simulation that the four peptides 252DLKN255, 283DLSN286, 314DLSK317, and 345DKFD348, embedded in a fragment that only comprised the repeat domain of tau (K18), showed high propensities to form turns [19]. In htau40, the peak overlap was strongly increased (Figure 2), and we could analyze reliably only K347, which showed a H-N RDC value of −4.6 Hz (Figure 4B). However, the Cα secondary chemical shifts observed for 252DLKN255, 283DLSN286, 314DLSK317 in K32 and htau40 were nearly identical (Figure 5A), suggesting that the four peptides also populate turn conformations in htau40.
Figure 5
Figure 5
Robustness of Secondary Structure Propensity of Tau
Temperature and Construct Dependence of Secondary Chemical Shifts
Comparison of the H-N HSQC spectra showed that most cross peaks of the three fragments (K25, K32, and K10; Figure 1) were found in very similar positions as in the htau40 spectrum. To further probe the influence of flanking domains on the structural propensities of different regions of htau40, we compared secondary chemical shifts of the three fragments with values observed in full-length tau. Close agreement was found between Cα secondary chemical shifts of K25, K32, K10, and htau40 (Figure 5A and S2). The largest deviation between any of these fragments and htau40 was 0.32 ppm. The rmsd values were 0.04, 0.05, 0.04 ppm for the comparisons K32-htau40 (171 residues), K25-htau40 (168 residues), and K10-htau40 (150 residues), respectively.
To further probe the robustness of the local conformational properties of htau40, we determined the sequence-specific assignment of backbone resonances at 20 °C. Backbone resonance assignment at 20 °C was achieved by following the shifts of cross peaks in H-N HSQCs of htau40 supported by 3D (HA)CANNH spectra of K25, K32, K10, and htau40 at 20 °C. The Cα secondary chemical shifts observed at 20 °C in htau40 were highly similar to the values observed at 5 °C (Figure 5B). We conclude that the structural propensities of monomeric tau are highly specific conformational fingerprints.
Flexibility of the Backbone of htau40
To probe the dynamics of the backbone of htau40, we measured spin relaxation rates. 15N R spin relaxation rates allow quantification of motions that occur on timescales of pico- to nanoseconds and micro- to milliseconds and reflect the flexibility of the protein in solution [37]. For the N-terminal domain up to residue 170, R rates were below 4.7 Hz with an average value of 3.8 Hz (Figure 6), indicating that this part of tau is highly mobile. In the proline-rich region P2, the R values increased and reached a maximum of 5.2 Hz for S235, indicating increased rigidity. The observed maximum is part of the 232Pro-Ala239 stretch that transiently populates polyproline II helical conformations, suggesting this forms a more periodic and less flexible structure. Similarly, we also observed R spin relaxation rates above 4.7 Hz for many residues that belong to elements of transient β-structure in the repeat domain (see above): Ile277, Asn279, and Leu282 (in R2), 309Val-Lys321 (in R3), and 343Lys-Lys353 (in R4). The largest R rates were detected for 370Lys-Lys395 in the region downstream of the repeat domain.
Figure 6
Figure 6
Intrinsic Flexibility of Tau
RDCs are not only excellent probes for structure, but are also sensitive to motions from picoseconds to milliseconds [38]. The large H-N RDC values, which were observed in the repeat domain (Figure 4), arise from locally more extended conformations that at the same time increase the rigidity of this region. Interestingly, R relaxation rates are not elevated in the proline-rich domain P1, whereas H-N RDCs in this region are a factor of two or more larger than for residues 1–140. This suggests that proline residues restrict the mobility of the backbone in the time window between the global correlation time of the protein and 50 μs that is invisible to NMR relaxation measurements.
Global Folding of Soluble Tau
To study the global folding of htau40, we employed paramagnetic relaxation enhancement (PRE) of NMR signals [39]. The primary sequence of htau40 contains two cysteines (C291 and C322) that provide convenient attachment points for the nitroxide radical (1-oxy-2,2,5,5-tetramethyl-D-pyrroline-3-methyl)-methanethiosulfonate (MTSL). In addition, five different mutants, which harbour only a single cysteine in the projection domain (A15C or A72C), in the proline-rich region (A239C), or near the C terminus (A384C or A416C), were constructed and labelled with MTSL. Figure 7 shows the PRE profile (ratio of NMR signal intensities in the paramagnetic and diamagnetic state) of the amide protons of htau40 for the six different MTSL-labelled htau40 samples. For a fully extended chain, the NMR signal intensities in the paramagnetic and diamagnetic state should be identical for residues that are more than ten to 15 residues away from the position of the spin label. Thus, if htau40 would be a fully extended chain most residues will show a PRE intensity ratio of one. In clear contrast, we observed for many residues far from the site of spin-labelling PRE intensity ratios below one, indicative of transient long-range contacts between the spin-label and distant areas of sequence (Figure 7). When the spin label was attached to position 15, intensity ratios of approximately 0.8 were observed in the proline-rich domain P2 (residues 200–240) as well as residues downstream of P2 (residues 130–200) and residues in the repeat region up to residue 330 (Figure 7A). A transient long-range interaction between the N-terminal region and its central domain is further supported by the PRE profile of C239-MTSL labelled htau40, for which intensity ratios of 0.6–0.8 were observed for residues 1–80 and weaker broadening extended up to residue 150 (Figure 7C). On the other hand, signal intensity ratios in case of C291/C322- and C384-MTSL labelled htau40 indicate that the C-terminal domain (residues 360–441) transiently contacts the repeat region and the 40 N-terminal residues (Figure 7D and and77E).
Figure 7
Figure 7
PRE of Amide Protons in Spin-Labelled Tau
Ensemble of Structures Populated by htau40 in Solution
To obtain more direct insight into the ensemble of structures populated by htau40 in solution, we converted all NMR signal intensity ratios <0.9 into distance restraints using the r−6 dependence of the PRE effect on the electron–proton distance [39]. In this way we obtained—from the six PRE profiles shown in Figure 3—1,224 intramolecular long-range contacts ranging between 0 Å and 25 Å. In addition, PRE intensity ratios close to one (here >0.9) indicate that the corresponding amide proton is on average more than 25 Å away from the spin label, allowing lower distance boundaries of 25 Å for these residues. The total of 2,288 PRE distance restraints was subjected to a structure calculation using simulated annealing [40]. Structure calculations of proteins are challenging when the protein exchanges rapidly between different conformations, such that only a single NMR signal is observed. Rapid exchange between multiple conformations is clearly present in the highly dynamic tau protein and the PRE intensity ratios shown in Figure 7 are values averaged over a large ensemble of conformations. Moreover, due to the r−6 dependence, conformations with short intramolecular distances contribute more strongly to the PRE broadening than more extended structures. To take into account the dynamic nature of htau40 we performed both single molecule calculations, in which all distance restraints were enforced simultaneously onto a single molecule, as well as ensemble calculations in which the PRE distance restraints had to be fulfilled not by single structure but collectively by an ensemble of 30 conformations, respectively [41,42]. Clearly, even the 30 conformer calculations are a compromise and underestimate the number of conformations htau40 can assume in solution. Nevertheless they better reflect the ensemble nature of the PRE distances, i.e., the observed PRE broadening arises from a very large number of different conformations and each conformation only fulfils a small subset of distance restraints at any given time. On the other hand, single molecule calculations can allow direct access to the more compact conformations that htau40 could potentially assume in solution.
The structure calculations in which all distance restraints were enforced onto a single molecule resulted in an ensemble of compact conformations (Figure 8A). The shown structures fulfil all 2,288 experimental distance constraints within 0.5 Å. It is readily apparent that many different conformations are in agreement with the experimental PRE distance restraints (structures shown in light grey). The conformation highlighted as ribbon diagram in Figure 8A has a radius of gyration Rg of 48 Å and is therefore at the lower end of the distribution of Rg values obtained from SAXS [26]. Subsequently, this conformation was used in the ensemble calculations. In agreement with the fact that a single compact structure could fulfil all distance restraints, the same was true for structures calculated by ensemble averaging. However, as distance restraints only had to be fulfilled by an ensemble of structures, more expanded conformations were obtained. The average radius of gyration of the ensemble of calculated structures was approximately 65 Å, in agreement with the average value obtained for htau40 from SAXS [26].
Figure 8
Figure 8
Native-State Conformations of 441-Residue Tau in Solution
htau40 is a highly dynamic protein and many conformations are in agreement with the experimental PRE profiles (Figure 8A). To extract long-range interactions that occur in many of the calculated structures, we determined all Cα distances in each structure and averaged the resulting contact map over the ensemble of structures (Figure 8B). Thus, dark spots in the contact map shown in Figure 8B indicate conserved intramolecular interactions. In detail, the following structural properties of htau40 were revealed: (i) The N-terminal 50 residues favour a compact conformation, as indicated by strong contacts within the residue stretch 1–20 and from this region to residues 30–50 (lower left corner in Figure 8B). (ii) The N-terminal 50 residues contact the C-terminal domain, as indicated by the contacts between residues 1–50 and residues 380–400 and seen in the PRE profile of C384-MTSL htau40 (Figure 7E). (iii) Inserts I1 and I2 fold back onto each other, as indicated by the short antidiagonal crossing the diagonal of the contact map at approximately the boundary between I1 and I2. (iv) Residues 113–124 interact with the N-terminal end of I1 (residues E45–D74). (v) The region separating I2 from the proline-rich region has a high propensity for compaction. (vi) Large regions of the N-terminal domain interact with the proline-rich region P2 and repeats R1 to R3. (vii) Residue stretches in the proline-rich domains, which transiently assume polyproline II helical conformations, are in contact. (viii) The proline-rich regions P1 and P2 interact with the hexapeptide in repeat R3. (iv) Repeats R1 and R2 assume compact conformations, favoured by the presence of turns in this region [19]. (v) Repeats R3 and R4 contact the C-terminal domain.
Dependence of Long-Range Contacts on Ionic Strength and Urea
To obtain insight into the nature of the long-range interactions observed in htau40, we performed NMR diffusion measurements in 50 mM phosphate buffer as well as in the presence of 600 mM NaCl and 8 M urea. NMR diffusion experiments allow estimation of the hydrodynamic radius of a protein in solution and therefore allow a global assessment of intramolecular long-range interactions [43]. For htau40 in buffer, the diffusion properties indicate a hydrodynamic radius of 54 Å (Figure 9A). Taking into account that for natively unfolded proteins the radius of gyration (Rg) is about 1.2 to 1.5 times larger than the hydrodynamic radius [44], this is consistent with an Rg value of 65 Å of htau40 as determined by SAXS [26]. In the presence of 600 mM NaCl the hydrodynamic radius of htau40 was increased to 57 Å, and in the presence of 8 M urea further to 64 Å (Figure 9A). In agreement with the increased hydrodynamic radius values, PRE-broadening between the spin label attached to residue 239 and the N-terminal domain was strongly reduced when the ionic strength was raised to 600 mM NaCl (Figure 9B).
Figure 9
Figure 9
Influence of Ionic Strength and Urea on Intramolecular Long-Range Interactions
Interaction of Tau with MTs
The binding of htau40 to MTs was characterized using NMR chemical shift perturbation in 2-D 1H-15N HSQC spectra. As shown previously, taxol-stabilized MTs are stable at 5 °C for several hours, sufficient for the time course of the NMR experiment [16]. Upon addition of taxol stabilized MTs to monomeric htau40 a nonuniform reduction of signal intensities in a 1H-15N HSQC of htau40 was observed (Figure 10A and and10B).10B). The broadening is caused by an exchange of tau molecules between the free and the MT-bound state that is intermediate on the NMR time scale. Strong signal broadening was observed for residues in the proline-rich region P2 and in repeats R1 to R3. For residues 214 to 241 in P2, signal intensities were reduced to below 70%. Within this region two minima were present, comprising residues Leu215 and 225KVAVVRT231 (Figure 10D). In the unbound state, 224KKVAVVR230 preferentially populate β-structure, whereas the two neighbouring residue stretches (216PTPPTREP223 and 232PPKSPSSA239) have a propensity for polyproline II helix (Figure 3G). In P1, 170RIPAKTPPAPKT181 showed a more pronounced broadening than other residues (Figure 10D). Interestingly, part of this stretch preferentially populates polyproline II helical conformations in the free state (Figure 3). In the repeat domain, strong signal broadening was observed for 13 residues in the beginning of repeats R2 and R3 with the minima located at I278 and V309 (Figure 10E). In addition, residues in the homologues region of R1 were strongly attenuated in the presence of MTs, although the signal reduction was more restricted and not as pronounced as in R2 and R3. Significantly less chemical exchange broadening was present in the N-terminal parts of R4 and R' (Figure 10B and and1010E).
Figure 10
Figure 10
Interaction of Tau with MTs
In agreement with the observed chemical exchange broadening, the presence of MTs induced chemical shift changes for residues in P2, R1–R4, and R' (Figure 10C). In addition, 15N chemical shift changes exceeding 0.025 ppm were observed for V75, T76, V80, V88, V122, A125, and I151 (Figure 3). As the N-terminal 150 residues of htau40 do not strongly contribute to binding and assembly of MTs [22,45], the chemical shift changes might be attributed to weak transient contacts with MTs or to changes—as a result of MT binding—in intramolecular long-range interactions in htau40.
Intrinsically Disordered Tau
Tau is important for neuronal cell biology because it stabilizes MTs and promotes axonal outgrowth, and for neurodegeneration because it undergoes abnormal aggregation in AD and other brain disorders [3,5,6]. However, the mode of action of tau is still enigmatic. As soon as the protein was discovered [1], its unusual behaviour became apparent because it retained its function even after heat denaturation. Subsequent biophysical characterization revealed that tau was highly soluble and almost devoid of secondary structure [12]. Cloning of the protein confirmed a high fraction of hydrophilic amino acid residues and an overall basic character, complementary to the acidic surface of MTs [7,8]. It also revealed three or four semiconserved repeats of ~31 residues that were involved both in the interactions with MTs and in the assembly of Alzheimer PHFs. However, the function of MTs was curiously distributed over many residues, each contributing only weakly [9]. Electron microscopy studies of tau showed that the molecule had very little contrast, and only special surface-rendering methods revealed tau as irregular elongated molecules [46,47]. X-ray scattering, circular dichroism (CD), and Fourier transform infrared (FTIR) studies all pointed to a seemingly random structure in solution that was termed “natively unfolded” [13].
NMR spectroscopy provided now a detailed view of the natively unfolded nature of 441-residue tau at single residue resolution. 343 out of 441 residues of the htau40 monomer are in a nonperiodic, disordered conformation (Figure 3). Transient elements of secondary structure were restricted to small regions (Figure 3G): (i) 274Lys-Leu284, 305Ser-Asp315, and 336Gln-Asp345 transiently populate β-structure, in agreement with previous studies on fragments covering the repeat region [17,18]. The high propensity of β-structure in the parts that are essential for PHF formation (274Lys-Leu284 and 305Ser-Asp315) underpins the fact that these residues serve as seeds of aggregation. (ii) 175TPPAPKTPPS184, 216PTPPTREP223, and 232PPKSPSSA239 in the proline-rich regions P1 and P2, transiently assume polyproline II helical conformation. These are the only short residue stretches in htau40 that comprise at least three prolines of which two are sequential (Figure 1B). Within the motifs 175TPPAPKTPPS184, 216PTPPTREP223, and 232PPKSPSSA239 in the proline-rich region, there are several phosphorylation sites that are characteristically elevated in AD, that is, Thr175, Thr181, Thr231, Ser235 [48]. Moreover, several antibodies against phosphorylated tau require dual phosphorylation, separated by three to four residues, e.g., 202 + 205, or 231 + 235. When a polyproline II helical structure is formed the two phosphorylated residues would be facing the same way on the helix. However, in the free state of htau40 the site that is recognized by the AT8 antibody (residues 199SPGSPGT205) does not show a clear propensity to adopt polyproline II helix. (iii) 114Leu-Thr123 transiently populates α-helical structure. The helical structure might promote intramolecular long-range interactions in tau and might be important for interaction with the dynein-activator complex dynactin [49]. Notably, Thr123 is one of only a few residues in the N-terminal domain that is phosphorylated in PHF tau [50]. (iv) 428Leu-Ala437 have a high propensity to form α-helical structure. Notably, truncation of the C terminus behind Asp421 was suggested to be an early molecular event in tau aggregation [51], suggesting that the conformational properties of 428Leu-Ala437 can influence proteolytic cleavage.
Based on their functional differences three different domains of tau were defined: (i) the projection domain comprising residues 1–200, i.e., the N-terminal part of tau up to the proline-rich region P2, (ii) the central region comprising the repeat domain and its flanking regions P2 and R', and (iii) the 40–50 C-terminal residues [52]. NMR dipolar couplings (Figure 4) demonstrated that the functional differences are associated with strong differences in the intrinsic flexibility of the three domains. Whereas the repeat domain and its flanking proline rich regions have a lower intrinsic flexibility on a time scale in the nanosecond to microsecond range detected by dipolar couplings, which is important for formation of secondary structure, in agreement with an increased propensity to populate polyproline II or β-structure, residues in the projection domain as well as in the C-terminal region more rapidly interconvert between different conformations, which is detected by relaxation measurements. The differences in intrinsic mobility are associated with a decreased number of hydrophobic and increased number of negatively charged residues in the projection domain (Figure 1B). Moreover, using NMR dipolar couplings a reduced mobility in regions that harbour many proline residues had been previously observed in the C-terminal tails of α- and β-synuclein [53,54]. Importantly, Figure 4 is also very suggestive of a possible folding of the N-tail and C-tail over the middle domain of tau, previously termed the paperclip model [27].
Intricate Network of Long-Range Interactions in Soluble Tau
The appearance of tau as an elongated molecule by some EM methods [46,47] suggested that the conformation in solution was extended in agreement with the natively unfolded state of tau and the accessibility to multiple kinases throughout the chain. However, evidence for global folding began to emerge from several antibodies that had discontinuous epitopes comprising residues near the N terminus and within the repeat domain (e.g., Alz-50, MC-1 [55]). This evidence was further confirmed by FRET studies showing that tau was able to form a double hairpin, leading to a “paperclip” structure whereby both N- and C terminus were folded into the vicinity of the repeat domain [27]. This concept is now substantially expanded and refined by the NMR analysis (Figure 8). Whereas fluorescence resonance energy transfer combined with electron paramagnetic resonance requires two labels, one label acting as donor and the other one as acceptor, PRE monitored by NMR requires only a single paramagnetic centre such as MTSL attached to a free cysteine. Even more important, whereas in FRET only a single distance can be measured for each donor-acceptor pair, all nuclei in the protein serve as acceptor. Thus, a large number of intraresidual distances (>100) can be measured from a single MTSL-labelled sample. In this study, six uniformly distributed MTSL positions provided a total of 2,288 distance restraints. The distance restraints were integrated into a structural model (Figure 8A), which shows tau in a much more compact form than previously expected from the EM images. However, the molecule is still loosely packed, highly flexible, and exchanges between a large number of conformations, consistent with large average values of the hydrodynamic radius (Figure 9A).
The Cα contact map, which reports on transient interactions that are found in many of the calculated structures, reveals an intricate network of long-range interactions in soluble tau (Figure 8B). In particular the hexapeptide in R3—a residue stretch that is essential for aggregation of tau into PHFs—is strongly involved into intramolecular contacts with both the C-terminal and N-terminal domain of tau. This includes a transient interaction with the amphiphatic C-terminal helix (Figures 7F and and8B).8B). The second residue stretch with increased propensity for formation of an amphiphatic helix, 114Leu-Thr123, interacts with the N-terminal end of I1 as well as with repeats R1 and R2. Striking is also the compaction in the N-terminal 50 residues, in the region between I2 and P1, in P1 and P2, in repeats R1 and R2, and in repeats R3 and R4. In agreement with the paper clip model proposed by FRET measurements [27], the N terminus weakly interacts with the C terminus (Figures 7E, E,7F,7F, F,8B,8B, and S4). Interestingly, the proline rich region P2 has many contacts with distant areas of the sequence, such as R1, R2, R4, and the amphiphatic helix at the C terminus, suggesting that phosphorylation of residues in P2 may modify the ensemble of tau conformations, thereby promoting or delaying aggregation into PHFs. Moreover, the intramolecular interactions between its repeat and proline-rich regions might prime the tau protein for MT binding [45].
Single-Residue Definition of the Hot Spots of the Tau-MT Interaction
The MAP tau is a critical regulator of diverse MT functions [9,45,52]. The repeat domain with its four repeats is essential for MT assembly, however, in the absence of the two flanking regions P2 and R', the repeat domain binds only weakly to MTs. The flanking domains, on the other hand, bind to MTs even in the absence of the repeats. This has led to the proposition of the “jaws” model of tau whereby the regions flanking the repeats are considered as targeting domains, responsible for positioning tau on the MT surface, and the repeats that act as catalytic domains for MT assembly [52,56].
Here we probed the interaction of htau40 with MTs using solution-state NMR spectroscopy. Despite the fact that tau molecules are invisible to solution-state NMR when they are bound to MTs owing to the high molecular weight of the complex, information about the residues of tau that are important for binding to MTs can be obtained when the exchange between the fully bound form and the free state is sufficiently fast. In this case, the observed NMR signals will be an average of the resonances originating from the unbound and bound forms of tau, causing changes in NMR signal position and intensity. The strength of these changes will depend on the conformation and chemical environment in the bound state and the concentration of the bound species. Particularly striking was the pattern of NMR signal intensity ratios in the presence and absence of MTs (Figures 10B and and11A).11A). Four highly localized regions were revealed, in which, because of chemical exchange between the unbound and MT-bound state, HN signal intensities were broadened below 40% of their value in the unbound state: 225KVAVVRT231, 245TAPVPMPDL253, 275VQIINKKLDLSNV287, and 306VQIVYKPVDLSKV318. In these regions, intensity minima were found for V228, M250, I278, and V309, respectively (Figure 11B and and11C).11C). In agreement with the NMR data, deletion analysis mapped the MT-binding activity of the proline-rich region to residues K224-N255 and in particular to the stretch 225KVAVVRT231. Moreover, site-directed mutagenesis indicated that K224, K225, and R230 are important for MT-binding and -assembly [45]. It is noteworthy that the region 225KVAVVRT231 is conserved between tau and two other MAPs, MAP-2 and MAP-4 [45]. The importance of 275VQIINKKLDLSNV287 for MT-binding is supported by biochemical analyses that reported a strong reduction of MT binding affinity upon mutation of K274, K280, K281 to alanine [57]. Importantly, the hot spots of interaction are separated by residue stretches that show smaller chemical shift changes and less signal broadening (Figures 10 and and11).11). These residue stretches might act as flexible linker sequences and suggest that tau protein can assume multiple conformations on the surface of MTs. Moreover, the flexible structure may allow tau to be easily displaced from the MT lattice, consistent with the rapid diffusion of tau in neurons [58].
Figure 11
Figure 11
Comparison of MT-Binding Profile with Hydrophobicity Pattern of Tau
In addition to the hot spots of MT-interaction, signal broadening and chemical shift changes were observed for most residues in R4 and R' and extended to a weaker degree even up to the C terminus, indicative of transient interactions of these regions with MTs. However, in contrast to the proline-rich region P1 and repeats R1, R2, and R3, no clear minimum in NMR signal intensities was observed in repeat R4, indicating that R4 may not be very important to the interaction between tau and MTs. This observation is in agreement with biochemical studies that suggested a core-MT binding domain comprising the N-terminal side of the repeat region [59]. On the other hand, it is in contrast to the view that tau possesses multiple independent tubulin-binding sites [7]. As far as the projection domain is concerned, significant chemical shift changes were observed for residues in insert I2 and in the region with helical propensity (114Leu-Thr123), consistent with the finding that the projection domain regulates the spacing of MTs [60].
Little is known about the nature of the cognate tau binding sites in tubulin. Based on digestion experiments it is believed that tau binds to the acidic carboxyl tail of tubulin, which is supposed to be exposed on the surface of MTs [61]. Moreover, mutational analysis pointed to the importance of positively charged lysine residues (K274, K280, K281) in tau for MT-binding, suggesting that electrostatic interactions are important for the tau-MT binding [45]. On the other hand, the homologues region in R1 contains only a single positively charged residue (K254) at the edge of the most affected residue stretch, but does contain a negatively charged residue (D252). Similarly, the hot spots of MT-binding in R2 and R3 also contain a negative charge (Figures 1B and and11),11), suggesting that the tau-MT interaction might be more complex. Indeed, there is a striking correlation between the NMR-based MT-interaction profile and the hydrophobicity pattern of tau (Figure 11A). The 13-residue stretch in R3 (306VQIVYKPVDLSKV318), which shows the strongest chemical exchange contribution in the presence of MTs, is the most hydrophobic region of 441-residue tau (Figure 11C). Consistent with the importance of hydrophobic interactions, substitution of the tyrosine residue in this residue stretch by an asparagine (Y310->N) reduced the MT affinity of tau [59]. Maxima are also found in the hydrophobicity profile for the other three hot spots of MT-interaction (225KVAVVRT231, 245TAPVPMPDL253, 275VQIINKKLDLSNV287), whereas the homologous region in repeat R4 only has a few hydrophobic residues but has many charged residues, and is less affected by the presence of MTs (Figure 11C). Further support for the importance of hydrophobic interactions for formation of the tau-MT complex comes from the MT-binding site in P2: 225KVAVVRT231 is the most hydrophobic residue stretch in the proline-rich regions P1 and P2 (Figure 11B). Taken together, the NMR and biochemical data suggest a complex mechanism of tau-MT interaction involving both electrostatic and hydrophobic contacts.
Over 30 phosphorylation sites have been identified in tau, many of which are elevated in AD [62,63]. Prominent sites are located in the flanking domains, e.g., S199, S202, T205, T212, S214, T231, S235 before the repeats, S396, S404, S422, and others after the repeats (Figures 11B and and11C).11C). The major sites within the repeats are located in the KXGS motifs, i.e., S262, S293, S324, and S356. These sites are phosphorylated by the kinase MARK, which results in the detachment of tau from MTs [64]. Interestingly, S262, S293, S324, and S356 are not part of the hot spots of MT-interaction, suggesting that phosphorylation at these sites might inhibit MT-binding by long-range electrostatic interactions. Alternatively, or in combination, phosphorylation can induce conformational changes that are incompatible with MT-binding. On the other hand, T231 is right in the middle of the MT-binding region in the proline-rich domain P2, making even a steric inhibition of MT-binding possible. Another potential mechanism is stabilization of a polyproline II helix by phosphorylation. Phosphorylation of residues within the fragments 216PTPPTREP223 and 232PPKSPSSA239 might stabilize their nascent polyproline II helical propensity, such that the resulting conformation is no longer able to efficiently bind to MTs.
Why would a natively unfolded protein evolve to stabilize axonal MTs? To consider this, we note that MTs bind a variety of proteins, some of which are natively unfolded (e.g., the tau-MAP2-MAP4 family), but others are typical well-folded molecules. Two cases in point are kinesin and doublecortin, both of which bind to the MT surface in a periodic and well-defined fashion. By contrast, the MAPs are rather diffusely distributed over the MT surface, there is little detectable periodicity, the MAPs bind relatively weakly, and they diffuse rapidly off MTs and along them [65]. To complicate matters further, binding of MAPs involves the C-terminal tails of tubulin subunits, which are themselves natively unfolded [61]. Not surprisingly, much of the N- and C-terminal domains of tau are highly flexible even when the repeat domain is attached to MTs [22]. Several explanations have been invoked to explain the functions of unfolded proteins such as MAPs [66]: They can act as entropic bristles to keep the spacing between MTs and other cell components (and indeed large MAPs keep larger spacings than small MAPs [60]); serve as assemblers for multisubunit structures (e.g., to pre-assemble tubulin into oligomers for incorporation into MT); serve as docking sites for enzymes (e.g., kinases and phosphatases for the case of MT-bound MAPs); and may even have regulatory functions for MT-related functions (e.g., interaction with motor proteins of axonal transport [67]). Thus, the multiplicity of functions would correspond to a multiplicity of conformations. Many of the above functions are inferred from biochemical evidence without detailed knowledge of the responsible residues and conformations of tau. The identification of residues reported here will provide a basis for future experiments to clarify the interactions of tau with interaction partners in cells, and hopefully the changes that occur during neurodegenerative tauopathies.
Disordered protein statistics.
Assignments of disordered proteins (IDPs) resolved by other groups with the corresponding protein size were found either in the BMRB databank (http://www.bmrb.wisc.edu/) or in publications listed in the PubMed (http://www.ncbi.nlm.nih.gov/pubmed/) databank.
Assembly of MTs.
Porcine brain tubulin was purified and incubated at concentrations higher than 200 M in MT assembly buffer (100 mM Pipes, [pH 6.9], 1 mM EDTA, 1 mM MgSO4, 1 mM dithiothreitol) in the presence of 1 mM GTP at 37 °C for 5 min. After addition of 100 M Paclitaxel (Sigma-Aldrich) the polymerization was performed for 20 min at 37 °C. Analysis of MTs showed that MTs remained stable over the entire duration of the NMR experiments.
Spin labelling of tau.
To label htau40 cysteine-containing mutants with the nitroxide spin label MTSL (Toronto Research Chemicals), DTT was removed before labelling from the buffer by using size exclusion chromatography (PD-10 columns, GE Healthcare), and the proteins were equilibrated in PBS buffer (pH 7.4). Free sulfhydryl groups were reacted with a 5-fold molar excess of the MTSL solubilized in ethyl acetate, at 21 °C for 2.5 h. Unreacted spin label was removed by using PD-10 columns equilibrated in 50 mM Na phosphate buffer (pH 6.8), and spin-labelled proteins were concentrated by using Amicon Ultra-15 (molecular weight cutoff, 3,000) (Millipore).
NMR spectroscopy.
Protein concentrations were between 0.2–0.9 mM of htau40. NMR spectra were acquired on a Bruker Avance 900 spectrometer equipped with a cryogenic probe. Aggregation did not occur under these low temperature conditions. 3-D (HA)CANNH [30] (100 [F1] × 72 [F2] × 1 K [F3] complex data points) and HNN [31] (100 [F1] × 100 [F2] × 1 K [F3] complex data points) experiments (four scans, 1.2-s recovery delay, for each experiment roughly one day of measurement time) were collected. To enable and validate assignment of MT-bound htau40, a 3-D (HA)CANNH experiment was measured at 20 °C (total experiment time: ~1.5 d). NMR data were processed and analyzed using NMRPipe [68] and Sparky 3 (T. D. Goddard and D. G. Kneller, http://www.cgl.ucsf.edu/home/sparky).
Secondary shift values were calculated as the differences between measured Cα chemical shifts and the empirical random coil value for the appropriate amino acid type [69]. Random coil values for histidines, glutamates, and aspartates were taken from Wishart and Sykes [70], as the chemical shifts of these residues are particularly sensitive to pH. To estimate the secondary structure propensity in contiguous segments of htau40, the observed Cα chemical shifts were normalized by the empirically determined secondary shift expected for that residue type in a regular secondary structure (β-sheet or α-helix) conformation [70], summed and normalized by the number of residues in the segment.
3J(HNHα) scalar couplings were measured using an intensity modulated HSQC [71] on a Bruker 900 Avance spectrometer (32 scans, relaxation delay 1.2 ms, 2τ = time for evolution of 3JHNHα: 18 ms). Coupling values were calculated from the intensity ratios using the relation Scross/Sdiag = cos(π3JHNHα2τ). Secondary 3J(HNHα) scalar couplings were calculated as the difference between experimental 3J(HNHα) scalar couplings and random coil values [72].
15N R relaxation rates were measured at 5 °C on a Bruker Avance 700-MHz spectrometer using a spinlock frequency of 2 kHz and relaxation periods of 20, 100, 220, and 300 ms. Relaxation times were calculated by fitting an exponential function to the decaying signal integrals.
One-bond N-H RDCs (DNH) were determined by using an inphase-antiphase (IPAP)-HSQC [73]. DNH values were calculated as the difference between splittings measured in the isotropic phase and in a sample, in which htau40 had been aligned in 5 mg/ml Pf1 bacteriophage (Asla). Errors estimated on the basis of the signal-to-noise ratio are 0.2 Hz for D(NH) and 0.4 Hz for 3J(HNHα) couplings, respectively.
NMR diffusion measurements.
For determining the hydrodynamic radius, htau40 was dissolved in 99.9 % D2O, 50 mM phosphate buffer (pH 6.9). The samples contained dioxane (concentration ~2%) as an internal radius standard and viscosity probe [43]. 1-D 1H spectra were collected employing the standard Bruker pulse program ledbpgp2s. The gradient strength was linearly increased from 2% to 95% of the maximum gradient strength in 16 steps, with 100% gradient strength corresponding to 56.9 G/cm. For each 1H spectrum 128 scans and 16 K complex with a spectral width of 7,200 Hz were acquired. Signal intensities corresponding to the aliphatic region of the 1H spectra (3.3–0.5 ppm) were readout with the TOPSPIN T1/T2 Relaxation module (Bruker Instruments). The diffusion data (signal intensity versus gradient strength) were fitted to exponential functions using Igor Pro 5.01 (WaveMetrics). From the apparent diffusion coefficients of htau40 and dioxane and the known Stokes radius of dioxane (2.12 Å), Stokes radii of monomeric htau40 were calculated [74].
Measurement of PRE.
PRE broadening was investigated using 15N-labelled htau40 at a concentration of 15 μM (MTSL at A15C, A239C, and C291 + C322) and 50 μM (MTSL at A416C, A384C, and A72C) in 50 mM phosphate buffer at (pH 6.8). PRE effects were measured from the peak intensity ratios between two 2D 15N-1H HSQC NMR spectra acquired in the presence of the nitroxide radical and after addition of 4 mM DTT (heated to 45 °C for 15 min before measurement) to the same sample. Addition of DTT will cleave the MTSL tag from the cysteine residue, such that the spin label is no longer attached to the protein and the protein is in the diamagnetic state. Oxidation of the MTSL tag with ascorbic acid, gave very similar results (Figure S3B).
To exclude intermolecular contacts as cause for PRE line-broadening, a mixture of 15 μM 15N-labelled htau40-(C291A/C322A) mutant (without any cysteine residues) and 15 μM 14N-labelled htau40-A239C-MTSL was measured at 5 °C (Figure S3A). In addition, intramolecular PRE broadening at htau40 concentrations of 15 and 50 μM were compared (Figure S3C).
Calculation of distance restraints.
Distance restraints were calculated as described from the intensity ratio between two 2D 15N-1H HSQC NMR spectra, in the diamagnetic and paramagnetic states of the protein [39]. To reduce the impact of peak overlap, for each residue the average of its own intensity ratio Ipara /Idia and that of the preceding and following residue was calculated. These smoothed intensity ratios were linearly fit for the enhancement of the transverse relaxation rate by the unpaired electron [75]. For calculation of distance restraints, amide proton R2 values were approximated by experimental amide nitrogen R values [76]. The correlation time for the electron-nuclear interaction was set to 4 ns, in agreement with previous studies [39]. For peaks broadened beyond detection, distances were set to 7 ± 5 Å. Peaks with intensity ratios below 0.95 were restrained to the calculated distance ±5 Å by using a harmonic square well potential. For residues that were not broadened in the paramagnetic state, a lower distance bound of 25 Å was used. All distances were imposed as restraints between the Cα atom of the residue with the cysteine-MTSL group and residue-specific amide protons.
Structure calculation and analysis.
Structure calculations were performed using XPLOR-NIH, version 2.9.7 [40]. An all-atom representation of htau40 was used. Structural energy terms from steric repulsion, bond length, bond angles, dihedral angles, and favoured regions of the Ramachandran map were employed.
For restraining a single molecule simultaneously by all PRE distance restraints, torsion angle dynamics were started at 3,000 K with the temperature reduced to 20 K, followed by a short energy minimization. 50 structures were calculated starting from a random coil. The seven lowest-energy structures that satisfied the 2,288 distance restraints with no violations greater than 1 Å were used for calculation of the average contact map shown in Figure 8B.
Single molecule calculations were followed by ensemble calculations, in which distance restraints do not have to be fulfilled by a single molecule, but collectively by the ensemble of molecules. Ensemble calculations were started from the lowest-energy structure obtained in the single molecule calculation (see above) and performed in two rounds. Initially, distance restraints were enforced onto an ensemble of 30 molecules [41]. Torsion angle dynamics were used with the temperature reduced from 10,000 K to 5,000 K. The lowest energy structure obtained from this first round of ensemble averaging was subjected to another round of structure calculation using an ensemble size of 5. Torsion angle dynamics was used and the temperature was reduced from 3,000 K to 1,000 K. A total of 100 structures were calculated. The average contact map obtained from the seven lowest-energy structures of the ensemble was very similar to the one obtained from single molecule calculations (Figure S4).
Figure S1: Degeneracy of Cα Connectivity
(A) Residue Y18 points back to T17, but for 18 other residues the Cα frequency also lies within a range of +/− 0.2 ppm (red strip) of the Cα(i-1)peak of Y18.
(B) Comparison of Cα chemical shift degeneracy observed in htau40 (441 residues) with values reported for the globular proteins calmodulin (148 residues), N-terminal domain enzyme I (259 residues), maltose-binding protein (370 residues), and malate synthase G (723 residues).
(411 KB PDF)
Figure S2: Robustness of Residual Secondary Structure of htau40 in Solution
(A),(C),(D), Comparison of Cα secondary chemical shifts, ΔCα, observed in htau40 and the tau fragments K32, K25, and K10, respectively (all at 5 °C).
(B) Comparison of Cα secondary chemical shifts, ΔCα, observed in htau40 at 278K and 293K. Red lines indicate y = x.
(337 KB PDF)
Figure S3: Control PRE Measurements
(A) Intensity ratio (Iparam/Idiam) obtained from HSQC spectra recorded on 15 μM 15N-labelled htau40-(C291A/C322A) (i.e., htau40 without any cysteine residues) in the presence of 15 μM 14N-labelled htau40-A239C-MTSL in the paramagnetic and diamagnetic state. No significant deviations from unity were observed, indicating that decreases in intensity ratios are not due to aggregation, and can be attributed instead to intramolecular long-range contacts in native htau40.
(B) Comparison of PRE intensity ratios for different diamagnetic reference states. For the data shown by blue bars, DTT was added to the sample, whereas for the red line ascorbic acid was used to quench the spin label.
(C) Comparison of intensity ratios obtained from HSQC spectra recorded on 15 μM (red line) and 50 μM (blue bars) 15N-labelled htau40-A239C-MTSL in the presence and absence of DTT.
(202 KB PDF)
Figure S4: Average Contact Map for the Seven Lowest-Energy Structures Obtained from Ensemble Calculations with an Ensemble Size of 5 (see Materials and Methods)
A continuous grey scale from 3 Å (black) to 22 Å (white) is used.
(2.14 MB PDF)
Acknowledgments
The authors wish to thank Lukasz Skora and Ilka Lindner for excellent technical assistance and Eva-Maria Mandelkow for discussions.
Abbreviations
ADAlzheimer disease
FRETFörster resonance energy transfer
HSQCheteronuclear single quantum coherence
htau40longest tau isoform found in the human central nervous system
MAPmicrotubule-associated protein
MTmicrotubule
MTSL(1-oxy-2,2,5,5-tetramethyl-D-pyrroline-3-methyl)-methanethiosulfonate
NMRnuclear magnetic resonance
PHFpaired helical filament
PREparamagnetic relaxation enhancement
RDCresidual dipolar coupling
SAXSsmall angle x-ray scattering

Footnotes
Author contributions. EM and MZ conceived the project. MDM, SB, JK, SJ, and JB performed the experiments. MDM, CG, EM, and MZ wrote the paper.
Funding. This work was supported by the Max Planck Society, the Fonds der Chemischen Industrie, the Boehringer Ingelheim Fonds, and the DFG through GRK 782 and through a DFG Heisenberg scholarship (ZW 71/2–1 and 3–1). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests. The authors have declared that no competing interests exist.
  • Weingarten MD, Lockwood AH, Hwo SY, Kirschner MW. A protein factor essential for microtubule assembly. Proc Natl Acad Sci U S A. 1975;72:1858–1862. [PubMed]
  • Drubin DG, Kirschner MW. Tau protein function in living cells. J Cell Biol. 1986;103:2739–2746. [PMC free article] [PubMed]
  • Garcia ML, Cleveland DW. Going new places using an old MAP: tau, microtubules and human neurodegenerative disease. Curr Opin Cell Biol. 2001;13:41–48. [PubMed]
  • Cassimeris L, Spittle C. Regulation of microtubule-associated proteins. Int Rev Cytol. 2001;210:163–226. [PubMed]
  • Goedert M, Spillantini MG. A century of Alzheimer's disease. Science. 2006;314:777–781. [PubMed]
  • Ballatore C, Lee VM, Trojanowski JQ. Tau-mediated neurodegeneration in Alzheimer's disease and related disorders. Nat Rev Neurosci. 2007;8:663–672. [PubMed]
  • Lee G, Cowan N, Kirschner M. The primary structure and heterogeneity of tau protein from mouse brain. Science. 1988;239:285–288. [PubMed]
  • Goedert M, Spillantini MG, Potier MC, Ulrich J, Crowther RA. Cloning and sequencing of the cDNA encoding an isoform of microtubule-associated protein tau containing four tandem repeats: differential expression of tau protein mRNAs in human brain. Embo J. 1989;8:393–399. [PubMed]
  • Butner KA, Kirschner MW. Tau protein binds to microtubules through a flexible array of distributed weak sites. J Cell Biol. 1991;115:717–730. [PMC free article] [PubMed]
  • Wischik CM, Novak M, Edwards PC, Klug A, Tichelaar W, et al. Structural characterization of the core of the paired helical filament of Alzheimer-disease. Proc Natl Acad Sci U S A. 1988;85:4884–4888. [PubMed]
  • von Bergen M, Friedhoff P, Biernat J, Heberle J, Mandelkow EM, et al. Assembly of tau protein into Alzheimer paired helical filaments depends on a local sequence motif ((306)VQIVYK(311)) forming beta structure. Proc Natl Acad Sci U S A. 2000;97:5129–5134. [PubMed]
  • Cleveland DW, Hwo SY, Kirschner MW. Physical and chemical properties of purified tau factor and the role of tau in microtubule assembly. J Mol Biol. 1977;116:227–247. [PubMed]
  • Schweers O, Schonbrunnhanebeck E, Marx A, Mandelkow E. Structural studies of tau-protein and Alzheimer paired helical filaments show no evidence for beta-structure. J Biol Chem. 1994;269:24290–24297. [PubMed]
  • Dunker AK, Silman I, Uversky VN, Sussman JL. Function and structure of inherently disordered proteins. Curr Opin Struct Biol. 2008;18:756–764. [PubMed]
  • Dyson HJ, Wright PE. Intrinsically unstructured proteins and their functions. Nat Rev Mol Cell Biol. 2005;6:197–208. [PubMed]
  • Mukrasch MD, von Bergen M, Biernat J, Fischer D, Griesinger C, et al. The “jaws” of the tau-microtubule interaction. J Biol Chem. 2007;282:12230–12239. [PubMed]
  • Mukrasch MD, Biernat J, von Bergen M, Griesinger C, Mandelkow E, et al. Sites of tau important for aggregation populate {beta}-structure and bind to microtubules and polyanions. J Biol Chem. 2005;280:24978–24986. [PubMed]
  • Eliezer D, Barre P, Kobaslija M, Chan D, Li X, et al. Residual structure in the repeat domain of tau: echoes of microtubule binding and paired helical filament formation. Biochemistry. 2005;44:1026–1036. [PubMed]
  • Mukrasch MD, Markwick P, Biernat J, Bergen M, Bernado P, et al. Highly populated turn conformations in natively unfolded tau protein identified from residual dipolar couplings and molecular simulation. J Am Chem Soc. 2007;129:5235–5243. [PubMed]
  • Barre P, Eliezer D. Folding of the repeat domain of tau upon binding to lipid surfaces. J Mol Biol. 2006;362:312–326. [PubMed]
  • Andronesi OC, von Bergen M, Biernat J, Seidel K, Griesinger C, et al. Characterization of Alzheimer's-like paired helical filaments from the core domain of tau protein using solid-state NMR spectroscopy. J Am Chem Soc. 2008;130:5922–5928. [PubMed]
  • Sillen A, Barbier P, Landrieu I, Lefebvre S, Wieruszeski JM, et al. NMR investigation of the interaction between the neuronal protein Tau and the microtubules. Biochemistry. 2007;46:3055–3064. [PubMed]
  • Landrieu I, Lacosse L, Leroy A, Wieruszeski JM, Trivelli X, et al. NMR analysis of a Tau phosphorylation pattern. J Am Chem Soc. 2006;128:3575–3583. [PubMed]
  • Sillen A, Leroy A, Wieruszeski JM, Loyens A, Beauvillain JC, et al. Regions of tau implicated in the paired helical fragment core as defined by NMR. Chembiochem. 2005;6:1849–1856. [PubMed]
  • Sibille N, Sillen A, Leroy A, Wieruszeski JM, Mulloy B, et al. Structural impact of heparin binding to full-length Tau as studied by NMR spectroscopy. Biochemistry. 2006;45:12560–12572. [PubMed]
  • Mylonas E, Hascher A, Bernado P, Blackledge M, Mandelkow E, et al. Domain conformation of tau protein studied by solution small-angle X-ray scattering. Biochemistry. 2008;47:10345–10353. [PubMed]
  • Jeganathan S, von Bergen M, Brutlach H, Steinhoff HJ, Mandelkow E. Global hairpin folding of tau in solution. Biochemistry. 2006;45:2283–2293. [PubMed]
  • Sevcik J, Skrabana R, Dvorsky R, Csokova N, Iqbal K, et al. X-ray structure of the PHF core C-terminus: insight into the folding of the intrinsically disordered protein tau in Alzheimer's disease. FEBS Lett. 2007;581:5872–5878. [PubMed]
  • Hiller S, Wasmer C, Wider G, Wuthrich K. Sequence-specific resonance assignment of soluble nonglobular proteins by 7D APSY-NMR spectroscopy. J Am Chem Soc. 2007;129:10823–10828. [PubMed]
  • Zweckstetter M, Bax A. Single-step determination of protein substructures using dipolar couplings: aid to structural genomics. J Am Chem Soc. 2001;123:9490–9491. [PubMed]
  • Panchal SC, Bhavesh NS, Hosur RV. Improved 3D triple resonance experiments, HNN and HN(C)N, for H-N and N-15 sequential correlations in (C-13, N-15) labeled proteins: application to unfolded proteins. J Biomol NMR. 2001;20:135–147. [PubMed]
  • Dyson HJ, Wright PE. Nuclear magnetic resonance methods for elucidation of structure and dynamics in disordered states. Methods Enzymol. 2001;339:258–270. [PubMed]
  • Tjandra N, Bax A. Direct measurement of distances and angles in biomolecules by NMR in a dilute liquid crystalline medium. Science. 1997;278:1111–1114. [PubMed]
  • Prestegard JH, al-Hashimi HM, Tolman JR. NMR structures of biomolecules using field oriented media and residual dipolar couplings. Q Rev Biophys. 2000;33:371–424. [PubMed]
  • Blackledge M. Recent progress in the study of biomolecular structure and dynamics in solution from residual dipolar couplings. Prog NMR Spectrosc. 2005;46:23–61.
  • Mohana-Borges R, Goto NK, Kroon GJ, Dyson HJ, Wright PE. Structural characterization of unfolded states of apomyoglobin using residual dipolar couplings. J Mol Biol. 2004;340:1131–1142. [PubMed]
  • Wang CY, Palmer AG. Solution NMR methods for quantitative identification of chemical exchange in N-15-labeled proteins. Magn Reson Chem. 2003;41:866–876.
  • Lange OF, Lakomek NA, Fares C, Schroder GF, Walter KF, et al. Recognition dynamics up to microseconds revealed from an RDC-derived ubiquitin ensemble in solution. Science. 2008;320:1471–1475. [PubMed]
  • Gillespie JR, Shortle D. Characterization of long-range structure in the denatured state of staphylococcal nuclease. I. Paramagnetic relaxation enhancement by nitroxide spin labels. J Mol Biol. 1997;268:158–169. [PubMed]
  • Schwieters CD, Kuszewski JJ, Tjandra N, Clore GM. The Xplor-NIH NMR molecular structure determination package. J Magn Reson. 2003;160:65–73. [PubMed]
  • Iwahara J, Schwieters CD, Clore GM. Ensemble approach for NMR structure refinement against (1)H paramagnetic relaxation enhancement data arising from a flexible paramagnetic group attached to a macromolecule. J Am Chem Soc. 2004;126:5879–5896. [PubMed]
  • Lindorff-Larsen K, Kristjansdottir S, Teilum K, Fieber W, Dobson CM, et al. Determination of an ensemble of structures representing the denatured state of the bovine acyl-coenzyme a binding protein. J Am Chem Soc. 2004;126:3291–3299. [PubMed]
  • Wilkins DK, Grimshaw SB, Receveur V, Dobson CM, Jones JA, et al. Hydrodynamic radii of native and denatured proteins measured by pulse field gradient NMR techniques. Biochemistry. 1999;38:16424–16431. [PubMed]
  • Sherman E, Haran G. Coil-globule transition in the denatured state of a small protein. Proc Natl Acad Sci U S A. 2006;103:11539–11543. [PubMed]
  • Goode BL, Denis PE, Panda D, Radeke MJ, Miller HP, et al. Functional interactions between the proline-rich and repeat regions of tau enhance microtubule binding and assembly. Mol Biol Cell. 1997;8:353–365. [PMC free article] [PubMed]
  • Hirokawa N, Shiomura Y, Okabe S. Tau proteins: the molecular structure and mode of binding on microtubules. J Cell Biol. 1988;107:1449–1459. [PMC free article] [PubMed]
  • Wille H, Drewes G, Biernat J, Mandelkow EM, Mandelkow E. Alzheimer-like paired helical filaments and antiparallel dimers formed from microtubule-associated protein tau in vitro. J Cell Biol. 1992;118:573–584. [PMC free article] [PubMed]
  • Mi K, Johnson GV. The role of tau phosphorylation in the pathogenesis of Alzheimer's disease. Curr Alzheimer Res. 2006;3:449–463. [PubMed]
  • Magnani E, Fan J, Gasparini L, Golding M, Williams M, et al. Interaction of tau protein with the dynactin complex. Embo J. 2007;26:4546–4554. [PubMed]
  • Gong CX, Liu F, Grundke-Iqbal I, Iqbal K. Post-translational modifications of tau protein in Alzheimer's disease. J Neural Transm. 2005;112:813–838. [PubMed]
  • Gamblin TC, Chen F, Zambrano A, Abraha A, Lagalwar S, et al. Caspase cleavage of tau: linking amyloid and neurofibrillary tangles in Alzheimer's disease. Proc Natl Acad Sci U S A. 2003;100:10032–10037. [PubMed]
  • Gustke N, Trinczek B, Biernat J, Mandelkow EM, Mandelkow E. Domains of tau protein and interactions with microtubules. Biochemistry. 1994;33:9511–9522. [PubMed]
  • Bertoncini CW, Jung YS, Fernandez CO, Hoyer W, Griesinger C, et al. Release of long-range tertiary interactions potentiates aggregation of natively unstructured alpha-synuclein. Proc Natl Acad Sci U S A. 2005;102:1430–1435. [PubMed]
  • Sung YH, Eliezer D. Residual structure, backbone dynamics, and interactions within the synuclein family. J Mol Biol. 2007;372:689–707. [PMC free article] [PubMed]
  • Jicha GA, Bowser R, Kazam IG, Davies P. Alz-50 and MC-1, a new monoclonal antibody raised to paired helical filaments, recognize conformational epitopes on recombinant tau. J Neurosci Res. 1997;48:128–132. [PubMed]
  • Preuss U, Biernat J, Mandelkow EM, Mandelkow E. The ‘jaws' model of tau-microtubule interaction examined in CHO cells. J Cell Sci. 1997;110(Pt 6):789–800. [PubMed]
  • Panda D, Goode BL, Feinstein SC, Wilson L. Kinetic stabilization of microtubule dynamics at steady state by tau and microtubule-binding domains of tau. Biochemistry. 1995;34:11117–11127. [PubMed]
  • Konzack S, Thies E, Marx A, Mandelkow EM, Mandelkow E. Swimming against the tide: mobility of the microtubule-associated protein tau in neurons. J Neurosci. 2007;27:9916–9927. [PubMed]
  • Goode BL, Chau M, Denis PE, Feinstein SC. Structural and functional differences between 3-repeat and 4-repeat tau isoforms. Implications for normal tau function and the onset of neurodegenetative disease. J Biol Chem. 2000;275:38182–38189. [PubMed]
  • Chen J, Kanai Y, Cowan NJ, Hirokawa N. Projection domains of MAP2 and tau determine spacings between microtubules in dendrites and axons. Nature. 1992;360:674–677. [PubMed]
  • Nogales E, Whittaker M, Milligan RA, Downing KH. High-resolution model of the microtubule. Cell. 1999;96:79–88. [PubMed]
  • Morishima-Kawashima M, Hasegawa M, Takio K, Suzuki M, Yoshida H, et al. Proline-directed and non-proline-directed phosphorylation of PHF-tau. J Biol Chem. 1995;270:823–829. [PubMed]
  • Hanger DP, Betts JC, Loviny TL, Blackstock WP, Anderton BH. New phosphorylation sites identified in hyperphosphorylated tau (paired helical filament-tau) from Alzheimer's disease brain using nanoelectrospray mass spectrometry. J Neurochem. 1998;71:2465–2476. [PubMed]
  • Drewes G, Ebneth A, Preuss U, Mandelkow EM, Mandelkow E. MARK, a novel family of protein kinases that phosphorylate microtubule-associated proteins and trigger microtubule disruption. Cell. 1997;89:297–308. [PubMed]
  • Santarella RA, Skiniotis G, Goldie KN, Tittmann P, Gross H, et al. Surface-decoration of microtubules by human tau. J Mol Biol. 2004;339:539–553. [PubMed]
  • Tompa P. Intrinsically unstructured proteins. Trends Biochem Sci. 2002;27:527–533. [PubMed]
  • Stamer K, Vogel R, Thies E, Mandelkow E, Mandelkow EM. Tau blocks traffic of organelles, neurofilaments, and APP vesicles in neurons and enhances oxidative stress. J Cell Biol. 2002;156:1051–1063. [PMC free article] [PubMed]
  • Delaglio F, Grzesiek S, Vuister GW, Zhu G, Pfeifer J, et al. Nmrpipe - a multidimensional spectral processing system based on unix pipes. J Biomol NMR. 1995;6:277–293. [PubMed]
  • Schwarzinger S, Kroon GJA, Foss TR, Chung J, Wright PE, et al. Sequence-dependent correction of random coil NMR chemical shifts. J Am Chem Soc. 2001;123:2970–2978. [PubMed]
  • Wishart DS, Sykes BD. The C-13 chemical-shift index - a simple method for the identification of protein secondary structure using C-13 chemical-shift data. J Biomol NMR. 1994;4:171–180. [PubMed]
  • Permi P, Kilpelainen I, Annila A, Heikkinen S. Intensity modulated HSQC and HMQC: two simple methods to measure 3J(HNH)alpha in proteins. J Biomol NMR. 2000;16:29–37. [PubMed]
  • Plaxco KW, Morton CJ, Grimshaw SB, Jones JA, Pitkeathly M, et al. The effects of guanidine hydrochloride on the ‘random coil' conformations and NMR chemical shifts of the peptide series GGXGG. J Biomol NMR. 1997;10:221–230. [PubMed]
  • Ottiger M, Delaglio F, Bax A. Measurement of J and dipolar couplings from simplified two-dimensional NMR spectra. J Magn Reson. 1998;131:373–378. [PubMed]
  • Jones JA, Wilkins DK, Smith LJ, Dobson CM. Characterisation of protein unfolding by NMR diffusion measurements. J Biomol NMR. 1997;10:199–203.
  • Battiste JL, Wagner G. Utilization of site-directed spin labeling and high-resolution heteronuclear nuclear magnetic resonance for global fold determination of large proteins with limited nuclear overhauser effect data. Biochemistry. 2000;39:5355–5365. [PubMed]
  • Ishima R, Torchia DA. Extending the range of amide proton relaxation dispersion experiments in proteins using a constant-time relaxation-compensated CPMG approach. J Biomol NMR. 2003;25:243–248. [PubMed]
Articles from PLoS Biology are provided here courtesy of
Public Library of Science