The length dependence of the polyQ-mediated protein aggregation

Molecular Dynamics Simulations Applied to Structural and Dynamical Transitions of the Huntingtin Protein: A Review

Over the recent years, Huntington's disease (HD) has become widely discussed in the scientific literature especially because at the mutant level there are several contradictions regarding the aggregation mechanism. The specific role of the physiological huntingtin protein remains unknown, due to the lack of characterization of its entire crystallographic structure, making the experimental and theoretical research even harder when taking into consideration its involvement in multiple biological functions and its high affinity for different interacting partners. Different types of models, containing fewer (not more than 35 Qs) polyglutamine residues for the WT structure and above 35 Qs for the mutants, were subjected to classical or advanced MD simulations to establish the proteins' structural stability by evaluating their conformational changes. Outside the polyQ tract, there are two other regions of interest (the N17 domain and the polyP rich domain) considered to be essential for the aggregation kinetics at the mutant level. The polymerization process is considered to be dependent on the polyQ length. As the polyQ tract's dimension increases, the structures present more β-sheet conformations. Contrarily, it is also considered that the aggregation stability is not necessarily dependent on the number of Qs, while the initial stage of the aggregation seed might play the decisive role. A general assumption regarding the polyP domain is that it might preserve the polyQ structures soluble by acting as an antagonist for β-sheet formation.

Is there codon usage bias for poly-Q stretches in the human proteome?

We have analyzed codon usage for poly-Q stretches of different lengths for the human proteome. First, we have obtained that all long poly-Q stretches in Protein Data Bank (PDB) belong to the disordered regions. Second, we have found the bias for codon usage for glutamine homo-repeats in the human proteome. In the cases when the same codon is used for poly-Q stretches only CAG triplets are found. Similar results are obtained for human proteins with glutamine homo-repeats associated with diseases. Moreover, for proteins associated with diseases (from the HraDis database), the fraction of proteins for which the same codon is used for glutamine homo-repeats is less (22%) than for proteins from the human proteome (26%). We have demonstrated for poly-Q stretches in the human proteome that in some cases (28) the splicing sites correspond to the homo-repeats and in 11 cases, these sites appear at the C -terminal part of the homo-repeats with statistical significance 10 -8 .

Structural insights into the aggregation mechanism of huntingtin exon 1 protein fragment with different polyQ-lengths

Huntington disease is a neurodegenerative disorder caused by the expansion of polyglutamine (polyQ) at the N-terminal of the huntingtin exon 1 protein. The detailed structure and the mechanism behind this aggregation remain unclear and it is assumed that the polyQ undergoes a conformational transition to the β-sheet structure when it aggregates. Investigating the misfolding of polyQ facilitates the determination of the molecular mechanism of aggregation and can potentially help in developing a novel approach to inhibit polyQ aggregation. Moreover, the flanking sequences of the polyQ region play a vital role in structural changes and the aggregation mechanism. We performed all-atom molecular dynamics simulations to gain structural insights into the aggregation mechanism using eight different models with glutamine repeat lengths Q₂₇ , Q₂₇ P₁₁ , Q₃₄ , Q₃₅ , Q₃₆ , Q₄₀ , Q₅₀ , and Q₅₀ P₁₁ . In the models without flanking polyPs, we noticed that the transformation of a random coil to β-sheet occurs when the number of Q increases. We also found that the flanking polyPs prevent aggregation by decreasing the probability of forming a β-sheet structure. When polyQ length increases, the 17 N-terminal flanking residues are more likely to adopt a β-sheet conformation from α-helix and coil. From our simulations, we suggest that at least 34 glutamines are required for initiating aggregation and 40 residues length is critical for the aggregation of huntingtin exon 1 protein for disease onset. This study provides structural insights into misfolding and the role of flanking sequences in huntingtin aggregation which will further help in developing therapeutic strategies for Huntington's disease.

Glutamine Codon Usage and polyQ Evolution in Primates Depend on the Q Stretch Length

Amino acid usage in a proteome depends mostly on its taxonomy, as it does the codon usage in transcriptomes. Here, we explore the level of variation in the codon usage of a specific amino acid, glutamine, in relation to the number of consecutive glutamine residues. We show that CAG triplets are consistently more abundant in short glutamine homorepeats (polyQ, four to eight residues) than in shorter glutamine stretches (one to three residues), leading to the evolutionary growth of the repeat region in a CAG-dependent manner. The length of orthologous polyQ regions is mostly stable in primates, particularly the short ones. Interestingly, given a short polyQ the CAG usage is higher in unstable-in-length orthologous polyQ regions. This indicates that CAG triplets produce the necessary instability for a glutamine stretch to grow. Proteins related to polyQ-associated diseases behave in a more extreme way, with longer glutamine stretches in human and evolutionarily closer nonhuman primates, and an overall higher CAG usage. In the light of our results, we suggest an evolutionary model to explain the glutamine codon usage in polyQ regions.

