Next Article in Journal
Conformational Dynamics of Herpesviral NEC Proteins in Different Oligomerization States
Previous Article in Journal
Ubiquitylation at the Fork: Making and Breaking Chains to Complete DNA Replication
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Filamentous Aggregates of Tau Proteins Fulfil Standard Amyloid Criteria Provided by the Fuzzy Oil Drop (FOD) Model

1
ABB Business Services Sp. z o.o. ul. Żegańska 1, 04-713 Warszawa, Poland
2
ACK—Cyfronet AGH, Nawojki 11, 30-950 Kraków, Poland
3
Department of Bioinformatics and Telemedicine, Medical College, Jagiellonian University, Łazarza 16, 31-530 Kraków, Poland
4
Faculty of Physics, Astronomy and Applied Computer Science, Jagiellonian University, Łojasiewicza 11, 30-348 Kraków, Poland
5
Chair of Medical Biochemistry, Medical College, Jagiellonian University, Kopernika 7, 31-034 Kraków, Poland
*
Author to whom correspondence should be addressed.
Currently address: Schibsted Tech Polska Sp. z o. o. ul. Armii Krajowej 28, 30-150 Kraków, Poland.
Int. J. Mol. Sci. 2018, 19(10), 2910; https://doi.org/10.3390/ijms19102910
Submission received: 14 July 2018 / Revised: 12 September 2018 / Accepted: 20 September 2018 / Published: 25 September 2018
(This article belongs to the Section Molecular Biophysics)

Abstract

:
Abnormal filamentous aggregates that are formed by tangled tau protein turn out to be classic amyloid fibrils, meeting all the criteria defined under the fuzzy oil drop model in the context of amyloid characterization. The model recognizes amyloids as linear structures where local hydrophobicity minima and maxima propagate in an alternating manner along the fibril’s long axis. This distribution of hydrophobicity differs greatly from the classic monocentric hydrophobic core observed in globular proteins. Rather than becoming a globule, the amyloid instead forms a ribbonlike (or cylindrical) structure.

1. Introduction

The origin of amyloid transformation has attracted scientific attention for more than 35 years—at least since being acknowledged as the cause of various neurodegenerative disorders [1]. The coexistence and mutual relations between Aβ amyloids and tau tangles, resulting in the damage and destruction of synapses, is believed to provoke behavioral changes that are associated with cognitive impairment [2,3]. The appearance of amyloid fibrils is the consequence of plasticity of proteins, which can adopt different conformational states [4]. The proteins of high content of intrinsically disordered structural forms seem to be the candidates ready for partially folded state which may transform to disordered aggregates with low packing [5,6].
Emergence of fibrillary structures is also thought of as the result of involvement of intrinsically disordered proteins, especially at early phases of the folding process [7].
Reaching the form of highly packed structuralised aggregates that are based mainly on β-structural forms opens the possibility for the unlimited elongation of highly packed ordered amyloid form [8]. The presence of Beta-structural form (cross Beta) allows for the propagation due to the possible H-bonds system to be organised on both sites.
The specificity of tau amyloidosis, as evidenced by abnormal phosphorylation, results not only in disorganization of microtubulin, but also in the appearance of intracellular tau filaments, referred to as neurofibrillary tangles [9]. Tauopathies are defined as clinically, morphologically, and biochemically heterogeneous neurodegenerative diseases that are characterized by the deposition of abnormal tau protein in the brain [10].
Tau amyloid fibrils are regarded as peculiar due to the existence of two distinct superfibrillary forms: straight filament and paired helical filament [11]. Two individual protofilaments may form different structures depending on their mutual arrangement in the dimer. The authors of [11] refer to the conformation of the protofilament as “C-shaped”. The dimer (superfibril), when analyzed under cryo-electron microscopy (cryo-EM) imaging [12], resembles two arched C-shaped structures that are bound back to back. In one form, the structure is symmetrical (with the same residues in both unit molecules involved in complexation), while the other form lacks such symmetry (in this case, the central fragment of one chain makes contact with a fragment which is somewhat closer to the N terminus of the adjacent chain).
According to the analysis shown in [11], the dimerization of protofilaments occurs by way of hydrogen bonds forming between adjacent fragments. In its native form, tau in complex with a microtubule adopts a conformation referred to as “natively unfolded”. This conformation is highly resistant to spikes in temperature and acidity [13].
The presented here analysis focuses on microtubule-associated tau neurofibrillary tangle protein, paired helical filaments (consisting of two individual fibrils), individual protofibrils, as well as individual chains.
Identification and characterization of amyloid structures, as shown in this work, bases on comparative analysis of the structure of globular and fibrillary proteins. The proposed fibrillarization model for globular structures is also discussed in [14,15]. At the core of the fuzzy oil drop (FOD) model lies the concept of applying a three-dimensional (3D) Gaussian form to express the idealized distribution of hydrophobicity in a globular protein. Such distribution has a distinct peak at the geometric center of the globule and then falls off along with distance from the center, reaching almost zero at the surface (i.e., at a distance of 3σ from the center). We have identified proteins whose structure closely corresponds to this theoretical distribution [16]. Any local deviations are usually associated with the protein’s intended function: local excess of hydrophobicity, if occurring on the surface, usually marks a complexation interface, while the hydrophobicity deficiencies tend to correspond to binding cavities, which are capable of housing ligands (or substrates, in the case of enzymes) [17,18]. The universality and ubiquity of hydrophobic cores is attested to by a study of a large number of proteins with varying secondary and supersecondary structural characteristics [19]. A specific type of discordance vs. the idealized (monocentric) distribution of hydrophobicity is observed in the case of proteins which contain so-called solenoid fragments, including some antifreeze proteins [20]. Such fragments deviate from a centralized core in favor of a distribution comprising alternating bands of high and low hydrophobicity, propagating along solenoids long axis [20]. However, in addition to solenoids, biologically active proteins also include fragments whose purpose is to restrict unchecked propagation of such linear sequences (thereby preventing dimerization or polymerization), as well as to ensure solubility. The special “caps” are identified in such proteins [20]. Structures that strongly deviate from the monocentric distribution and lack suitable “caps” are prone to amyloid filaments [14,15]. In the context of the FOD model, such structures can be likened to ribbonlike micelles, capable of unrestricted propagation [21].
A thorough description of the FOD model can be found in [22].
In the FOD model, the emergent structure of the hydrophobic core is though to result from interactions between the protein and the aqueous solvent, which directs hydrophobic residues towards the center of the globule and favors exposure of hydrophilic residues on the surface. As already indicated, local discordances are often associated with the presence of external factors. The fact that amyloid forms do not require mutations to emerge suggests that misfolding is not caused by factors internal to the polypeptide itself. One shall mention the mutation-related amyloidosis [23], however the prion-based amyloid transformation does not require the presence of mutation, as it is the discussed case.
Shaking is known to promote amyloidogenesis—and it can hardly be called a chemical factor. Many factors—including chemical ones—were identified to support the amyloidosis transformation [24]. However other factors than environmental (shaking in particular) are not the object of our analysis. One shall also take into account that the folding as well the misfolding processes take place in macromolecular crowding conditions, however the immanent presence of water makes the water environment of high importance [25].
Perhaps shaking disrupts the structure of the solvent in such a way as to prevent it from guiding “natural” conformational changes within the protein chain. Alternatively, shaking is notable for aerating the solvent. The resulting increase in the area of the liquid/gas interface may produce structural changes within the solvent itself.
In addition to analysis of the tau amyloid, as listed in Protein Data Bank [26], this work proposes an in silico experiment, which involves determining alternative structures that the tau amyloid sequence may attain (using specialized protein folding software, such as Robetta [27,28] and I-Tasser [29,30]), and performing folding simulations based on the FOD model. It turns out that the sequence is indeed capable of producing a globular form with a single, monocentric hydrophobic core. Subjecting globular structures to FOD characterization enables us to track changes that result in amyloidogenesis. The work focuses on three distinct structures: (1) the superfibril (seeking the causes behind its structural variability); (2) the protofibril (identifying the characteristic properties of amyloid structures); and (3) a single chain participating in the fibril. Our research is based on observations rooted in the FOD model, specifically, the linear propagation of hydrophobicity in amyloids (which prevents a shared hydrophobic core from forming). As discussed in [14,15], the presence of alternating bands of high and low hydrophobicity can be regarded as one of the principal indicators of amyloid transformation.

2. Results

2.1. Abbreviations Used

  • FOD—Fuzzy Oil Drop model
  • RD—Relative Distance—The divergence entropy introduced by Kullback and Leibler (described in Methods) used to express the distance between two compared profiles is of entropy category thus it requires an introduction of reference distribution. This is why the distance between T-O (T theoretical—idealized distribution and O-observed distribution) measured by divergence entropy is compared with the O-R (O-observed, R-uniform distribution deprived on any form of hydrophobicity concentration), also measured by divergence entropy. The parameter expressing the relative distance O|T versus (O|T + O|R) measures the closeness of O distribution versus T distribution in respect to O versus R distribution. The RD parameter becomes polypeptide chain length independent. It makes possible comparison of different proteins.
  • RD (T-O-R)—RD parameter calculated for two reference distributions T-theoretical and R-uniform
  • RD (T-O-H)—RD parameter calculated in respect to reference distribution called H-distribution based on intrinsic hydrophobicity of amino acids present in particular polypeptide chain fragment
  • HvT—correlation coefficient expressing the relation between H-intrinsic hydrophobicity of amino acids versus the T-theoretical (expected) hydrophobicity for the idealized status of the residue
  • HvO—correlation coefficient expressing relation between H-intrinsic hydrophobicity of amino acid versus is status as observed in particular protein
  • TvO—correlation coefficient expressing relation between T-idealized hydrophobicity and O-observed in protein under consideration
  • phf-tau—paired helical filament-tau
  • phf-tauO—paired helical filament-tau—as it is available in 5O3O phf-tau in symmetrical form of superfibril
  • phf-tauT—paired helical filament-tau—as it is available in 5O3T phf-tau in asymmetrical form of superfibril
  • phf-tauL—paired helical filament-tau—as it is available in 5O3L phf-tau in form of superfibril similar to phf-tauO
  • IT-#—identification of the model constructed using I-Tasser program with number 1–5 since 5 models were constructed using this program
  • ROB-#—identification of the model constructed using Robetta program with number 1–5 since 5 models were constructed using this program
  • FOD-#—identification of the model constructed using FOD model with number 1–5 since 5 models were constructed using this program
  • Tau (267–312)—fragment of tau peptide—protein under PDB ID 2MZ7
  • Tpp—tau phosphothreonine peptide—protein under PDB ID 1I8H
  • Tau (306–311A)—fragment of tau to identify the structure available in PDB as 2ON9
  • Tau (306–310)—fragment of tau to identify the structure available in PDB as 3Q9G
  • Tau (306–311B)—fragment of tau to identify the structure available in PDB as 3OVL
  • Tau (305–311)—fragment of tau to identify the structure available in PDB as 4E0M
  • Tau (623–628)—fragment of tau to identify the structure available in PDB as 4NP8
  • Tau (306–311C)—fragment of tau to identify the structure available in PDB as 5K7N
  • F-actin—actin, alpha skeletal muscle as available in PDB as 3J8I
  • PDB—Protein Data Bank
  • CASP—Critical Assessment Protein Structure Prediction
  • BLAST—Basic Local Alignment Search Tool
  • PSI-BLAST—Position-Specific Iterated BLAST
  • MSA—Multiple Sequence Alignment

