Molecular dynamics simulations on the oligomer-formation process of the GNNQQNY peptide from yeast prion protein Sup35

Oligomeric intermediates are possible cytotoxic species in diseases associated with amyloid deposits. Understanding the early steps of fibril formation at atomic details may provide useful information for the rational therapeutic design. In this study, using the heptapeptide GNNQQNY from the yeast prion-like protein Sup35 as a model system, for which a detailed atomic structure of the fibril formed has been determined by x-ray microcrystallography, we investigated its oligomer-formation process from monomer to tetramer at the atomistic level by means of a molecular dynamics simulation with explicit water. Although the number of simulations was limited, the qualitative statistical data gave some interesting results, which indicated that the oligomer formation might start from antiparallel beta-sheet-like dimers. When a new single peptide strand was added to the preformed dimers to form trimers and then tetramers, the transition time from disorder aggregates to regular ones for the parallel alignment was found to be obviously much less than for the antiparallel one. Moreover, the parallel pattern also statistically stayed longer, providing more chances for oligomer extending, although the number of parallel stack events was almost equal to antiparallel ones. Therefore, our simulations showed that new strands might prefer to extend in a parallel arrangement to form oligomers, which agrees with the microcrystal structure of the amyloid fibril formed by this peptide. In addition, analysis of the pi-pi stacking of aromatic residues showed that this type of interaction did not play an important role in giving directionality for beta-strand alignment but played a great influence on stabilizing the structures formed in the oligomer-formation process.     

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.     

Protein misfolding and amyloid formation for the peptide GNNQQNY from yeast prion protein Sup35: simulation by reaction path annealing

We study the early steps of amyloid formation of the seven residue peptide GNNQQNY from yeast prion-like protein Sup35 by simulating the random coil to beta-sheet and alpha-helix to beta-sheet transition both in the absence and presence of a cross-beta amyloid nucleus. The simulation method at atomic resolution employs a new implementation of a Langevin dynamics "reaction path annealing" algorithm. The results indicate that the presence of amyloid-like cross-beta-sheet strands both facilitates the transition into the cross-beta conformation and substantially lowers the free energy barrier for this transition. This model systems allows us to investigate the energetic and kinetic details of this transition, which is consistent with an auto-catalyzed, nucleation-like mechanism for the formation of beta-amyloid. In particular, we find that electrostatic interactions of peptide backbone dipoles contribute significantly to the stability of the beta-amyloid state. Furthermore, we find water exclusion and interactions of polar side-chains to be driving forces of amyloid formation: the cross-beta conformation is stabilized by burial of polar side-chains and inter-residue hydrogen bonds in the presence of an amyloid-like "seed". These findings are in support of a "dry, polar zipper model" of amyloid formation.     

Molecular dynamics simulations to investigate the stability and aggregation behaviour of the amyloid-forming peptide VQIVYK from tau protein

The formation of paired helical filaments arising from the short hexapeptide in the third repeat of tau protein, ³⁰⁶VQIVYK³¹¹, is critical for tau polymerisation. The atomic structure of the VQIVYK oligomer has revealed a dry, tightly self-complementing structure between the neighbouring β-sheet layers, termed as ‘steric zipper’. In this study, several molecular dynamics simulations with all-atom explicit water were conducted to investigate the structural stability and aggregation behaviour of the VQIVYK peptide with various sizes and its single alanine replacement mutations. Our results indicate that the van der Waals interaction between side chains of Q2, the π–π stacking interaction between aromatic rings of Y5, and the electrostatic interaction between K6 and the C-terminus play an important role in stabilising the VQIVYK oligomers within the same β-sheet layer, while hydrophobic steric zipper involving V1, I3 and Y5 is responsible for holding the neighbouring β-sheet layers together. The twisted angles of the VQIVYK oligomers were also analysed and shown to be size dependent. The present results not only provide atomic insights into amyloid formation, but are also helpful for designing new or modified capping peptides and inhibitors to prevent fibril formation of the VQIVYK peptide from tau protein.     

Flexible docking of an amyloid-forming peptide from beta(2)-microglobulin

