Molecular origin of polyglutamine aggregation in neurodegenerative diseases

Expansion of polyglutamine (polyQ) tracts in proteins results in protein aggregation and is associated with cell death in at least nine neurodegenerative diseases. Disease age of onset is correlated with the polyQ insert length above a critical value of 35-40 glutamines. The aggregation kinetics of isolated polyQ peptides in vitro also shows a similar critical-length dependence. While recent experimental work has provided considerable insights into polyQ aggregation, the molecular mechanism of aggregation is not well understood. Here, using computer simulations of isolated polyQ peptides, we show that a mechanism of aggregation is the conformational transition in a single polyQ peptide chain from random coil to a parallel beta-helix. This transition occurs selectively in peptides longer than 37 glutamines. In the beta-helices observed in simulations, all residues adopt beta-strand backbone dihedral angles, and the polypeptide chain coils around a central helical axis with 18.5 +/- 2 residues per turn. We also find that mutant polyQ peptides with proline-glycine inserts show formation of antiparallel beta-hairpins in their ground state, in agreement with experiments. The lower stability of mutant beta-helices explains their lower aggregation rates compared to wild type. Our results provide a molecular mechanism for polyQ-mediated aggregation.     

Slow amyloid nucleation via α-helix-rich oligomeric intermediates in short polyglutamine-containing huntingtin fragments

The 17-amino-acid N-terminal segment (htt(NT)) that leads into the polyglutamine (polyQ) segment in the Huntington's disease protein huntingtin (htt) dramatically increases aggregation rates and changes the aggregation mechanism, compared to a simple polyQ peptide of similar length. With polyQ segments near or above the pathological repeat length threshold of about 37, aggregation of htt N-terminal fragments is so rapid that it is difficult to tease out mechanistic details. We describe here the use of very short polyQ repeat lengths in htt N-terminal fragments to slow this disease-associated aggregation. Although all of these peptides, in addition to htt(NT) itself, form α-helix-rich oligomeric intermediates, only peptides with Q(N) of eight or longer mature into amyloid-like aggregates, doing so by a slow increase in β-structure. Concentration-dependent circular dichroism and analytical ultracentrifugation suggest that the htt(NT) sequence, with or without added glutamine residues, exists in solution as an equilibrium between disordered monomer and α-helical tetramer. Higher order, α-helix rich oligomers appear to be built up via these tetramers. However, only htt(NT)Q(N) peptides with N=8 or more undergo conversion into polyQ β-sheet aggregates. These final amyloid-like aggregates not only feature the expected high β-sheet content but also retain an element of solvent-exposed α-helix. The α-helix-rich oligomeric intermediates appear to be both on- and off-pathway, with some oligomers serving as the pool from within which nuclei emerge, while those that fail to undergo amyloid nucleation serve as a reservoir for release of monomers to support fibril elongation. Based on a regular pattern of multimers observed in analytical ultracentrifugation, and a concentration dependence of α-helix formation in CD spectroscopy, it is likely that these oligomers assemble via a four-helix assembly unit. PolyQ expansion in these peptides appears to enhance the rates of both oligomer formation and nucleation from within the oligomer population, by structural mechanisms that remain unclear.     

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 β-rich misfolded state and their expansion length. We also find that the flanking regions of XN1 affect its probability to^x_page_count=28 form a β-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 β-rich state. Additionally, by lengthening polyQ, the first N-terminal 17 residues are more likely to adopt a β-sheet conformation rather than an α-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.     

A variational model for oligomer-formation process of GNNQQNY peptide from yeast prion protein Sup35

Many human neurodegenerative diseases are associated with the aggregation of insoluble amyloid-like fibrous proteins. However, the processes by which the randomly diffused monomer peptides aggregate into the highly regulated amyloid fibril structures are largely unknown. We proposed a residue-level coarse-grained variational model for the investigation of the aggregation pathway for a small assembly of amyloid proteins, the peptide GNNQQNY from yeast prion protein Sup35. By examining the free energy surface, we identified the residue-level sequential pathways for double parallel and antiparallel β-peptides, which show that the central dry polar zipper structure is the major folding core in both cases. The critical nucleus size is determined to be three peptides for the homogeneous nucleation process, whereas the zig-zag growth pattern appears most favorably for heterogeneous nucleation. Consistent with the dock-and-lock mechanism, the aggregation process of free peptides to the fibril core was found to be highly cooperative. The quantitative validation with the computational simulations and experimental data demonstrated the usefulness of the proposed model in understanding the general mechanism of the amyloid fibril system.     

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.     