Graphene oxide inhibits hIAPP amyloid fibrillation and toxicity in insulin-producing NIT-1 cells

Human islet amyloid polypeptide (hIAPP or amylin) aggregation is directly associated with pancreatic β-cell death and subsequent insulin deficiency in type 2 diabetes (T2D). Since no cure is currently available for T2D, it is of great benefit to devise new anti-aggregation molecules, which protect β-cells against hIAPP aggregation-induced toxicity. Engineered nanoparticles have been recently exploited as anti-aggregation nanomedicines. In this work, we studied graphene oxide (GO) nanosheets for their potential for hIAPP aggregation inhibition by combining computational modeling, biophysical characterization and cell toxicity measurements. Using discrete molecular dynamics (DMD) simulations and in vitro studies, we showed that GO exhibited an inhibitory effect on hIAPP aggregation. DMD simulations indicated that the strong binding of hIAPP to GO nanosheets was driven by hydrogen bonding and aromatic stacking and that the strong peptide-GO binding efficiently inhibited hIAPP self-association and aggregation on the nanosheet surface. Secondary structural changes of hIAPP upon GO binding derived from DMD simulations were consistent with circular dichroism (CD) spectroscopy measurements. Transmission electron microscopy (TEM) images confirmed the reduction of hIAPP aggregation in the presence of GO. Furthermore, we carried out a cell toxicity assay and found that these nanosheets protected insulin-secreting NIT-1 pancreatic β-cells against hIAPP-induced toxicity. Our multidisciplinary study suggests that GO nanosheets have the potential to be utilized as an anti-aggregation nanomedicine itself in addition to a biosensor or delivery vehicle for the mitigation of T2D progression.

A plant cell model of polyglutamine aggregation: Identification and characterisation of macromolecular and small-molecule anti-protein aggregation activity in vivo

In vitro studies have shown that LEA proteins from plants and invertebrates protect and stabilise other proteins under conditions of water stress, suggesting a role in stress tolerance. However, there is little information on LEA protein function in whole plants or plant cells, particularly with respect to their anti-aggregation activity. To address this, we expressed in tobacco BY-2 suspension cells an aggregation-prone protein based on that responsible for Huntington's disease (HD). In HD, abnormally long stretches of polyglutamine (polyQ) in huntingtin (Htt) protein cause aggregation of Htt fragments within cells. We constructed stably transformed BY-2 cell lines expressing enhanced green fluorescent protein (EGFP)-HttQ23 or EGFP-HttQ52 fusion proteins (encoding 23 or 52 glutamine residues, pertaining to the normal and disease states, respectively), as well as an EGFP control. EGFP-HttQ52 protein aggregated in the cytoplasm of transformed tobacco cells, which showed slow growth kinetics; in contrast, EGFP-HttQ23 or EGFP did not form aggregates and cells expressing these constructs grew normally. To test the effect of LEA proteins on protein aggregation in plant cells, we constructed cell lines expressing both EGFP-HttQ52 and LEA proteins (PM1, PM18, ZLDE-2 or AavLEA1) or a sHSP (PM31). Of these, AavLEA1 and PM31 reduced intracellular EGFP-HttQ52 aggregation and alleviated the associated growth inhibition, while PM18 and ZLDE-2 partially restored growth rates. Treatment of EGFP-HttQ52-expressing BY2 cells with the polyphenol epigallocatechin-3-gallate (EGCG) also reduced EGFP-HttQ52 aggregation and improved cell growth rate. The EGFP-HttQ52 cell line therefore has potential for characterising both macromolecular and small molecule inhibitors of protein aggregation in plant cells.

Influence of the protein context on the polyglutamine length-dependent elongation of amyloid fibrils

Polyglutamine (polyQ) diseases, including Huntington's disease, are neurodegenerative disorders associated with the abnormal expansion of a polyQ tract within nine proteins. The polyQ expansion is thought to be a major determinant in the development of neurotoxicity, triggering protein aggregation into amyloid fibrils, although non-polyQ regions play a modulating role. In this work, we investigate the relative importance of the polyQ length, its location within a host protein, and the conformational state of the latter in the amyloid fibril elongation. Model polyQ proteins made of the β-lactamase BlaP containing up to 79Q inserted at two different positions, and quartz crystal microbalance and atomic force microscopy were used for this purpose. We demonstrate that, independently of the polyQ tract location and the conformational state of the host protein, the relative elongation rate of fibrils increases linearly with the polyQ length. The slope of the linear fit is similar for both sets of chimeras (i.e., the elongation rate increases by ~1.9% for each additional glutamine), and is also similar to that previously observed for polyQ peptides. The elongation rate is, however, strongly influenced by the location of the polyQ tract within BlaP and the conformational state of BlaP. Moreover, comparison of our results with those previously reported for aggregation in solution indicates that these two parameters also modulate the ability of BlaP-polyQ chimeras to form the aggregation nucleus. Altogether our results suggest that non-polyQ regions are valuable targets in order to interfere with the process of amyloid fibril formation associated with polyQ diseases.