Using an all-atom, molecular dynamics-based, flexible docking method, the tertiary and quaternary structures of protofilaments of the "K3" fragment from beta(2)-microglobulin (residues Ser20-Lys41) were predicted at low pH in a continuous mixture of water and 2,2,2-trifluoroethanol (TFE). Tetramers with energies very close to the global minimum were produced with C(alpha) root-mean square deviation values under 4A over 88 residues compared to a recently solved SSNMR structure. The most accurate model distinguishes itself from other low-energy solutions in that it shows high structural similarity to another known fold, the parallel beta-helix, in agreement with models proposed previously by several other groups. The method achieves efficiency without loss of generality or atomic detail by enforcing local symmetry on the individual peptides, rewarding intermolecular contacts, and iteratively building up the protofilaments by successively doubling the number of chains. Solvent effects were included in the model by treating the dielectric constant and surface tension as functions of the TFE concentration. In order to understand the physical basis for the stabilizing effects of TFE, the TFE concentration was varied from 0% to 50% (v/v) and a peak in stability was observed at 16%, where the polar and hydrophobic terms cancel out and close to the experimentally determined value of 20%.     

Cross-beta order and diversity in nanocrystals of an amyloid-forming peptide

The seven-residue peptide GNNQQNY from the N-terminal region of the yeast prion protein Sup35, which forms amyloid fibers, colloidal aggregates and highly ordered nanocrystals, provides a model system for characterizing the elusively protean cross-beta conformation. Depending on preparative conditions, orthorhombic and monoclinic crystals with similar lath-shaped morphology have been obtained. Ultra high-resolution (<0.5A spacing) electron diffraction patterns from single nanocrystals show that the peptide chains pack in parallel cross-beta columns with approximately 4.86A axial spacing. Mosaic striations 20-50 nm wide observed by electron microscopy indicate lateral size-limiting crystal growth related to amyloid fiber formation. Frequently obtained orthorhombic forms, with apparent space group symmetry P2(1)2(1)2(1), have cell dimensions ranging from /a/=22.7-21.2A, /b/=39.9-39.3A, /c/=4.89-4.86A for wet to dried states. Electron diffraction data from single nanocrystals, recorded in tilt series of still frames, have been mapped in reciprocal space. However, reliable integrated intensities cannot be obtained from these series, and dynamical electron diffraction effects present problems in data analysis. The diversity of ordered structures formed under similar conditions has made it difficult to obtain reproducible X-ray diffraction data from powder specimens; and overlapping Bragg reflections in the powder patterns preclude separated structure factor measurements for these data. Model protofilaments, consisting of tightly paired, half-staggered beta strands related by a screw axis, can be fit in the crystal lattices, but model refinement will require accurate structure factor measurements. Nearly anhydrous packing of this hydrophilic peptide can account for the insolubility of the crystals, since the activation energy for rehydration may be extremely high. Water-excluding packing of paired cross-beta peptide segments in thin protofilaments may be characteristic of the wide variety of anomalously stable amyloid aggregates.     

Alternative assembly pathways of the amyloidogenic yeast prion determinant Sup35-NM

The self-perpetuating conformational change of the translation termination factor Sup35 is associated with a prion phenomenon of Saccharomyces cerevisiae. In vitro, the prion-determining region (NM) of Sup35 assembles into amyloid-like fibres through a mechanism of nucleated conformational conversion. Here, we describe an alternative assembly pathway of NM that produces filaments that are composed of beta-strands and random coiled regions with several-fold smaller diameters than the amyloid fibres. NM filaments are not detectable with either thioflavin T or Congo Red and do not show SDS or protease resistance. As filaments do not self-convert into fibres and do not act as seed, they are not intermediates of amyloid fibre formation. Instead, they represent a stable off-pathway form. Similar to mammalian prion proteins, Sup35 contains oligopeptide repeats located in the NM region. We found that the number of repeats determines the partitioning of the protein between filaments and amyloid-like fibres. Low numbers of repeats favour the formation of the filamentous structure, whereas high numbers of repeats favour the formation of amyloid-like fibres.     

Characterizing the Assembly of the Sup35 Yeast Prion Fragment, GNNQQNY: Structural Changes Accompany a Fiber-to-Crystal Switch