Toxic prefibrillar α-synuclein amyloid oligomers adopt a distinctive antiparallel β-sheet structure

Parkinson's disease is an age-related movement disorder characterized by the presence in the mid-brain of amyloid deposits of the 140-amino-acid protein AS (α-synuclein). AS fibrillation follows a nucleation polymerization pathway involving diverse transient prefibrillar species varying in size and morphology. Similar to other neurodegenerative diseases, cytotoxicity is currently attributed to these prefibrillar species rather than to the insoluble aggregates. Nevertheless, the underlying molecular mechanisms responsible for cytotoxicity remain elusive and structural studies may contribute to the understanding of both the amyloid aggregation mechanism and oligomer-induced toxicity. It is already recognized that soluble oligomeric AS species adopt β-sheet structures that differ from those characterizing the fibrillar structure. In the present study we used ATR (attenuated total reflection)-FTIR (Fourier-transform infrared) spectroscopy, a technique especially sensitive to β-sheet structure, to get a deeper insight into the β-sheet organization within oligomers and fibrils. Careful spectral analysis revealed that AS oligomers adopt an antiparallel β-sheet structure, whereas fibrils adopt a parallel arrangement. The results are discussed in terms of regions of the protein involved in the early β-sheet interactions and the implications of such conformational arrangement for the pathogenicity associated with AS oligomers.     

Structures and dynamics of β-barrel oligomer intermediates of amyloid-beta16-22 aggregation

Accumulating evidence suggests that soluble oligomers are more toxic than final fibrils of amyloid aggregations. Among the mixture of inter-converting intermediates with continuous distribution of sizes and secondary structures, oligomers in the β-barrel conformation – a common class of protein folds with a closed β-sheet – have been postulated as the toxic species with well-defined three-dimensional structures to perform pathological functions. A common mechanism for amyloid toxicity, therefore, implies that all amyloid peptides should be able to form β-barrel oligomers as the aggregation intermediates. Here, we applied all-atom discrete molecular dynamics (DMD) simulations to evaluate the formation of β-barrel oligomers and characterize their structures and dynamics in the aggregation of a seven-residue amyloid peptide, corresponding to the amyloid core of amyloid-β with a sequence of ¹⁶KLVFFAE²² (Aβ16-22). We carried out aggregation simulations with various numbers of peptides to study the size dependence of aggregation dynamics and assembly structures. Consistent with previous computational studies, we observed the formation of β-barrel oligomers in all-atom DMD simulations. Using a network-based approach to automatically identify β-barrel conformations, we systematically characterized β-barrels of various sizes. Our simulations revealed the conformational inter-conversion between β-barrels and double-layer β-sheets due to increased structural strains upon forming a closed β-barrel while maximizing backbone hydrogen bonds. The potential of mean force analysis further characterized the free energy barriers between these two states. The obtained structural and dynamic insights of β-barrel oligomers may help better understand the molecular mechanism of oligomer toxicities and design novel therapeutics targeting the toxic β-barrel oligomers. This article is part of a Special Issue entitled: Protein Aggregation and Misfolding at the Cell Membrane Interface edited by Ayyalusamy Ramamoorthy     

d-Polyglutamine Amyloid Recruits l-Polyglutamine Monomers and Kills Cells

