Assessing polyglutamine conformation in the nucleating event by molecular dynamics simulations

Advances in the understanding of protein misfolding and aggregation through molecular dynamics simulation

Aberrant protein folding known as protein misfolding is counted as one of the striking factors of neurodegenerative diseases. The extensive range of pathologies caused by protein misfolding, aggregation and subsequent accumulation are mainly classified into either gain of function diseases or loss of function diseases. In order to seek for novel strategies for treatment and diagnosis of neurodegenerative diseases, insights into the mechanism of misfolding and aggregation is essential. A comprehensive knowledge on the factors influencing misfolding and aggregation is required as well. An extensive experimental study on protein aggregation is somewhat challenging due to the insoluble and noncrystalline nature of amyloid fibrils. Thus there has been a growing use of computational approaches including Monte Carlo simulation, docking simulation, molecular dynamics simulation in the study of protein misfolding and aggregation. The review presents a discussion on molecular dynamics simulation alone as to how it has emerged as a promising tool in the understanding of protein misfolding and aggregation in general, detailing upon three different aspects considering four misfold prone proteins in particular. It is noticeable that all four proteins considered in this review i.e prion, superoxide dismutase1, huntingtin and amyloid β are linked to chronic neurodegenerative diseases with debilitating effects. Initially the review elaborates on the factors influencing the misfolding and aggregation. Next, it addresses our current understanding of the amyloid structures and the associated aggregation mechanisms, finally, summarizing the contribution of this computational tool in the search for therapeutic strategies against the respective protein-deposition diseases.

Membrane Interactions Accelerate the Self-Aggregation of Huntingtin Exon 1 Fragments in a Polyglutamine Length-Dependent Manner

The accumulation of aggregated protein is a typical hallmark of many human neurodegenerative disorders, including polyglutamine-related diseases such as chorea Huntington. Misfolding of the amyloidogenic proteins gives rise to self-assembled complexes and fibres. The huntingtin protein is characterised by a segment of consecutive glutamines which, when exceeding ~ 37 residues, results in the occurrence of the disease. Furthermore, it has also been demonstrated that the 17-residue amino-terminal domain of the protein (htt17), located upstream of this polyglutamine tract, strongly correlates with aggregate formation and pathology. Here, we demonstrate that membrane interactions strongly accelerate the oligomerisation and β-amyloid fibril formation of htt17-polyglutamine segments. By using a combination of biophysical approaches, the kinetics of fibre formation is investigated and found to be strongly dependent on the presence of lipids, the length of the polyQ expansion, and the polypeptide-to-lipid ratio. Finally, the implications for therapeutic approaches are discussed.

A Perspective on Biophysical Studies of Crystallin Aggregation and Implications for Cataract Formation

Lens crystallins are subject to various types of damage during their lifetime which triggers protein misfolding and aggregation, ultimately causing cataracts. There are several models for crystallin aggregation, but a comprehensive picture of the mechanism of cataract is still underway. The complex biomolecular interactions underlying crystallin aggregation have motivated major efforts to resolve the structural details and mechanism of aggregation using multiple biophysical techniques at different resolutions. Together, experimental and computational approaches identify and characterize both amyloidogenic and amorphous aggregates leading to an improved understanding of crystallin aggregation. A rigorous characterization of the aggregation-prone intermediates is crucial in cataract-mediated drug discovery. This Perspective summarizes recent biophysical studies on lens crystallin aggregation. We evaluate the outstanding challenges, future outlook, and rewards in this fertile field of research. With lessons learned from protein folding and multiple pathways of aggregation, we highlight the differences in the overall mechanisms of age-related and congenital cataracts. We expect that a correlation between the existing and developing biophysical techniques would provide a platform to study amyloid architecture in the eye lens and reduce the existing gaps in our understanding of crystallin biophysics.

Length-Dependent Structural Transformations of Huntingtin PolyQ Domain Upon Binding to 2D-Nanomaterials

