Prion proteins are known to misfold into a range of different aggregated forms, showing different phenotypic and pathological states. Understanding strain specificities is an important problem in the field of prion disease. Little is known about which PrPSc structural properties and molecular mechanisms determine prion replication, disease progression and strain phenotype. The aim of this work is to investigate, through a mathematical model, how the structural stability of different aggregated forms can influence the kinetics of prion replication. The model-based results suggest that prion strains with different conformational stability undergoing in vivo replication are characterizable in primis by means of different rates of breakage. A further role seems to be played by the aggregation rate (i.e. the rate at which a prion fibril grows). The kinetic variability introduced in the model by these two parameters allows us to reproduce the different characteristic features of the various strains (e.g., fibrils' mean length) and is coherent with all experimental observations concerning strain-specific behavior.
Prion diseases are caused by the accumulation of a cellular prion protein with an altered conformation, which acts as a catalyst for the further recruitment and the modification of the normal form of the protein. Protein polymerization appears to have a central role in the progression of the disease, an aspect shared with several other neurodegenerative diseases. The aim of this work is to investigate at the kinetic level the “prion strain phenomenon”, i.e., the ability of prion proteins to misfold into a range of different aggregated forms exhibiting different replication and propagation properties. The dynamics of prion replication is investigated with the help of a mathematical model. We relate a measurement accessible in vitro (prion structural stability) to a mathematical description of the fibrils' kinetics in vivo. The analysis of the model suggests that the replication kinetics of the different prion strains is characterizable by means of two parameters, representing the rates of breakage and aggregation. This result is coherent with various experimental findings concerning strain-specific behavior, such as, for example, the observation of the fibril mean length of the various strains.
Citation: Zampieri M, Legname G, Altafini C (2009) Investigating the Conformational Stability of Prion Strains through a Kinetic Replication Model. PLoS Comput Biol 5(7): e1000420. doi:10.1371/journal.pcbi.1000420
Editor: Eugene I. Shakhnovich, Harvard University, United States of America
Received: August 6, 2008; Accepted: May 20, 2009; Published: July 3, 2009
Copyright: © 2009 Zampieri et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work is supported by a grant from Illy Caffè, Trieste, Italy. The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Prions are infectious agents composed solely of proteins, whose replication does not rely upon the presence of nucleic acids . Although the molecular mechanisms of prion replication are poorly understood, the current working hypothesis is based on the assumption that prions replicate by means of an autocatalytic process which converts cellular prion protein () to the disease-associated misfolded PrP isoform (). This process of replication of a prion depends upon the capacity of the pathogenic protein form to bind to and to catalyze the conversion of existing intermediate molecules. Recent studies  have observed that the prion protein can misfold into a range of different aggregated forms derived from a continuum of structural conformation templates  from which different phenotypic and pathological states derive. The ability of the same encoded protein to encipher a multitude of phenotypic states is known as the “prion strain phenomenon” . Prion strains are defined as infectious isolates that, when transmitted to identical hosts, exhibit the following distinct prion disease phenotypes:
- i) Proteinase K (PK) digestion profile;
- ii) Incubation time;
- iii) Histopathological lesion profiles;
- iv) Specific neuronal target areas.
A reason for the strain phenomenon can be the association of to several disease conformations, characterizable by means of a different stability against denaturation, different post-translational modifications (e.g. glycosylation) and distinct cleavage sites. These observations are reinforced by , where it is reported that the amyloid fibrils (formed by the 40-residues -amyloid peptide) with different morphologies have significantly different molecular structures. These differences are shown to be self-propagating and to be associated with different toxicities, suggesting the possibility for a structural origin of prion strains. Moreover, recent studies on prion disease have confirmed that the incubation time is related not only to the inoculum dosage and the prion protein expression, but also to the resistance of prion strains against denaturation  in terms of the concentration of guanidine hydrochloride (Gdn-HCl) required to denaturate 50% of the disease-causing protein (see Text S1 for further discussions). Other studies have highlighted a strong relationship between the stability of the prion protein against denaturation and neuropathological lesion profiles ,. Lesions due to stable prions tend to show large vacuolations localized in specific small brain regions, whilst lesions due to unstable prion strains show a less intense vacuolation and are more widely distributed in the brain. Apart from these properties, crucial details of the molecular mechanisms enabling the characterization of different prion strains are still missing. For example, neither structural characterizations of , nor maps of protein-protein interactions have so far been provided, and even the biological function is unclear. Hence, in order to use the existing data to gain some insight into the properties of the different prion strains, we decided to follow a model-based approach.
In this paper, using a well established model for the kinetics of the in vivo prion replication , we relate the evidence about conformational stability to the parameters of the model describing the evolution in time of the fibril length. The main points we deduce from our analysis are:
- i) In terms of the model, the key parameter describing strain-dependent replication kinetics is the fibril breakage rate.
- ii) A precise fitting of the model prediction to the experimental data is obtained assuming that also the aggregation rate changes with the strain. In particular, a functional dependence on the breakage rate is assumed.
- iii) The prediction of the model is that the stability against denaturation is inversely correlated to both breakage and aggregation rates and directly correlated with the mean length of fibrils.
- iv) By fitting experimental data, we can quantitatively predict the fibril length distributions associated to different prion strains.
Protein polymerization seems to have a central role in the progression of the prion pathology, an aspect shared with several other neurodegenerative diseases associated with different aggregating proteins, such as Alzheimer's (A ), Parkinson's (-synuclein) and Huntington's (huntingtin) diseases. The aggregation kinetics of amyloid peptides has been studied extensively (see ,), and has shed light on the wide range of amyloid aggregation mechanisms observed. Many modeling approaches have been introduced for this purpose in recent years, e.g. theoretical models consisting of nonlinear ordinary differential equations (ODEs), two-dimensional lattice-based statistical models and molecular dynamics simulations , –. In this paper we explore a mathematical description of the prion replication dynamics through nonlinear ODEs. This class of models explain the appearance of the disease by means of a bistability induced by a quadratic term, as in classical epidemic models . The model we used is drawn from , and is based on a nucleated polymerization mechanism  (see Materials and Methods). This approach has been shown to overcome the limitations of the “heterodimer model”  and to be a reasonable simplification of the “cooperative autocatalysis” approach . Furthermore, it is able to explain the kinetics of spontaneous generation , the association between infectivity and aggregated PrP, the linear appearance of the fibrils and to take into account fundamental processes of an in vivo replication (i.e. fibrils splitting), all while remaining relatively mathematically tractable. Moreover, its dynamical behavior has been extensively studied ,, and experimental measurements were used in  to provide an estimation of the full set of parameters for a particular prion strain. The model has three state variables (Eq. 10) describing the amount of monomer (), polymer () and the mass of polymer (), and it involves 6 parameters (see Table 1). We reproduce here only the features essential to discuss the strain dependence of its parameters; the details are covered in Materials and Methods.
In  it has been shown that for any prion strain two parameters, the rate of growth () and the reproductive ratio (), can be estimated from experimental data. The former (Eq. 2) represents the exponential growth of the number of infectious particles. The latter (Eq. 3) is defined as the average number of prion fibrils that a single infectious particle can give rise to, before splitting into fibrils smaller than a critical size or being degraded. In other words, represents the ability of the fibrils to survive to critical breakage and degradation events. The equations for the and parameters obtained from the kinetic model of  can be reparametrized in terms of the mean length of the fibrils (1)obtaining:(2)(3)
In order to estimate from experimental measures both parameters ( and ) certain assumptions are necessary (see Materials and Methods for full details). An estimation of and from in vivo experiments and for different prion strains characterized by different values of stability against denaturation () is listed in Table 2. The dataset currently available is limited (as not many prion strains can be fully characterized) and many error sources are potentially affecting the estimation of the parameters. Nevertheless, Figure 1 shows the existence of a negative trend between these two empirical parameters (Pearson correlation = −0.91, p-value = 0.01). If we now turn to the kinetic model and look at the corresponding expressions (Eq. 2, 3) the interesting question is whether such a behavior is predicted by the model itself, and is explainable in terms of some of its parameters, in a way that is both mathematically and biologically plausible. Otherwise stated, we investigate which, if any, among the model parameters best describe the strain variability. The critical size of the nucleus (parameter in the model) plays a marginal role in our analysis and is likely to be a fixed integer, in between 2 and 4, across different strains . Even though it has been argued that a hexamer is the minimum infectious unit , it can be shown that the model-based conclusions are not conditioned by the value of . In addition is clearly independent of the prion strains, so we remain with three possible choices: , and . From Eq. 1, increasing means incrementing and this affects and in a similar manner, so that this parameter alone cannot explain the inverse relationship derived in Figure 1. The same can be said for and which, if increased/decreased, would induce changes of equal sign in and . Different conclusions can be drawn when considering as the only strain-varying parameter. This dependence becomes clearer assuming that fibrils cannot be degraded in the exponential phase (, identical results can be obtained supposing that the degradation of the fibrils scales as the fibrils breakage rate, , see Text S2). Such assumption leads to the following expressions:(4)(5)(6)
The reproductive ratio is plotted against the rate of growth. The downward trend is not well described by the linear model with negative angular coefficient () and an intercept () (dotted blue line). In addition, the model prediction with fixed (dashed-dot black line) fails to precisely represent the data, even if it provides a more reasonable relationship (notice that high stable prions, such as MK4985, would always be associated to positive values). Introducing one more degree of freedom (exponent ) yields a higher value (red line, ). This result corresponds to a prediction of . In addition, we tested a further simplified model version (where is considered to be much smaller than ) according to which (i.e. , shown in green). Similar conclusions could be drawn.
From these simplified formulas it is clear that an increase in the frangibility of the fibers (i.e., in ) produces an increment of (Eq. 7) and a decrement of (Eq. 8) in agreement with the trend in Figure 1. Therefore, from the model we expect to give the best fitting result. As a matter of fact, this relationship (black dash-dotted line in Figure 1) does not provide the optimal fit, although it reproduces the qualitative observed behavior (). The fittings of Figure 1 (see Table 3) suggest that, approximately, (red line) implying that we are observing proportional to and to (see Materials and Methods, Eq. 12). This means that the estimated exponents for are somewhat different from the expected values of () predicted in Eq. 7 and 8. In order to improve the model prediction, we introduce a strain-dependence on a second parameter. The simplest solution suggested by the model for this scope (deducible from Eq. 4 and 5) points to the aggregation rate . By linking to , we are still left with a one-parameter family of models describing the strain-dependence. In doing so, we obtain the estimate (see again Materials and Methods, Eq. 13). This correction yields and , this time respecting the predictions of Eq. 4 and 5. Therefore, on the one hand we can show that at a qualitative level is the only parameter that alone can explain the inverse relationship between and . On the other hand, the variation of by itself is not able to quantitatively describe the experimental data in a precise way. An additional correction, obtained relating to , leads to a substantially improved fitting. Apart from Eq. 4 and 5, our choice of alongside as strain-dependent parameter is suggested by the structure of the model of Eq. 10, in which, of all parameters, those describing the kinetics of fibril aggregation/breakage are the most likely to vary across strains. Both the fitting and the model structure suggest an interplay between and , with partially balancing the effect of .
In the following, we will describe how the previous results can be extended to the stability to denaturation of the prion strains, providing experimental observations in support of our claims. From Figure 2A a direct linear proportionality between and is inferred. Therefore, combining the fitting between and and and , a similar inverse relationship (see Figure 2B) relates and (see Table 3). A point of note is that a linear model (i.e., ) is not only associated to a low coefficient of determination but is also implausible, as it predicts negative values of in correspondence of very stable prion strains (such as MK4985, see Table 4). Owing to the linear proportionality () of Figure 2A, the inferred functional dependencies from extend to (i.e., ). This result, in light of the experimental observations in , contributes to validate the results of the kinetic model and provides us with a simple practical tool to interpret prion strain stability. As a matter of fact, the experimental data in  report a relationship between the chemical stability of yeast prion strains and their structural properties, hence reinforcing our conclusions. In particular, the frangibility of different Sup35NM amyloid conformations was measured and shown to be consistent with an increase in sensitivity to denaturants and proteases. Thus, confirming the main role of the breakage rate, as predicted here by the model. Furthermore, the authors observed also a variation in the aggregation rate (parameter in the model), which was however overcome by the stronger effect of the division rate; an additional observation in agreement with our results, where the best match with the experimental data is obtained for a variation of that only partially compensates for that of . The importance of breakage events for the in vivo prion propagation is also underlined in , where the authors observed that membrane-anchored PrP is necessary for the exponential growth of prion aggregates. In transgenic mice, expressing anchorless prion protein inoculated with different prion strains, the aggregates seem to grow quadratically in time . This feature is explainable by a linear aggregation model (i.e setting equal to 0). Moreover, in , different prion strains show a common inability to induce the disease. The absence of fibrils disruption can prevent the formation of oligomeric species, thus hiding the difference between prion strains. Our model-based analysis suggests that an experiment monitoring the propagation of prion strains lacking the GPI anchor would be useful to characterize in more depth the strain phenomenon.
In (A) and (B) the stability against denaturation is plotted against the reproductive ratio and the rate of growth. A direct proportionality links to . As expected, an inverse proportionality emerges between and , reinforcing the previous results.
In the last part of the Section, we investigate how the prion stability () is reflected in the mean length of the fibrils (). Combining the fitting of Figure 2 with Eq. 6, (and consequently , from Eq. 1) can be inferred directly from and :(9)In Table 4, we compare the approach of Eq. 9 with the results obtained in , where the authors give a complete estimation (including a range of uncertainty) of all the parameters for the RML prion strain (, highlighted in bold in Table 4). The comparison between these two approaches shows that the predictions obtained through Eq. 9 are similar to the values reported in  for the RML strain. In addition, we can compare the values of for the strains inferred from Eq. 9, with the ones computed using Eq. 11 and then imposing equal to the values of  for the RML strain (see Figure S1). Our predictions are approximately within the range of values computed considering constant among strains. This result reinforces the major role of in explaining strain variability. Owing to the fact that is now strain-dependent (), we can also predict the mean length of the fibrils (Eq. 1) for each considered strain (see Table 4, ). For instance the mean length of the fibrils population for two prion strains with different stabilities (e.g. RecMoPrP (89–230) and Sc237) can be compared. For the unstable prion strain (Sc237) this is approximately 7 monomer units, while for the stable prion strain (RecMoPrP (89–230)) it is approximately 14 monomer units. This theoretical approach provides a valuable method to simplify the model characterization. Furthermore, it contributes to understanding the properties associated to prion strains with different stability against guanidine denaturation.
While it is reasonable that the parameters of the kinetic model might all be affected by strain specificities (i.e. stability against denaturation), the dominant contribution seems to be due to the susceptibility to frangibility (i.e. ), with only a minor correction due to . The inverse relationship between and shown in Figure 1 is the main argument in the identification of as the key physical aspect differentiating prion strains. In addition, is suggested as the most plausible and parsimonious correcting factor, in order to improve the data fitting.
Several aspects can influence the estimation of the parameters and the model predictions. For example, the uncertainty affecting the estimation of and (respectively inferred from an exponential curve and from a ratio of exponentials); or the possibility that the breakage rate is not equal across all the different polymer lengths (e.g. mechanical stress can differently affects longer fibers); or even the impact of the mouse age on the model parameters (affecting e.g. the production rate). In spite of these (and potentially many other) disregarded aspects characterizing an in vivo system, this simple model is able to capture and explain the observed data dependencies through arguments supported by multiple independent experimental observations. Our analysis reveals that stable prion strains can be characterized by a “stronger” aggregated structure which is less prone to breakage events. This will further imply a longer mean length of the fibrils. Instead, unstable prion strains are subject to a higher fragmentation rate. The role of is essentially to partially balance the increased breakage and is coherent with the experimental observations in yeast. Furthermore, the increased number of catalytic sites may be also responsible for the shorter incubation time.
As already mentioned, such phenomenon was observed in yeast prions . The yeast prion proteins, although fundamentally different from the mammalian prion proteins, show the same ability to convert into aggregated forms, propagate and be infectious. This simpler unicellular system is a valuable model as it enables a deeper analysis of the fibril formation process , not possible to the same extent in higher organisms.
The framework proposed allows for a model-based analysis of these properties in mammalian prions in vivo. In the context of mammals, our results are consistent with , where fibrils with different conformational stability are generated in vitro from full length mammalian PrP. In that paper, the authors relate the stability to the size of the smallest possible fibrillar fragment without taking into account the kinetics of the replication (reproducing the in vivo behavior). We draw similar conclusions from a different point of view. As a matter of fact, we investigate the dynamic evolution of prion propagation in a multicellular in vivo system, in which molecular and cellular mechanisms are present as well. Our model-based conclusions provide further evidence that in vitro systems and yeast prion propagation mechanisms can be transposed in mammals. Moreover, linking the strain phenomena to dynamical features leads to a characterization of the evolution of the length of the fibrils in vivo.
We can, in fact, speculate (in agreement with ) that stable prion strains exhibit a proliferation of longer fibrils that, upon splitting, still manifest the same stability properties (Figure 3), giving rise to a preferential proliferation of relatively long fibrils with a low toxic effect. On the other hand, less stable prion strains tend to form shorter fibrils, to proliferate faster and to be more neurotoxic.
The cartoon describes the pathways of kinetic replication of two prion strains with a different stability against denaturation: a stable one (high ) and an unstable one (low ) are drawn. These act as templates bringing the same cellular prion protein (triangle) to the two different strain conformations ( ▴→ ▪, ♦). The model assumes that the aggregation of monomers to polymers produces a very fast change of conformation and that this aggregation is unfavorable below a critical size (), which is assumed to be independent of the prion strain in our model. The experimental data suggest that stable prions are characterized by a higher and a corresponding lower . In the model, this is translated into strain-specificity of the rates of breakage and of aggregation (which are both lower for stable prions). This implies that stable fibrils are longer and prefer to proliferate while maintaining themselves as fibrils larger than the nucleus size (pathway on the left). On the contrary, unstable prions are more frangible (i.e. more sensitive to breakage), implying a shorter mean length. This means that breakage events are more likely to be associated with the formation of very short fibrils, even under the critical size. The increase in the aggregation rate is not enough to avoid an increased growth in the number of fibrils. We can therefore hypothesize that an apoptotic pathway is most likely for these last strains (pathway on the right). These conclusions are in agreement with the working hypothesis of oligomer toxicity .
It is worth noting the connection with , where the kinetics of aggregation of amyloid peptides is studied by means of coarse-grained molecular dynamics. The authors showed how the relative stability of -prone states of a polypeptide can influence the pathway of aggregation. Their results suggest that the -stable amyloids follow an aggregation pathway without intermediates, while -unstable amyloids seem to involve on-pathway oligomers.
The characterization of prion strains in terms of polymer mean size is per se a significant observation. It provides a new possible explanation of the observation that stability is correlated with lesion profiles and vacuolation areas. Several hypotheses have been made, such as the existence of a co-factor that supports the conversion of distinct prion strains in precise brain regions. Here, another possibility emerges: the increased size associated to stable prions can decrease their ability to diffuse, and can circumscribe them to small brain regions. On the contrary, oligomers can spread around the brain more easily, causing a more homogeneous damage.
In conclusion, we show that linking the conformational stability property of prions, acquired during in vivo propagation in mammals, to their replication kinetic properties is achievable through a rather simple model. For a wide range of parameters, the model predicts that a higher breakage rate implies shorter and shorter incubation time (in Figure S2 two simulations are compared). Our model-based approach suggests that the amount of information that can be extrapolated from the knowledge of goes beyond the expected incubation time.
Materials and Methods
In vitro prion propagation is characterized phenomenologically by the following properties: (i) a critical concentration threshold below which fibrils cannot form; (ii) a delay before their propagation (which can be eliminated by the addition of seeds of preformed fibrils); (iii) a direct proportionality between the initial rate of fiber growth and the monomer concentration . The overall behavior resembles a sigmoidal growth curve : an exponential growth of infectious particles followed by a plateau. The simplest description of the underlying observed mechanism of protein aggregation consists of a slow continuous nucleation followed by a fast autocatalytic growth. Therefore a simple two-step model is able to reproduce the dynamics of the in vitro prion propagation . An in vivo prion propagation model should explain the fact that the spontaneous prion-induced disease is rare but progresses inevitably after infection, that the incubation period is long and followed by a brief fatal clinical disease and that prions undergo several molecular processes within the cell (e.g. fibrils breakage, degradation, endogenous production). The model derived in  is obtained as a closed form of an infinite set of differential equations describing the variation in time of the monomer and fibrils of each possible length (from to ). The biological mechanisms taken into account are the lengthening at the fiber end by the addition of monomers, the degradation of polymers, and their splitting into smaller polymers. Only if several monomeric molecules are mounted into a highly ordered seed, can further monomeric be recruited and form amyloid aggregates. If, after the breakage, the fibril has a length under the critical size, it degrades instantaneously into normal monomers. The model in Eq. 10 has three state variables, describing the amount of monomer (), polymer () and the mass of polymer (), and it comprises of 6 independent parameters: nucleus size (), rates of production (), degradation (), aggregation (), clearance () and breakage ():(10)
The assumption that is negligible, made in in order to simplify the parameters equations, changes the qualitative behavior of the model, that no longer has two stable steady states but only one, which is unstable. This means that the exponential growth will never reach a plateau. As mentioned in the text, this does not affect our previous considerations, especially in light of the fact that in vivo death occurs during the exponential growth phase (see also the Text S2 for similar conclusions on the full model).
Measuring the parameters and
In this section we summarize the procedures mentioned in  and adopted here to derive a measure for and . The assumptions deemed, in order to measure and from the observed effect of different levels of PrP expression and inoculum dosage, are as follows:
- i) The linear relationship relating the incubation time to the inoculum log dilution reflects the exponential growth of the infectious units.
- ii) The only parameter that varies between two transgenic mice with an altered level of expression is considered to be the production rate ().
- iii) The termination stage (animal death) occurs during the exponential growth phase.
- iv) The level of in the brain at the termination stage can be considered to be the same in all experiments.
Of all assumptions, the last one is the most important. It is considered valid also for transgenic mice expressing different quantities of cellular prion protein. Currently there is wide debate about the cause of cell death in prion neurodegeneration. From knockout mutants, it seems that loss of function is not sufficient to cause cell death. What has been observed is that the conversion of to the isoform has a key role in the disease. In spite of their apparent low neurotoxic effect , fibrils have been proven to be the main ingredient in catalyzing variations of protein conformation . Therefore, it is reasonable to assume that even if toxicity is not directly associated to fibrils aggregates, it has to be closely related to their amount, implying that a critical concentration of is required to provoke cell death. The current working hypothesis is that oligomeric species are the most infectious  and a substantial body of evidence suggests that they are also highly cytotoxic . According to the previous observations, a possible explanation is that an equal mass of prion fibrils with smaller mean size provides a larger number of active sites for catalysis, hence inducing a higher lethality.
In order to extrapolate a measure for we follow the method described in  based on relating the incubation time to the inoculum dose and implying an exponential growth in the number of infectious particles. Taking advantage of these data (e.g. incubation time inoculum dosage), we can infer the parameter just by fitting an exponential growth curve. More precisely, before inoculation of prions, PrP () can be considered at steady state (). After inoculation, it is reasonable to assume that it remains almost constant for a while. According to the model equations, the steady state of the mean polymers distribution length ( in Eq. 1), is typically reached before the exponential phase. Immediately after reaching , the polymer amount () and the polymer mass () start to grow exponentially. Thus, is defined as the dominant mode of this exponential growth (i.e., ) (Eq. 2).
To have an indirect measurement of , the inverse relationship between incubation time and the PrP expression is exploited. We take into account the previous assumptions reporting that the number of infectious units in two inoculated mice expressing different level of PrP (, ) at the times of death (, ) can be considered almost equal. Thus imposing we can derive :(11)where and .
Computing the and exponents
Rather that using Eq. 7 and 8, the exponents and such that , can be computed from the fitted curve in Figure 1. We approximate the numerical value 0.38 of Table 3 with 0.4 (i.e. ). From the above expressions, , which yields , i.e., or . Examples of values on this line are:From Eq. 6 it is clear that the only admissible pair of values is(12)If, following Eq. 4 and 5, we add the extra functional dependence of from as , we can look for a value of that satisfies simultaneouslyyielding:(13)
Disease evolution for different values of b
(0.05 MB PDF)
Incubation time (t_G) vs stability (G) and rate of growth (r)
(0.12 MB PDF)
Full model considerations
(0.04 MB PDF)
The authors would like to gratefully acknowledge the contribution of our reviewers. In particular we want to thank Dr. Joanna Masel for her many incisive comments and useful suggestions. In addition, we thank Diane Latawiec for a careful reading of the manuscript.
Conceived and designed the experiments: MZ GL CA. Analyzed the data: MZ CA. Wrote the paper: MZ CA.
- 1. Prusiner SB (1982) Novel proteinaceous infectious particles cause scrapie. Science 216: 136–144.
- 2. Safar J, Wille H, Itri V, Groth D, Serban H, et al. (1998) Eight prion strains have prp(sc) molecules with different conformations. Nat Med 4: 1157–1165.
- 3. Legname G, Nguyen HO, Peretz D, Cohen FE, DeArmond SJ, et al. (2006) Continuum of prion protein structures enciphers a multitude of prion isolate-specified phenotypes. Proc Natl Acad Sci U S A 103: 19105–19110.
- 4. Peretz D, Williamson RA, Legname G, Matsunaga Y, Vergara J, et al. (2002) A change in the conformation of prions accompanies the emergence of a new prion strain. Neuron 34: 921–932.
- 5. Petkova AT, Leapman RD, Guo Z, Yau WM, Mattson MP, et al. (2005) Self-propagating, molecular-level polymorphism in alzheimer's beta-amyloid fibrils. Science 307: 262–265.
- 6. Legname G, Nguyen HO, Baskakov IV, Cohen FE, Dearmond SJ, et al. (2005) Strain-specified characteristics of mouse synthetic prions. Proc Natl Acad Sci U S A 102: 2168–2173.
- 7. Sigurdson CJ, Nilsson KP, Hornemann S, Manco G, Polymenidou M, et al. (2007) Prion strain discrimination using luminescent conjugated polymers. Nat Methods 4: 1023–1030.
- 8. Masel J, Jansen VA, Nowak MA (1999) Quantifying the kinetic parameters of prion replication. Biophys Chem 77: 139–152.
- 9. Sun Y, Makarava N, Lee CI, Laksanalamai P, Robb FT, et al. (2008) Conformational stability of prp amyloid fibrils controls their smallest possible fragment size. J Mol Biol 376: 1155–1167.
- 10. Tanaka M, Collins SR, Toyama BH, Weissman JS (2006) The physical basis of how prion conformations determine strain phenotypes. Nature 442: 585–589.
- 11. Kryndushkin DS, Alexandrov IM, Ter-Avanesyan MD, Kushnirov VV (2003) Yeast [psi+] prion aggregates are formed by small sup35 polymers fragmented by hsp104. J Biol Chem 278: 49636–49643.
- 12. Morris AM, Watzky MA, Agar JN, Finke RG (2008) Fitting neurological protein aggregation kinetic data via a 2-step, minimal/“Ockham's razor” model: the Finke-Watzky mechanism of nucleation followed by autocatalytic surface growth. Biochemistry 47: 2413–2427.
- 13. Pellarin R, Caflisch A (2006) Interpreting the aggregation kinetics of amyloid peptides. J Mol Biol 360: 882–892.
- 14. Masel J, Genoud N, Aguzzi A (2005) Efficient inhibition of prion replication by prp-fc(2) suggests that the prion is a prp(sc) oligomer. J Mol Biol 345: 1243–1251.
- 15. Nowak M, Krakauer D, Klug A, May RM (1999) Prion infection dynamics. Integrative Biology: Issues, News, and Reviews 1: 3–15.
- 16. Slepoy A, Singh RR, Pázmándi F, Kulkarni RV, Cox DL (2001) Statistical mechanics of prion diseases. Phys Rev Lett 87: 058101.
- 17. Mobley DL, Cox DL, Singh RR, Kulkarni RV, Slepoy A (2003) Simulations of oligomeric intermediates in prion diseases. Biophys J 85: 2213–2223.
- 18. Eigen M (1996) Prionics or the kinetic basis of prion diseases. Biophys Chem 63: 1–18.
- 19. Nowak MA, Bangham CR (1996) Population dynamics of immune responses to persistent viruses. Science 272: 74–79.
- 20. Jarrett JT, Lansbury PT (1993) Seeding “one-dimensional crystallization” of amyloid: a pathogenic mechanism in alzheimer's disease and scrapie? Cell 73: 1055–1058.
- 21. Greer ML, Pujo-Menjouet L, Webb GF (2006) A mathematical analysis of the dynamics of prion proliferation. J Theor Biol 242: 598–606.
- 22. Rubenstein R, Gray PC, Cleland TJ, Piltch MS, Hlavacek WS, et al. (2007) Dynamics of the nucleated polymerization model of prion replication. Biophys Chem 125: 360–367.
- 23. Govaerts C, Wille H, Prusiner SB, Cohen FE (2004) Evidence for assembly of prions with lefthanded beta-helices into trimers. Proc Natl Acad Sci U S A 101: 8342–8347.
- 24. Nakamura HK, Takano M, Kuwata K (2007) Modeling of a propagation mechanism of infectious prion protein; a hexamer as the minimum infectious unit. Biochem Biophys Res Commun 361: 789–793.
- 25. Cox DL, Sing RR, Yang S (2006) Prion disease: exponential growth requires membrane binding. Biophys J 90: 77–79.
- 26. Chesebro B, Trifilo M, Race R, Meade-White K, Teng C, et al. (2005) Anchorless prion protein results in infectious amyloid disease without clinical scrapie. Science 308: 1435–1439.
- 27. Osherovich LZ, Cox BS, Tuite MF, Weissman JS (2004) Dissection and design of yeast prions. PLoS Biol 2: e86. doi:10.1371/journal.pbio.0020086.
- 28. Toyama BH, Kelly MJ, Gross JD, Weissman JS (2007) The structural basis of yeast prion strain variants. Nature 449: 233–237.
- 29. Collins SR, Douglass A, Vale RD, Weissman JS (2004) Mechanism of prion propagation: amyloid growth occurs by monomer addition. PLoS Biol 2: e321. doi:10.1371/journal.pbio.0020321.
- 30. Stohr J, Weinmann N, Wille H, Kaimann T, Nagel-Steger L, et al. (2008) Mechanisms of prion protein assembly into amyloid. Proc Natl Acad Sci U S A 105: 2409–2414.
- 31. Jones EM, Surewicz WK (2005) Fibril conformation as the basis of species- and strain-dependent seeding specificity of mammalian prion amyloids. Cell 121: 63–72.
- 32. Silveira JR, Raymond GJ, Hughson AG, Race RE, Sim VL, et al. (2005) The most infectious prion protein particles. Nature 437: 257–261.
- 33. Bucciantini M, Giannoni E, Chiti F, Baroni F, Formigli L, et al. (2002) Inherent toxicity of aggregates implies a common mechanism for protein misfolding diseases. Nature 416: 507–511.
- 34. Somerville RA, Carp RI (1983) Altered scrapie infectivity estimates by titration and incubation period in the presence of detergents. J Gen Virol 64(Pt 9): 2045–2050.
- 35. Castilla J, Gutiérrez Adán A, Brun A, Pintado B, Ramírez MA, et al. (2003) Early detection of prpres in bse-infected bovine prp transgenic mice. Arch Virol 148: 677–691.
- 36. Prusiner SB, Scott M, Foster D, Pan KM, Groth D, et al. (1990) Transgenetic studies implicate interactions between homologous prp isoforms in scrapie prion replication. Cell 63: 673–686.
- 37. Prusiner SB, Groth DF, Cochran SP, McKinley MP, Masiarz FR (1980) Gel electrophoresis and glass permeation chromatography of the hamster scrapie agent after enzymatic digestion and detergent extraction. Biochemistry 19: 4892–4898.
- 38. Thackray AM, Klein MA, Aguzzi A, Bujdoso R (2002) Chronic subclinical prion disease induced by low-dose inoculum. J Virol 76: 2510–2517.
- 39. Manuelidis L, Sklaviadis T, Manuelidis EE (1987) Evidence suggesting that prp is not the infectious agent in Creutzfeldt-Jakob disease. EMBO J 6: 341–347.
- 40. Lau AL, Yam AY, Michelitsch MM, Wang X, Gao C, et al. (2007) Characterization of prion protein (prp)-derived peptides that discriminate full-length prpsc from prpc. Proc Natl Acad Sci U S A 104: 11551–11556.
- 41. Prusiner SB, Groth D, Serban A, Koehler R, Foster D, et al. (1993) Ablation of the prion protein (prp) gene in mice prevents scrapie and facilitates production of anti-prp antibodies. Proc Natl Acad Sci U S A 90: 10608–10612.
- 42. Sakaguchi S, Katamine S, Shigematsu K, Nakatani A, Moriuchi R, et al. (1995) Accumulation of proteinase k-resistant prion protein (prp) is restricted by the expression level of normal prp in mice inoculated with a mouse-adapted strain of the Creutzfeldt-Jakob disease agent. J Virol 69: 7586–7592.
- 43. Taylor DM, Fernie K, Steele PJ, McConnell I, Somerville RA (2002) Thermostability of mousepassaged bse and scrapie is independent of host prp genotype: implications for the nature of the causal agents. J Gen Virol 83: 3199–3204.
- 44. Simoneau S, Rezaei H, Salès N, Kaiser-Schulz G, Lefebvre-Roque M, et al. (2007) In vitro and in vivo neurotoxicity of prion protein oligomers. PLoS Pathog 3: e125. doi:10.1371/journal.ppat.0030125.