Amyloid-like fibrils can be formed by many different proteins and peptides. The structural characteristics of these fibers are very similar to those of amyloid fibrils that are deposited in a number of protein misfolding diseases, including Alzheimer's disease and the transmissible spongiform encephalopathies. The elucidation of two crystal structures from an amyloid-like fibril-forming fragment of the yeast prion, Sup35, with sequence GNNQQNY, has contributed to knowledge regarding side-chain packing of amyloid-forming peptides. Both structures share a cross-β steric zipper arrangement but vary in the packing of the peptide, particularly in terms of the tyrosine residue. We investigated the fibrillar and crystalline structure and assembly of the GNNQQNY peptide using x-ray fiber diffraction, electron microscopy, intrinsic and quenched tyrosine fluorescence, and linear dichroism. Electron micrographs reveal that at concentrations between 0.5 and 10 mg/mL, fibers form initially, followed by crystals. Fluorescence studies suggest that the environment of the tyrosine residue changes as crystals form. This is corroborated by linear dichroism experiments that indicate a change in the orientation of the tyrosine residue over time, which suggests that a structural rearrangement occurs as the crystals form. Experimental x-ray diffraction patterns from fibers and crystals also suggest that these species are structurally distinct. A comparison of experimental and calculated diffraction patterns contributes to an understanding of the different arrangements accessed by the peptide.     

Origins and kinetic consequences of diversity in Sup35 yeast prion fibers

A remarkable feature of prions is that infectious particles composed of the same prion protein can give rise to different phenotypes. This strain phenomenon suggests that a single prion protein can adopt multiple infectious conformations. Here we use a novel single fiber growth assay to examine the heterogeneity of amyloid fibers formed by the yeast Sup35 prion protein. Sup35 spontaneously forms multiple, distinct and faithfully propagating fiber types, which differ dramatically both in their degrees of polarity and overall growth rates. Both in terms of the number of distinct self-propagating fiber types, as well as the ability of these differences to dictate the rate of prion growth, this diversity is well suited to account for the range of prion strain phenotypes observed in vivo.     

Structure analysis of an amyloid-forming model peptide by a systematic glycine and proline scan

The ability to adopt at least two different stable conformations is a common feature of proteins involved in many neurodegenerative diseases. The involved molecules undergo a conformational transition from native, mainly helical states to insoluble amyloid structures that have high β-sheet content. A detailed characterization of the molecular architecture of highly ordered amyloid structures, however, is still challenging. Their intrinsically low solubility and high tendency to aggregate often considerably limits the application of established high-resolution techniques such as NMR and X-ray crystallography. An alternative approach to elucidating the tertiary and quaternary organization within an amyloid fibril is the systematic replacement of residues with amino acids that exhibit special conformational characteristics, such as glycine and proline. Substitutions within the β-sheet-prone sequences of the molecules usually severely affect their ability to form fibrils, whereas incorporation at external loop- and bend-like positions often has only marginal effects. Here we present the characterization of the internal architecture of a de novo designed coiled-coil-based amyloid-forming model peptide by means of a series of systematic single glycine and proline replacements in combination with a set of simple low-resolution methods. The folding and assembly behavior of the substituted peptides was monitored simultaneously using circular dichroism spectroscopy, Thioflavin T fluorescence staining, and transmission electron microscopy. On the basis of the obtained data, we successfully identify characteristic bend and core positions within the peptide sequence and propose a detailed structural model of the internal fibrillar arrangement.     

Suppression of polyglutamine toxicity by the yeast Sup35 prion domain in Drosophila

The propensity of proteins to form beta-sheet-rich amyloid fibrils is related to a variety of biological phenomena, including a number of human neurodegenerative diseases and prions. A subset of amyloidogenic proteins forms amyloid fibrils through glutamine/asparagine (Q/N)-rich domains, such as pathogenic polyglutamine (poly(Q)) proteins involved in neurodegenerative disease, as well as yeast prions. In the former, the propensity of an expanded poly(Q) tract to abnormally fold confers toxicity on the respective protein, leading to neuronal dysfunction. In the latter, Q/N-rich prion domains mediate protein aggregation important for epigenetic regulation. Here, we investigated the relationship between the pathogenic ataxin-3 protein of the human disease spinocerebellar ataxia type 3 (SCA3) and the yeast prion Sup35, using Drosophila as a model system. We found that the capacity of the Sup35 prion domain to mediate protein aggregation is conserved in Drosophila. Although select yeast prions enhance poly(Q) toxicity in yeast, the Sup35N prion domain suppressed poly(Q) toxicity in the fly. Suppression required the oligopeptide repeat of the Sup35N prion domain, which is critical for prion properties in yeast. These results suggest a trans effect of prion domains on pathogenic poly(Q) disease proteins in a multicellular environment and raise the possibility that Drosophila may allow studies of prion mechanisms.     