There is a strong negative correlation between the polyglutamine (polyQ) domain length (Q-length) in the intrinsically disordered Huntingtin protein (Htt) exon-1 and the age of onset of Huntington's disease (HD). PolyQ of Q-length longer than 40 has the propensity of forming very compact aggregate structures, leading to HD at full penetrance. Recent advances in nanobiotechnology provided a new platform for the development of novel diagnosis and therapeutics. Here, we explore the possibility of utilizing 2D-nanomaterials to inhibit the formation of supercompact polyQ structures through the so-called “folding-upon-binding” where the protein structure is dependent on the binding substrate. Using molecular dynamics simulations, we characterize two polyQ peptides with Q-length of 22 (Q22, normal length) and 46 (Q46, typical length causing HD) binding to both graphene and molybdenum disulfide (MoS₂) nanosheets, which have been applied as antibacterial or anticancer agents. Upon binding, Q22 unfolds and elongates on both grapheme and MoS₂ surfaces, regardless of its initial conformation, with graphene showing slightly stronger effect. In contrast, initially collapsed Q46 remains mostly collapsed within our simulation time on both nanosheets even though they do provide some “stretching” to Q46 as well. Further analyses indicate that the hydrophobic nature of graphene/MoS₂ promotes the stretching of polyQ on nanosheets. However, there is strong competition with the intra-polyQ interactions (mainly internal hydrogen bonds) leading to the disparate folding/binding behaviors of Q22 and Q46. Our results present distinct Q-length specific behavior of the polyQ domain upon binding to two types of 2D-nanomaterials which holds clinical relevance for Huntington's disease.

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.

Conformational studies of pathogenic expanded polyglutamine protein deposits from Huntington's disease

Huntington’s disease, like other neurodegenerative diseases, continues to lack an effective cure. Current treatments that address early symptoms ultimately fail Huntington’s disease patients and their families, with the disease typically being fatal within 10–15 years from onset. Huntington’s disease is an inherited disorder with motor and mental impairment, and is associated with the genetic expansion of a CAG codon repeat encoding a polyglutamine-segment-containing protein called huntingtin. These Huntington’s disease mutations cause misfolding and aggregation of fragments of the mutant huntingtin protein, thereby likely contributing to disease toxicity through a combination of gain-of-toxic-function for the misfolded aggregates and a loss of function from sequestration of huntingtin and other proteins. As with other amyloid diseases, the mutant protein forms non-native fibrillar structures, which in Huntington’s disease are found within patients’ neurons. The intracellular deposits are associated with dysregulation of vital processes, and inter-neuronal transport of aggregates may contribute to disease progression. However, a molecular understanding of these aggregates and their detrimental effects has been frustrated by insufficient structural data on the misfolded protein state. In this review, we examine recent developments in the structural biology of polyglutamine-expanded huntingtin fragments, and especially the contributions enabled by advances in solid-state nuclear magnetic resonance spectroscopy. We summarize and discuss our current structural understanding of the huntingtin deposits and how this information furthers our understanding of the misfolding mechanism and disease toxicity mechanisms.

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.

Two zinc-binding domains in the transporter AdcA from facilitate high-affinity binding and fast transport of zinc

Zinc is an essential metal in bacteria. One important bacterial zinc transporter is AdcA, and most bacteria possess AdcA homologs that are single-domain small proteins due to better efficiency of protein biogenesis. However, a double-domain AdcA with two zinc-binding sites is significantly overrepresented in Streptococcus species, many of which are major human pathogens. Using molecular simulation and experimental validations of AdcA from Streptococcus pyogenes, we found here that the two AdcA domains sequentially stabilize the structure upon zinc binding, indicating an organization required for both increased zinc affinity and transfer speed. This structural organization appears to endow Streptococcus species with distinct advantages in zinc-depleted environments, which would not be achieved by each single AdcA domain alone. This enhanced zinc transport mechanism sheds light on the significance of the evolution of the AdcA domain fusion, provides new insights into double-domain transporter proteins with two binding sites for the same ion, and indicates a potential target of antimicrobial drugs against pathogenic Streptococcus species.

Emerging β-Sheet Rich Conformations in Supercompact Huntingtin Exon-1 Mutant Structures

There exists strong correlation between the extended polyglutamines (polyQ) within exon-1 of Huntingtin protein (Htt) and age onset of Huntington's disease (HD); however, the underlying molecular mechanism is still poorly understood. Here we apply extensive molecular dynamics simulations to study the folding of Htt-exon-1 across five different polyQ-lengths. We find an increase in secondary structure motifs at longer Q-lengths, including β-sheet content that seems to contribute to the formation of increasingly compact structures. More strikingly, these longer Q-lengths adopt supercompact structures as evidenced by a surprisingly small power-law scaling exponent (0.22) between the radius-of-gyration and Q-length that is substantially below expected values for compact globule structures (∼0.33) and unstructured proteins (∼0.50). Hydrogen bond analyses further revealed that the supercompact behavior of polyQ is mainly due to the "glue-like" behavior of glutamine's side chains with significantly more side chain-side chain H-bonds than regular proteins in the Protein Data Bank (PDB). The orientation of the glutamine side chains also tend to be "buried" inside, explaining why polyQ domains are insoluble on their own.