2.2. Superfibril

This analysis concerns the amyloid form that is listed in PDB as 5O3O, 5O3L, and 5O3T (pronase-treated paired helical filament in Alzheimer’s disease brain neurofibrillary tangle protein, paired helical filament-tau, phf-tau, Homo sapiens). Fragment: residues 623–695 of tau protein (306–378 according to PDB numbering) Chains A, C, E, G, I, along with their counterparts (B, D, F, H, J) make up the proto-fibrils [11]. In order to characterize individual chains in the context of the superfibril and proto-fibrils, we have singled out chains E and F. These two chains are located in the central part of the fibril and can be regarded as representative of an arbitrarily long structure. This selection also minimizes edge effects caused by the finite width of the complex.

Properties of Superfibrils and Interfaces—What Is the Source of Different Isoforms of Tau Filaments?

Table 1 presents the status of tau amyloid structures in terms of RD values, revealing large discordances between T and O profiles in both models (T-O-R and T-O-H). This means that the distribution of hydrophobicity does not involve a central hydrophobic core. Further analysis will reveal that the amyloids are dominated by a pattern that consists of alternating bands of high and low hydrophobicity. High values of RD further indicate that the folding process is driven by the intrinsic properties of each residue rather than by a global force field—this is also typical for amyloids [14,15]. Regarding the hydrophobicity profile correlation coefficients, HvT and TvO lag behind HvO. This is also due to the absence of a central hydrophobic core, which is replaced by linear propagation of narrow “bands”. In further sections we will specifically describe locations that exhibit these properties.
Visual comparison of T and O (Figure 1) highlights the major differences between these distributions. It should be noted that the chart consists of many overlapping profiles, which means that the distribution of local minima and maxima is replicated in each adjacent chain, resulting in a set of narrow bands, as suggested above.
The FOD model may also be used to predict the properties of shared hydrophobic cores in protein complexes [31]. In order to properly characterize a given complex, it is important to assess the status of its interface. With regard to proto-fibrils, the distribution that was observed in phf-tauT differs from those exhibited by the remaining structures. However the difference is limited only to the structure of interface, which is discussed in this paper. In phf-tauO and phf-tauL the status is similar and it suggests that the superfibril emerges as a result of factors consistent with the FOD model, i.e., under the influence of the aqueous solvent. This interpretation is supported by the high values of all correlation coefficients. We may conclude that the interface is shaped by all factors which determine the structure of the complex itself, with major involvement of water.
The picture changes, however, when dealing with phf-tauT. Its high value of RD (T-O-R), coupled with negative values of HvT and TvO coefficients and a high value of the HvO coefficient, suggest that, in this case, the solvent does not play a significant role in complexation.
It should be noted that the status of the interface is computed by taking into account all interface residues in the entire fibril (following protein-protein contact distance criteria of PDBsum [32]). When all three correlation coefficients adopt strongly positive values, we may assume that the structure of the interface represents a compromise between all three hydrophobicity profiles (observed, intrinsic and theoretical). In contrast, negative values of HvT and TvO are understood to mean that the interface folds “in spite of” the FOD model and in consequence in spite of environmental effects that act upon the protofibril complex.
The characterization concerns chains E and F, which are located centrally and therefore representative of an unrestricted fibril. The calculated values are typical for amyloid forms, and include high values of RD and HvO, along with very low (sometimes even negative) values of HvT and TvO. In phf-tauT, both individual chains, as well as the interface fragment, are shaped by intrinsic hydrophobicity rather than by the external environment, which would favor the formation of a monocentric hydrophobic core.
As the status of individual chains (in the context of the superfibril) is largely similar in all structures, we will limit their presentation to Phf-tau in form, as observed in phf-tauO and phf-tauT (Table 2).
When analyzing individual chains as components of the superfibril, we arrive at similar RD values. In contrast, when the same chains are analyzed as components of individual proto-fibrils (phf-tauT), their values differ due to differences in the orientation of each proto-fibril. Negative values of correlation coefficients for HvT and TvO with a high value of correlation coefficient for HvO relation suggest that phf-tauT is a typical amyloid form.
Figure 1 provides a graphical representation of the superfibril and both chains (E and F) treated as components of the superfibril.
The status of phf-tauO superfibril is visualized in Figure 1A, which shows that the theoretical distribution involves two local maxima, along with hydrophilic fragments that are exposed on the surface. Neither maximum is evident in the observed distribution, however O includes other local maxima, located in areas where low hydrophobicity is expected. It should be noted that each of these local maxima (as well as minima) represents an entire band stretching along the fibril’s long axis. The overlap is due to the repeating pattern that is present in each individual chain, with only the outlying chains exhibiting slightly lower hydrophobicity. On the other hand, the differences between the theoretical distributions are readily apparent since this distribution predicts that hydrophobicity should decrease along with distance from the center. The degree of discordance between T and O can be analyzed by comparing theoretical charts with the observed distributions for chains E and F (which are centrally located and therefore representative for the entire fibril—see Figure 1B,C).
The interface fragment appears to be consistent with the FOD model. Given the central location of the interface, a hydrophobicity peak is expected and—to a certain extent—present in the actual complex. Comparing O with T reveals that two outlying residues exhibit relatively low hydrophobicity, while the central section corresponds to a major spike. Consequently, we rate the interface fragment as being accordant with the model.
A characteristic feature of amyloids is the presence of numerous local maxima in areas where low hydrophobicity (and vice versa) is predicted by the theoretical model. However, it is important to remember that, unlike globular proteins (which may also exhibit this phenomenon), the complexed chains form here bands which stretch along the entire long axis of the fibril. These observations are confirmed by analysis of T and O for chains E and F (treated as components of the superfibril). The discordance between T and O distribution in most of proteins is of local character.

2.3. Properties of Proto-Fibrils

Table 3 presents the hydrophobicity parameters for each proto-fibril. In this case, each proto-fibril is treated as a distinct structural unit. This means that a separate Gaussian is constructed for each proto-fibril (in the preceding section, a shared Gaussian form was computed for the entire superfibril). Results are indicative of an amyloid form: high values of RD and HvO along with very low (even negative) values of HvT and TvO. The correlation coefficients reveal that the structure is dominated by the intrinsic properties of its component residues—an observation that is supported by the observed high values of RD in both models (T-O-R and T-O-H). Thus, the observed distribution is more closely aligned with R (or H) rather than T.
The status of chains E and F, treated as components of their respective proto-fibrils, confirms that they adopt amyloid-like forms, although this effect is less pronounced than in the case of the superfibril (lower values of T-O-H RD and HvO—see Table 3).
Figure 2 provides a visual representation of these results, showing T and O distributions for one proto-fibril (chains A, C, E, G and I) and the status of the E chain within this structural unit.
An asymmetrical distribution of local maxima is observed in the proto-fibril as a result of significant displacement of the system’s central point as compared to the superfibril. Numerous local maxima are present in areas where low hydrophobicity is expected. The involvement of a local maximum in the interface fragment indicates that complexation of protofibrils is generated as the effect of the influence of environment (according to FOD model).

2.4. Properties of Individual Chains Treated as Distinct Structural Units

Our analysis also covers individual chains that are treated as distinct structural units, with a separate Gaussian being plotted for each chain (under the assumption that each chain folds in separation from other chains). To determine the causes of the discordance between the observed and theoretical distributions, we have singled out fragments for which this discordance is particularly evident. Note that we are not dealing with isolated deviations—in many areas, both distributions strongly oppose each other, indicating that the chain does not produce a globule and is likely insoluble due to the lack of a polar surface.
Figure 3 illustrates the status of chains E and F treated as individual structural units.
It is clear that even when analyzed as distinct units, the discussed chains still diverge from the theoretical distribution of hydrophobicity (Figure 3). No C-terminal maxima (predicted by T) are present in the observed distributions.
Regarding phf-tauL, both of the distributions are similar to those calculated for phf-tauO, with only the interface being somewhat different. The observed distribution, while discordant, does not resemble an amyloid (which would appear as a sinusoidal pattern consisting of similar local maxima).
Taking into account the discussed distributions, it is easy to pinpoint fragments where O deviates from T (see Figure 4).
Figure 4 provides a sample set of distributions (T and O) for a single chain—E from phf-tauO. As already noted, the chains differ in detail, while the overall pattern remains largely identical, regardless of the structural unit in question (superfibril, protofibril, or individual chain). The highlighted fragments have been singled out on the basis of visual inspection, supplemented with correlation coefficients computed for successive five-residue segments. Fragments for which HvT and TvO are negative while HvO assumes that a large value will be subjected to further analysis.
Table 4 summarizes the results obtained for all structural units in phf-tauO.
As shown in Table 4, the status of selected fragments is quite similar, regardless of the structural unit in question—in all cases, these fragments are strongly discordant vs. the theoretical distribution.

2.5. Comparative Analysis Involving Theoretical Models

As previously noted, we have carried out an in silico experiment that consisted of predicting the conformation of a tau protein whose sequence matches the discussed amyloids. Our analysis concerned the entire molecule as well as fragments that are highlighted in Figure 4.
Figure 5 presents 3D models of tau polypeptides obtained using software packaged described in the Materials and Methods section. Visual inspection reveals the possible emergence of globular forms: I-Tasser produces four such structures (out of five input cases), while Robetta produces one (out of five). The tendency of the FOD model to produce globular forms should come as no surprise given the model’s propensity to direct hydrophobic residues towards the center of the molecule (due to interactions with the aqueous solvent).
One of the models that was produced by I-Tasser appears to involve a hydrophobic core (in the sense of the FOD model—cf. underlined structures in Figure 5). None of the models produced by Robetta satisfies this criterion. Regarding the FOD model, despite its natural tendency to generate hydrophobic cores, only two among 500 structures analyzed in the course of the study contain a hydrophobic core (i.e., satisfy the RD < 0.5 condition).
Analysis of numerical values that are listed in Table 5, along with visual inspection of 3D forms reveals that some of these fragments adopt structures consistent with the Gaussian distribution. Of particular note is the fragment marked in red on Figure 5 (residues 50–61, numbered 356–357 according to PDB), which does not conform to theoretical predictions in any model. In interpreting this fact, we may refer to the dominant role of the fragment whose sequence does not adapt to the centralized distribution of hydrophobicity. We can speculate that this fragment (GSLDNITHVPGG) is therefore the most amyloidogenic sequence in the set under consideration. This suggestion is supported by the data shown in Table 5, particularly RD values that are either highest or second highest in the entire set.

2.6. Other Fragments of Tau Proteins

