- Research article
- Open Access
- Published:

# Nucleation of protein aggregation kinetics as a basis for genotype-phenotype correlations in polyglutamine diseases

*Molecular Neurodegeneration*
**volume 4**, Article number: 29 (2009)

## Abstract

Recent studies of inherited neurodegenerative disorders have suggested a linkage between the propensity toward aggregation of mutant protein and disease onset. This is particularly apparent for polyglutamine (polyQ) diseases caused by expansion of CAG-trinucleotide repeats. However, a quantitative framework for relating aggregation kinetics with molecular mechanisms of neurodegeneration initiation is lacking. Here, using the repeat-length-dependent age-of-onset in polyQ diseases, we derived a mathematical model based on nucleation of aggregation kinetics to describe genotype-phenotype correlations, and validated the model using both *in vitro* data and clinical data. Instead of describing polyQ aggregation kinetics with a derivative equation, our model divided age-of-onset (equivalent to the time required for aggregation) into two processes: nucleation lag time (a first-order exponential function of CAG-repeat length) and elongation time. With the exception of spinocerebellar ataxia (SCA) 3, the relation between CAG-repeat length and age-of-onset in all examined polyQ diseases, including Huntington's disease, dentatorubral-pallidoluysian atrophy and SCA1, -2, -6 and -7, could be well explained by three parameters derived from linear regression analysis based on the nucleated growth polymerization model. These parameters composed of probability of nucleation at an individual repeat, a protein concentration associated factor, and elongation time predict the overall features of neurodegeneration initiation, including constant risk for cell death, toxic polyQ species, main pathological subcellular site and the contribution of cellular factors. Our model also presents an alternative therapeutic strategy according to the distinct subcellular loci by the finding that nuclear localization of soluble mutant protein monomers itself has great impact on disease onset.

## Background

Accumulation of misfolded proteins into protein aggregates is a hallmark of various aging-associated neurodegenerative diseases, including Alzheimer's, Parkinson's and polyglutamine (polyQ) diseases [1, 2]. The biochemical properties of the affected proteins dictate their propensity to aggregate as well as the age-of-onset of these diseases. To date, nine inherited neurodegenerative disorders known as polyQ disease have been identified, including Huntington's disease (HD), spinal and bulbar muscular atrophy (SBMA), dentatorubral-pallidoluysian atrophy (DRPLA) and spinocerebellar ataxia (SCA) 1–3, -6, -7 and -17 [3]. These diseases have little in common at the genetic level other than the presence of polyQ sequence. However, there is a strong and consistent inverse relation between the length of the expansion and the age-of-onset in these disorders [3, 4]. Despite the fact that an understanding of genotype-phenotype relationships might offer insights into the intrinsic toxicity of polyQ peptides and the contribution of tissue context factors, such correlational analyses have rarely been attempted for the various polyQ diseases [4–6]. Those studies that have been performed have provided limited information because the parameters derived from a simple exponential regression analysis are highly variable. Rather than providing evidence that the genotype-phenotype relationship is modified by the nature of the protein encoded by each disease gene, this variability suggests deficiencies in the basic model describing the relationship.

Another common feature of polyQ diseases is the neuronal accumulation of the mutant protein in nuclear or cytoplasmic inclusion [2, 7]. In addition to the length dependence of disease onset, the length of polyQ sequence also predicted the propensity toward aggregation of polyQ-containing peptides [8, 9]. The aggregation of polyQ peptides *in vitro* follows a simple nucleated growth polymerization pathway, implying crystallization or, in some cases, amyloid fibril formation [9, 10]. Nucleated growth polymerization is a two-stage process consisting of the energetically unfavorable formation of a nucleus (i.e., nucleation), followed by efficient elongation of the nucleus via sequential additions of monomers [10]. Its kinetics is exemplified by long lag time followed by rapid aggregate growth, with a strong dependence of aggregation lag time on monomer concentration [10]. Nucleated growth polymerization has been proposed to govern disease progression kinetics in Alzheimer's and prion-related diseases [11]. We and others have previously suggested a linkage between the biophysics of polyQ aggregation nucleation and HD onset [9, 12]. The actual mechanism of the generation of nuclei based on polyQ sequences will be structurally complex, but a kinetic parameter of nucleation is expected to be an exponential function of repeat length [13]. The polyQ length dependence of disease onset correlates strongly with the tendency of expanded polyQ proteins to aggregate in disease models [14, 15]. Accordingly, we have focused on this length dependence of age-of-onset and nucleation kinetics to derive a stochastic mathematical model describing genotype-phenotype correlations in polyQ diseases.

Accumulating evidence strongly suggests that the cell nucleus is the main pathological subcellular site for SCA1, -7 and HD [16–18], whereas the cytoplasm is thought to be the site for SCA2 and SCA6 [19, 20]. Our mathematical model clearly subdivided polyQ diseases into two groups in accordance with the presumed main pathological subcellular site. From a comparison of the parameters based on the nucleated growth polymerization model, our present study leads us to propose alternative therapeutic targets according to the distinct subcellular loci.

## Methods

### Data collection

Clinical data from affected patients with a mutation in the responsible gene were obtained from previous studies [5, 6, 21–37]. Patients with a homozygous mutation were excluded. An insufficient number of patients with SBMA and SCA17 were available, so these diseases were not included in the analysis. A total of 1398 patients were analyzed: 221 with SCA1, 141 with SCA2, 308 with SCA3, 132 with SCA6, 188 with SCA7, 317 with HD and 91 with DRPLA.

### Modeling

The derivation of equation describing polyQ peptide aggregation kinetics in a nucleated growth polymerization mechanism has been described by Chen *et al*. [9] as follows:

where Δ is the concentration of monomer that has been incorporated into polymers, *k*_{+} is the forward elongation rate constant, *k*_{n*} is the equilibrium constant describing the monomer-nucleus equilibrium, *c* is the bulk concentration of monomers, *n** is the critical nucleus, and *t* is time. This equation represents the overall pathway of nucleated growth polymerization. Further, kinetic studies also suggest that, in this mathematical model, the only factor related to CAG-repeat length is the nucleation constant, *k*_{n*}[9, 38, 39]. However, nucleation kinetics cannot be determined directly through physical measurement of nuclei because nucleation is a very rare event and nuclei, once formed, quickly either collapse to bulk phase monomer or proceed along the productive aggregation pathway. Aggregation occurs with a lag phase and a growth phase that reflect an underlying nucleation-polymerization mechanism. The lag time of aggregation (aggregation lag time) as a kinetic parameter of nucleation was calculated by extrapolation of the linear region of the growth phase to the base line of the lag phase by experimental observation using Equation 1 (Figure 1A). However, polyQ aggregation kinetics also feature lag phases that can be abbreviated by seeding [8, 40]. Instead of using polyQ aggregation kinetics to describe the overall pathway of nucleated growth polymerization, we hypothesized a mathematical model that divides the time required for aggregation into two processes: CAG-repeat-length-dependent nucleation lag time, and repeat-length-independent elongation time from elongation kinetics using aggregates of the polyQ peptides as a seed [8, 38, 39]. In this mathematical model, nucleation lag time (a proxy for aggregation lag time) is determined by the time required for aggregation and elongation time, as a function of repeat length.

In the nucleated growth polymerization pathway, aggregate growth is exponential. A previous study that employed a mutagenic analysis to the nucleation propensity of Alzheimer's β-amyloid showed that nucleation lag time of mutant β-amyloid could be related directly to the rate of aggregate growth [41]. Aggregate growth is largely dependent on nucleation lag time. In polyQ diseases, by contrast, the relationship between age-of-onset and the number of repeats is typically characterized by an exponential curve in which the change in the age-of-onset as a result of additional inherited repeats decreases with the number of repeats [4]. Moreover, as noted by Clarke *et al*., in clinically affected HD patients there is a first-order exponential decline in the number of surviving caudate nucleus neurons with time (based on metabolic activity), regardless of CAG-repeat length [42]. These findings, taken together with the linkage between aggregation nucleation of polyQ peptides and HD age-of-onset [9, 12], suggest that the nucleation lag time for polyQ aggregation could be represented by a first-order exponential function of CAG-repeat lengths. *In vitro* studies also showed that the elongation rates of polyQ peptides changed little as a function of repeat-length [8, 38, 39]. Assuming such an identity in the elongation rate constant, and given the fact that the aggregation rate is *t*^{2} dependent [9, 39], the age-of-onset in polyQ diseases could be represented as a function of CAG-repeat lengths as follows:

where *t*_{A} is age-of-onset of the disease (equivalent to the time required for aggregation), *t*_{E} is elongation time, *t*_{N} is nucleation lag time and *x* is the repeat-length number.