Computational approaches to understanding protein aggregation in neurodegeneration

The generation of toxic non-native protein conformers has emerged as a unifying thread among disorders such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. Atomic-level detail regarding dynamical changes that facilitate protein aggregation, as well as the structural features of large-scale ordered aggregates and soluble non-native oligomers, would contribute significantly to current understanding of these complex phenomena and offer potential strategies for inhibiting formation of cytotoxic species. However, experimental limitations often preclude the acquisition of high-resolution structural and mechanistic information for aggregating systems. Computational methods, particularly those combine both all-atom and coarse-grained simulations to cover a wide range of time and length scales, have thus emerged as crucial tools for investigating protein aggregation. Here we review the current state of computational methodology for the study of protein self-assembly, with a focus on the application of these methods toward understanding of protein aggregates in human neurodegenerative disorders.

Atomistic mechanisms of huntingtin N-terminal fragment insertion on a phospholipid bilayer revealed by molecular dynamics simulations

The huntingtin protein is characterized by a segment of consecutive glutamines (Q(N)) that is responsible for its fibrillation. As with other amyloid proteins, misfolding of huntingtin is related to Huntington's disease through pathways that can involve interactions with phospholipid membranes. Experimental results suggest that the N-terminal 17-amino-acid sequence (htt(NT)) positioned just before the Q(N) region is important for the binding of huntingtin to membranes. Through all-atom explicit solvent molecular dynamics simulations, we unveil the structure and dynamics of the htt(NT)Q(N) fragment on a phospholipid membrane at the atomic level. We observe that the insertion dynamics of this peptide can be described by four main steps-approach, reorganization, anchoring, and insertion-that are very diverse at the atomic level. On the membrane, the htt(NT) peptide forms a stable α-helix essentially parallel to the membrane with its nonpolar side-chains-mainly Leu-4, Leu-7, Phe-11 and Leu-14-positioned in the hydrophobic core of the membrane. Salt-bridges involving Glu-5, Glu-12, Lys-6, and Lys-15, as well as hydrogen bonds involving Thr-3 and Ser-13 with the phospholipids also stabilize the structure and orientation of the htt(NT) peptide. These observations do not significantly change upon adding the Q(N) region whose role is rather to provide, through its hydrogen bonds with the phospholipids' head group, a stable scaffold facilitating the partitioning of the htt(NT) region in the membrane. Moreover, by staying accessible to the solvent, the amyloidogenic Q(N) region could also play a key role for the oligomerization of htt(NT)Q(N) on phospholipid membranes.

Biomineral proteins from Mytilus edulis mantle tissue transcriptome

The common blue mussel, Mytilus edulis, has a bimineralic shell composed of approximately equal proportions of the two major polymorphs of calcium carbonate: calcite and aragonite. The exquisite biological control of polymorph production is the focus of research interest in terms of understanding the details of biomineralisation and the proteins involved in the process of complex shell formation. Recent advances in ease and availability of pyrosequencing and assembly have resulted in a sharp increase in transcriptome data for invertebrate biominerals. We have applied Roche 454 pyrosequencing technology to profile the transcriptome for the mantle tissue of the bivalve M. edulis. A comparison was made between the results of several assembly programs: Roche Newbler assembler versions 2.3, 2.5.2 and 2.6 and MIRA 3.2.1 and 3.4.0. The Newbler and MIRA assemblies were subsequently merged using the CAP3 assembler to give a higher consensus in alignments and a more accurate estimate of the true size of the M. edulis transcriptome. Comparison sequence searches show that the mantle transcripts for M. edulis encode putative proteins exhibiting sequence similarities with previously characterised shell proteins of other species of Mytilus, the Bivalvia Pinctada and haliotid gastropods. Importantly, this enhanced transcriptome has detected several transcripts that encode proteins with sequence similarity with previously described shell biomineral proteins including Shematrins and lysine-rich matrix proteins (KRMPs) not previously found in Mytilus.

The effects of flanking sequences in the interaction of polyglutamine peptides with a membrane bilayer

Huntington's disease (HD) is caused by the presence of an extended polyglutamine (polyQ) region at the N-terminus of the huntingtin (htt) protein. The presence of flanking sequences adjacent to the polyQ region has been reported to modulate the effects of potentially toxic protein-membrane interactions. In this study, we consider four peptide systems with various combinations of flanking sequences (KKQ35KK, KKQ35P11KK, N17Q35KK, N17Q35P11KK) and use atomistic molecular dynamics simulations to study the interactions with a DOPC lipid bilayer. We observe significant membrane thinning, disorderliness of lipid molecules, and compensation effects between the top and the bottom leaflets of the bilayer depending on the presence of particular flanking sequences. Overall, we find that the presence of the N-17 flanking sequence is crucial for membrane interactions. Polyproline decreases the interaction with the membrane in the absence of N-17, but enhances it when present along N-17.