In order for our comparative analysis to be as comprehensive as possible, we also include tau proteins (or fragments of tau amyloids) that PDB lists as being capable of adopting non-amyloid conformations. The Tpp (1I8H) represents the 541–553 fragment of the previously described tau protein, in complex with a microtubule—specifically, ww domain complexed with human tau phosphothreonine peptide microtubule-associated protein tau. In this complex, the A chain comprises residues 541–553 (so-called phf-tau), while chain B represents the ww domain (6–44) [33].
The Tpp sequence does not fully match phf-tau (and the others), however we have included it in our analysis due to functional similarities.
Results shown in Table 6 and on Figure 6 reveal that the Tpp complex does not conform to the FOD model. The status of chain A (tau), when analyzed on its own, is also discordant. On the other hand, the same chain conforms to the model when analyzed as part of the complex. This means that chain B creates suitable conditions for chain A to produce a shared hydrophobic core that is consistent with the 3D Gaussian.
Considering the individual fragments of Tpp, it turns out that the fragment 27–29 of chain B is the most discordant, and that eliminating this fragment from calculations lowers the RD value of the complex. This may indicate the conformational alignment between chain B and chain A (given that the presence of chain A disrupts the distribution of hydrophobicity in chain B).
Analysis of our results indicates the need of a “chaperone”, which chain A requires to reach together a conformation consistent with the FOD model.
Tau (267–312) is another protein related to the discussed tau structure. Its sequence matches the short fragment at 306–312 in tau as it is present in phf-tauO. According to [34] this fragment (267–312) of tau protein is bound to microtubules.
The status of the tau chain in tau (267–312) reveals strong discordance with regard to the theoretical model (see Figure 7), with RD = 0.680 (T-O-R) and 0.527 (T-O-H). Correlation coefficients are −0.049, −0.042, and 0.673 for HvT, TvO, and HvO, respectively. These values indicate that the structure of the chain is dominated by the conformational tendencies of individual residues rather than by the external hydrophobic force field. Interaction with microtubules is likely to be the driving force behind conformational adaptation. The structure of the entire complex is not known, however information regarding the interaction of individual residues with the microtubule might explain the discordance that was observed throughout the chain. Only the helical fragment at 295–299 appears accordant with the FOD model (with RD values and correlation coefficients of 0.264, 0.170, 0.876, 0.921, and 0.965, respectively, showing good alignment between the theoretical and observed distribution).

2.7. Peptides

To complete our study of tau-derived structures that are listed in PDB we also need to consider peptides capable of amyloid transformation. The possible mechanism driving this process, discussed in [14,15], remains applicable in the case under consideration.
Peptides that match the tau protein sequence are mostly related to the fragment at 306–311—the short N-terminal fragment of the tau protein present in phf-tauO, phf-tauL, and phf-tauT. A short peptide which does not produce a globular form should not, in principle, be analyzed using the FOD model. Nevertheless, for the sake of completeness, we will present RD (T-O-R) values calculated for such peptides—see Table 7.
The values shown in Table 7 only reveal the type of hydrophobicity distribution, with no assessment of the hydrophobic core structure. Low values indicate that hydrophobic residues are located in the central part of the chain, surrounded by N- and C-terminal hydrophilic residues. As shown, the short VQIVYK fragment, despite including an outlying Val residue, is a good match for the centralized hydrophobic core structure. When additional neighboring residues (adjacent to the 306–311 fragment) are included in analysis, the value of RD increases significantly.
Peptides that are identified as capable of amyloidogenesis appear to adopt amyloid-like conformations themselves. As discussed in [14,15], the distribution of hydrophobicity in a peptide may—regardless of its accordance with the theoretical distribution—give rise to amyloid formation, as long as the environment favors linear complexation of additional peptides, with alternating bands of high and low hydrophobicity emerging along the axis of the fibril. This is visualized in Figure 8, which compares two fringe cases (in terms of RD values).
Figure 8 evidences the appearance of structures that are characterized by linear propagation of hydrophobicity peaks/troughs, which is a precondition of amyloid formation. As highlighted by to-date observations and interpretations, the environment must “support” the creation of such forms. It is thought that under natural conditions the structure of water does not favor the formation of amyloid fibrils.
The value of RD computed for a short peptide implies how many local maxima are present. A low value indicates that hydrophobicity is concentrated in the central part of the peptide (e.g., Tau (306–311B)), while a high value suggests the presence of numerous local maxima (e.g., Tau (305–311)).

2.8. Is It Possible to Differentiate between the Amyloid Fibril and the Fibrillary Structure Present in the Microfilament?

The defining property of amyloids is their fibrillary nature. This phenomenon, however, is not restricted to amyloids. Many biologically active fibrillary proteins exist, often serving as biological scaffolds—this includes polymer microfilaments, such as F-actin (filamentous actin). An example structure of this protein is listed in PDB under ID 3J8I [35].
Analysis of protein 3J8I reveals an alternative approach to a fibrillary structural formation that relies on different mechanisms than those, which drive amyloidogenesis. The PDB structure comprises five monomeric units arranged into a linear complex. Each monomer is a single-domain chain, 375 aa in length, with varying secondary folds: three beta sheets (11 beta strands in total) and 17 helices. The monomers are spatially arranged in a shape of a helix, forming an elongated fibril. Our analysis focuses on the F chain, which is placed in the center of this fibril. We believe that this chain best represents any inner subunit of a long fibril, with adjacent neighbors on either side.
In analyzing the F chain we apply a twofold approach: first, we treat the chain as an independent structure (constructing a 3D Gaussian capsule calculated specifically for that chain), and subsequently we analyze it as part of the complex (with a broader Gaussian encapsulating the entire complex). The former approach enables us to determine the status of the chain itself, while the latter provides clues regarding its role in the formation of a fibril.
The same observations that we relied on when an analyzing the tau amyloid (high values of RD in both configurations; negative HvT and TvO correlation coefficients and strongly positive HvO) will be sought in our study of F-actin to determine whether the conditions that give rise to amyloid fibrils also apply in the presented case.
The above mentioned set of parameters shows that T not only deviates from O, but can, in some respects, be viewed as its polar opposite. This is taken as evidence that the given conformation is driven by intrinsic hydrophobicity of individual residues.
Table 8 presents FOD parameters that describe F-actin (as listed in PDB under ID 3J8I), which were derived from T and O distributions plotted in Figure 9. It appears that the entire complex, as well as both versions of the F chain, deviate from the theoretical distribution, with no monocentric hydrophobic core being observed in either case. Note that Table 8 only lists the status of selected fragments—those, whose amyloid-like conformation may be important in light of the current discussion regarding identification of amyloid forms.
Results show that local amyloid-like properties may be attributed (in varying degrees) to the beta sheet. Such localized amyloid-like folds can indeed be found in many biologically active proteins (e.g., antifreeze proteins that contain solenoid fragments) [36]. Likewise, the beta sheet found in the lysozyme may also be regarded as amyloid-like [37]. (Note that this particular beta sheet plays an important role given its proximity to the active site—it even contributes one of the catalytic residues of the lysozyme). It seems that the presence of a similar structure in actin is not a unique phenomenon, especially given the structure of its immediate neighborhood. As it turns out, local amyloid-like folds in biologically proteins are typically bracketed by “stop” signals (or “caps”), which prevent unchecked linear propagation. They do so by ensuring that the structure, as a whole, conforms to the theoretical distribution of hydrophobicity and mediating entropically advantageous contact with water. This is highlighted in Table 8 with “stop” annotations.
In the scope of our analysis we also computed FOD correlation coefficients for successive fragments of the input chain while using a 5 aa moving frame. This reveals the exact placement of residues which exhibit amyloid-like characteristics. Eliminating such residues lowers the value of RD (although not below 0.5—see the “No neg CC” annotations in Table 8). Visual inspection of both profiles (theoretical—T and observed—O) reveals residues that contribute to the discordance (these are highlighted in Figure 9 and marked in Figure 10, which presents the protein’s 3D structure). It is worth noting that these residues together comprise only 10% of the chain. Eliminating visually inspected residues brings RD down below 0.5, which means that the remainder of the chain conforms to the monocentric distribution of hydrophobicity.
The amyloid-like beta sheet-1 is characterized by the linear propagation of alternating bands that differ in terms of hydrophobicity. This effect manifests itself as a strong discordance between T and O profiles, where—in some cases—the observed distribution appears to be a polar opposite of the theoretical distribution.
Linear propagation can be observed by studying the status of successive fragments that comprise the beta sheet. It is therefore interesting to speculate about the participation of such beta sheets in formation of a complex with a clearly fibrillary nature. The sheets in question are dispersed and do not form a continuous band of alternating hydrophobicity. Consequently, they cannot be regarded as a structural scaffold for the complex. What is more, the beta sheets that are contributed by different chains are not in contact (as shown in Figure 11). The presence of “stop” fragments, also shown in Figure 11, which arrest linear propagation, suggests that amyloid-like conditions are intended to remain local and not dominate the structure. Similar “caps” can be found in many other proteins, which include amyloid-like fragments [38] and prevent the unrestricted elongation of such structures.
In summarizing our comparative analysis of fibrillary structures, it should be noted that these structures owe their existence to different mechanisms. An amyloid emerges as a consequence of linear propagation of alternating bands of high and low hydrophobicity, whereas globular proteins form complexes via nonbinding interactions (including salt bridges and hydrogen bonds). In the latter case, even when amyloid-like fragments can be found in the proteins’ structures, they are dispersed and protected by “stoppers”, which prevent them from interacting with one another to form complexes.
In effect, we can state that the structure referred to as a “fibril” might be produced in various ways. Linear propagation of hydrophobicity bands is the prerequisite of amyloid formation (as well as a useful criterion for identifying amyloids), whereas other fibrillary structures (such as F-actin) are formed through nonbinding interactions (including salt bridges and hydrogen bonds). Thus, even though the end result (elongated fibril) is similar, the underlying mechanisms differ. Against this background we propose that the criteria listed in this book differentiate amyloids and enable their identification. We also show that the presence of beta folds is not required (e.g., as evidenced by the tau amyloid). Instead, amyloids may form whenever the folding process is driven by the intrinsic properties of individual residues, as confirmed by the parameters that are studied in this work.

3. Discussion

The comparative analysis of proteins associated with amyloid tau confirms the previously stated hypothesis concerning the structural properties of the amyloid. According to this hypothesis, the amyloid is characterized by the presence of alternating bands of variable hydrophobicity. It seems that linear propagation—which can be regarded as contrary to the emergence of a centralized hydrophobic core (as seen in globular proteins)—is a characteristic property of amyloids. A similar phenomenon can be observed in Aβ amyloids [39,40]. The network of hydrogen bonds that is discussed in numerous studies [11,41] favors this type of conformation and is thought to be associated with the linear properties of beta folds. In the tau amyloid, however, β-strands play a much smaller role than in other known amyloids. This suggests that while hydrogen bonds are important, their role is not necessarily linked to β-structures.
Hydrophobicity is capable of binding together proximate charged residues, however, electrostatic interactions should, in principle, prevent such clustering. Under such conditions only hydrophobic forces can result in the observed arrangement. Thus, a conformation that is driven by intrinsic hydrophobicity (and does not generate a central hydrophobic core) may be regarded as both the cause and the mechanism of amyloid transformation.
The FOD model recognizes several possible forms for the tau superfibril. This diversity is likely caused by interactions between the solvent and the emerging amyloid. We suggest that, while phf-tauO and phf-tauL emerge as the effect of the influence of surrounding water, in phf-tauT, the structure is driven by the specific band-like arrangement of hydrophobicity in the amyloid itself.
The tau protein, whose task is to mediate interaction with microtubules, must align itself to the complexed object. When the protein is subjected to folding on its own, in an independent manner, it may adopt a globular conformation and remain soluble. An open question is why the same protein undergoes complexation in a form which does not resemble a globule. As shown, a chain that is sequentially identical to the amyloid fragment of the tau chain cannot produce a globular structure. In this context, microtubules may be viewed as a “chaperone”, which ensures that the protein adopts its intended conformation, required for biological activity.
Conclusions that are related to the process of amyloidogenesis and the role of the FOD model in explaining this process, all point to the need for further research into the properties of the aqueous solvent. While we possess good knowledge of the properties of ice, the corresponding “normal” (or physiological) condition of liquid water is poorly understood—for example, we are still unsure of why the density of water peaks at 4 degrees C. This may explain the recent uptake in investigations that aim to explain such phenomena [41,42,43,44,45,46]. We believe that these studies may also cast a new light on the process of amyloidogenesis, which—in all likelihood—is associated with the (heretofore unknown) influence of the force field exerted by the surrounding water. This field should be modeled as a continuum rather than (as is common practice in modern molecular dynamics packages) as a collection of distinct molecules. The FOD model provides a good baseline for such research.
The analysis allows for distinguishing of critical short sequences especially resistant to adopt the conformation accordant with the expected uni-centric hydrophobic distribution. This phenomenon is also observed in other amyloids, especially Aβ(1–42) amyloid [39].