Polyglutamine (polyQ) amyloid fibrils are observed in disease tissue and have been implicated as toxic agents responsible for neurodegeneration in expanded CAG repeat diseases such as Huntington's disease. Despite intensive efforts, the mechanism of amyloid toxicity remains unknown. As a novel approach to probing polyQ toxicity, we investigate here how some cellular and physical properties of polyQ amyloid vary with the chirality of the glutamine residues in the polyQ. We challenged PC12 cells with small amyloid fibrils composed of either l- or d-polyQ peptides and found that d-fibrils are as cytotoxic as l-fibrils. We also found using fluorescence microscopy that both aggregates effectively seed the aggregation of cell-produced l-polyQ proteins, suggesting a surprising lack of stereochemical restriction in seeded elongation of polyQ amyloid. To investigate this effect further, we studied chemically synthesized d- and l-polyQ in vitro. We found that, as expected, d-polyQ monomers are not recognized by proteins that recognize l-polyQ monomers. However, amyloid fibrils prepared from d-polyQ peptides can efficiently seed the aggregation of l-polyQ monomers in vitro, and vice versa. This result is consistent with our cell results on polyQ recruitment but is inconsistent with previous literature reports on the chiral specificity of amyloid seeding. This chiral cross-seeding can be rationalized by a model for seeded elongation featuring a “rippled β-sheet” interface between seed fibril and docked monomers of opposite chirality. The lack of chiral discrimination in polyQ amyloid cytotoxicity is consistent with several toxicity mechanisms, including recruitment of cellular polyQ proteins.     

Sequence-independent control of peptide conformation in liposomal vaccines for targeting protein misfolding diseases

Synthetic peptide immunogens that mimic the conformation of a target epitope of pathological relevance offer the possibility to precisely control the immune response specificity. Here, we performed conformational analyses using a panel of peptides in order to investigate the key parameters controlling their conformation upon integration into liposomal bilayers. These revealed that the peptide lipidation pattern, the lipid anchor chain length, and the liposome surface charge all significantly alter peptide conformation. Peptide aggregation could also be modulated post-liposome assembly by the addition of distinct small molecule β-sheet breakers. Immunization of both mice and monkeys with a model liposomal vaccine containing β-sheet aggregated lipopeptide (Palm1-15) induced polyclonal IgG antibodies that specifically recognized β-sheet multimers over monomer or non-pathological native protein. The rational design of liposome-bound peptide immunogens with defined conformation opens up the possibility to generate vaccines against a range of protein misfolding diseases, such as Alzheimer disease.     

Amyloid Evolution: Antiparallel Replaced by Parallel

Several atomic structures have now been found for micrometer-scale amyloid fibrils or elongated microcrystals using a range of methods, including NMR, electron microscopy, and X-ray crystallography, with parallel β-sheet appearing as the most common secondary structure. The etiology of amyloid disease, however, indicates nanometer-scale assemblies of only tens of peptides as significant agents of cytotoxicity and contagion. By combining solution X-ray with molecular dynamics, we show that antiparallel structure dominates at the first stages of aggregation for a specific set of peptides, being replaced by parallel at large length scales only. This divergence in structure between small and large amyloid aggregates should inform future design of molecular therapeutics against nucleation or intercellular transmission of amyloid. Calculations and an overview from the literature argue that antiparallel order should be the first appearance of structure in many or most amyloid aggregation processes, regardless of the endpoint. Exceptions to this finding should exist, depending inevitably on the sequence and on solution conditions.     

Dimers initiate and propagate serine protease inhibitor polymerisation

The serine protease inhibitor (serpin) family can readily form long-chain polymers by a process that underlies a variety of diseases. We show here that monomers of plasma serpins α₁-antitrypsin and antithrombin are stable on incubation with the rate-limiting step in their polymerisation being the formation of the initial dimer. Once formed, the dimers readily interlink to form tetramers and can bind monomers to form trimers and longer oligomers. Cleavage of the only exposed reactive loop, in unit I of the dimers, prevents their interlinkage, but these cleaved dimers can still link to monomers. The rapid binding by the cleaved dimers of a peptide specific to the lower half of β-sheet A of the molecule indicates the ready opening of this β-sheet in unit II of the dimers. The failure of the cleaved dimers to bind peptide-complexed monomers, together with the relative inaccessibility of the P14 hinge residue in the oligomers, is evidence that partial insertion of the reactive loop into its own A-sheet is required for polymer formation. We propose that serpin dimers initiate and propagate polymerisation by having one exposed loop with an optimal conformation as a β-strand donor and a readily opened β-sheet as an acceptor. The sequential reformation of these activated β-interfaces as the oligomer extends, molecule by molecule, provides a model for the fibril and amyloid formation of conformational diseases in general as well as for the infectivity of prion encephalopathies.     