Are long-range structural correlations behind the aggregration phenomena of polyglutamine diseases?

We have characterized the conformational ensembles of polyglutamine peptides of various lengths (ranging from to ), both with and without the presence of a C-terminal polyproline hexapeptide. For this, we used state-of-the-art molecular dynamics simulations combined with a novel statistical analysis to characterize the various properties of the backbone dihedral angles and secondary structural motifs of the glutamine residues. For (i.e., just above the pathological length for Huntington's disease), the equilibrium conformations of the monomer consist primarily of disordered, compact structures with non-negligible -helical and turn content. We also observed a relatively small population of extended structures suitable for forming aggregates including - and -strands, and - and -hairpins. Most importantly, for we find that there exists a long-range correlation (ranging for at least residues) among the backbone dihedral angles of the Q residues. For polyglutamine peptides below the pathological length, the population of the extended strands and hairpins is considerably smaller, and the correlations are short-range (at most residues apart). Adding a C-terminal hexaproline to suppresses both the population of these rare motifs and the long-range correlation of the dihedral angles. We argue that the long-range correlation of the polyglutamine homopeptide, along with the presence of these rare motifs, could be responsible for its aggregation phenomena.

Nine neurodegenerative diseases are caused by polyglutamine (polyQ) expansions greater than a given threshold in proteins with little or no homology except for the polyQ regions. The diseases all share a common feature: the formation of polyQ aggregates and eventual neuronal death. Using molecular dynamics simulations, we have explored the conformations of polyQ peptides. Results indicate that for peptides (i.e., just above the pathological length for Hungtington's disease), the equilibrium conformations were found to consist primarily of disordered, compact structures with a non-negligible -helical and turn content. We also observed a small population of extended structures suitable for forming aggregates. For peptides below the pathological length, the population of these structures was found to be considerably lower. For longer peptides, we found evidence for long-range correlations among the dihedral angles. This correlation turns out to be short-range for the smaller polyQ peptides, and is suppressed (along with the extended structural motifs) when a C-terminal polyproline tail is added to the peptides. We believe that the existence of these long-range correlations in above-threshold polyQ peptides, along with the presence of rare motifs, could be responsible for the experimentally observed aggregation phenomena associated with polyQ diseases.

Communication: Conformation state diagram of polypeptides: a chain length induced α-β transition

By using a generic coarse grained polypeptide model, we perform multicanonical molecular dynamics simulations for determining the equilibrium conformation state diagram of a single homopolypeptide chain as a function of the chain length and temperature. The state diagram highlights the thermal regimes of stability for various conformational patterns in polypeptides, including swollen, random and collapsed coils, globular structures, extended and bended α helices, and compact β bundles. Remarkably, at low temperatures we observe a sharp transition from extended α helix to compact β bundles as the chain length increases. This finding indicates that the chain length is one of the intrisic factors that can trigger α-β transformations in a broad class of polypeptides.

Neuronal aggregates are associated with phenotypic onset in the R6/2 Huntington's disease transgenic mouse

Huntington's disease (HD) is caused by the expansion of the polyglutamine tract expressed in the huntingtin protein. Data from patients show a strong negative correlation between CAG repeat size and age of disease onset. Recent studies in mixed background C57×CBA R6/2 mice suggest the inverse correlation observed in the human disease may not be replicated in some animal models of HD. To further clarify the relationship between repeat length and age of onset, congenic C57BL6/J R6/2 transgenic mice expressing 110, 260 or 310 CAG were tested in a comprehensive behavioral battery at multiple ages. Data confirmed the findings of earlier studies and indicate that on a pure C57BL6/J genetic background, R6/2 mice with larger repeats exhibit a delay in phenotypic onset with increasing polyglutamine size (6 weeks in 110 CAG and 17 weeks in 310 CAG mice). Further analysis confirmed a decrease in transgene transcript expression in 310 CAG mice as well as differential aggregated protein localization in association with repeat length. Mice expressing 110 CAG developed aggregates that localized almost exclusively to the nucleus of neuronal cells in the striatum and cortex. In contrast, tissue from 310 CAG mice exhibited predominantly extranuclear inclusions. Novel mutant protein analysis obtained using time-resolved fluorescence resonance energy transfer (FRET) revealed that soluble protein levels decreased with disease onset in R6/2 mice while aggregated protein levels increased. We believe that these data suggest a role for aggregation and inclusion localization in HD pathogenesis and propose a mechanism for the age of onset delay observed in R6/2 mice.

Location trumps length: polyglutamine-mediated changes in folding and aggregation of a host protein

Expanded CAG diseases are progressive neurodegenerative disorders in which specific proteins have an unusually long polyglutamine stretch. Although these proteins share no other sequence or structural homologies, they all aggregate into intracellular inclusions that are believed to be pathological. We sought to determine what impact the position and number of glutamines have on the structure and aggregation of the host protein, apomyoglobin. Variable-length polyQ tracts were inserted either into the loop between the C- and D-helices (Q(n)CD) or at the N-terminus (Q(n)NT). The Q(n)CD mutants lost some α-helix and gained unordered and/or β-sheet in a length-dependent manner. These mutants were partially unfolded and rapidly assembled into soluble chain-like oligomers. In sharp contrast, the Q(n)NT mutants largely retained wild-type tertiary structure but associated into long, fibrillar aggregates. Control proteins with glycine-serine repeats (GS(8)CD and GS(8)NT) were produced. GS(8)CD exhibited similar structural perturbations and aggregation characteristics to an analogously sized Q(16)CD, indicating that the observed effects are independent of amino acid composition. In contrast to Q(16)NT, GS(8)NT did not form fibrillar aggregates. Thus, soluble oligomers are produced through structural perturbation and do not require polyQ, whereas classic fibrils arise from specific polyQ intermolecular interactions in the absence of misfolding.

Proteomic analysis of the organic matrix of the abalone Haliotis asinina calcified shell

Background

The formation of the molluscan shell is regulated to a large extent by a matrix of extracellular macromolecules that are secreted by the shell forming tissue, the mantle. This so called "calcifying matrix" is a complex mixture of proteins and glycoproteins that is assembled and occluded within the mineral phase during the calcification process. While the importance of the calcifying matrix to shell formation has long been appreciated, most of its protein components remain uncharacterised.

Results

Recent expressed sequence tag (EST) investigations of the mantle tissue from the tropical abalone (Haliotis asinina) provide an opportunity to further characterise the proteins in the shell by a proteomic approach. In this study, we have identified a total of 14 proteins from distinct calcified layers of the shell. Only two of these proteins have been previously characterised from abalone shells. Among the novel proteins are several glutamine- and methionine-rich motifs and hydrophobic glycine-, alanine- and acidic aspartate-rich domains. In addition, two of the new proteins contained Kunitz-like and WAP (whey acidic protein) protease inhibitor domains.

Conclusion

This is one of the first comprehensive proteomic study of a molluscan shell, and should provide a platform for further characterization of matrix protein functions and interactions.

Neurotoxic protein oligomerisation associated with polyglutamine diseases

Polyglutamine (polyQ) diseases are associated with a CAG/polyQ expansion mutation in unrelated proteins. Upon elongation of the glutamine tract, disease proteins aggregate within cells, mainly in the central nervous system (CNS) and this aggregation process is associated with neurotoxicity. However, it remains unclear to what extent and how this aggregation causes neuronal dysfunction in the CNS. Aiming at preventing neuronal dysfunction, it will be crucial to determine the links between aggregation and cellular dysfunction, understand the folding pathway of polyQ proteins and discover the relative neurotoxicity of polyQ protein species formed along the aggregation pathway. Here, we review what is known about conformations of polyQ peptides and proteins in their monomeric state from experimental and modelling data, how conformational changes of polyQ proteins relate to their oligomerisation and morphology of aggregates and which cellular function are impaired by oligomers, in vitro and in vivo. We also summarise the key modulatory cellular mechanisms and co-factors, which could affect the folding pathway and kinetics of polyQ aggregation. Although many studies have investigated the relationship between polyQ aggregation and toxicity, these have mainly focussed on investigating changes in the formation of the classical hallmark of polyQ diseases, i.e. microscopically visible inclusion bodies. However, recent studies in which oligomeric species have been considered start to shed light on the identity of neurotoxic oligomeric species. Initial evidence suggests that conformational changes induced by polyQ expansions and their surrounding sequence lead to the formation of particular oligomeric intermediates that may differentially affect neurotoxicity.

Polyglutamine induced misfolding of huntingtin exon1 is modulated by the flanking sequences

Polyglutamine (polyQ) expansion in exon1 (XN1) of the huntingtin protein is linked to Huntington's disease. When the number of glutamines exceeds a threshold of approximately 36-40 repeats, XN1 can readily form amyloid aggregates similar to those associated with disease. Many experiments suggest that misfolding of monomeric XN1 plays an important role in the length-dependent aggregation. Elucidating the misfolding of a XN1 monomer can help determine the molecular mechanism of XN1 aggregation and potentially help develop strategies to inhibit XN1 aggregation. The flanking sequences surrounding the polyQ region can play a critical role in determining the structural rearrangement and aggregation mechanism of XN1. Few experiments have studied XN1 in its entirety, with all flanking regions. To obtain structural insights into the misfolding of XN1 toward amyloid aggregation, we perform molecular dynamics simulations on monomeric XN1 with full flanking regions, a variant missing the polyproline regions, which are hypothesized to prevent aggregation, and an isolated polyQ peptide (Q(n)). For each of these three constructs, we study glutamine repeat lengths of 23, 36, 40 and 47. We find that polyQ peptides have a positive correlation between their probability to form a beta-rich misfolded state and their expansion length. We also find that the flanking regions of XN1 affect its probability to form a beta-rich state compared to the isolated polyQ. Particularly, the polyproline regions form polyproline type II helices and decrease the probability of the polyQ region to form a beta-rich state. Additionally, by lengthening polyQ, the first N-terminal 17 residues are more likely to adopt a beta-sheet conformation rather than an alpha-helix conformation. Therefore, our molecular dynamics study provides a structural insight of XN1 misfolding and elucidates the possible role of the flanking sequences in XN1 aggregation.

Modulation of polyglutamine conformations and dimer formation by the N-terminus of huntingtin

Polyglutamine expansions within different proteins are associated with nine different neurodegenerative diseases. There is growing interest in understanding the roles of flanking sequences from disease-relevant proteins in the intrinsic conformational and aggregation properties of polyglutamine. We report results from atomistic simulations and circular dichroism experiments that quantify the effect of the N-terminal 17-residue (Nt17) segment of the huntingtin protein on polyglutamine conformations and intermolecular interactions. We show that the Nt17 segment and polyglutamine domains become increasingly disordered as polyglutamine length (N) increases in Nt17-Q(N) constructs. Hydrophobic groups within Nt17 become sequestered in intramolecular interdomain interfaces. We also show that the Nt17 segment suppresses the intrinsic propensity of polyglutamine aggregation. This inhibition arises from the incipient micellar structures adopted by monomeric forms of the peptides with Nt17 segments. The degree of intermolecular association increases with increasing polyglutamine length and is governed mainly by associations between polyglutamine domains. Comparative analysis of intermolecular associations for different polyglutamine-containing constructs leads to clearer interpretations of recently published experimental data. Our results suggest a framework for fibril formation and identify roles for flanking sequences in the modulation of polyglutamine aggregation.

Cerebellar disorders--at the crossroad of molecular pathways and diagnosis

Our understanding of the pathogenesis of cerebellar ataxias has started several decades ago and is continuously growing. The numerous mechanisms of cerebellar dysfunction are being discovered by numerous groups of researchers worldwide. Neuronal damage results from a complex interaction of metabolic pathways, which leads to symptoms observed in cerebellar disorders. The main mechanisms at the molecular level are the following: impairment of DNA repair and replication, deregulation of transcription/deficits of processing/transport of RNA, abnormal protein transport and misfolding, aggregates both at the nuclear and cytosolic level, activation of caspases, apoptosis, involvement of autophagic mechanisms, oxidative stress and mitochondrial dysfunction, excitotoxicity, abnormal lipid metabolism, impaired axonal transport and vesicle trafficking, and defects of the neurotransmission. The convergence of the current clinical classification with molecular findings is expected. This integration is a basic substrate for the rationale development of therapies.

Dynamics and stability of amyloid-like steric zipper assemblies with hydrophobic dry interfaces

Recent seminal investigations have suggested that the basic structural motif of amyloid fibers may be constituted by a tight association of two facing beta-sheets (steric zipper). Although this model has been derived from crystal structures of small peptide models, several theoretical investigations, essentially focused on steric zipper interface containing large polar and/or aromatic side chains, have confirmed the stability of this motif in a crystal-free context. To analyze the general validity of these findings, we carried out molecular dynamics (MD) simulations on aggregates stabilized by steric zipper interfaces made also of small or hydrophobic residues. In particular, we here characterized assemblies formed by the peptides SSTSAA and VQIVYK, whose structures have been recently solved at high resolution. In contrast to previous results obtained for polar/aromatic aggregates of the same size and with similar interface area, steric zipper assemblies composed of a pair of 10-stranded beta-sheets show high fluctuations and significant distortions in the simulation timescales (40-60 ns). Taking into account the crystal packing, the effect of the addition of an extra sheet to the assemblies was also evaluated. The MD results indicate that this addition does not provide extra-stabilization to the pair of sheet models. Although present data do not preclude the possibility that the steric zipper association identified in the crystal structure is the basic motif of SSTSAA and VQIVYK fibers, our findings highlight the importance of the nature of residues directly involved in the motif. Indeed, polar and aromatic residues that may form intrasheet and intersheet interactions likely provide a strong contribution to the steric zipper motif stability. Along this line, assemblies endowed with hydrophobic residues presumably require larger interfaces. In line with this suggestion, MD analysis of the HET-s(218-289) prion models composed of a similar number of strands shows that the assembly is endowed with a remarkable stability.

Thermodynamics of beta-sheet formation in polyglutamine

The role of beta-sheets in the early stages of protein aggregation, specifically amyloid formation, remains unclear. Interpretations of kinetic data have led to a specific model for the role of beta-sheets in polyglutamine aggregation. According to this model, monomeric polyglutamine, which is intrinsically disordered, goes through a rare conversion into an ordered, metastable, beta-sheeted state that nucleates aggregation. It has also been proposed that the probability of forming the critical nucleus, a specific beta-sheet conformation for the monomer, increases with increasing chain length. Here, we test this model using molecular simulations. We quantified free energy profiles in terms of beta-content for monomeric polyglutamine as a function of chain length. In accord with estimates from experimental data, the free energy penalties for forming beta-rich states are in the 10-20 kcal/mol range. However, the length dependence of these free energy penalties does not mirror interpretations of kinetic data. In addition, although homodimerization of disordered molecules is spontaneous, the imposition of conformational restraints on polyglutamine molecules does not enhance the spontaneity of intermolecular associations. Our data lead to the proposal that beta-sheet formation is an attribute of peptide-rich phases such as high molecular weight aggregates rather than monomers or oligomers.

Insights from coarse-grained Gō models for protein folding and dynamics

Exploring the landscape of large scale conformational changes such as protein folding at atomistic detail poses a considerable computational challenge. Coarse-grained representations of the peptide chain have therefore been developed and over the last decade have proved extremely valuable. These include topology-based Gō models, which constitute a smooth and funnel-like approximation to the folding landscape. We review the many variations of the Gō model that have been employed to yield insight into folding mechanisms. Their success has been interpreted as a consequence of the dominant role of the native topology in folding. The role of local contact density in determining protein dynamics is also discussed and is used to explain the ability of Gō-like models to capture sequence effects in folding and elucidate conformational transitions.

Atomistic simulations of the effects of polyglutamine chain length and solvent quality on conformational equilibria and spontaneous homodimerization

Aggregation of expanded polyglutamine tracts is associated with nine different neurodegenerative diseases, including Huntington's disease. Experiments and computer simulations have demonstrated that monomeric forms of polyglutamine molecules sample heterogeneous sets of collapsed structures in water. The current work focuses on a mechanistic characterization of polyglutamine homodimerization as a function of chain length and temperature. These studies were carried out using molecular simulations based on a recently developed continuum solvation model that was designed for studying conformational and binding equilibria of intrinsically disordered molecules such as polyglutamine systems. The main results are as follows: Polyglutamine molecules form disordered, collapsed globules in aqueous solution. These molecules spontaneously associate at conditions approaching those of typical in vitro experiments for chains of length N>/=15. The spontaneity of these homotypic associations increases with increasing chain length and decreases with increasing temperature. Similar and generic driving forces govern both collapse and spontaneous homodimerization of polyglutamine in aqueous milieus. Collapse and dimerization maximize self-interactions and reduce the interface between polyglutamine molecules and the surrounding solvent. Other than these generic considerations, there do not appear to be any specific structural requirements for either chain collapse or chain dimerization; that is, both collapse and dimerization are nonspecific in that disordered globules form disordered dimers. In fact, it is shown that the driving force for intermolecular associations is governed by spontaneous conformational fluctuations within monomeric polyglutamine. These results suggest that polyglutamine aggregation is unlikely to follow a homogeneous nucleation mechanism with the monomer as the critical nucleus. Instead, the results support the formation of disordered, non-beta-sheet-like soluble molten oligomers as early intermediates--a proposal that is congruent with recent experimental data.

Formation and growth of oligomers: a Monte Carlo study of an amyloid tau fragment

Small oligomers formed early in the process of amyloid fibril formation may be the major toxic species in Alzheimer's disease. We investigate the early stages of amyloid aggregation for the tau fragment AcPHF6 (Ac-VQIVYK-NH2) using an implicit solvent all-atom model and extensive Monte Carlo simulations of 12, 24, and 36 chains. A variety of small metastable aggregates form and dissolve until an aggregate of a critical size and conformation arises. However, the stable oligomers, which are β-sheet-rich and feature many hydrophobic contacts, are not always growth-ready. The simulations indicate instead that these supercritical oligomers spend a lengthy period in equilibrium in which considerable reorganization takes place accompanied by exchange of chains with the solution. Growth competence of the stable oligomers correlates with the alignment of the strands in the β-sheets. The larger aggregates seen in our simulations are all composed of two twisted β-sheets, packed against each other with hydrophobic side chains at the sheet–sheet interface. These β-sandwiches show similarities with the proposed steric zipper structure for PHF6 fibrils but have a mixed parallel/antiparallel β-strand organization as opposed to the parallel organization found in experiments on fibrils. Interestingly, we find that the fraction of parallel β-sheet structure increases with aggregate size. We speculate that the reorganization of the β-sheets into parallel ones is an important rate-limiting step in the formation of PHF6 fibrils.

It is believed that the self association of certain protein molecules into aggregated structures, known as amyloid fibrils, plays an important role in a variety of human diseases, such as Alzheimer's disease and Parkinson's disease. Although the ability to form such amyloid fibrils is a common property for proteins, the process leading to these fibrils is incompletely understood. The early stages of the process involve small transient heterogeneous structures made of a few protein chains and are especially difficult to characterize. Here we use atomic-level simulations to explore the early part of the aggregation process for a fibril-forming fragment of the protein tau associated with Alzheimer's disease. We find that a multitude of small aggregates, rich in sheetlike structures, form through a nucleation process. Interestingly, a statistically preferred type of aggregate, consisting of two tightly packed sheets, emerges with increasing aggregate size. Growth of these larger aggregates seems to be a slow process that correlates with the emergence of more uniformly ordered sheets. We speculate that reorganization of the protein chains leading to that ordered arrangement is an important bottleneck to amyloid fibril formation for this peptide.

CAG repeat lengths > or =335 attenuate the phenotype in the R6/2 Huntington's disease transgenic mouse

With spontaneous elongation of the CAG repeat in the R6/2 transgene to > or =335, resulting in a transgene protein too large for passive entry into nuclei via the nuclear pore, we observed an abrupt increase in lifespan to >20 weeks, compared to the 12 weeks common in R6/2 mice with 150 repeats. In the > or =335 CAG mice, large ubiquitinated aggregates of mutant protein were common in neuronal dendrites and perikaryal cytoplasm, but intranuclear aggregates were small and infrequent. Message and protein for the > or =335 CAG transgene were reduced to one-third that in 150 CAG R6/2 mice. Neurological and neurochemical abnormalities were delayed in onset and less severe than in 150 CAG R6/2 mice. These findings suggest that polyQ length and pathogenicity in Huntington's disease may not be linearly related, and pathogenicity may be less severe with extreme repeats. Both diminished mutant protein and reduced nuclear entry may contribute to phenotype attenuation.

Expansion of the octarepeat domain alters the misfolding pathway but not the folding pathway of the prion protein

A misfolded conformation of the prion protein (PrP), PrP (Sc), is the essential component of prions, the infectious agents that cause transmissible neurodegenerative diseases. Insertional mutations that lead to an increase in the number of octarepeats (ORs) in PrP are linked to familial human prion disease. In this study, we investigated how expansion of the OR domain causes PrP to favor a prion-like conformation. Therefore, we compared the conformational and aggregation modulating properties of wild-type versus expanded OR domains, either as a fusion construct with the protein G B1 domain (GB1-OR) or as an integral part of full-length mouse PrP (MoPrP). Using circular dichroism spectroscopy, we first demonstrated that ORs are not unfolded but exist as an ensemble of three distinct conformers: polyproline helix-like, beta-turn, and "Trp-related". Domain expansion had little effect on the conformation of GB1-OR fusion proteins. When part of MoPrP however, OR domain expansion changed PrP's folding landscape, not by hampering the production of native alpha-helical monomers but by greatly reducing the propensity to form amyloid and by altering the assembly of misfolded, beta-rich aggregates. These features may relate to subtle pH-dependent conformational differences between wild-type and mutant monomers. In conclusion, we propose that PrP insertional mutations are pathogenic because they enhance specific misfolding pathways of PrP rather than by undermining native folding. This idea was supported by a trial bioassay in transgenic mice overexpressing wild-type MoPrP, where intracerebral injection of recombinant MoPrP with an expanded OR domain but not wild-type MoPrP caused prion disease.

Probing protein aggregation using discrete molecular dynamics

Understanding the role of biomolecular dynamics in cellular processes leading to human diseases and the ability to rationally manipulate these processes is of fundamental importance in scientific research. The last decade has witnessed significant progress in probing biophysical behavior of proteins. However, we are still limited in understanding how changes in protein dynamics and inter-protein interactions occurring in short length- and time-scales lead to aberrations in their biological function. Bridging this gap in biology probed using computer simulations marks a challenging frontier in computational biology. Here we examine hypothesis-driven simplified protein models in conjunction with discrete molecular dynamics in the study of protein aggregation, implicated in series of neurodegenerative diseases, such as Alzheimer's and Huntington's diseases. Discrete molecular dynamics simulations of simplified protein models have emerged as a powerful methodology with its ability to bridge the gap in time and length scales from protein dynamics to aggregation, and provide an indispensable tool for probing protein aggregation.