4. Materials and Methods

4.1. Data

The analysis concerns tau protein amyloids listed in PDB as capable of forming highly ordered superfibrils. In addition, we also consider selected fragments of the tau protein, including short peptides. Table 9 gives the full list of structures subjected to analysis.
Table 8 includes tau superfibrils (phf-tauL, phf-tauO, phf-tauT), smaller structural units (including individual chains—tau (267–312)), as well as complexes with other proteins (Tpp). We also consider individual peptides that are widely characterized as capable of forming amyloid structures.
All of the above structures are subjected to FOD characterization in the context of the superfibril, the protofibril and the individual chain. Our analysis further extends to peptides whose composition is similar or identical to PDB sequences. The status of such molecules is determined by computing their RD coefficients. It should be noted that seeking proper hydrophobic cores in very short peptides (<15 aa) makes little sense—such peptides are characterized while using FOD criteria only in order to provide a coherent platform for comparative studies. The FOD model provides useful information regarding the relationship between each residue’s intrinsic hydrophobicity and its placement in a fully folded chain.

4.2. Folding of Peptides—Components of Amyloid Structures

Peptide sequences which form parts of the tau amyloid (e.g., 306–378, as listed under phf-tauO) have been subjected to folding simulations while using Robetta [26,28] and I-Tasser [29,30], as well as to simulations based on the FOD model [53]. This operation can be regarded as an in silico experiment whose aim is to provide alternatives to structures generated by specialized 3D structure prediction software. Our goal is to identify theoretical opportunities for alternative folds (unlike those listed under phf-tauO and similar entries). The globular forms that are generated by the FOD model may provide clues regarding the discordance between the theoretical distribution of hydrophobicity and the actual location of hydrophobicity maxima/minima. A ranking list of the resulting structures may be composed in order to identify factors that increase similarities between the theorized conformation and the corresponding amyloid form.
Robetta is a software package that is aimed at the modeling and analysis of protein structures [27,28]. It is a strong performer in successive editions of the CASP challenge, which focuses on predicting the 3D conformations of input residue sequences [54]. Robetta works in the following manner: the user is asked to input a sequence of amino acids comprising a given protein chain. This sequence is then subdivided into fragments (called domains) while using the “Ginzu” hierarchical scanning algorithm. The algorithm recognises fragments homologous to sequences for which the preferred secondary conformation has been established on the basis of experimental studies. Such homologous areas are detected by (in the order of accuracy) BLAST, PSI-BLAST [55], FFAS03 [56], and 3D-Jury [57] taking as input the sequences produced in the preceding step. The identified domains are modeled by applying a comparative modeling protocol, while all other chain fragments are treated as linkers (if they consist of fewer than 50 residues) or are assigned to structural families as defined in the Pfam-A database [58] using HMMER [59]. Fragments and sequences that have not been recognized as putative domains are analyzed via MSA of the full-length target derived from a PSI-BLAST search against the NCBI non-redundant (NR) protein sequence database [60]. Putative domains identified through Pfam-A and MSA are modeled using de novo structure prediction. Finally, following assembly, side chains are modeled by applying Monte Carlo algorithms [61]. The description is based on [60].
I-Tasser (Iterative Threading ASSEmbly Refinement) is a software package which can predict the structure of a protein given its sequence. In this application, prediction bases on querying PDB for templates using the multiple threading approach. I-Tasser is a strong contender in CASP challenges, topping the ranking in editions 7 through 12 [62,63,64].
The user submits a sequence of amino acids, which is then compared (by LOMETS [65]) to template proteins with similar structural characteristics. An optimal template is then selected and overlapping fragments are assembled into an output model while using replica-exchange Monte Carlo simulations [66], while differing fragments are modeled ab initio. If LOMETS is unable to identify a suitable template, the entire structure is subjected to ab initio modeling. The next step involves a search for low energy states (using SPICKER [67]) in the resulting chain via clustering simulation decoys. This is followed by the reassembly of the template protein starting with SPICER cluster centroid, however this time the simulation is guided by spatial constraints that are provided by TM-align on the basis of LOMETS templates and PDB data. The purpose of the second iteration is to remove steric clashes as well as to refine the global topology of the cluster centroids. The decoys generated in the second simulations are clustered and the lowest energy structures are selected [68]. The final step involves the construction of a detailed model from the available structures via optimization of the hydrogen bond network using REMO. Further information can be found in [68].
Robetta computations were carried out while using the publicly available service [27]. I-Tasser computations were carried out using [68]. FOD model computations were carried out using PL-Grid platform on “Cyfronet” Computer Center AGH Krakow infrastructure [69], a detailed description of which can be found in [70].
The FOD model involves two intermediate folding stages: the early-stage intermediate [71,72,73,74] and the late-stage intermediate [72,73,74]. The initial step, which is meant to generate a starting structure for further optimization, is omitted since the conformation listed under 5O3O is taken as the starting structure (treated as early-stage in this case). This chain is then immersed in an aqueous solvent, whose effects are modeled using a 3D Gaussian (as an external force field). In line with the FOD model, hydrophobic residues tend to congregate at the center of the protein body while hydrophilic residues are exposed on its surface. The process produces a prominent hydrophobic core that is encapsulated by a hydrophilic “shell” (with near-zero values of hydrophobicity on the surface). Optimization of hydrophobic interactions and optimization of nonbinding internal interactions is carried out in an alternating fashion, with each step being repeated several times.
Nonbonding interactions are optimized using Gromacs 4.6.5 software suite (Groningen, The Netherlands) [75], available on the PL-Grid infrastructure at ACK Cyfronet AGH Kraków [69]. FOD-based optimization aims to minimize differences between the idealized (3D Gauss function) and observed distribution of hydrophobicity in the target protein. The workflow interleaves both procedures in order to converge on the final conformation.
Folding simulations that rely on the FOD model are relevant since they acknowledge the effects exerted by the aqueous solvent, and treat them as a global phenomenon (i.e., external force field producing a molecule-wide hydrophobic core).

4.3. Comparative Analysis

All tertiary conformations that were produced by the modeling algorithms, as well as structures that are listed in PDB, were analyzed with regard to the status of their hydrophobic cores, which is described by the RD (relative distance) coefficient. Comparing RD values brings the information about the degree of disorder in respect to ideal distribution. In consequence, the approach to amyloid form can be assessed. RD expresses the degree of order present in the protein’s hydrophobic core and indirectly indicates whether the protein is globular or not. Generating a globular structure with a prominent hydrophobic core (hydrophobicity peaking at the center of the molecule and decreasing along with distance from the center, becoming very low on the surface) suggests that the given amyloid peptide may, under certain circumstances, adopt a globular conformation. The ranking of protein structures, sorted in the order of decreasing globularity, reveals changes which cause proteins to forfeit their centralized hydrophobic cores and that may—in extreme cases—produce amyloid forms.
The RD coefficient can be computed for two independent cases: T-O-R and T-O-H. The former case expressing the relative distance between the observed distribution (O) and two boundary distributions: theoretical (T) which is given by the 3D Gaussian, and uniform (R, random) where each residue is ascribed a hydrophobicity value of 1/N (N being the number of residues in the input chain). R-distribution represents the case of uniform (absence of any local hydrophobicity concentration) distribution, which is the opposite one versus the centralized distribution. In the latter case the uniform distribution is replaced by a distribution corresponding to the intrinsic hydrophobicity of each residue in the input chain (H). Comparing both values reveals factors that guide the folding process (this is particularly true in the T-O-H case). A high value of RD (T-O-H) indicates that folding is dominated by the intrinsic properties of each residue with no regard to cooperative generation of a shared hydrophobic core. When this type of distribution is repeated in successive fragment of the polypeptide, the result is a linear sequence of alternating bands of high and low hydrophobicity. This, in turn, enables the unrestricted elongation of the fibril. An interpretation of this phenomenon (referred to as “ladders”) can also be found in [76].
A comparative assessment of T-O-R and T-O-H coefficients in fibrillary/amyloid structures as well as the structures that are produced by various folding algorithms may enable us to identify the “seeds” of linear propagation. FOD criteria have previously been used to assess the distribution of hydrophobicity in structures published by the CASP project [77].
It is hard to compare the interpretation based on FOD with other methods due to the fact that the hydrophobic interaction is underestimated in the discussion concerning amyloid transformation. However, some aspects of intrinsically disordered proteins that were extensively investigated [78] remain in agreement with the results of the analysis of these proteins in respect to the FOD model [79]. The development of techniques as cryo-electron microscopy [80] as well as solid state NMR [81] makes the availability of amyloid structures possible. Structuralization of water is recently in focus of attention, especially the ordering of water on surface [82] what remains in close relation to our model

Author Contributions

Conceptualization, I.R. and L.K.; Methodology, D.D. and M.G.; Software, M.B. and M.G.; Validation, D.D. and M.G.; Formal Analysis, I.R.; Investigation, I.R.; Resources, I.R. and M.B.; Writing-Original Draft Preparation, I.R.; Writing-Review & Editing, I.R. and M.B.; Visualization, M.B.; Supervision, L.K.; Project Administration, Z.W; Funding Acquisition, M.P.

Funding

This research was funded by Jagiellonian University—Medical College grant number [K/ZDS/006363 and K/ZDS/006366] and PLGrid Infrastructure—Cyfronet AGH, University of Science and Technology, 30-059 Krakow, Al. Mickiewicza 30, Poland.

Acknowledgments