Molecular dynamics analysis of the aggregation propensity of polyglutamine segments

Protein misfolding and aggregation is a pathogenic feature shared among at least ten polyglutamine (polyQ) neurodegenerative diseases. While solvent-solution interaction is a key factor driving protein folding and aggregation, the solvation properties of expanded polyQ tracts are not well understood. By using GPU-enabled all-atom molecular dynamics simulations of polyQ monomers in an explicit solvent environment, this study shows that solvent-polyQ interaction propensity decreases as the lengths of polyQ tract increases. This study finds a predominance in long-distance interactions between residues far apart in polyQ sequences with longer polyQ segments, that leads to significant conformational differences. This study also indicates that large loops, comprised of parallel β-structures, appear in long polyQ tracts and present new aggregation building blocks with aggregation driven by long-distance intra-polyQ interactions. Finally, consistent with previous observations using coarse-grain simulations, this study demonstrates that there is a gain in the aggregation propensity with increased polyQ length, and that this gain is correlated with decreasing ability of solvent-polyQ interaction. These results suggest the modulation of solvent-polyQ interactions as a possible therapeutic strategy for treating polyQ diseases.

Cataract-associated P23T γD-crystallin retains a native-like fold in amorphous-looking aggregates formed at physiological pH

Cataracts cause vision loss through the large-scale aggregation of eye lens proteins as a result of ageing or congenital mutations. The development of new treatments is hindered by uncertainty about the nature of the aggregates and their mechanism of formation. We describe the structure and morphology of aggregates formed by the P23T human γD-crystallin mutant associated with congenital cataracts. At physiological pH, the protein forms aggregates that look amorphous and disordered by electron microscopy, reminiscent of the reported formation of amorphous deposits by other crystallin mutants. Surprisingly, solid-state NMR reveals that these amorphous deposits have a high degree of structural homogeneity at the atomic level and that the aggregated protein retains a native-like conformation, with no evidence for large-scale misfolding. Non-physiological destabilizing conditions used in many in vitro aggregation studies are shown to yield qualitatively different, highly misfolded amyloid-like fibrils.

Aggregation of eye lens proteins leads to cataracts, a major cause of blindness. Here the authors use solid state NMR to probe the structure of γD-crystallin eye lens proteins aggregates, which are found to retain a native-like conformation.

The Aggregation Free Energy Landscapes of Polyglutamine Repeats

Aggregates of proteins containing polyglutamine (polyQ) repeats are strongly associated with several neurodegenerative diseases. The length of the repeats correlates with the severity of the disease. Previous studies have shown that pure polyQ peptides aggregate by nucleated growth polymerization and that the size of the critical nucleus (n*) decreases from tetrameric to dimeric and monomeric as length increases from Q₁₈ to Q₂₆. Why the critical nucleus size changes with repeat-length has been unclear. Using the associative memory, water-mediated, structure and energy model, we construct the aggregation free energy landscapes for polyQ peptides of different repeat-lengths. These studies show that the monomer of the shorter repeat-length (Q₂₀) prefers an extended conformation and that its aggregation indeed has a trimeric nucleus (n* ∼ 3), while a longer repeat-length monomer (Q₃₀) prefers a β-hairpin conformation which then aggregates in a downhill fashion at 0.1 mM. For an intermediate length peptide (Q₂₆), there is an equal preference for hairpin and extended forms in the monomer which leads to a mixed inhomogeneous nucleation mechanism for fibrils. The predicted changes of monomeric structure and nucleation mechanism are confirmed by studying the aggregation free energy profile for a polyglutamine repeat with site-specific PG mutations that favor the hairpin form, giving results in harmony with experiments on this system.

Effects of the enlargement of polyglutamine segments on the structure and folding of ataxin-2 and ataxin-3 proteins

Spinocerebellar ataxia type 2 (SCA2) and type 3 (SCA3) are two common autosomal-dominant inherited ataxia syndromes, both of which are related to the unstable expansion of trinucleotide CAG repeats in the coding region of the related ATXN2 and ATXN3 genes, respectively. The poly-glutamine (poly-Q) tract encoded by the CAG repeats has long been recognized as an important factor in disease pathogenesis and progress. In this study, using the I-TASSER method for 3D structure prediction, we investigated the effect of poly-Q tract enlargement on the structure and folding of ataxin-2 and ataxin-3 proteins. Our results show good agreement with the known experimental structures of the Josephin and UIM domains providing credence to the simulation results presented here, which show that the enlargement of the poly-Q region not only affects the local structure of these regions but also affects the structures of functional domains as well as the whole protein. The changes observed in the predicted models of the UIM domains in ataxin-3 when the poly-Q track is enlarged provide new insights on possible pathogenic mechanisms.