The role of the N-terminal oligopeptide repeats of the yeast Sup35 prion protein in propagation and transmission of prion variants

The cytoplasmic [PSI+] determinant of Saccharomyces cerevisiae is the prion form of the Sup35 protein. Oligopeptide repeats within the Sup35 N-terminal domain (PrD) presumably are required for the stable [PSI+] inheritance that in turn involves fragmentation of Sup35 polymers by the chaperone Hsp104. The nonsense suppressor [PSI+] phenotype can vary in efficiency probably due to different inheritable Sup35 polymer structures. Here we study the ability of Sup35 mutants with various deletions of the oligopeptide repeats to support [PSI+] propagation. We define the minimal region of the Sup35-PrD necessary to support [PSI+] as amino acids 1-64, which include the first two repeats, although a longer fragment, 1-83, is required to maintain weak [PSI+] variants. Replacement of wild-type Sup35 with deletion mutants decreases the strength of the [PSI+] phenotype. However, with one exception, reintroducing the wild-type Sup35 restores the original phenotype. Thus, the specific prion fold defining the [PSI+] variant can be preserved by the mutant Sup35 protein despite the change of phenotype. Coexpression of wild-type and mutant Sup35 containing three, two, one, or no oligopeptide repeats causes variant-specific [PSI+] elimination. These data suggest that [PSI+] variability is primarily defined by differential folding of the Sup35-PrD oligopeptide-repeat region.     

Oligopeptide repeats in the yeast protein Sup35p stabilize intermolecular prion interactions

The nuclear-encoded Sup35p protein is responsible for the prion-like [PSI(+)] determinant of yeast, with Sup35p existing largely as a high molecular weight aggregate in [PSI(+)] strains. Here we show that the five oligopeptide repeats present at the N-terminus of Sup35p are responsible for stabilizing aggregation of Sup35p in vivo. Sequential deletion of the oligopeptide repeats prevented the maintenance of [PSI(+)] by the truncated Sup35p, although deletants containing only two repeats could be incorporated into pre-existing aggregates of wild-type Sup35p. The mammalian prion protein PrP also contains similar oligopeptide repeats and we show here that a human PrP repeat (PHGGGWGQ) is able functionally to replace a Sup35p oligopeptide repeat to allow stable [PSI(+)] propagation in vivo. Our data suggest a model in which the oligopeptide repeats in Sup35p stabilize intermolecular interactions between Sup35p proteins that initiate establishment of the aggregated state. Modulating repeat number therefore alters the rate of yeast prion conversion in vivo. Furthermore, there appears to be evolutionary conservation of function of the N-terminally located oligopeptide repeats in prion propagation.     

The actin cytoskeletal network plays a role in yeast prion transmission and contributes to prion stability

Chaperone networks are required for the shearing and generation of transmissible propagons from pre-existing prion aggregates. However, other cellular networks needed for maintaining yeast prions are largely uncharacterized. Here, we establish a novel role for actin networks in prion maintenance. The [PIN⁺] prion, also known as [RNQ⁺], exists as stable variants dependent upon the chaperone machinery for the transmission of propagons to daughter cells during cell division and cytoplasmic transfer. Loss of the Hsp104 molecular chaperone leads to the growth of prion particles until they are too large to be transmitted. Here, we isolated a unique [PIN⁺] variant, which is unstable in actin mutants. This prion loss is observed over many generations, and coincides with the detection of both high molecular weight species of Rnq1 and large visible aggregates that are asymmetrically retained during cell division. Our data suggest that the irregular actin networks found in these mutants may influence propagon number by slowly permitting aggregate growth over time, resulting in the generation of nontransmissible large aggregates. Thus, we show the potential contribution of cytoskeletal networks in the transmission of prion propagons, which parallels models that have been proposed for cell-to-cell transmission of small amyloids in neurodegenerative protein aggregation diseases.     