β-Barrel topology of Alzheimer's β-amyloid ion channels

Emerging evidence supports the ion channel mechanism for Alzheimer's disease pathophysiology wherein small β-amyloid (Aβ) oligomers insert into the cell membrane, forming toxic ion channels and destabilizing the cellular ionic homeostasis. Solid-state NMR-based data of amyloid oligomers in solution indicate that they consist of a double-layered β-sheets where each monomer folds into β-strand-turn-β-strand and the monomers are stacked atop each other. In the membrane, Aβ peptides are proposed to be β-type structures. Experimental structural data available from atomic force microscopy (AFM) imaging of Aβ oligomers in membranes reveal heterogeneous channel morphologies. Previously, we modeled the channels in a non-tilted organization, parallel with the cross-membrane normal. Here, we modeled a β-barrel-like organization. β-Barrels are common in transmembrane toxin pores, typically consisting of a monomeric chain forming a pore, organized in a single-layered β-sheet with antiparallel β-strands and a right-handed twist. Our explicit solvent molecular dynamics simulations of a range of channel sizes and polymorphic turns and comparisons of these with AFM image dimensions support a β-barrel channel organization. Different from the transmembrane β-barrels where the monomers are folded into a circular β-sheet with antiparallel β-strands stabilized by the connecting loops, these Aβ barrels consist of multimeric chains forming double β-sheets with parallel β-strands, where the strands of each monomer are connected by a turn. Although the Aβ barrels adopt the right-handed β-sheet twist, the barrels still break into heterogeneous, loosely attached subunits, in good agreement with AFM images and previous modeling. The subunits appear mobile, allowing unregulated, hence toxic, ion flux.     

A hIAPP-derived all-d-amino-acid inhibits hIAPP fibrillation efficiently at membrane surface by targeting α-helical oligomeric intermediates

A variety of peptides and peptide derivatives have been constructed using the “β-sheet core segment” of amyloid proteins as inhibitors of amyloidogenic fibrillation. A novel all-d-amino-acid from hIAPP β-sheet core segment (hIAPP 22–27) is demonstrated to inhibit hIAPP fibril formation efficiently both at the phospholipid membrane and in bulk solution. The inhibitor terminates hIAPP aggregation to the α-helical oligomeric intermediates at the membrane surface, whereas it stops the aggregation at the stage of β-sheet oligomeric intermediates in bulk solution. This is the first evidence that the inhibition mechanism of the inhibitor at membrane surface is significantly different from that in bulk solution.     

The possible structural models for polyglutamine aggregation: a molecular dynamics simulations study

Huntington's disease is a neurodegenerative disorder caused by a polyglutamine (polyQ) expansion near the N-terminus of huntingtin. Previous studies have suggested that polyQ aggregation occurs only when the number of glutamine (Q) residues is more than 36-40, the disease threshold. However, the structural characteristics of polyQ nucleation in the very early stage of aggregation still remain elusive. In this study, we designed 18 simulation trials to determine the possible structural models for polyQ nucleation and aggregation with various shapes and sizes of initial β-helical structures, such as left-handed circular, right-handed rectangular, and left- and right-handed triangular. Our results show that the stability of these models significantly increases with increasing the number of rungs, while it is rather insensitive to the number of Qs in each rung. In particular, the 3-rung β-helical models are stable when they adopt the left-handed triangular and right-handed rectangular conformations due to the fact that they preserve high β-turn and β-sheet contents, respectively, during the simulation courses. Thus, we suggested that these two stable β-helical structures with at least 3 rungs might serve as the possible nucleation seeds for polyQ depending on how the structural elements of β-turn and β-sheet are sampled and preserved during the very early stage of aggregation.     

An Analysis of Biomolecular Force Fields for Simulations of Polyglutamine in Solution