Polyglutamine Fibrils: New Insights into Antiparallel β-Sheet Conformational Preference and Side Chain Structure

Understanding the structure of polyglutamine (polyQ) amyloid-like fibril aggregates is crucial to gaining insights into the etiology of at least ten neurodegenerative disorders, including Huntington's disease. Here, we determine the structure of D2Q10K2 (Q10) fibrils using ultraviolet resonance Raman (UVRR) spectroscopy and molecular dynamics (MD). Using UVRR, we determine the fibril peptide backbone Ψ and glutamine (Gln) side chain χ3 dihedral angles. We find that most of the fibril peptide bonds adopt antiparallel β-sheet conformations; however, a small population of peptide bonds exist in parallel β-sheet structures. Using MD, we simulate three different potential fibril structural models that consist of either β-strands or β-hairpins. Comparing the experimentally measured Ψ and χ3 angle distributions to those obtained from the MD simulated models, we conclude that the basic structural motif of Q10 fibrils is an extended β-strand structure. Importantly, we determine from our MD simulations that Q10 fibril antiparallel β-sheets are thermodynamically more stable than parallel β-sheets. This accounts for why polyQ fibrils preferentially adopt antiparallel β-sheet conformations instead of in-register parallel β-sheets like most amyloidogenic peptides. In addition, we directly determine, for the first time, the structures of Gln side chains. Our structural data give new insights into the role that the Gln side chains play in the stabilization of polyQ fibrils. Finally, our work demonstrates the synergistic power and utility of combining UVRR measurements and MD modeling to determine the structure of amyloid-like fibrils.

Huntingtin exon 1 fibrils feature an interdigitated β-hairpin-based polyglutamine core

Polyglutamine expansion within the exon1 of huntingtin leads to protein misfolding, aggregation, and cytotoxicity in Huntington's disease. This incurable neurodegenerative disease is the most prevalent member of a family of CAG repeat expansion disorders. Although mature exon1 fibrils are viable candidates for the toxic species, their molecular structure and how they form have remained poorly understood. Using advanced magic angle spinning solid-state NMR, we directly probe the structure of the rigid core that is at the heart of huntingtin exon1 fibrils and other polyglutamine aggregates, via measurements of long-range intramolecular and intermolecular contacts, backbone and side-chain torsion angles, relaxation measurements, and calculations of chemical shifts. These experiments reveal the presence of β-hairpin-containing β-sheets that are connected through interdigitating extended side chains. Despite dramatic differences in aggregation behavior, huntingtin exon1 fibrils and other polyglutamine-based aggregates contain identical β-strand-based cores. Prior structural models, derived from X-ray fiber diffraction and computational analyses, are shown to be inconsistent with the solid-state NMR results. Internally, the polyglutamine amyloid fibrils are coassembled from differently structured monomers, which we describe as a type of "intrinsic" polymorphism. A stochastic polyglutamine-specific aggregation mechanism is introduced to explain this phenomenon. We show that the aggregation of mutant huntingtin exon1 proceeds via an intramolecular collapse of the expanded polyglutamine domain and discuss the implications of this observation for our understanding of its misfolding and aggregation mechanisms.

An Exploration of the Universe of Polyglutamine Structures

Deposits of misfolded proteins in the human brain are associated with the development of many neurodegenerative diseases. Recent studies show that these proteins have common traits even at the monomer level. Among them, a polyglutamine region that is present in huntingtin is known to exhibit a correlation between the length of the chain and the severity as well as the earliness of the onset of Huntington disease. Here, we apply bias exchange molecular dynamics to generate structures of polyglutamine expansions of several lengths and characterize the resulting independent conformations. We compare the properties of these conformations to those of the standard proteins, as well as to other homopolymeric tracts. We find that, similar to the previously studied polyvaline chains, the set of possible transient folds is much broader than the set of known-to-date folds, although the conformations have different structures. We show that the mechanical stability is not related to any simple geometrical characteristics of the structures. We demonstrate that long polyglutamine expansions result in higher mechanical stability than the shorter ones. They also have a longer life span and are substantially more prone to form knotted structures. The knotted region has an average length of 35 residues, similar to the typical threshold for most polyglutamine-related diseases. Similarly, changes in shape and mechanical stability appear once the total length of the peptide exceeds this threshold of 35 glutamine residues. We suggest that knotted conformers may also harm the cellular machinery and thus lead to disease.