beta-Helix is a likely core structure of yeast prion Sup35 amyloid fibers

Helicobacter hepaticus, a causal agent of hepatocarcinoma in mice, exhibits a cytolethal distending toxin activity. The three subunits of this holotoxin, CdtA, CdtB, and CdtC, and three CdtB mutants were produced as recombinant histidine-tagged proteins by using an in vitro cell-free protein expression system. We found that the presence of the three H. hepaticus Cdt subunits is required for cellular toxicity and that only a C-terminal CdtB mutation abolishes the activity of the complex. In vitro, H. hepaticus CdtB exhibits a DNase activity which is also abolished by this C-terminal CdtB mutation. These results suggest that the effect of H. hepaticus CDT probably involves the DNase activity of CdtB.     

Qualitative and quantitative multiplexed proteomic analysis of complex yeast protein fractions that modulate the assembly of the yeast prion Sup35p

Background

The aggregation of the baker''s yeast prion Sup35p is at the origin of the transmissible [PSI⁺] trait. We and others have shown that molecular chaperones modulate Sup35p aggregation. However, other protein classes might be involved in [PSI⁺] formation.

Results

We designed a functional proteomic study that combines two techniques to identify modulators of Sup35p aggregation and describe the changes associated to [PSI⁺] formation. The first allows measuring the effect of fractionated Saccharomyces cerevisiae cytosolic extracts from [PSI⁺] and [psi⁻] yeast cells on Sup35p assembly. The second is a multiplex qualitative and quantitative comparison of protein composition of active and inactive fractions using a gel-free and label-free LC-MS approach. We identify changes in proteins involved in translation, folding, degradation, oxido-reduction and metabolic processes.

Conclusion

Our functional proteomic study provides the first inventory list of over 300 proteins that directly or indirectly affect Sup35p aggregation and [PSI⁺] formation. Our results highlight the complexity of the cellular changes accompanying [PSI⁺] formation and pave the way for in vitro studies aimed to document the effect of individual and/or combinations of proteins identified here, susceptible of affecting Sup35p assembly.     

Molecular dynamics studies on the structural stability of wild-type dog prion protein

Prion diseases such as Creutzfeldt-Jakob disease, variant Creutzfeldt-Jakob diseases, Gerstmann-Str?ussler-Scheinker syndrome, Fatal Familial Insomnia, Kuru in humans, scrapie in sheep, bovine spongiform encephalopathy (or 'mad-cow' disease) and chronic wasting disease in cattle are invariably fatal and highly infectious neurodegenerative diseases affecting humans and animals. However, by now there have not been some effective therapeutic approaches to treat all these prion diseases. In 2008, canine mammals including dogs (canis familials) were the first time academically reported to be resistant to prion diseases (Vaccine 26: 2601-2614 (2008)). Thus, it is very worth studying the molecular structures of dog prion protein to obtain insights into the immunity of dogs to prion diseases. This paper studies the molecular structural dynamics of wild-type dog prion protein. The comparison analyses with rabbit prion protein show that the dog prion protein has stable molecular structures whether under neutral or low pH environments. We also find that the salt bridges such as D177-R163 contribute to the structural stability of wild-type rabbit prion protein under neutral pH environment.     

Structure, dynamics, and stability of assemblies of the human prion fragment SNQNNF

Misfolding of the prion protein (PrP) is associated with the development of Transmissible Spongiform Encephalopathies. The recent crystal structure of ‘steric zipper’ aggregates of the peptide SNQNNF (human PrP fragment 170-175) has highlighted its potential involvement in the misfolding process. A detailed molecular dynamics investigation on SNQNNF aggregates has been performed to analyze the behavior of the assemblies in a non-crystalline context. Stability, dynamics, and structural features suggest that SNQNNF assemblies are very good candidates to be involved in the structure of PrP fibrils. In addition, the analysis of small aggregates shows that steric zipper interfaces are able to stabilize assemblies composed of four strands per sheet. Altogether, the present findings indicate that steric zipper may play a key role in prion diseases. This suggestion is also corroborated by MD analyses of point mutations within the region 170-175.