The identical relationship between nucleation lag time and CAG-repeat length was also proposed by Perutz et al. based on chemicalthermodynamics [13]. According to the theory of nucleation of aggregates, the probability of nucleation would be an exponential function of the free energy of formation of the nucleus, and is proportional to exp [-Δ*G*_{crit}/*kT*] [13]. Δ*G*_{crit} is the critical free energy required to create a spherical nucleus, *k* is Boltzmann's constant, and *T* is absolute temperature. Since the addition of each glutamine stabilizes the helix structure by the formation of another three or four hydrogen bonds [43], alternation of free energy by the addition of each repeat to be expected constant. Therefore, probability of nucleation is likely to rise exponentially with the number of repeats. *P*_{n+1}/*P*_{
n
}= exp [Δ*G*_{+1}/*kT*], where *P*_{n+1 }is the probability of nucleation at repeat number *n* + 1, *P*_{
n
}is the probability of nucleation at repeat number n, Δ*G*_{+1} is the change of free energy by one additional repeat. Nucleation lag time is defined as that lag time which is required for the formation of a critical number of stable nuclei leading to polymerization. *In vitro* aggregation kinetics of polyQ peptide also showed that the critical nucleus – the number of monomeric units comprising the nucleus – is equal to 1 [9]. Thus, under the condition of monomeric nucleus, nucleation lag time could be represented by an exponential function of repeat length as follows:

where *N*_{cri} is a critical number of nucleus leading to polymerization at a given space, *N*_{A} is Avogadro's number, *C*_{mo} is the bulk concentration of mutant protein monomers, *R*_{0} is the rate of formation of stable nuclei at the first structure, *m* is the number of repeat length at the first structure stable enough to form a nucleus, and e^{b}is the probability of nucleation at an individual repeat.

### Statistical analysis

For statistical analysis, we used the UNISTAT 5.6 statistical package for Windows (Unistat). To examine the association between age-of-onset and CAG size, we used linear regression with logarithmic transformation, allowing an exponential function to be treated as an intrinsically linear model as follows:

where In(*a*) is the intercept, *b* is the slope, and *ε* represents the residual error. Of course, the slope factor represents the probability of nucleation at an individual repeat and is determined by a nucleation constant, *k*_{n*}, because *k*_{n*} = *c**/*c*^{n}* (*c** is the concentration of nuclei, *c* is the bulk concentration of monomers, *n** is the critical nucleus) [10]. The intercept is composed of the complex factors (Equation 3). However, it is expected that its variation is mainly dependent on the concentration of soluble mutant protein monomers. A linear regression analysis was then applied to determine *t*_{E}^{2} for the best fit to a linear relationship between the logarithm of *t*_{N} and CAG size. This was achieved by comparing the *R*^{2} for a quadratic curve and a linear model. Then, to confirm that our mathematical model is a significantly better fit than a simple exponential regression model, we compared the *R*^{2} values for linear regression analyses derived from simple natural log-transformed age-of-onset data without *t*_{E}^{2}. We used the residual error as another measure of goodness-of-fit. The small number of individuals with shortest and longest CAG-repeat sizes precluded a rigorous statistical analysis of these groups because they would have overruled the bulk of the other data and thus had too great an impact on the model. Applying this procedure, 94% to 98% of patients in each disease category were included in the analysis.

## Results

### Validated the model using *in vitro*data

In previous study, using the *in vitro* aggregation lag time of a series of polyQ peptides as a function of CAG size, described by Chen *et al*. [9], we compared the relationship between median age-at-onset of HD and particular repeat sizes, demonstrating a significant linear correlation [12]. Instead of using polyQ peptide aggregation kinetics to describe the overall pathway of nucleated growth polymerization [9], we successfully developed a mathematical model that divided the time required for aggregation (*t*_{A}) into two processes governed by nucleation lag time (*t*_{N}), a first-order exponential function of CAG size, and elongation time (*t*_{E}), to yield the following relationship: (*t*_{A}^{2} - *t*_{E}^{2})^{1/2} = *t*_{N} = *a*e^{-bx}(see Methods). Then, we tested the linear correlation with logarithmic transformation, allowing an exponential function to be treated as an intrinsically linear model. A simple logarithmic transformation of the *in vitro* aggregation lag times against the polyQ tract number did not show a linear relationship (data not shown). These data represents the time required for 0.9% aggregation of polyQ peptides [9]. In contrast, when *t*_{E}^{2} = 27, the application of our model to *in vitro* nucleation lag times versus CAG size yields a plot that is perfectly fit with an exponential function (Figure 1A).

### Validated the model using clinical data

Based on the nucleated growth polymerization model of polyQ peptides aggregation, we used the model for describing the relationship between CAG-repeat length and age-of-onset in polyQ diseases as follows: In(*t*_{A}^{2} - *t*_{E}^{2})^{1/2} = In(*a*) - *bx* + *ε* (see Methods), where *t*_{A} is age-of-onset of the disease (equivalent to the time required for aggregation), *ε* represents the residual error and *x* is the repeat-length number. We subsequently validated the model using available clinical data from patients with HD, DRPLA and SCA1–3, -6 and -7 [5, 6, 21–37]. A recent study of the relationship between age-of-onset in HD and CAG size in the HD gene has shown that a two-segment exponential regression curve provides a significantly better fit than a simple exponential regression model [44]. A similar tendency was also observed for SCA1, -2, -6, and -7. In fact, in these disorders, a linear regression analysis of natural log-transformed age-of-onset versus CAG size revealed that plots of residual error versus CAG size were consistently better fit to a U-shaped curve (Figure 1B). These findings strongly support our model, since a regression analysis of natural log-transformed age-of-onset, *t*_{A} = (*t*_{N}^{2} + *t*_{E}^{2})^{1/2}, versus CAG size is better fit to a U-shaped curve than a simple linear model, except for the case where elongation time is 0 (*t*_{E} = 0). Therefore, elongation time is determined for each polyQ disease by identifying when the coefficients of determination (*R*^{2}) for a quadratic curve and a linear model are identical (Figures 2, 3 and 4). This procedure restores the uncorrelation of the residual error, and consistently obtained higher *R*^{2} values than linear regression analyses derived from simple natural log-transformed age-of-onset data without *t*_{E}^{2} in demonstrating a significant better fit than a simple exponential regression model (Table 1). We also analyzed the relationship using a logarithmic transformation of (*t*_{A} - *t*) against CAG size, corresponding to an off-pathway of nucleated growth polymerization. Then, *t* is determined for each polyQ disease by identifying when the *R*^{2} for a quadratic curve and linear model are identical (data not shown). However, a higher *R*^{2} value was consistently obtained with (*t*_{A}^{2} - *t*_{E}^{2})^{1/2} versus CAG size, indicating an on-pathway of nucleated growth polymerization.

In contrast to the linear regression analysis of genotype-phenotype correlation in other polyQ diseases we examined, SCA3 clearly did not conform to the nucleated growth polymerization model. A natural log-transformed age-of-onset in a total of 308 SCA3 patients versus CAG size in the ataxin-3 gene was the best fit to a negative quadratic function (Figure 5). This suggests that, despite shared pathological findings, including nuclear inclusion, a different mechanism underlies SCA3 pathogenesis.

The validity of the model is also supported by clinical data on disease progression. A statistical analysis of the detailed information provided by activity of daily living (ADL) milestones in SBMA patients and their relationship to the CAG-repeat length of the androgen receptor gene has demonstrated that the ADL milestone is significantly related to CAG size [45]. In contrast, the rate of disease progression showed no significant correlation to CAG size [45]. A similar lack of correlation between the rate of disease progression and CAG size was observed for HD patients [46]. Interestingly, the slopes of the regression curves for age-at-onset of each ADL milestone versus CAG size likely parallel one another [45], suggesting that the repeat-length-dependent parameter is conserved with respect to disease progression. In our model, the slope (*b*) factor is only dependent on CAG-repeat length, which is determined by a nucleation constant, *k*_{n*} (see Methods). Therefore, the constancy of the slope could account for the fact that CAG-repeat length has little effect on the rate of disease progression, but is related to age-at-onset for each ADL milestone. In fact, the relationship between the age-at-onset of each ADL milestone in SBMA patients and CAG-repeat length also conformed to the nucleated growth polymerization model (data not shown).

### Linear regression analysis of HD, DRPLA, SCA1, -2, -6 and -7 based on the model

HD: A total of 317 HD patients with expanded CAG repeats in the huntingtin gene was analyzed. A regression analysis of natural log-transformed (*t*_{A}^{2} - *t*_{E}^{2})^{1/2} against CAG size provided the best fit to a linear model (-0.0662*x* + 6.657) when *t*_{E}^{2} = 215 (Figure 2A). DRPLA: A total of 91 DRPLA patients with expanded CAG repeats in the atrophin-1 gene was analyzed. A regression analysis of natural log-transformed (*t*_{A}^{2} - *t*_{E}^{2})^{1/2} against CAG size provided the best fit to a linear model (-0.1247*x* + 11.011) when *t*_{E}^{2} = 0 (Figure 2B). SCA1: A total of 221 SCA1 patients with expanded CAG repeats in the ataxin-1 gene was analyzed. A regression analysis of natural log-transformed (*t*_{A}^{2} - *t*_{E}^{2})^{1/2} against CAG size provided the best fit to a linear model (-0.0569*x* + 6.249) when *t*_{E}^{2} = 170 (Figure 3A). SCA2: A total of 141 SCA2 patients with expanded CAG repeats in the ataxin-2 gene was analyzed. A regression analysis of natural log-transformed (*t*_{A}^{2} - *t*_{E}^{2})^{1/2} against CAG size provided the best fit to a linear model (-0.1217*x* + 8.377) when *t*_{E}^{2} = 80 (Figure 3B). SCA6: A total of 132 SCA6 patients with expanded CAG repeats in the alpha-1A-calcium channel gene was analyzed. A regression analysis of natural log-transformed (*t*_{A}^{2} - *t*_{E}^{2})^{1/2} against CAG size provided the best fit to a linear model (-0.2567*x* + 9.553) when *t*_{E}^{2} = 700 (Figure 4A). SCA7: A total of 188 SCA7 patients with expanded CAG repeats in the ataxin-7 gene was analyzed. A regression analysis of natural log-transformed (*t*_{A}^{2} - *t*_{E}^{2})^{1/2} against CAG size provided the best fit to a linear model (-0.0658*x* + 6.513) when *t*_{E}^{2} = 30 (Figure 4B).

### Parameters derived from linear regression analysis

Three parameters, represented by the slope and intercept of the natural log-transformed linear regression curve, and elongation time, as well as descriptive statistics for each disease are summarized in Table 1. From Equations 3 (see Methods), the slope factor of the regression curve, e^{b}, represents the probability of nucleation at an individual repeat (governed by the nucleation constant, *k*_{n*}), and the intercept (In(*a*)) is interpreted as a variable that is inversely related to the concentration of soluble mutant protein monomers.

In genotype-phenotype correlations, each polyQ disease shows a characteristic threshold of CAG-repeat length [3, 4]. SCA6 arises from a relatively small expansion with as few as 21 repeats. This contrasts with SCA1, -2, -7 and HD, where 35–40 repeats cause disease, and DRPLA, which arises from an even larger expansion (> 50 repeats). Three parameters derived from statistical analysis can well explain the relationship between CAG-repeat length and age-of-onset in any of these disorders. For example, these parameters indicate that a relatively large expansion required for disease onset in DRPLA is mainly due to the lower concentration of soluble mutant protein monomers (Figure 2B). In addition, a comparison of these parameters should reveal which factors strongly influence disease onset. In polyQ diseases, symptoms typically appear in adulthood; in DRPLA, infant-onset of the disease is occasionally seen, but this is never the case for SCA6. Our model clearly pointed to elongation time as a contributor to infant-onset, while adult-onset was mostly attributable to the nucleation lag time (Table 1).

### Variance in nucleation lag times

It is now believed that the identification of non-CAG size-dependent factors that explain residual onset age variance would possibly allow treatments that retard disease onset of polyQ disease. Regression analysis of HD, DRPLA and SCA1, -2, -6, and -7 showed that 57 to 68% variance of nucleation lag times is explained by the number of repeat units on the mutant allele (Table 1). However, using median nucleation lag time with a particular repeat size, there was a high significant correlation (*R* = 0.96 ~ 0.99) between CAG size and median nucleation lag times (data not shown), implying that the residual variance in nucleation lag times is mainly due to the variance at a particular repeat size. Our mathematical model (Equation 3) suggests that the concentration of soluble mutant protein monomers accounts for a major contributor to the residual variance in age-of-onset. Importantly, nucleation is a very rare event [13]. Therefore, only a small difference in concentration of soluble mutant protein monomers could result in a large difference in nucleation lag times, which directly reflects to onset age variance. In addition to the expression levels of a gene, the complexity of the cellular environment, including degradation and transport processes capable of partitioning proteins into different molecular forms and compartments, and the presence of chaperones will contribute to the differences in intracellular concentration of mutant protein.

## Discussion

Three parameters derived from linear regression analysis based on the nucleated growth polymerization model predict the overall features of neurodegeneration initiation of polyQ diseases, including toxic polyQ species, constant risk for cell death, main pathological subcellular site and the contribution of cellular factors, and provide an explanation for the aggregation of relatively short expanded polyQ tracts in SCA6.

Unexpectedly, a linear regression analysis of each disease showed that *t*_{E}^{2} varied from 0 to 700 (Table 1). This result, which is in contrast to the findings of a previous study [47], suggests that the specific oligomeric conformation of expanded polyQ is not uniformly toxic to neuronal cells in polyQ diseases. Instead, the validity of the nucleated growth polymerization model in any of these diseases indicates that toxic polyQ species accrues uniformly on-pathway of nucleated growth polymerization. Even though the disease-causing polyQ proteins are widely expressed, specific collections of neurons are more susceptible in each polyQ disease, resulting in characteristic patterns of pathology and clinical symptoms. Recent studies have suggested that altered protein function is fundamental to pathogenesis, with protein context of the expanded polyQ having key roles in disease specific processes [48]. Therefore, varied *t*_{E} suggest that a growth phase of aggregation may represent disease specific processes affected by protein context of the expanded polyQ.

The clinical manifestations of inherited neurodegenerative diseases are often delayed for periods from years to decades. This observation has led to the idea that, in these disorders, neurons die from cumulative damage. Consistent with a cumulative damage model, a failure of protein folding quality control has been proposed for the pathogenesis of polyQ diseases [49, 50]. A critical prediction of the cumulative damage hypothesis is that the probability of neuronal death increase over time (sigmoidal kinetics). However, Clarke *et al*. demonstrated that, in many aging-associated neurodegenerative diseases, including HD, the kinetics of neuronal death appear to be exponential [42, 51]. Exponential kinetics, which also describe radioactive decay, indicates that in the neurodegenerations, the risk of cell death is constant (in some cases, as described by Clarke *et al*., an exponentially decreasing risk of death) over time. They accounted for this constancy by proposing a one-hit model in which the death of a neuron is initiated randomly in time by a single, rare catastrophic event rather than resulting from cumulative damage [42, 51]. Thereafter, Perutz *et al*. proposed that nucleation is responsible for this kind of rare event, based on the fact that nucleation occurs randomly in time and space [13]. We have successfully established the common quantitative connection among polyQ diseases between the nucleation kinetics based on polyQ sequence and the repeat-length-dependent age-of-onset, implicating nucleation of protein aggregation kinetics as the basis for the genotype-phenotype correlations. In our model, the slope of the regression curve, which is a constant factor related only to CAG-repeat length, represents the probability of nucleation at an individual repeat and is equivalent to the constant risk of neuronal cells death at an individual repeat, consistent with Perutz's hypothesis [13]. Remarkably, nucleation and nucleation dependent polymerization exhibits kinetics consistent with both a constant and an exponentially decreasing risk of neuronal death in the neurodegenerations. Taken together, these findings provide indirect evidence that the genetic gain-of-function mechanism of polyQ pathogenesis is attributable to a critical nucleation event; thus, cytotoxicity accrues through the process of nucleation and nucleation-dependent polymerization.

In contrast to other polyQ diseases we examined, SCA3 does not fit our derived model (Figure 5), suggesting that another factor, dependent on the repeat-length of polyQ expansions, plays a critical role in determining disease onset of SCA3. Recent studies have suggested that ataxin-3 normally participates in protein quality control pathways in the cell [52, 53]. Endoplasmic reticulum (ER)-associated degradation (ERAD) is a quality control system in the secretory pathway responsible for degrading misfolded proteins [54]. ERAD involves a series of steps to extract proteins from the ER and deliver them to proteasomes. The key protein essential for extracting substrates from the ER to the cytosol is valosin-containing protein [55]. It is of interest that ataxin-3 binds valosin-containing protein and regulates retrotranslocation of ERAD substrates in a repeat-length-dependent manner [53]. These findings underscore the critical role of protein quality control in SCA3 pathogenesis and raise a possibility that neuronal death in SCA3 exhibits sigmoidal kinetics.

Further study of disease progression as a function of CAG-repeat length for each polyQ disease will be needed to confirm the constancy of the slope and toxicity relationships. These analyses will more precisely characterize rates of disease progression whether it consists with exponential kinetics or sigmoidal kinetics.

One of the parameters, the slope of the regression curve, clearly subdivided polyQ diseases into two groups. One group included SCA1, -7 and HD; the other included SCA2, -6 and DRPLA (Table 1). Accumulating evidence strongly suggests that the cell nucleus is the main pathological subcellular site for SCA1, -7 and HD [16–18], whereas the cytoplasm is thought to be the site for SCA2 and SCA6 [19, 20]. Clearly, the distribution into two groups is consistent with the presumed main pathological subcellular site in each disease. Further support for this idea is provided by the observation that a smaller intercept is consistently observed in the case of SCA1, -7 and HD (Table 1). Given the nucleus-to-cytoplasm volume ratio, it is expected that the concentration of degradation-resistant mutant proteins would be higher in the nucleus than in the cytoplasm. Another explanation for the smaller intercept in SCA1, -7 and HD might be related to the finding that the autophagy-lysosomal system is also involved in cytoplasmic degradation of aggregate-prone protein with polyQ, but the nucleus lacked this activity [56].

The slope of *in vitro* aggregation of polyQ peptides is 0.1512 (Figure 1A), which is close to the slope for SCA2 (0.1217). In contrast to conventional models of nucleated growth polymerization, the critical nucleus (the number of monomeric unit comprising the nucleus) for polyQ peptides aggregation is a monomer (*n* = 1) [9]. Assuming that the cytoplasmic pathological sites are identical, compared to the slope of SCA2, the slope of SCA6 (0.2567) corresponds to a dimeric nucleus (*n* = 2). To date, there are no data capable of providing a mechanism to explain the aggregation of such a short polyQ expansion in SCA6. Our model predicts that the aggregation in SCA6 is attributable to the dimeric nucleus, which promotes disease-onset about 20 times more effectively than a monomeric nucleus at a CAG size of 22. In contrast, the slopes of SCA1, -7 and HD are substantially lower than the value of 0.1512 (Table 1), which corresponds to the smallest critical nucleus size. Therefore, a certain inhibitory factor such as molecular chaperone or ubiquitin-proteasome system is inferred to exist in the cell nucleus that reduces the probability of nucleation to less than half the value of the slope. This is compatible with the idea that nuclear inclusion serves to protect against neurodegeneration in polyQ diseases [57]. Taken together with the finding that SCA3 operates by a different mechanism, this suggests that nuclear inclusion itself is an off-pathway product of nucleated growth polymerization.

Based on the nucleated growth polymerization model, the parameters derived from linear regression analysis have led us to propose alternative therapeutic targets according to the distinct subcellular loci (Figure 6). The slope (probability of nucleation at an individual repeat) and intercept (inversely related to the concentration of soluble mutant protein monomers) determines a critical number of the nucleus, which initiates polymerization pathway, considered as a threshold of polyQ cytotoxicity. Our results suggest that, in the cell nucleus, host defense factors considerably inhibit the probability of nucleation at an individual mutant protein. By contrast, it is estimated that the concentration of soluble mutant protein monomers in the cell nucleus is quite higher than the values in the cytoplasm (e^{8.377~11.011} versus e^{6.249~6.657}, the mean ratio is about 40:1). Therefore, nuclear localization of the mutant protein itself has great impact on the age-of-onset in SCA1, -7 and HD. Theoretically, an attractive therapeutic target is to reduce the concentration of soluble mutant protein monomers in the cell nucleus, particularly by inhibition of translocation of the mutant proteins from the cytoplasm to the cell nucleus, and thereby delay polyQ nucleation aggregation.

## Conclusion

One of the striking findings of neurodegeneration research is the observation that most of the proteins implicated in disease have a strong propensity to aggregate. Aggregation is a central aspect of the biology of many neurodegenerative diseases. However, the role of aggregates in neurodegeneration is unclear. In this article, we present the first development of a mathematical model, which describes the basis for genotype-phenotype correlations in inherited neurodegenerative disorders known as polyQ disease and have successfully established the common quantitative connection among polyQ diseases between the repeat-length-dependent age-of-onset and aggregation kinetics based on polyQ sequence. Our results have clear implications for polyQ disease pathogenesis and therapy.

## References

- 1.
Chiti F, Dobson CM: Protein misfolding, functional amyloid, and human disease. Annu Rev Biochem. 2006, 75: 333-366.

- 2.
Ross CA, Poirier MA: Protein aggregation and neurodegenerative disease. Nat Med. 2004, 10: S10-17.

- 3.
Gatchel JR, Zoghbi HY: Disease of unstable repeatexpansion: mechanisms and common principles. Nat Rev Genet. 2005, 6: 743-755.

- 4.
Gusella JF, MacDonald ME: Molecular genetics: unmasking polyglutamine triggers in neurodegenerative disease. Nat Rev Neurosci. 2000, 1: 109-115.

- 5.
van de Warrenburg BP, Sinke RJ, Verschuuren-Bemelmans CC, Scheffer H, Brunt ER, Ippel PF, Maat-Kievit JA, Dooijes D, Notermans NC, Lindhout D, Knoers NV, Kremer HP: Spinocerebellar ataxia in the Netherlands: prevalence and age at onset variance analysis. Neurology. 2002, 58: 702-708.

- 6.
Warrenburg van de BP, Hendriks H, Durr A, van Zuijelen MC, Stevanin G, Camuzat A, Sinke RJ, Brice A, Kremer BPH: Age at onset variance analysis in spinocerebellar ataxias: a study in Dutch-French cohort. Ann Neurol. 2005, 57: 505-512.

- 7.
Michalik A, van Broeckhoven C: Pathogenesis of polyglutamine disorders: aggregation revisited. Hum Mol Genet. 2003, 12: R173-R186.

- 8.
Chen S, Berthelier V, Yang W, Wetzel R: Polyglutamine aggregation behavior in vitro supports a recruitment mechanism of cytotoxicity. J Mol Biol. 2001, 311: 173-182.

- 9.
Chen S, Ferrone FA, Wetzel R: Huntington's disease age-of-onset linked to polyglutamine aggregation nucleation. Proc Natl Acad Sci USA. 2002, 99: 11884-11889.

- 10.
Ferrone F: Analysis of protein aggregation kinetics. Methods Enzymol. 1999, 309: 256-274.

- 11.
Jarrelt JT, Lansbury PT: Seeding "one-dimensional crystallization" of amyloid: a pathogenic mechanism in Alzheimer's disease and scrapie?. Cell. 1993, 73: 1055-1058.

- 12.
Sugaya K, Matsubara S, Kagamihara Y, Kawata A, Hayashi H: Polyglutamine expansion mutation yield a pathological epitope linked to nucleation of protein aggregate: determinant of Huntington's disease onset. PLoS ONE. 2007, 2: e635-

- 13.
Perutz MF, Windle AH: Cause of neural death in neurodegenerative diseases attributable to expansion of glutamine repeats. Nature. 2001, 412: 143-144.

- 14.
Scherzinger E, Lurz R, Turmaine M, Mangiarini L, Hollenbach B, Hasenbank R, Bates GP, Davies SW, Lehrach H, Wanker EE: Huntingtin-encoded polyglutamine expansions form amyloid-like protein aggregates in vitro and in vivo. Cell. 1997, 90: 549-558.

- 15.
Jackson GR, Salecker I, Dong X, Yao X, Arnheim N, Faber PW, MacDonald ME, Zipursky SL: Polyglutamine-expanded human huntingtin transgenes induce degeneration of drosophila photoreceptor neuron. Neuron. 1998, 21: 633-642.

- 16.
Saudou F, Finkbeiner S, Devys D, Greenberg ME: Huntingtin acts in the nucleus to induce apoptosis but death does not correlate with the formation of intranuclear inclusion. Cell. 1998, 95: 55-66.

- 17.
Orr HT, Zoghbi HY: SCA1 molecular genetics: a history of a 13 year collaboration against glutamines. Hum Mol Genet. 2001, 10: 2307-2311.

- 18.
La Spada AR, Fu YH, Sopher BI, Libby RT, Wang X, Li LY, Einum DD, Huang J, Possin DE, Smith AC, Martinez RA, Koszdin KL, Treuting PM, Ware CB, Hurley JB, Ptacek LJ, Chen S: Polyglutamine-expanded ataxin-7 antagonizes CRX function and induces cone-rod dystrophy in a mouse model of SCA7. Neuron. 2001, 31: 913-927.

- 19.
Huynh DP, Figueroa K, Hoang N, Pulst SM: Nuclear localization or inclusion body formation of ataxin-2 are not necessary for SCA2 pathogenesis in mouse or human. Nat Genet. 2000, 26: 44-50.

- 20.
Watase K, Barrett CF, Miyazaki T, Ishiguro T, Ishikawa K, Hu Y, Unno T, Sun Y, Kasai S, Watanabe M, Gomez CM, Mizusawa H, Tsien RW, Zoghbi HY: Spinocerebellar ataxia type 6 knock in mice develop a progressive neuronal dysfunction with age-dependent accumulation of mutant Ca

_{v}2.1 channels. Proc Natl Acad Sci USA. 2008, 105: 11987-11992. - 21.
Andrew SE, Goldberg YP, Kremer B, Telenius H, Theilmann J, Adam S, Starr E, Squitieri F, Lin B, Kalchman MA, Graham RK, Hayden MR: The relationship between trinucleotide (CAG) repeat length and clinical features of Huntington's disease. Nat Genet. 1993, 4: 398-403.

- 22.
Komure O, Sano A, Nishino N, Yamauchi N, Ueno S, Kondoh K, Sano N, Takahashi M, Murayama N, Kondo I, Nagafuchi S, Yamada M, Kanazawa I: DNA analysis in hereditary dentatorubral-pallidolysian atrophy: correlation between CAG repeat length and phenotype variation and molecular basis of anticipation. Neurology. 1995, 45: 143-149.

- 23.
Ikeuchi T, Koide R, Tanaka H, Onodera O, Igarashi S, Takahashi S, Kondo R, Ishikawa A, Tomoda A, Miike T, Sato K, Ihara Y, Hayabara T, Isa F, Tanabe H, Tokiguchi S, Hayashi M, Shimizu N, Ikuta F, Naito H, Tsuji S: Dentatorubral-pallidolysian atrophy: clinical features are closely related to unstable expansions of trinucleotide (CAG) repeat. Ann Neurol. 1995, 37: 769-775.

- 24.
Ranum LP, Chung MY, Banfi S, Bryer A, Schut LJ, Ramesar R, Duvick LA, McCall A, Subramony SH, Goldfarb L, Gomez C, Sandkuijl LA, Orr HT, Zoghbi HY: Molecular and clinical correlations in spinocerebellar ataxia type 1: evidence for familial effects on the age at onset. Am J Hum Genet. 1994, 55: 244-252.

- 25.
Cancel G, Durr A, Didierjean O, Imbert G, Burk K, Lezin A, Belal S, Benomar A, Abada-Bendib M, Vial C, Guimaraes J, Chneiweiss H, Stevanin G, Yvert G, Abbas N, Saudou F, Lebre AS, Yahyaoui M, Hentati F, Vernant JC, Klockgether T, Mandel JL, Agid Y, Brice A: Molecular and clinical correlations in spinocerebellar ataxia2: a study of 32 families. Hum Mol Genet. 1997, 6: 709-715.

- 26.
Giunti P, Sabbadini G, Sweeney MG, Davis MB, Veneziano L, Mantuano E, Federico A, Plasmati P, Frontali M, Wood NW: The role of the SCA2 trinucleotide repeat expansion in 89 autosomal dominant cerebellar ataxia families. Frequency, clinical and genetic correlates. Brain. 1998, 121: 459-467.

- 27.
Maciel P, Gaspar C, DeStefano AL, Silveira I, Coutinho P, Anita L, Radvany J, Dawson DM, Sudarsky L, Guimaraes J, Loureiro JE, Nezarati MM, Corwin LI, Lopes-Cendes I, Rooke K, Rosenberg R, MacLeod P, Sequeiros LA, Rouleau GA: Correlation between CAG repeat length and clinical features in Machod-Joseph disease. Am J Hum Genet. 1995, 57: 54-61.

- 28.
Maruyama H, Nakamura S, Matsuyama Z, Sakai T, Doyu M, Sobue G, Seto M, Tsujihata M, Oh-i T, Nishio T, Sunohara N, Takahashi R, Hayashi M, Nishino I, Ohtake T, Oda T, Nishimura M, Saida T, Matsumoto H, Baba M, Kawaguchi Y, Kakizuka A, Kawakami H: Molecular features of the CAG repeats and clinical manifestation of Machado-Joseph disease. Hum Mol Genet. 1995, 4: 807-812.

- 29.
Matsumura R, Futamura N, Fujimoto Y, Yanagimoto S, Horikawa H, Suzumura A, Takayanagi T: Spinocerebellar ataxia type 6. Molecular and clinical features of Japanese patients including one homozygous for the CAG repeat expansion. Neurology. 1997, 49: 1238-1243.

- 30.
Ikeuchi T, Takano H, Koide R, Horikawa Y, Honma Y, Onishi Y, Igarashi S, Tanaka H, Nakao N, Sahashi K, Tsukagoshi H, Inoue K, Takahashi H, Tsuji S: Spinocerebellar ataxia type 6: CAG repeat expansion in α

_{1A}voltage-dependent calcium channel gene and clinical variations in Japanese population. Ann Neurol. 1997, 42: 879-884. - 31.
Geschwind DH, Perlman S, Figueroa KP, Karrim J, Baloh RW, Pulst SM: Spinocerebellar ataxia type 6. Frequency of the mutation and genotype-phenotype correlations. Neurology. 1997, 49: 1247-1251.

- 32.
Schols L, Kruger R, Amoiridis G, Przuntek H, Epplen JT, Riess O: Spinocerebellar ataxia type 6: genotype and phenotype in German kindreds. J Neurol Neurosurg Psychiatry. 1998, 64: 67-73.

- 33.
Sinke RJ, Ippel EF, Diepstraten CM, Beemer FA, Wokke JHJ, van Hilten BJ, Knoers NVAM, van Amstel HKP: Clinical and molecular correlations in spinocerebellar ataxia type 6. Arch Neurol. 2001, 58: 1839-1844.

- 34.
David G, Durr A, Stevanin G, Cancel G, Abbas N, Benomar A, Belal S, Lebre AS, Abada-Bendib M, Grid D, Holmberg M, Yahyaoui M, Hentati F, Chkili T, Agid Y, Brice A: Molecular and clinical correlations in autosomal dominant cerebellar ataxia with progressive macular dystrophy (SCA7). Hum Mol Genet. 1998, 7: 165-170.

- 35.
Gouw LG, Castaneda MA, McKenna CK, Digre KB, Pulst SM, Perlman S, Lee MS, Gomez C, Fischbeck K, Gagnon D, Storey E, Bird T, Jeri FR, Ptacek LJ: Analysis of the dynamic mutation in the SCA7 gene shows marked paternal effects on CAG repeat transmission. Hum Mol Genet. 1998, 7: 525-532.

- 36.
Johansson J, Forsgren L, Sandgren O, Brice A, Holmgren G, Holmberg M: Expanded CAG repeats in Swedish spinocerebellar ataxia type 7 (SCA7) patients: effect of CAG repeat length on the clinical manifestation. Hum Mol Genet. 1998, 7: 171-176.

- 37.
Giunti P, Stevanin G, Worth PF, Brice A, Wood NW: Molecular and clinical study of 18 families with ADCA type II: evidence for genetic heterogeneity and de novo mutation. Am J Hum Genet. 1999, 64: 1594-1603.

- 38.
Bhattacharyya AM, Thanker AK, Wetzel R: Polyglutamine aggregation nucleation: thermodynamics of a highly unfavorable protein folding reaction. Proc Natl Acad Sci USA. 2005, 102: 15400-15405.

- 39.
O'Nuallain B, Thanker AK, Williams AD, Bhattacharyya AM, Chen S, Thiagarajan G, Wetzel R: Kinetics and thermodynamics of amyloid assembly using a high-performance liquid chromatography-based sedimentation assay. Methods Enzymol. 2006, 413: 34-74.

- 40.
Thakur AK, Wetzel R: Mutational analysis of the structural organization of polyglutamine aggregates. Proc Natl Acad Sci USA. 2002, 99: 17014-17019.

- 41.
Christopeit T, Hortschansky P, Schroeckh V, Guhrs K, Zandomeneghi G, Fandrich M: Mutagenic analysis of the nucleation propensity of oxidized Alzheimer's β-amyloid peptide. Protein Sci. 2005, 14: 2125-2131.

- 42.
Clarke G, Collins RA, Leavitt BR, Andrews DF, Hayden MR, Charles J, Lumsden CJ, McInnes RR: A one-hit model of cell death in inherited neuronal degeneration. Nature. 2000, 406: 195-199.

- 43.
Perutz MF, Finch JT, Berriman J, Lesk A: Amyloid fibers are water-filled nanotubes. Proc Natl Acad Sci USA. 2002, 99: 5591-5595.

- 44.
Andresen JM, Gayan J, Djousse L, Roberts S, Brocklebank D, Cherny SS, Cardon LR, Gusella JF, Macdonald ME, Myers RH, Housman DE, Wexler NS: The relationship between CAG repeat length and age of onset differs for Huntington's disease patients with juvenile onset or adult onset. Ann Hum Genet. 2007, 71: 295-301.

- 45.
Atsuta N, Watanabe H, Ito M, Banno H, Suzuki K, Katsuno M, Tanaka F, Tamakoshi A, Sobue G: Natural history of spinal and bulbar muscular atrophy (SBMA): a study of 223 Japanese patients. Brain. 2006, 129: 1446-1455.

- 46.
Andrews TC, Weeks RA, Turjanski N, Gunn RN, Watkins LHA, Sahakian B, Hodges JR, Rosser AE, Wood NW, Brooks DJ: Huntington's disease progression. PET and clinical observation. Brain. 1999, 122: 2353-2363.

- 47.
Kayed R, Head E, Thompson JL, McIntire TM, Milton SC, Cotman CW, Glabe CG: Common structure of soluble amyloid oligomers implies common mechanism of pathogenesis. Science. 2003, 300: 486-489.

- 48.
Riley BE, Orr HT: Polyglutamine neurodegenerative diseases and regulation of transcription: assembling the puzzle. Genes & Dev. 2006, 20: 2183-2192.

- 49.
Gidalevitz T, Ben-Zvi AB, Ho KH, Brignull HR, Morimoto RI: Progressive disruption of cellular protein folding in Models of polyglutamine diseases. Science. 2006, 311: 1471-1474.

- 50.
Shao J, Diamond MI: Polyglutamine diseases: emerging concepts in pathogenesis and therapy. Hum Mol Genet. 2007, 16: R115-R123.

- 51.
Clarke G, Lumsden CJ, Mclnnes RR: Inherited neurodegenerative diseases: the one-hit model of neurodegeneration. Hum Mol Genet. 2001, 10: 2269-2275.

- 52.
Burnett BG, Pittman RN: The polyglutamine neurodegenerative protein ataxin 3 regulates aggresome formation. Proc Natl Acad Sci USA. 2005, 102: 4330-4335.

- 53.
Zhong X, Pittman RN: Ataxin-3 binds VCP/p97 and regulates retrotranslocation of ERAD substrates. Hum Mol Genet. 2006, 15: 2409-2420.

- 54.
Hampton RY: ER-associated degradation in protein quality control and cellular regulation. Curr Opin Cell Biol. 2002, 14: 476-482.

- 55.
Ye Y, Meyer HH, Rapoport TA: The AAA ATPase Cdc48/p97 and its partners transport proteins from the ER into the cytosol. Nature. 2001, 414: 652-656.

- 56.
Ravikumar B, Duden R, Rubinsztein DC: Aggregate-prone proteins with polyglutamine and polyalanine expansions are degraded by autophagy. Hum Mol Genet. 2002, 11: 1107-1117.

- 57.
Arrasate M, Mitra S, Schweitzer ES, Segal MR, Finkbeiner S: Inclusion body formation reduces levels of mutant huntingtin and the risk of neuronal death. Nature. 2004, 431: 805-810.

## Acknowledgements

This study was funded by the Tokyo Metropolitan Government.

## Author information

### Affiliations

### Corresponding author

## Additional information

### Competing interests

The authors declare that they have no competing interests.

### Authors' contributions

Conceived and designed the experiments: KS. Analyzed the data: KS. Wrote the paper: KS SM. All authors read and approved the final manuscript.

## Authors’ original submitted files for images

Below are the links to the authors’ original submitted files for images.

## Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

## About this article

### Cite this article

Sugaya, K., Matsubara, S. Nucleation of protein aggregation kinetics as a basis for genotype-phenotype correlations in polyglutamine diseases.
*Mol Neurodegeneration* **4, **29 (2009). https://doi.org/10.1186/1750-1326-4-29

Received:

Accepted:

Published:

### Keywords

- Aggregation Kinetic
- polyQ Disease
- Elongation Time
- Individual Repeat
- polyQ Protein