Polyglutamine (polyQ) peptides are a useful model system for biophysical studies of protein folding and aggregation, both for their intriguing aggregation properties and their own relevance to human disease. The genetic expansion of a polyQ tract triggers the formation of amyloid aggregates associated with nine neurodegenerative diseases. Several clearly identifiable and separable factors, notably the length of the polyQ tract, influence the mechanism of aggregation, its associated kinetics, and the ensemble of structures formed. Atomistic simulations are well positioned to answer open questions regarding the thermodynamics and kinetics of polyQ folding and aggregation. The additional, explicit representation of water permits deeper investigation of the role of solvent dynamics, and it permits a direct comparison of simulation results with infrared spectroscopy experiments. The generation of meaningful simulation results hinges on satisfying two essential criteria: achieving sufficient conformational sampling to draw statistically valid conclusions, and accurately reproducing the intermolecular forces that govern system structure and dynamics. In this work, we examine the ability of 12 biomolecular force fields to reproduce the properties of a simple, 30-residue polyQ peptide (Q₃₀) in explicit water. In addition to secondary and tertiary structure, we consider generic structural properties of polymers that provide additional dimensions for analysis of the highly degenerate disordered states of the molecule. We find that the 12 force fields produce a wide range of predictions. We identify AMBER ff99SB, AMBER ff99SB^∗, and OPLS-AA/L to be most suitable for studies of polyQ folding and aggregation.     

Elongation kinetics of polyglutamine peptide fibrils: a quartz crystal microbalance with dissipation study

Abnormally expanded polyglutamine domains in proteins are associated with several neurodegenerative diseases, including Huntington's disease. Expansion of the polyglutamine (polyQ) domain facilitates aggregation of the affected protein, and several studies directly link aggregation to neurotoxicity. Studies of synthetic polyQ peptides have contributed substantially to our understanding of the mechanism of aggregation. In this report, polyQ fibrils were immobilized onto a sensor, and their elongation by polyQ peptides of various length and conformation was examined using quartz crystal microbalance with dissipation monitoring (QCM-D). The rate of elongation increased as the peptide length increased from 8 to 24 glutamines (Q8, Q20, and Q24). Monomer conformation affected elongation rates: insertion of a β-turn template d-Pro-Gly in the center of the peptide increased elongation rates several-fold, while insertion of Pro-Pro dramatically slowed elongation. Dissipation measurements of the QCM-D provided qualitative information about mechanical properties of the elongating fibrils. These data showed clear differences in the characteristics of the elongating aggregates, depending on the specific identity of the associating polyQ peptide. Elongation rates were sensitive to the pH and ionic strength of the buffer. Comparison of QCM-D data with those obtained by optical waveguide lightmode spectroscopy revealed that very little water was associated with the elongation of fibrils by the peptide containing d-Pro-Gly, but a significant amount of water was associated when the fibrils were elongated by Q20. Together, the data indicate that elongation of polyQ fibrils can occur without full consolidation to the fibril structure, resulting in variations to the aggregate structure during elongation.     

Nucleation-dependent Aggregation Kinetics of Yeast Sup35 Fragment GNNQQNY

An N-terminal hepta-peptide sequence of yeast prion protein Sup35 with the sequence GNNQQNY is widely used as a model system for amyloid fibril formation. In this study, we used a reproducible solubilisation protocol that allows the generation of a homogenous monomeric solution of GNNQQNY to uncover the molecular details of its self-assembly mechanism. The aggregation kinetics data show that the GNNQQNY sequence follows nucleation-dependent aggregation kinetics with a critical nucleus of size ~7 monomers and that the efficiency of nucleation were found to be inversely related to the reaction temperature. The nucleus reduces the thermodynamic energy barrier by acting as a template for further self-assembly and results in highly ordered amyloid fibrils. The fibers grown at different temperatures showed similar Thioflavin T fluorescence, Congo-red binding and β-sheet rich structures displaying a characteristic cross-β diffraction pattern. These aggregates also share morphological and structural identity with those reported earlier. The mature GNNQQNY fibers did not exert significant oxidative stress or cytotoxicity upon incubating with differentiated SHSY5Y cells. To our knowledge, this is the first study to experimentally validate previous nucleus size predictions based on theoretical and molecular dynamics simulations. These findings provide the basis for understanding the kinetics and thermodynamics of amyloid nucleation and elongation of amyloidogenic proteins/peptides associated with many systemic and neurodegenerative diseases.     