Misfolding and aggregation of several proteins are known to be related to neurodegenerative diseases. Among them, polyglutamine expansions are known to be responsible for at least 9 diseases, including Huntington. Nonetheless, the structural properties of these intrinsically disordered proteins are difficult to study using classical techniques because of their rapid fluctuations that result in high conformational polymorphism. Here, we use molecular dynamics simulations to study polyglutamines of different chain lengths, starting with short non-pathogenic ones, and study the independent structures they are able to form. We characterize all structures by their geometrical properties, connectivity, putative mechanical stability and residence time (life span). Similar to the findings of a previous study with polyvalines, only some of the conformers are similar to those found in natural globular proteins. Moreover, we find structures that contain knots in both polyglutamine and polyvaline 60-mers, although the former contains many more knotted conformers than the latter. We suggest that these knotted conformers may impair the cell machinery for degradation and eventually lead to toxicity.

Structural determinants of polyglutamine protofibrils and crystallites

Nine inherited neurodegenerative diseases are associated with the expansion of the CAG codon. Once the translated polyglutamine expansion becomes longer than ~36 residues, it triggers the formation of intraneural protein aggregates that often display the signature of cross-β amyloid fibrils. Here, we use fully atomistic molecular dynamics simulations to probe the structural stability and conformational dynamics of both previously proposed and new polyglutamine aggregate models. We test the relative stability of parallel and antiparallel β sheets, and characterize possible steric interfaces between neighboring sheets and the effects of different alignments of the side-chain carboxamide dipoles. Results indicate that (i) different initial oligomer structures converge to crystals consistent with available diffraction data, after undergoing cooperative side-chain rotational transitions and quarter-stagger displacements on a microsecond time scale, (ii) structures previously deemed stable on a hundred nanosecond time scale are unstable over the microsecond time scale, and (iii) conversely, structures previously deemed unstable did not account for the correct side-chain packing and once the correct symmetry is considered the structures become stable for over a microsecond, due to tightly interdigitated side chains, which lock into highly regular polar zippers with inter-side-chain and backbone-side-chain hydrogen bonds. With these insights, we built Q40 monomeric models with different combinations of arc and hairpin turns and tested them for stability. The stable monomers were further probed as a function of repeat length. Our results are consistent with the aggregation threshold. These results explain and reconcile previously reported experimental and model discrepancies about polyglutamine aggregate structures.

Unfolding and folding internal friction of β-hairpins is smaller than that of α-helices

By the forced unfolding of polyglutamine and polyalanine homopeptides in competing α-helix and β-hairpin secondary structures, we disentangle equilibrium free energetics from nonequilibrium dissipative effects. We find that α-helices are characterized by larger friction or dissipation upon unfolding, regardless of whether they are free energetically preferred over β-hairpins or not. Our analysis, based on MD simulations for atomistic peptide models with explicit water, suggests that this difference is related to the internal friction and mostly caused by the different number of intrapeptide hydrogen bonds in the α-helix and β-hairpin states.

Architecture of polyglutamine-containing fibrils from time-resolved fluorescence decay

The disease risk and age of onset of Huntington disease (HD) and nine other repeat disorders strongly depend on the expansion of CAG repeats encoding consecutive polyglutamines (polyQ) in the corresponding disease protein. PolyQ length-dependent misfolding and aggregation are the hallmarks of CAG pathologies. Despite intense effort, the overall structure of these aggregates remains poorly understood. Here, we used sensitive time-dependent fluorescent decay measurements to assess the architecture of mature fibrils of huntingtin (Htt) exon 1 implicated in HD pathology. Varying the position of the fluorescent labels in the Htt monomer with expanded 51Q (Htt51Q) and using structural models of putative fibril structures, we generated distance distributions between donors and acceptors covering all possible distances between the monomers or monomer dimensions within the polyQ amyloid fibril. Using Monte Carlo simulations, we systematically scanned all possible monomer conformations that fit the experimentally measured decay times. Monomers with four-stranded 51Q stretches organized into five-layered β-sheets with alternating N termini of the monomers perpendicular to the fibril axis gave the best fit to our data. Alternatively, the core structure of the polyQ fibrils might also be a zipper layer with antiparallel four-stranded stretches as this structure showed the next best fit. All other remaining arrangements are clearly excluded by the data. Furthermore, the assessed dimensions of the polyQ stretch of each monomer provide structural evidence for the observed polyQ length threshold in HD pathology. Our approach can be used to validate the effect of pharmacological substances that inhibit or alter amyloid growth and structure.