We are indebted to Anna Śmietańska for technical work and to Piotr Nowakowski for translation services.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Wilcock, G.K.; Esiri, M.M. Plaques, tangles and dementia. A quantitative study. J. Neurol. Sci. 1982, 56, 343–356. [Google Scholar] [CrossRef]
  2. Bloom, G.S. Amyloid-β and tau: The trigger and bullet in Alzheimer disease pathogenesis. JAMA Neurol. 2014, 71, 505–508. [Google Scholar] [CrossRef] [PubMed]
  3. Rafii, M.S.; Lukic, A.S.; Andrews, R.D.; Brewer, J.; Rissman, R.A.; Strother, S.C.; Wernick, M.N.; Pennington, C.; Mobley, W.C.; Ness, S.; et al. PET Imaging of Tau Pathology and Relationship to Amyloid, Longitudinal MRI, and Cognitive Change in Down Syndrome: Results from the Down Syndrome Biomarker Initiative (DSBI). J. Alzheimers Dis. 2017, 60, 439–450. [Google Scholar] [CrossRef] [PubMed]
  4. Chiti, F.; Dobson, C.M. Protein Misfolding, Functional Amyloid, and Human Disease. Ann. Rev. Biochem. 2006, 75, 333–366. [Google Scholar] [CrossRef] [PubMed]
  5. Uversky, V.N. Unusual biophysics of intrinsically disordered proteins. Biochim. Biophys. Acta 2013, 1834, 932–951. [Google Scholar] [CrossRef] [PubMed]
  6. Dobson, C.M. Protein folding and misfolding. Nature 2003, 426, 884–890. [Google Scholar] [CrossRef] [PubMed]
  7. Fichou, Y.; Eschmann, N.A.; Keller, T.J.; Han, S. Conformation-based assay of tau protein aggregation. Methods Cell Biol. 2017, 141, 89–112. [Google Scholar] [CrossRef] [PubMed]
  8. Eakin, C.M.; Berman, A.J.; Miranker, A.D. A native to amyloidogenic transition regulated by a backbone trigger. Nat. Struct. Mol. Biol. 2006, 13, 202–208. [Google Scholar] [CrossRef] [PubMed]
  9. Maccioni, R.B.; Muñoz, J.P.; Barbeito, L. The Molecular Bases of Alzheimer’s Disease and Other Neurodegenerative Disorders. Arch. Med. Res. 2001, 32, 367–381. [Google Scholar] [CrossRef]
  10. Kovacs, G.G. Neuropathology of tauopathies: Principles and practice. Neuropathol. Appl. Neurobiol. 2015, 41, 3–23. [Google Scholar] [CrossRef] [PubMed]
  11. Fitzpatrick, A.W.P.; Falcon, B.; He, S.; Murzin, A.G.; Murshudov, G.; Garringer, H.J.; Crowther, R.A.; Ghetti, B.; Goedert, M.; Scheres, S.H.W. Cryo-EM structures of tau filaments from Alzheimer’s disease. Nature 2017, 547, 185–190. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  12. Bai, X.C.; McMullan, G.; Scheres, S.H. How cryo-EM is revolutionizing structural biology. Trends Biochem. Sci. 2015, 40, 49–57. [Google Scholar] [CrossRef] [PubMed]
  13. Cleveland, D.W.; Hwo, S.Y.; Kirschner, M.W. Physical and chemical properties of purified tau factor and the role of tau in microtubule assembly. J. Mol. Biol. 1977, 116, 227–247. [Google Scholar] [CrossRef]
  14. Roterman, I.; Banach, M.; Kalinowska, B.; Konieczny, L. Influence of the Aqueous Environment on Protein Structure—A Plausible Hypothesis Concerning the Mechanism of Amyloidogenesis. Entropy 2016, 18, 351. [Google Scholar] [CrossRef]
  15. Roterman, I.; Banach, M.; Konieczny, L. Application of the Fuzzy Oil Drop Model Describes Amyloid as a Ribbonlike Micelle. Entropy 2017, 19, 167. [Google Scholar] [CrossRef]
  16. Banach, M.; Konieczny, L.; Roterman, I. The fuzzy oil drop model, based on hydrophobicity density distribution, generalizes the influence of water environment on protein structure and function. J. Theor. Biol. 2014, 359, 6–17. [Google Scholar] [CrossRef] [PubMed]
  17. Banach, M.; Konieczny, L.; Roterman, I. Ligand-binding-site recognition. In Protein Folding In Silico; Roterman-Konieczna, I., Ed.; Woodhead Publishing (Currently Elsevier): Oxford/Cambridge, UK; Philadelphia, PA, USA; New Dehli, India, 2012; pp. 79–93. ISBN 978-1-907568-17-6. [Google Scholar]
  18. Banach, M.; Konieczny, L.; Roterman, I. Use of the fuzzy oil drop model to identify the complexation area in protein homodimers. In Protein Folding In Silico; Roterman-Konieczna, I., Ed.; Woodhead Publishing (Currently Elsevier): Oxford/Cambridge, UK; Philadelphia, PA, USA; New Dehli, India, 2012; pp. 95–122. ISBN 978-1-907568-17-6. [Google Scholar]
  19. Kalinowska, B.; Banach, M.; Wiśniowski, Z.; Konieczny, L.; Roterman, I. Is the hydrophobic core a universal structural element in proteins? J. Mol. Model. 2017, 23, 205. [Google Scholar] [CrossRef] [PubMed]
  20. Banach, M.; Konieczny, L.; Roterman, I. Why do antifreeze proteins require a solenoid? Biochimie 2018, 144, 74–84. [Google Scholar] [CrossRef] [PubMed]
  21. Banach, M.; Konieczny, L.; Roterman, I. Secondary and super-secondary structure of proteins in light of the structure of hydrophbic core. In Supersecondary Structure of Proteins; Kister, A., Ed.; Springer: Cham, Switzerland, 2018; in press. [Google Scholar]
  22. Kalinowska, B.; Banach, M.; Konieczny, L.; Roterman, I. Hydrophobic Core in DNA Interacting Proteins. Entropy 2015, 17, 1477–1507. [Google Scholar] [CrossRef]
  23. Condello, C.; Lemmin, T.; Stöhr, J.; Nick, M.; Wu, Y.; Maxwell, A.M.; Watts, J.C.; Caro, C.D.; Oehler, A.; Keene, C.D.; et al. Structural heterogeneity and intersubject variability of Aβ in familial and sporadic Alzheimer’s disease. Proc. Natl. Acad. Sci. USA 2018, 115, E782–E791. [Google Scholar] [CrossRef] [PubMed]
  24. Serpell, L.C. Alzheimer’s amyloid fibrils: Structure and assembly. Biochim. Biophys. Acta 2000, 1502, 16–30. [Google Scholar] [CrossRef]
  25. Trovato, F.; Fumagalli, G. Molecular simulations of cellular processes. Biophys. Rev. 2017, 9, 941–958. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Berman, H.M.; Westbrook, J.; Feng, Z.; Gilliland, G.; Bhat, T.N.; Weissig, H.; Shindyalov, I.N.; Bourne, P.E. The Protein Data Bank. Nucleic Acids Res. 2000, 28, 235–242. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Ovchinnikov, S.; Park, H.; Kim, D.E.; DiMaio, F.; Baker, D. Protein structure prediction using Rosetta in CASP12. Proteins 2018, 86, 113–121. [Google Scholar] [CrossRef] [PubMed]
  28. Robetta: Full-Chain Protein Structure Prediction Server. Available online: http://robetta.bakerlab.org (accessed on 5 July 2018).
  29. Zhang, Y. I-TASSER: Fully automated protein structure prediction in CASP8. Proteins 2009, 77, 100–113. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  30. Yang, J.; Zhang, Y. I-TASSER server: New development for protein structure and function predictions. Nucleic Acids Res. 2015, 43, W174–W181. [Google Scholar] [CrossRef] [PubMed]
  31. Dygut, J.; Kalinowska, B.; Banach, M.; Piwowar, M.; Konieczny, L.; Roterman, I. Structural Interface Forms and Their Involvement in Stabilization of Multidomain Proteins or Protein Complexes. Int. J. Mol. Sci. 2016, 17, 1741. [Google Scholar] [CrossRef] [PubMed]
  32. De Beer, T.A.P.; Berka, K.; Thornton, J.M.; Laskowski, R.A. PDBsum additions. Nucleic Acids Res. 2014, 42, D292–D296. [Google Scholar] [CrossRef] [PubMed]
  33. Wintjens, R.; Wieruszeski, J.M.; Drobecq, H.; Rousselot-Pailley, P.; Buée, L.; Lippens, G.; Landrieu, I. 1H NMR study on the binding of Pin1 Trp-Trp domain with phosphothreonine peptides. J. Biol. Chem. 2001, 276, 25150–25156. [Google Scholar] [CrossRef] [PubMed]
  34. Kadavath, H.; Jaremko, M.; Jaremko, Ł.; Biernat, J.; Mandelkow, E.; Zweckstetter, M. Folding of the Tau Protein on Microtubules. Angew. Chem. Int. Ed. Engl. 2015, 54, 10347–10351. [Google Scholar] [CrossRef] [PubMed]
  35. Galkin, V.E.; Orlova, A.; Vos, M.R.; Schröder, G.F.; Egelman, E.H. Near-atomic resolution for one state of F-actin. Structure 2015, 23, 173–182. [Google Scholar] [CrossRef] [PubMed]
  36. Roterman, I.; Banach, M.; Konieczny, L. Antifreeze proteins. Bioinformation 2017, 13, 400–401. [Google Scholar] [CrossRef] [PubMed]
  37. Banach, M.; Kalinowska, B.; Konieczny, L.; Roterman, I. Possible Mechanism of Amyloidogenesis of V Domains. In Self-Assembled Molecules—New Kind of Protein Ligands—Supramolecular Ligands; Roterman, I., Konieczny, L., Eds.; Springer: Cham, Switzerland, 2018; pp. 77–100. ISBN 978-3-319-65638-0. [Google Scholar]
  38. Roterman, I.; Banach, M.; Konieczny, L. Propagation of Fibrillar Structural Forms in Proteins Stopped by Naturally Occurring Short Polypeptide Chain Fragments. Pharmaceuticals 2017, 10, 89. [Google Scholar] [CrossRef] [PubMed]
  39. Roterman, I.; Dułak, D.; Gadzała, M.; Banach, M.; Konieczny, L. Structural analysis of the Aβ(11-42) amyloid fibril based hydrophobicity distribution. Unpublished work. 2018. [Google Scholar]
  40. Dułak, D.; Gadzała, M.; Banach, M.; Konieczny, L.; Roterman, I. Analysis of the structure of Aβ(15-40) amyloid. Acta Biochim. Pol. 2018, submitted. [Google Scholar]
  41. Chiti, F.; Dobson, C.M. Protein misfolding, amyloid formation and human disease: A summary of progress over the last decade. Annu. Rev. Biochem. 2017, 86, 27–68. [Google Scholar] [CrossRef] [PubMed]
  42. Kim, K.H.; Späh, A.; Pathak, H.; Perakis, F.; Mariedahl, D.; Amann-Winkel, K.; Sellberg, J.A.; Lee, J.H.; Kim, S.; Park, J.; et al. Maxima in the thermodynamic response and correlation functions of deeply supercooled water. Science 2017, 358, 1589–1593. [Google Scholar] [CrossRef] [PubMed]
  43. Gallo, P.; Stanley, H.E. Supercooled water reveals its secrets. Science 2017, 358, 1543–1544. [Google Scholar] [CrossRef] [PubMed]
  44. Palmer, J.C.; Martelli, F.; Liu, Y.; Car, R.; Panagiotopoulos, A.Z.; Debenedetti, P.G. Metastable liquid-liquid transition in a molecular model of water. Nature 2014, 510, 385–388. [Google Scholar] [CrossRef] [PubMed]
  45. Tanaka, H. A self-consistent phase diagram for supercooled water. Nature 1996, 380, 328–330. [Google Scholar] [CrossRef]
  46. Poole, P.H.; Sciortino, F.; Essmann, U.; Stanley, H.E. Phase behaviour of metastable water. Nature 1992, 360, 324–328. [Google Scholar] [CrossRef]
  47. Sawaya, M.R.; Sambashivan, S.; Nelson, R.; Ivanova, M.I.; Sievers, S.A.; Apostol, M.I.; Thompson, M.J.; Balbirnie, M.; Wiltzius, J.J.; McFarlane, H.T.; et al. Atomic structures of amyloid cross-β spines reveal varied steric zippers. Nature 2007, 447, 453–457. [Google Scholar] [CrossRef] [PubMed]
  48. Liu, C.; Sawaya, M.R.; Cheng, P.N.; Zheng, J.; Nowick, J.S.; Eisenberg, D. Characteristics of amyloid-related oligomers revealed by crystal structures of macrocyclic β-sheet mimics. J. Am. Chem. Soc. 2011, 133, 6736–6744. [Google Scholar] [CrossRef] [PubMed]
  49. Landau, M.; Sawaya, M.R.; Faull, K.F.; Laganowsky, A.; Jiang, L.; Sievers, S.A.; Liu, J.; Barrio, J.R.; Eisenberg, D. Towards a pharmacophore for amyloid. PLoS Biol. 2011, 9, e1001080. [Google Scholar] [CrossRef] [PubMed]
  50. Liu, C.; Zhao, M.; Jiang, L.; Cheng, P.N.; Park, J.; Sawaya, M.R.; Pensalfini, A.; Gou, D.; Berk, A.J.; Glabe, C.G.; et al. Out-of-register β-sheets suggest a pathway to toxic amyloid aggregates. Proc. Natl. Acad. Sci. USA 2012, 109, 20913–20918. [Google Scholar] [CrossRef] [PubMed]
  51. Wiltzius, J.J.; Landau, M.; Nelson, R.; Sawaya, M.R.; Apostol, M.I.; Goldschmidt, L.; Soriaga, A.B.; Cascio, D.; Rajashankar, K.; Eisenberg, D. Molecular mechanisms for protein-encoded inheritance. Nat. Struct. Mol. Biol. 2009, 16, 973–978. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. De la Cruz, M.J.; Hattne, J.; Shi, D.; Seidler, P.; Rodriguez, J.; Reyes, F.E.; Sawaya, M.R.; Cascio, D.; Weiss, S.C.; Kim, S.K.; et al. Atomic-resolution structures from fragmented protein crystals with the cryoEM method MicroED. Nat. Methods 2017, 14, 399–402. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Konieczny, L.; Bryliński, M.; Roterman, I. Gauss-function-Based model of hydrophobicity density in proteins. In Silico Biol. 2006, 6, 15–22. [Google Scholar] [PubMed]
  54. Simons, K.T.; Bonneau, R.; Ruczinski, I.; Baker, D. Ab initio protein structure prediction of CASP III targets using ROSETTA. Proteins 1999, 37 (Suppl. 3), 171–176. [Google Scholar] [CrossRef]
  55. Altschul, S.F.; Madden, T.L.; Schaffer, A.A.; Zhang, J.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [PubMed]
  56. Rychlewski, L.; Jaroszewski, L.; Li, W.; Godzik, A. Comparison of sequence profiles. Strategies for structural predictions using sequence information. Protein Sci. 2000, 9, 232–241. [Google Scholar] [CrossRef] [PubMed]
  57. Ginalski, K.; Elofsson, A.; Fischer, D.; Rychlewski, L. 3D-Jury: A simple approach to improve protein structure predictions. Bioinformatics 2003, 19, 1015–1018. [Google Scholar] [CrossRef] [PubMed]
  58. Bateman, A.; Birney, E.; Cerruti, L.; Durbin, R.; Etwiller, L.; Eddy, S.R.; Griffiths-Jones, S.; Howe, K.L.; Marshall, M.; Sonnhammer, E.L. The Pfam protein families database. Nucleic Acids Res. 2002, 30, 276–280. [Google Scholar] [CrossRef] [PubMed]
  59. Eddy, S.R. Profile hidden Markov models. Bioinformatics 1998, 14, 755–763. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  60. Kim, D.E.; Chivian, D.; Baker, D. Protein structure prediction and analysis using the Robetta server. Nucleic Acids Res. 2004, 32, W526–W531. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  61. Kuhlman, B.; Baker, D. Native protein sequences are close to optimal for their structures. Proc. Natl. Acad. Sci. USA 2000, 97, 10383–10388. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  62. Zhang, Y. I-TASSER server for protein 3D structure prediction. BMC Bioinform. 2008, 9, 40. [Google Scholar] [CrossRef] [PubMed]
  63. Roy, A.; Kucukural, A.; Zhang, Y. I-TASSER: A unified platform for automated protein structure and function prediction. Nat. Protoc. 2010, 5, 725–738. [Google Scholar] [CrossRef] [PubMed]
  64. Yang, J.; Yan, R.; Roy, A.; Xu, D.; Poisson, J.; Zhang, Y. The I-TASSER Suite: Protein structure and function prediction. Nat. Methods 2015, 12, 7–8. [Google Scholar] [CrossRef] [PubMed]
  65. Wu, S.; Zhang, Y. LOMETS: A local meta-threading-server for protein structure prediction. Nucleic Acids Res. 2007, 35, 3375–3382. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  66. Swendsen, R.H.; Wang, J.S. Replica Monte Carlo simulation of spin glasses. Phys. Rev. Lett. 1986, 57, 2607–2609. [Google Scholar] [CrossRef] [PubMed]
  67. Zhang, Y.; Skolnick, J. SPICKER: A clustering approach to identify near-native protein folds. J. Comput. Chem. 2004, 25, 865–871. [Google Scholar] [CrossRef] [PubMed]
  68. About I-TASSER Server. Available online: https://zhanglab.ccmb.med.umich.edu/I-TASSER/about.html (accessed on 5 July 2018).
  69. PLGrid User Manual. Available online: https://docs.cyfronet.pl/pages/viewpage.action?pageId=23626756https://docs.cyfronet.pl/display/PLGDoc/User+manual (accessed on 5 July 2018).
  70. Roterman-Konieczna, I. Protein Folding—Protein Folding Versus Protein Structure Prediction; Woodehead Publishing (Currently Elsevier): Oxford/Cambridge, UK; Philadelphia, PA, USA; New Dehli, India, 2012; ISBN 978-1-907568-17-6. [Google Scholar]
  71. Roterman, I.; Konieczny, L.; Banach, M.; Jurkowski, W. Intermediates in the protein folding process: A computational model. Int. J. Mol. Sci. 2011, 12, 4850–4860. [Google Scholar] [CrossRef] [PubMed]
  72. Kalinowska, B.; Krzykalski, A.; Roterman, I. Contingency Table Browser—Prediction of early stage protein structure. Bioinformation 2015, 11, 486–488. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  73. Jurkowski, W.; Bryliński, M.; Konieczny, L.; Wiíniowski, Z.; Roterman, I. Conformational subspace in simulation of early-stage protein folding. Proteins 2004, 55, 115–127. [Google Scholar] [CrossRef] [PubMed]
  74. Kalinowska, B.; Fabian, P.; Stąpor, K.; Roterman, I. Statistical dictionaries for hypothetical in silico model of the early-stage intermediate in protein folding. J. Comput. Aided Mol. Des. 2015, 29, 609–618. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  75. Van Der Spoel, D.; Lindahl, E.; Hess, B.; Groenhof, G.; Mark, A.E.; Berendsen, H.J. GROMACS: Fast, flexible, and free. J. Comput. Chem. 2005, 26, 1701–1718. [Google Scholar] [CrossRef] [PubMed]
  76. Biancalana, M.; Makabe, K.; Koide, S. Minimalist design of water-soluble cross-β architecture. Proc. Natl. Acad. Sci. USA 2010, 107, 3469–3474. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  77. Gadzała, M.; Kalinowska, B.; Banach, M.; Konieczny, L.; Roterman, I. Determining protein similarity by comparing hydrophobic core structure. Heliyon 2017, 3, e00235. [Google Scholar] [CrossRef] [PubMed]
  78. Uversky, V.N. Amyloidogenesis of natively unfolded proteins. Curr. Alzheimer Res. 2008, 5, 260–287. [Google Scholar] [CrossRef] [PubMed]
  79. Kalinowska, B.; Banach, M.; Konieczny, L.; Marchewka, D.; Roterman, I. Intrinsically disordered proteins—Relation to general model expressing the active role of the water environment. Adv. Protein Chem. Struct. Biol. 2014, 94, 315–346. [Google Scholar] [CrossRef] [PubMed]
  80. Gremer, L.; Schölzel, D.; Schenk, C.; Reinartz, E.; Labahn, J.; Ravelli, R.B.G.; Tusche, M.; Lopez-Iglesias, C.; Hoyer, W.; Heise, H.; et al. Fibril structure of amyloid-β(1-42) by cryo-electron microscopy. Science 2017, 358, 116–119. [Google Scholar] [CrossRef] [PubMed]
  81. Tycko, R. Indirect detection in solid state NMR: An illustrious history and a bright future. J. Magn. Reson. 2018, 288, 122–123. [Google Scholar] [CrossRef] [PubMed]
  82. Macias-Romero, C.; Nahalka, I.; Okur, H.I.; Roke, S. Optical imaging of surface chemistry and dynamics in confinement. Science 2017, 357, 784–788. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Theoretical (T—blue) and observed (O—red) distributions of hydrophobicity in: (A) the phf-tauO superfibril; (B,C) chains E and F as components of the phf-tauO superfibril. Orange highlighting marks residues which comprise the inter-fibril interface area (331–338).