Secondary structure propensity and chirality of the amyloidophilic peptide p5 and its analogues impacts ligand binding - In vitro characterization

BackgroundPolybasic helical peptides, such as peptide p5, bind human amyloid extracts and synthetic amyloid fibrils. When radiolabeled, peptide p5 has been shown to specifically bind amyloid in vivo thereby allowing imaging of the disease. Structural requirements for heparin and amyloid binding have been studied using analogues of p5 that modify helicity and chirality.MethodsPeptide-ligand interactions were studied using CD spectroscopy and solution-phase binding assays with radiolabeled p5 analogues. The interaction of a subset of peptides was further studied by using molecular dynamics simulations.ResultsDisruption of the peptide helical structure reduced peptide binding to heparin and human amyloid extracts. The all-D enantiomer and the β-sheet-structured peptide bound all substrates as well as, or better than, p5. The interaction of helical and β-sheet structured peptides with Aβ fibrils was modeled and shown to involve both ionic and non-ionic interactions.ConclusionsThe α-helical secondary structure of peptide p5 is important for heparin and amyloid binding; however, helicity is not an absolute requirement as evidenced by the superior reactivity of a β-sheet peptide. The differential binding of the peptides with heparin and amyloid fibrils suggests that these molecular interactions are different. The all-D enantiomer of p5 and the β-sheet peptide are candidates for amyloid targeting reagents in vivo.General SignificanceEfficient binding of polybasic peptides with amyloid is dependent on the linearity of charge spacing in the context of an α-helical secondary structure. Peptides with an α-helix or β-sheet propensity and with similar alignment of basic residues is optimal.     

Comparison of molecular dynamics simulation methods for amyloid β(1-42) monomers containing D-aspartic acid residues for predicting retention times in chromatography

Molecular dynamics simulations of amyloid β(1-42) containing D-aspartic acid residues were performed using several continuous solvent models to investigate the usefulness of simulation methods for D-amino acid-containing proteins and peptides. Normal molecular dynamics simulations and replica exchange molecular dynamics simulations, which are one of the generalized-ensemble algorithms, were performed. Because the β-structure contents of amyloid β(1-42) peptides obtained by replica exchange molecular dynamics simulations with Onufriev-Bashford-Case generalized Born implicit solvent were qualitatively consistent with experimental data, replica exchange molecular dynamics rather than other methods appeared to be more reasonable for calculations of amyloid β(1-42) containing D-aspartic acid residues. Computational results revealed that peptides with stereoinversion of Asp23 tend to form β-sheet structures by themselves, in contrast to the wild-type peptides that form β-sheet structures only after aggregation. These results are expected to be useful for computational investigations of proteins and peptides such as prediction of retention time of peptides and proteins containing D-aspartic acid residues.     

The Possible Structural Models for Polyglutamine Aggregation: A Molecular Dynamics Simulations Study

Abstract

Huntington's disease is a neurodegenerative disorder caused by a polyglutamine (polyQ) expansion near the N-terminus of huntingtin. Previous studies have suggested that polyQ aggregation occurs only when the number of glutamine (Q) residues is more than 36-40, the disease threshold. However, the structural characteristics of polyQ nucleation in the very early stage of aggregation still remain elusive. In this study, we designed 18 simulation trials to determine the possible structural models for polyQ nucleation and aggregation with various shapes and sizes of initial β-helical structures, such as left-handed circular, right-handed rectangular, and left- and right-handed triangular. Our results show that the stability of these models significantly increases with increasing the number of rungs, while it is rather insensitive to the number of Qs in each rung. In particular, the 3-rung β-helical models are stable when they adopt the left-handed triangular and right-handed rectangular conformations due to the fact that they preserve high β-turn and β-sheet contents, respectively, during the simulation courses. Thus, we suggested that these two stable β-helical structures with at least 3 rungs might serve as the possible nucleation seeds for polyQ depending on how the structural elements of β-turn and β-sheet are sampled and preserved during the very early stage of aggregation.