Structure prediction of polyglutamine disease proteins: comparison of methods

Background

The expansion of polyglutamine (poly-Q) repeats in several unrelated proteins is associated with at least ten neurodegenerative diseases. The length of the poly-Q regions plays an important role in the progression of the diseases. The number of glutamines (Q) is inversely related to the onset age of these polyglutamine diseases, and the expansion of poly-Q repeats has been associated with protein misfolding. However, very little is known about the structural changes induced by the expansion of the repeats. Computational methods can provide an alternative to determine the structure of these poly-Q proteins, but it is important to evaluate their performance before large scale prediction work is done.

Results

In this paper, two popular protein structure prediction programs, I-TASSER and Rosetta, have been used to predict the structure of the N-terminal fragment of a protein associated with Huntington's disease with 17 glutamines. Results show that both programs have the ability to find the native structures, but I-TASSER performs better for the overall task.

Conclusions

Both I-TASSER and Rosetta can be used for structure prediction of proteins with poly-Q repeats. Knowledge of poly-Q structure may significantly contribute to development of therapeutic strategies for poly-Q diseases.

Association of polyalanine and polyglutamine coiled coils mediates expansion disease-related protein aggregation and dysfunction

The expansion of homopolymeric glutamine (polyQ) or alanine (polyA) repeats in certain proteins owing to genetic mutations induces protein aggregation and toxicity, causing at least 18 human diseases. PolyQ and polyA repeats can also associate in the same proteins, but the general extent of their association in proteomes is unknown. Furthermore, the structural mechanisms by which their expansion causes disease are not well understood, and these repeats are generally thought to misfold upon expansion into aggregation-prone β-sheet structures like amyloids. However, recent evidence indicates a critical role for coiled-coil (CC) structures in triggering aggregation and toxicity of polyQ-expanded proteins, raising the possibility that polyA repeats may as well form these structures, by themselves or in association with polyQ. We found through bioinformatics screenings that polyA, polyQ and polyQA repeats have a phylogenetically graded association in human and non-human proteomes and associate/overlap with CC domains. Circular dichroism and cross-linking experiments revealed that polyA repeats can form—alone or with polyQ and polyQA—CC structures that increase in stability with polyA length, forming higher-order multimers and polymers in vitro. Using structure-guided mutagenesis, we studied the relevance of polyA CCs to the in vivo aggregation and toxicity of RUNX2—a polyQ/polyA protein associated with cleidocranial dysplasia upon polyA expansion—and found that the stability of its polyQ/polyA CC controls its aggregation, localization and toxicity. These findings indicate that, like polyQ, polyA repeats form CC structures that can trigger protein aggregation and toxicity upon expansion in human genetic diseases.

Stable polyglutamine dimers can contain β-hairpins with interdigitated side chains-but not α-helices, β-nanotubes, β-pseudohelices, or steric zippers

A common thread connecting nine fatal neurodegenerative protein aggregation diseases is an abnormally expanded polyglutamine tract found in the respective proteins. Although the structure of this tract in the large mature aggregates is increasingly well described, its structure in the small early aggregates remains largely unknown. As experimental evidence suggests that the most toxic species along the aggregation pathway are the small early ones, developing strategies to alleviate disease pathology calls for understanding the structure of polyglutamine peptides in the early stages of aggregation. Here, we present a criterion, grounded in available experimental data, that allows for using kinetic stability of dimers to assess whether a given polyglutamine conformer can be on the aggregation path. We then demonstrate that this criterion can be assessed using present-day molecular dynamics simulations. We find that although the α-helical conformer of polyglutamine is very stable, dimers of α-helices lack the kinetic stability necessary to support further oligomerization. Dimers of steric zipper, β-nanotube, and β-pseudohelix conformers are also too short-lived to initiate aggregation. The β-hairpin-containing conformers, instead, invariably form very stable dimers when their side chains are interdigitated. Combining these findings with the implications of recent solid-state NMR data on mature fibrils, we propose a possible pathway for the initial stages of polyglutamine aggregation, in which β-hairpin-containing conformers act as templates for fibril formation.

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.