Figure 1. Theoretical (T—blue) and observed (O—red) distributions of hydrophobicity in: (A) the phf-tauO superfibril; (B,C) chains E and F as components of the phf-tauO superfibril. Orange highlighting marks residues which comprise the inter-fibril interface area (331–338).
Ijms 19 02910 g001
Figure 2. Theoretical (T—blue) and observed (O—red) distributions: (A) as calculated for the proto-fibril (chains A, C, E, G and I) listed under phf-tauO; (B) for the E chain treated as a component of the proto-fibril. Orange highlighting marks residues which comprise the inter-fibril interface (331–338).
Figure 2. Theoretical (T—blue) and observed (O—red) distributions: (A) as calculated for the proto-fibril (chains A, C, E, G and I) listed under phf-tauO; (B) for the E chain treated as a component of the proto-fibril. Orange highlighting marks residues which comprise the inter-fibril interface (331–338).
Ijms 19 02910 g002
Figure 3. Theoretical (T—blue) and observed (O—red) distributions for phf-tauO: (A) chain E treated as a distinct structural unit; and (B) chain F treated as a distinct structural unit. Orange highlighting marks residues which comprise the inter-fibril interface (331–338).
Figure 3. Theoretical (T—blue) and observed (O—red) distributions for phf-tauO: (A) chain E treated as a distinct structural unit; and (B) chain F treated as a distinct structural unit. Orange highlighting marks residues which comprise the inter-fibril interface (331–338).
Ijms 19 02910 g003
Figure 4. Theoretical (T—blue) and observed (O—red) distribution for the chain E from phf-tauO (the distribution is similar in all structural units). Fragments where O significantly deviates from T are highlighted by different colors (see Table 4) that match the colors of three-dimensional (3D) presentations on Figure 5.
Figure 4. Theoretical (T—blue) and observed (O—red) distribution for the chain E from phf-tauO (the distribution is similar in all structural units). Fragments where O significantly deviates from T are highlighted by different colors (see Table 4) that match the colors of three-dimensional (3D) presentations on Figure 5.
Ijms 19 02910 g004
Figure 5. 3D presentation of structures obtained using I-Tasser upper row (BF), Robetta central row (GK) and the fuzzy oil drop (FOD) model bottom row (LP). The top structure (A) is the structure phf-tauO as appears in PDB database. The colors used correspond to highlights on Figure 4. Models whose status is consistent with FOD predictions (RD < 0.5) have been additionally underlined in blue (D,M,N).
Figure 5. 3D presentation of structures obtained using I-Tasser upper row (BF), Robetta central row (GK) and the fuzzy oil drop (FOD) model bottom row (LP). The top structure (A) is the structure phf-tauO as appears in PDB database. The colors used correspond to highlights on Figure 4. Models whose status is consistent with FOD predictions (RD < 0.5) have been additionally underlined in blue (D,M,N).
Ijms 19 02910 g005
Figure 6. Distribution of hydrophobicity in chain A in the Tpp in complex (T—blue, O—red, H—green): (A) when treated as a distinct structural unit; (B) when treated as a component of the complex. Good alignment between O, T, and H can be observed in the latter case.
Figure 6. Distribution of hydrophobicity in chain A in the Tpp in complex (T—blue, O—red, H—green): (A) when treated as a distinct structural unit; (B) when treated as a component of the complex. Good alignment between O, T, and H can be observed in the latter case.
Ijms 19 02910 g006
Figure 7. T (blue) and O (red) hydrophobicity distributions in Tau (267–312) (2MZ7), revealing overall strong discordance.
Figure 7. T (blue) and O (red) hydrophobicity distributions in Tau (267–312) (2MZ7), revealing overall strong discordance.
Ijms 19 02910 g007
Figure 8. H (green), T (blue) and O (red) distributions for selected peptides: (A) Tau (305–311), (B) Tau (306–311B); (C) for (A) and (D) for (B)—pseudo-3D view, presenting the observed hydrophobicity of a theorized structure of fibrillary tangles formed by linear propagation of the corresponding peptides. Dashed lines on (A) represent deleted residue number 11 Tau (305–311) according to data available in PDB.
Figure 8. H (green), T (blue) and O (red) distributions for selected peptides: (A) Tau (305–311), (B) Tau (306–311B); (C) for (A) and (D) for (B)—pseudo-3D view, presenting the observed hydrophobicity of a theorized structure of fibrillary tangles formed by linear propagation of the corresponding peptides. Dashed lines on (A) represent deleted residue number 11 Tau (305–311) according to data available in PDB.
Ijms 19 02910 g008
Figure 9. T (blue) and O (red) profiles for the F chain from F-actin (PDB ID 3J8I), divided into two parts for visibility: 1–200 (A) and 201–375 (B). Highlighted positions mark residues that cause discordance between those distributions (on the basis of visual inspection). The remainder of the chain is regarded as accordant with the theoretical distribution. It likely contributes to the protein’s structural stability—under the assumption that a well-ordered hydrophobic core and the presence of disulfide bonds both play a role in stabilizing tertiary conformations.
Figure 9. T (blue) and O (red) profiles for the F chain from F-actin (PDB ID 3J8I), divided into two parts for visibility: 1–200 (A) and 201–375 (B). Highlighted positions mark residues that cause discordance between those distributions (on the basis of visual inspection). The remainder of the chain is regarded as accordant with the theoretical distribution. It likely contributes to the protein’s structural stability—under the assumption that a well-ordered hydrophobic core and the presence of disulfide bonds both play a role in stabilizing tertiary conformations.
Ijms 19 02910 g009
Figure 10. 3D presentation of F chain from F-actin (PDB ID 3J8I). Beta sheets are displayed in different shades of yellow. Red fragments distinguish residues that cause discordance between T and O distributions. These fragments correspond to highlighted parts of hydrophobicity profiles presented in Figure 9.
Figure 10. 3D presentation of F chain from F-actin (PDB ID 3J8I). Beta sheets are displayed in different shades of yellow. Red fragments distinguish residues that cause discordance between T and O distributions. These fragments correspond to highlighted parts of hydrophobicity profiles presented in Figure 9.
Ijms 19 02910 g010
Figure 11. 3D presentation of F chain (blue) and H chain (red) from F-actin (PDB ID 3J8I). Beta sheets are displayed in different shades of yellow (same three shades in each chain). “Stop” fragments are shown in green.
Figure 11. 3D presentation of F chain (blue) and H chain (red) from F-actin (PDB ID 3J8I). Beta sheets are displayed in different shades of yellow (same three shades in each chain). “Stop” fragments are shown in green.
Ijms 19 02910 g011
Table 1. Relative Distance (RD) values and correlation coefficients computed for the superfibril, proto-fibrils and individual chains treated as components of the superfibril.
Table 1. Relative Distance (RD) values and correlation coefficients computed for the superfibril, proto-fibrils and individual chains treated as components of the superfibril.
Phf-tauRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
phf-tauO0.7450.6870.0120.2590.675
phf-tauL0.7310.6690.0130.2960.646
phf-tauT0.7240.6410.0080.3010.716
Inter-fibril interface
phf-tauO0.4010.5500.6660.7540.950
phf-tauL0.3680.5870.6960.7720.929
phf-tauT0.5270.386−0.164−0.0030.850
Chains in superfibril
Phf-tauO
Chain E0.7410.7000.0110.2370.731
Chain F0.7610.7220.0100.2240.731
Chains E and F0.7520.7110.0130.2300.731
Phf-tauL
Chain E0.7280.6750.0170.2740.696
Chain F0.7470.6980.0120.2650.696
Chains E and F0.7410.6900.0140.2680.606
Phf-tauT
Chain E0.7320.6660.0340.4000.769
Chain F0.7280.666−0.0180.1360.785
Chain E and F0.7300.6660.0080.2730.777
Table 2. Status of individual chains treated as components of the superfibril. The presentation of phf-tauL is omitted since due to its similarity to phf-tauO.
Table 2. Status of individual chains treated as components of the superfibril. The presentation of phf-tauL is omitted since due to its similarity to phf-tauO.
ChainRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
Phf-tauOPhf-tauTPhf-tauOPhf-tauTPhf-tauOPhf-tauTPhf-tauOPhf-tauTPhf-tauOPhf-tauT
A0.7800.8010.6960.6810.0100.0410.1540.4410.6690.710
B0.7870.7390.7110.6000.031−0.0220.3030.1150.6360.723
C0.7560.7660.7150.7050.0140.0380.2020.3940.7310.772
D0.7790.7280.7410.6650.022−0.0200.2480.1140.7300.788
E0.7410.7320.7000.6660.0160.0340.2370.4000.7310.769
F0.7610.7280.7220.6660.010−0.0100.2240.1360.7310.785
G0.7330.7000.6900.6280.0160.0280.2650.3940.7310.770
H0.7450.7330.7040.6700.003−0.0170.1920.1570.7320.786
I0.7270.6600.6130.5220.0150.0220.3430.4050.6520.686
J0.7500.7380.6340.618−0.019−0.0140.1370.1920.6910.707
Table 3. Status of individual chains treated as components of their respective proto-fibrils. The table lists only values obtained for the proto-fibril (chains A, C, E, G and I) (Differences with regard to the other proto-fibril are negligible.
Table 3. Status of individual chains treated as components of their respective proto-fibrils. The table lists only values obtained for the proto-fibril (chains A, C, E, G and I) (Differences with regard to the other proto-fibril are negligible.
PDB IDRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
Phf-tauO0.6610.607−0.0120.0890.772
Chain E0.6790.430−0.0270.0950.548
Chain F0.6790.439−0.0270.0950.548
Phf-tauL0.6640.595−0.0220.0820.767
Chain E0.6740.410−0.0390.0910.545
Chain F0.6740.410−0.0390.0910.545
Phf-tauT0.6730.602−0.0170.1180.773
Chain E0.6830.415−0.0330.0980.550
Chain F0.6830.415−0.0330.0960.551
Table 4. RD values in both models (T-O-R and T-O-H), along with HvT, TvO, and HvO correlation coefficients for chains E and F analyzed as part of the superfibril, as part of a protofibril and on their own. Figure 4 illustrates the division of the chain into individual fragments.
Table 4. RD values in both models (T-O-R and T-O-H), along with HvT, TvO, and HvO correlation coefficients for chains E and F analyzed as part of the superfibril, as part of a protofibril and on their own. Figure 4 illustrates the division of the chain into individual fragments.
Phf-tauO
Chain E—Superfibril
FragmentRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
1–60.5210.4050.2070.2700.855
16–260.8200.771−0.126−0.4910.874
23–300.4960.495−0.3500.4350.580
39–490.5750.5440.0400.1470.971
50–610.8510.802−0.106−0.7820.446
61–720.5630.491−0.157−0.0240.700
CHAIN0.7410.7000.0160.2370.731
Chain F—Superfibril
FragmentRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
1–60.6270.5120.0380.1170.875
16–260.8220.773−0.129−0.4950.875
23–300.4960.495−0.3180.4650.580
39–490.5630.5330.0540.1530.971
50–610.8610.814−0.112−0.7200.444
61–720.5830.512−0.1530.0060.699
CHAIN0.7610.7220.0100.2240.731
Chain E—Protofibril
FragmentRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
1–60.6680.553−0.0060.0640.875
16–260.6150.5570.288−0.1410.875
23–300.5710.5730.146−0.1510.810
39–490.6570.623−0.139−0.1100.971
50–610.7000.6220.016−0.5580.443
61–720.6130.540−0.221−0.2000.700
CHAIN0.6610.607−0.0120.0890.772
Chain F—Protofibril
FragmentRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
1–60.6680.553−0.0050.0640.875
16–260.6150.5580.286−0.1440.875
23–300.5720.5730.143−0.1550.809
39–490.6570.623−0.139−0.1100.971
50–610.7000.6210.016−0.5580.442
61–720.6120.539−0.221−0.2000.698
CHAIN0.6610.607−0.0120.0890.772
Chain E—Individual
FragmentRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
1–60.7410.372−0.093−0.0370.543
16–260.6760.4330.286−0.3870.696
23–300.5820.3870.117−0.2280.620
39–490.6870.388−0.1140.0170.888
50–610.7000.4890.057−0.5370.116
61–720.6560.398−0.237−0.2400.507
CHAIN0.6790.430−0.0270.0950.548
Chain F—Individual
FragmentRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
1–60.7410.372−0.093−0.0370.543
16–260.6760.4330.286−0.3870.696
23–300.5820.3870.117−0.2280.620
39–490.6870.388−0.1140.0170.888
50–610.7000.4890.057−0.5370.115
61–720.6560.397−0.233−0.2400.507
CHAIN0.6790.430−0.0270.0950.548
Table 5. RD values obtained using I-Tasser, Robetta and FOD. In addition to RD values (for T-O-R and T-O-H) the table also lists HvT, TvO and HvO correlation coefficients. The division into fragments is consistent with the one shown in Figure 4. Values in bold distinguish the status accordant with FOD model (RD < 0.5).
Table 5. RD values obtained using I-Tasser, Robetta and FOD. In addition to RD values (for T-O-R and T-O-H) the table also lists HvT, TvO and HvO correlation coefficients. The division into fragments is consistent with the one shown in Figure 4. Values in bold distinguish the status accordant with FOD model (RD < 0.5).
I-Tasser
IT-1RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.4110.1420.0320.6000.469
16–260.5090.2000.3080.6590.706
23–300.6410.3570.3080.2190.415
39–490.4940.2830.5010.3060.841
50–610.6140.6180.1690.1960.466
61–720.6270.2320.0890.1760.559
CHAIN0.5190.2660.2000.4730.574
IT-2RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.8250.138−0.0160.5510.242
16–260.6070.2790.3720.0390.730
23–300.5030.2360.3360.2260.651
39–490.7650.313−0.230−0.3580.849
50–610.6670.5480.027−0.4160.122
61–720.6050.367−0.486−0.1330.475
CHAIN0.6350.353−0.0280.1720.519
IT-3RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.7290.258−0.438−0.1090.416
16–260.5580.1790.3810.5670.830
23–300.4080.1730.6070.6350.805
39–490.3690.1910.2820.6590.828
50–610.7010.7020.102−0.1960.587
61–720.5450.2850.2120.1620.638
CHAIN0.4780.2620.2700.5240.684
IT-4RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.4070.1830.2160.7020.554
16–260.8170.2300.1910.1250.730
23–300.7460.4680.253−0.4530.345
39–490.4320.1400.1900.4940.818
50–610.6300.5850.4260.0940.367
61–720.6010.3500.4120.1930.781
CHAIN0.6050.3150.2290.2520.581
IT-5RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.5450.138−0.0240.1980.590
16–260.5370.2320.3210.3130.840
23–300.6630.3230.2540.0400.670
39–490.5040.1620.4440.4260.848
50–610.7230.8620.1940.2080.872
61–720.4770.2570.5010.5240.858
CHAIN0.5130.3270.2600.4180.659
Robetta
ROB-1RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.6270.4740.0100.1650.633
16–260.4730.3800.6870.4000.808
23–300.5660.6190.6190.1570.830
39–490.6060.3810.0420.2750.889
50–610.6500.4950.070−0.1050.376
61–720.6750.405−0.27−0.1820.594
CHAIN0.7450.568−0.0010.1090.683
ROB-2RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.5350.345−0.3000.2180.690
16–260.6010.5180.4480.0040.809
23–300.6210.6750.388−0.1470.824
39–490.8390.2920.0140.0360.904
50–610.6810.5310.083−0.2970.404
61–720.6990.364−0.2500.0170.552
CHAIN0.6560.4220.0320.2490.611
ROB-3RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.6000.382−0.495−0.0120.623
16–260.5230.4100.5760.1440.795
23–300.5720.5100.463−0.0380.821
39–490.8420.473−0.173−0.1890.955
50–610.6460.5120.127−0.1680.394
61–720.6980.396−0.1820.2260.576
CHAIN0.6820.4440.0740.3510.650
ROB-4RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.6240.519−0.0130.1410.766
16–260.4410.3290.7140.4640.806
23–300.5930.5570.6180.1710.831
39–490.7820.453−0.120−0.0720.939
50–610.6580.5750.160−0.1870.344
61–720.6780.346−0.1870.0160.568
CHAIN0.7480.5420.0220.1540.657
ROB-5RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.5450.111−0.0450.5710.558
16–260.7010.3410.2010.1660.916
23–300.7480.3140.251−0.2840.666
39–490.3720.2200.2550.6760.809
50–610.6430.7120.2580.3980.640
61–720.4320.1550.5130.5870.794
CHAIN0.5660.3510.2720.3750.733
FOD
FOD-1RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.7740.288−0.555−0.0400.023
16–260.3460.1160.2380.6430.503
23–300.8080.0410.510−0.3190.339
39–490.6020.187−0.2640.2670.606
50–610.7330.7140.369−0.2910.457
61–720.5680.217−0.2690.2790.200
CHAIN0.5030.2530.1260.5630.390
FOD-2RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.5020.017−0.6870.723−0.013
16–260.3150.1100.1650.7700.527
23–300.6750.2610.2490.1050.521
39–490.5430.271−0.1750.2580.693
50–610.5930.6980.3750.0920.706
61–720.5220.191−0.2110.2680.268
CHAIN0.3640.1910.0930.6650.410
FOD-3RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.6050.104−0.3730.4400.386
16–260.3250.1070.0540.7000.552
23–300.5800.1940.3280.1270.574
39–490.5750.187−0.1710.3130.689
50–610.6760.6820.331−0.2820.554
61–720.3380.0800.1350.6030.395
CHAIN0.3690.1750.1520.6350.414
FOD-4RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.4920.137−0.3490.4430.231
16–260.3870.1360.2240.6010.546
23–300.7050.4700.301−0.4140.290
39–490.5890.140−0.2110.4530.410
50–610.8400.744−0.151−0.3660.252
61–720.5560.3190.0380.0920.357
CHAIN0.6480.3460.0590.1500.395
FOD-5RDCorrelation Coefficient
FragmentT-O-RT-O-HHvTTvOHvO
1–60.6200.075−0.3630.4840.004
16–260.5130.1210.0250.5760.446
23–300.6680.247−0.1130.0310.070
39–490.8060.205−0.180.2200.646
50–610.8070.588−0.458−0.3260.033
61–720.6380.344−0.154−0.0220.348
CHAIN0.6440.313−0.0490.2000.357
Table 6. Parameters describing the 541–553 fragment (chain A) of the tau protein Tpp in complex with the ww domain (residues 6–44).
Table 6. Parameters describing the 541–553 fragment (chain A) of the tau protein Tpp in complex with the ww domain (residues 6–44).
TppRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
Complex0.6800.5460.1710.2350.700
Chain A (tau) in complex0.4520.1520.3650.7300.756
Chain A(tau) individual0.6300.2400.0800.3480.753
Table 7. Properties of peptides which match the N-terminal fragment of the tau protein. The table lists PDB ID, sequences (indicating which fragments are identical to the tau chain) and RD (T-O-R) values. Positions + aa inform about positions of amino acids present in the chain but these residues are not the object of analysis.
Table 7. Properties of peptides which match the N-terminal fragment of the tau protein. The table lists PDB ID, sequences (indicating which fragments are identical to the tau chain) and RD (T-O-R) values. Positions + aa inform about positions of amino acids present in the chain but these residues are not the object of analysis.
PeptideSequencePositionRD
Tau (306–311)VQIVYK306–3110.164
Tau (306–310)VQIVYK + LA306–311 + 2 aa0.332
Tau (306–311B)VQIVYK306–3110.160
Tau (305–311)AS + VQIVYK + AEFYK2 aa + 306–311 + 5 aa0.707
Tau (623–628)VQIVYK306–3110.174
Tau (306–311C)VQIVYK306–3110.170
Table 8. Parameters describing the microfilament structure as present in F-actin. β-sheet-1—sheet with starting numbers 8–11, β-sheet-2—sheet with starting numbers 150–154, β-sheet-3—sheet with starting numbers 70–72. Stop sign—fragment “stopping” linear propagation: *—fragment 351–374, **—fragment 325–331, ***—fragment 173–176, No neg CC—status of the chain with residues representing negative correlation coefficients (CC) for HvT and TvO with high HvO, No selected—status of the chain with residues identified as discordant under visual analysis of the T and O profile—shown in Figure 9 and Figure 10.
Table 8. Parameters describing the microfilament structure as present in F-actin. β-sheet-1—sheet with starting numbers 8–11, β-sheet-2—sheet with starting numbers 150–154, β-sheet-3—sheet with starting numbers 70–72. Stop sign—fragment “stopping” linear propagation: *—fragment 351–374, **—fragment 325–331, ***—fragment 173–176, No neg CC—status of the chain with residues representing negative correlation coefficients (CC) for HvT and TvO with high HvO, No selected—status of the chain with residues identified as discordant under visual analysis of the T and O profile—shown in Figure 9 and Figure 10.
F-actinRDCorrelation Coefficient
T-O-RT-O-HHvTTvOHvO
Complex0.7830.7040.0580.0810.722
Chain F in complex 0.6140.5120.1110.1550.730
β-sheet-1 *0.7430.465−0.163−0.4110.723
β-sheet-2 **0.5930.415−0.08−0.0360.635
β-sheet-3 ***0.5320.593−0.398−0.2530.939
Stop sign β10.4450.3460.5040.4200.867
No neg CC0.6050.4670.1210.1790.717
No selected0.4880.3910.3120.4250.748
C-term0.4450.3460.5040.4200.867
Chain F individual0.6410.5440.1400.3820.709
P-P0.7540.5940.0620.1360.687
No P-P0.6350.5430.1580.3570.737
β-sheet-1 *0.7740.503−0.161−0.3930.721
β-sheet-2 **0.5560.402−0.137−0.3670.489
β-sheet-3 ***0.4380.5070.2210.4600.940
Stop sign β10.4920.3760.6100.4140.866
No neg CC0.6150.5020.2090.4020.714
No selected0.4870.4000.3380.5280.729
C-term0.4920.3760.6100.4140.866
Table 9. Set of proteins subjected to analysis, along with an indication of chain length and complexation capabilities. The rightmost column provides references.
Table 9. Set of proteins subjected to analysis, along with an indication of chain length and complexation capabilities. The rightmost column provides references.
PDB IDCharacteristics LengthComplexReference
Tau—amyloid
5O3O phf-tauOMicrotubule-associated protein tau73 aa10 chains[11]
5O3L phf-tauLMicrotubule-associated protein tau73 aa10 chains[11]
5O3T phf-tauTMicrotubule-associated protein tau73 aa10 chains[11]
Tau—non-amyloid
2MZ7 Tau (267–312)Tau (267–312) bound to microtubules46 aa [34]
1I8H TppPin1 ww domain complexed with human tau phosphothreonine peptide14 aaComplex[33]
Peptides
2ON9 Tau (306–311A)Amyloid forming peptide VQIVYK from the repeat region of tau (in tau 306-311)6 aa [47]
3Q9G Tau (306–310)VQIVY segment from Alzheimer’s tau displayed on 42-membered macrocycle scaffold Cyclic pseudo-peptide vqiv(4bf)(orn)(hao)kl(orn)5 aa [48]
3OVL Tau (306–311B)Microtubule-associated protein. VQIVYK (residues 306–311)6 aa [49]
4E0M Tau (305–311)SVQIVYK segment from human tau (305–311) displayed on 54-membered macrocycle scaffold (form i)7 aa [50]
4NP8 Tau (623–628)Structure of an amyloid forming peptide VQIVYK from the second repeat region of tau (alternate polymorph) (in tau 623–628)7 aa [51]
5K7N Tau (306–311C)tau VQIVYK peptide6 aa [52]
Fibrilar form as appears in microfilament
F-actin Actin, alpha skeletal muscle375 aa [53]

Share and Cite

MDPI and ACS Style

Dułak, D.; Gadzała, M.; Banach, M.; Ptak, M.; Wiśniowski, Z.; Konieczny, L.; Roterman, I. Filamentous Aggregates of Tau Proteins Fulfil Standard Amyloid Criteria Provided by the Fuzzy Oil Drop (FOD) Model. Int. J. Mol. Sci. 2018, 19, 2910. https://doi.org/10.3390/ijms19102910

AMA Style

Dułak D, Gadzała M, Banach M, Ptak M, Wiśniowski Z, Konieczny L, Roterman I. Filamentous Aggregates of Tau Proteins Fulfil Standard Amyloid Criteria Provided by the Fuzzy Oil Drop (FOD) Model. International Journal of Molecular Sciences. 2018; 19(10):2910. https://doi.org/10.3390/ijms19102910

Chicago/Turabian Style

Dułak, Dawid, Małgorzata Gadzała, Mateusz Banach, Magdalena Ptak, Zdzisław Wiśniowski, Leszek Konieczny, and Irena Roterman. 2018. "Filamentous Aggregates of Tau Proteins Fulfil Standard Amyloid Criteria Provided by the Fuzzy Oil Drop (FOD) Model" International Journal of Molecular Sciences 19, no. 10: 2910. https://doi.org/10.3390/ijms19102910

APA Style

Dułak, D., Gadzała, M., Banach, M., Ptak, M., Wiśniowski, Z., Konieczny, L., & Roterman, I. (2018). Filamentous Aggregates of Tau Proteins Fulfil Standard Amyloid Criteria Provided by the Fuzzy Oil Drop (FOD) Model. International Journal of Molecular Sciences, 19(10), 2910. https://doi.org/10.3390/ijms19102910

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop