Stability of single sheet GNNQQNY aggregates analyzed by replica exchange molecular dynamics: antiparallel versus parallel association

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.

Assembly modes of hexaphenylalanine variants as function of the charge states of their terminal ends

The ability of peptides to self-assemble represents a valuable tool for the development of biomaterials of biotechnological and/or biomedical interest. Diphenylalanine homodimer (FF) and its analogues are among the most promising systems in this field. The longest Phe-based building block hitherto characterized is pentaphenylalanine (F5). We studied the aggregation propensity and the structural/morphological features of assemblies of zwitterionic hexaphenylalanine H+-F6-O- and of three variants characterized by different charged states of the terminal ends (Ac-F6-Amide, H+-F6-Amide and Ac-F6-O-). As previously observed for PEGylated hexaphenylalanine (PEG8-F6), all F6 variants show a strong tendency to form β-rich assemblies in which the structural motif is constituted by antiparallel β-strands in the cross-β framework. Extensive replica exchange molecular dynamics simulations carried out on a pairs of F6 peptides indicate that the antiparallel β-structure of the final assemblies is likely dictated by the preferred association modes of the individual chains in the very early stages of the aggregation process. Our data suggest that even very small F6 peptides are properly pre-organized and prone to the build-up of the final assembly.

The Associative Memory, Water Mediated, Structure and Energy Model (AWSEM)-Amylometer: Predicting Amyloid Propensity and Fibril Topology Using an Optimized Folding Landscape Model

Amyloids are fibrillar protein aggregates with simple repeated structural motifs in their cores, usually β-strands but sometimes α-helices. Identifying the amyloid-prone regions within protein sequences is important both for understanding the mechanisms of amyloid-associated diseases and for understanding functional amyloids. Based on the crystal structures of seven cross-β amyloidogenic peptides with different topologies and one recently solved cross-α fiber structure, we have developed a computational approach for identifying amyloidogenic segments in protein sequences using the Associative memory, Water mediated, Structure and Energy Model (AWSEM). The AWSEM-Amylometer performs favorably in comparison with other predictors in predicting aggregation-prone sequences in multiple data sets. The method also predicts well the specific topologies (the relative arrangement of β-strands in the core) of the amyloid fibrils. An important advantage of the AWSEM-Amylometer over other existing methods is its direct connection with an efficient, optimized protein folding simulation model, AWSEM. This connection allows one to combine efficient and accurate search of protein sequences for amyloidogenic segments with the detailed study of the thermodynamic and kinetic roles that these segments play in folding and aggregation in the context of the entire protein sequence. We present new simulation results that highlight the free energy landscapes of peptides that can take on multiple fibril topologies. We also demonstrate how the Amylometer methodology can be straightforwardly extended to the study of functional amyloids that have the recently discovered cross-α fibril architecture.

Intrinsic structural versatility of the highly conserved 412-423 epitope of the Hepatitis C Virus E2 protein

HCV infection is a major threaten for human health as it affects hundreds of million people worldwide. Here we investigated the conformational properties of the 412-423 fragment of the envelope E2 protein, one of the most immunogenic regions of the virus proteome whose characterization may provide interesting insights for anti-HCV vaccine development. The spectroscopic characterization of the polypeptide unravels its unexpected tendency to form amyloid-like aggregates. When kept in monomeric state, it shows a limited tendency to adopt regular secondary structure. Enhanced molecular dynamics simulations, starting from four distinct conformational states, highlight its structural versatility. Interestingly, all multiform conformational states of the polypeptide detected in crystallographic complexes with antibodies are present in the structural ensemble of all simulations. This observation corroborates the idea that known antibodies recognize this region through a conformational selection mechanism. Accordingly, the design of effective anti-HCV vaccines should consider the intrinsic flexibility of this region. The structural versatility of the 412-423 region is particularly puzzling if its remarkable sequence conservation is considered. It is likely that flexibility and sequence conservation are important features that endow this epitope with the ability to accomplish distinct functions such as immunity escape and interaction with host receptors.

Inhibition of GNNQQNY prion peptide aggregation by trehalose: a mechanistic view

Deposition of amyloid fibrils is the seminal event in the pathogenesis of numerous neurodegenerative diseases. The formation of this amyloid assembly is the manifestation of a cascade of structural transitions including toxic oligomer formation in the early stages of aggregation. Thus a viable therapeutic strategy involves the use of small molecular ligands to interfere with this assembly. In this perspective, we have explored the kinetics of aggregate formation of the fibril forming GNNQQNY peptide fragment from the yeast prion protein SUP35 using multiple all atom MD simulations with explicit solvent and provided mechanistic insights into the way trehalose, an experimentally known aggregation inhibitor, modulates the aggregation pathway. The results suggest that the assimilation process is impeded by different barriers at smaller and larger oligomeric sizes: the initial one being easily surpassed at higher temperatures and peptide concentrations. The kinetic profile demonstrates that trehalose delays the aggregation process by increasing both these activation barriers, specifically the latter one. It increases the sampling of small-sized aggregates that lack the beta sheet conformation. Analysis reveals that the barrier in the growth of larger stable oligomers causes the formation of multiple stable small oligomers which then fuse together bimolecularly. The PCA of 26 properties was carried out to deconvolute the events within the temporary lag phases, which suggested dynamism in lags involving an increase in interchain contacts and burial of SASA. The predominant growth route is monomer addition, which changes to condensation on account of a large number of depolymerisation events in the presence of trehalose. The favourable interaction of trehalose specifically with the sidechain of the peptide promotes crowding of trehalose molecules in its vicinity - the combination of both these factors imparts the observed behaviour. Furthermore, increasing trehalose concentration leads to faster expulsion of water molecules than interpeptide interactions. These expelled water molecules have larger translational movement, suggesting an entropy factor to favor the assembly process. Different conformations observed under this condition suggest the role of water molecules in guiding the morphology of the aggregates as well. A similar scenario exists on increasing peptide concentration.

Inhibitory effect of hydrophobic fullerenes on the β-sheet-rich oligomers of a hydrophilic GNNQQNY peptide revealed by atomistic simulations

Fullerenes suppress fibril-like β-sheet oligomers by interacting strongly with the nonpolar aliphatic groups of polar residues of GNNQQNY peptide, thus inhibit peptide aggregation.

Amyloid Properties of Asparagine and Glutamine in Prion-like Proteins

Sequences rich in glutamine (Q) and asparagine (N) are intrinsically disordered in monomeric form, but can aggregate into highly ordered amyloids, as seen in Q/N-rich prion domains (PrDs). Amyloids are fibrillar protein aggregates rich in β-sheet structures that can self-propagate through protein-conformational chain reactions. Here, we present a comprehensive theoretical study of N/Q-rich peptides, including sequences found in the yeast Sup35 PrD, in parallel and antiparallel β-sheet aggregates, and probe via fully atomistic molecular dynamics simulations all their possible steric-zipper interfaces in order to determine their protofibril structure and their relative stability. Our results show that polyglutamine aggregates are more stable than polyasparagine aggregates. Enthalpic contributions to the free energy favor the formation of polyQ protofibrils, while entropic contributions favor the formation of polyN protofibrils. The considerably larger phase space that disordered polyQ must sample on its way to aggregation probably is at the root of the associated slower kinetics observed experimentally. When other amino acids are present, such as in the Sup35 PrD, their shorter side chains favor steric-zipper formation for N but not Q, as they preclude the in-register association of the long Q side chains.

Polymorphism of amyloid fibrils formed by a peptide from the yeast prion protein Sup35: AFM and Tip-Enhanced Raman Scattering studies

Aggregation of prion proteins is the cause of various prion related diseases. The infectious form of prions, amyloid aggregates, exist as multiple strains. The strains are thought to represent structurally different prion protein molecules packed into amyloid aggregates, but the knowledge on the structure of different types of aggregates is limited. Here we report on the use of AFM (Atomic Force Microscopy) and TERS (Tip-Enhanced Raman Scattering) to study morphological heterogeneity and access underlying conformational features of individual amyloid aggregates. Using AFM we identified the morphology of amyloid fibrils formed by the peptide (CGNNQQNY) from the yeast prion protein Sup35 that is critically involved in the aggregation of the full protein. TERS results demonstrate that morphologically different amyloid fibrils are composed of a distinct set of conformations. Fibrils formed at pH 5.6 are composed of a mixture of peptide conformations (β-sheets, random coil and α-helix) while fibrils formed in pH~2 solution primarily have β-sheets. Additionally, peak positions in the amide III region of the TERS spectra suggested that peptides have parallel arrangement of β-sheets for pH~2 fibrils and antiparallel arrangement for fibrils formed at pH 5.6. We also developed a methodology for detailed analysis of the peptide secondary structure by correlating intensity changes of Raman bands in different regions of TERS spectra. Such correlation established that structural composition of peptides is highly localized with large contribution of unordered secondary structures on a fibrillar surface.

Role of monomer arrangement in the amyloid self-assembly

Assembly of amyloid proteins into aggregates requires the ordering of the monomers in oligomers and especially in such highly organized structures as fibrils. This ordering is accompanied by structural transitions leading to the formation of ordered β-structural motifs in proteins and peptides lacking secondary structures. To characterize the effect of the monomer arrangements on the aggregation process at various stages, we performed comparative studies of the yeast prion protein Sup35 heptapeptide (GNNQQNY) along with its dimeric form CGNNQQNY-(d-Pro)-G-GNNQQNY. The (d-Pro)-G linker in this construct is capable of adopting a β-turn, facilitating the assembly of the dimer into the dimeric antiparallel hairpin structure (AP-hairpin). We applied Atomic Force Microscopy (AFM) techniques to follow peptide-peptide interactions at the single molecule level, to visualize the morphology of aggregates formed by both constructs, thioflavin T (ThT) fluorescence to follow the aggregation kinetics, and circular dichroism (CD) spectroscopy to characterize the secondary structure of the constructs. The ThT fluorescence data showed that the AP-hairpin aggregation kinetics is insensitive to the external environment such as ionic strength and pH contrary to the monomers the kinetics of which depends dramatically on the ionic strength and pH. The AFM topographic imaging revealed that AP-hairpins primarily assemble into globular aggregates, whereas linear fibrils are primary assemblies of the monomers suggesting that both constructs follow different aggregation pathways during the self-assembly. These morphological differences are in line with the AFM force spectroscopy experiments and CD spectroscopy measurements, suggesting that the AP-hairpin is structurally rigid regardless of changes of environmental factors.

Interplay of sequence, topology and termini charge in determining the stability of the aggregates of GNNQQNY mutants: a molecular dynamics study

This study explores the stabilities of single sheet parallel systems of three sequence variants of ¹GNNQQNY⁷, N2D, N2S and N6D, with variations in aggregate size (5–8) and termini charge (charged or neutral). The aggregates were simulated at 300 and 330 K. These mutations decrease amyloid formation in the yeast prion protein Sup35. The present study finds that these mutations cause instability even in the peptide context. The protonation status of termini is found to be a key determinant of stabilities; other determinants are sequence, position of mutation and aggregate size. All systems with charged termini are unstable, whereas both stable and unstable systems are found when the termini are neutral. When termini are charged, the largest stable aggregate for the N2S and N6D systems has 3 to 4 peptides whereas N2D mutation supports oligomers of larger size (5-and 6-mers) as well. Mutation at 2^nd position (N2S and N2D) results in fewer H-bonds at the mutated as well as neighboring (Gly1/Gln4) positions. However, no such effect is found if mutation is at 6^th position (N6D). The effect of Asn→Asp mutation depends on the position and termini charge: it is more destabilizing at the 2^nd position than at the 6^th in case of neutral termini, however, the opposite is true in case of charged termini. Appearance of twist in stable systems and in smaller aggregates formed in unstable systems suggests that twist is integral to amyloid arrangement. Disorder, dissociation or rearrangement of peptides, disintegration or collapse of aggregates and formation of amorphous aggregates observed in these simulations are likely to occur during the early stages of aggregation also. The smaller aggregates formed due to such events have a variety of arrangements of peptides. This suggests polymorphic nature of oligomers and presence of a heterogeneous mixture of oligomers during early stages of aggregation.

A kinetic approach to the sequence-aggregation relationship in disease-related protein assembly

It is generally accepted that oligomers of aggregating proteins play an important role in the onset of neurodegenerative diseases. While in silico aggregation studies of full length amyloidogenic proteins are computationally expensive, the assembly of short protein fragments derived from these proteins with similar aggregating properties has been extensively studied. In the present work, molecular dynamics simulations are performed to follow peptide aggregation on the microsecond time scale. By defining aggregation states, we identify transition networks, disconnectivity graphs, and first passage time distributions to describe the kinetics of the assembly process. This approach unravels differences in the aggregation into hexamers of two peptides with different primary structures. The first is GNNQQNY, a hydrophilic fragment from the prion protein Sup35, and the second is KLVFFAE, a fragment from amyloid-β protein, with a hydrophobic core delimited by two charged amino acids. The assembly of GNNQQNY suggests a mechanism of monomer addition, with a bias toward parallel peptide pairs and a gradual increase in the amount of β-strand content. For KLVFFAE, a mechanism involving dimers rather than monomers is revealed, involving a generally higher β-strand content and a transition toward a larger number of antiparallel peptide pairs during the rearrangement of the hexamer. The differences observed for the aggregation of the two peptides suggests the existence of a sequence-aggregation relationship.

Factors that drive peptide assembly and fibril formation: experimental and theoretical analysis of Sup35 NNQQNY mutants

Residue mutations have substantial effects on aggregation kinetics and propensities of amyloid peptides and their aggregate morphologies. Such effects are attributed to conformational transitions accessed by various types of oligomers such as steric zipper or single β-sheet. We have studied the aggregation propensities of six NNQQNY mutants: NVVVVY, NNVVNV, NNVVNY, VIQVVY, NVVQIY, and NVQVVY in water using a combination of ion-mobility mass spectrometry, transmission electron microscopy, atomic force microscopy, and all-atom molecular dynamics simulations. Our data show a strong correlation between the tendency to form early β-sheet oligomers and the subsequent aggregation propensity. Our molecular dynamics simulations indicate that the stability of a steric zipper structure can enhance the propensity for fibril formation. Such stability can be attained by either hydrophobic interactions in the mutant peptide or polar side-chain interdigitations in the wild-type peptide. The overall results display only modest agreement with the aggregation propensity prediction methods such as PASTA, Zyggregator, and RosettaProfile, suggesting the need for better parametrization and model peptides for these algorithms.

Kinetics of amyloid aggregation: a study of the GNNQQNY prion sequence

The small amyloid-forming GNNQQNY fragment of the prion sequence has been the subject of extensive experimental and numerical studies over the last few years. Using unbiased molecular dynamics with the OPEP coarse-grained potential, we focus here on the onset of aggregation in a 20-mer system. With a total of 16.9 of simulations at 280 K and 300 K, we show that the GNNQQNY aggregation follows the classical nucleation theory (CNT) in that the number of monomers in the aggregate is a very reliable descriptor of aggregation. We find that the critical nucleus size in this finite-size system is between 4 and 5 monomers at 280 K and 5 and 6 at 300 K, in overall agreement with experiment. The kinetics of growth cannot be fully accounted for by the CNT, however. For example, we observe considerable rearrangements after the nucleus is formed, as the system attempts to optimize its organization. We also clearly identify two large families of structures that are selected at the onset of aggregation demonstrating the presence of well-defined polymorphism, a signature of amyloid growth, already in the 20-mer aggregate.

Protein aggregation plays an important pathological role in numerous neurodegenerative diseases such as Alzheimer's, Parkinson's, Creutzfeldt-Jakob, the Prion disease and diabetes mellitus. In most cases, misfolded proteins are involved and aggregate irreversibly to form highly ordered insoluble macrostructures, called amyloid fibrils, which deposit in the brain. Studies have revealed that all proteins are capable of forming amyloid fibrils that all share common structural features and therefore aggregation mechanisms. The toxicity of amyloid aggregates is however not attributed to the fibrils themselves but rather to smaller more disordered aggregates, oligomers, forming parallel to or prior to fibrils. Understanding the assembly process of these amyloid oligomers is key to understanding their toxicity mechanism in order to devise a possible treatment strategy targeting these toxic aggregates. Our approach here is to computationally study the aggregation dynamics of a 20-mer of an amyloid peptide GNNQQNY from a prion protein. Our findings suggest that the assembly is a spontaneous process that can be described as a complex nucleation and growth mechanism and which can lead to two classes of morphologies for the aggregates, one of which resembles a protofibril-like structure. Such numerical studies are crucial to understanding the details of fast biological processes and complement well experimental studies.

Intrinsic disorder modulates protein self-assembly and aggregation

Protein molecules have evolved to adopt distinctive and well-defined functional and soluble states under physiological conditions. In some circumstances, however, proteins can self-assemble into fibrillar aggregates designated as amyloid fibrils. In vivo these processes are normally associated with severe pathological conditions but can sometimes have functional relevance. One such example is the hydrophobins, whose aggregation at air-water interfaces serves to create robust protein coats that help fungal spores to resist wetting and thus facilitate their dispersal in the air. We have performed multiscale simulations to address the molecular determinants governing the formation of functional amyloids by the class I fungal hydrophobin EAS. Extensive samplings of full-atom replica-exchange molecular dynamics and coarse-grained simulations have allowed us to identify factors that distinguish aggregation-prone from highly soluble states of EAS. As a result of unfavourable entropic terms, highly dynamical regions are shown to exert a crucial influence on the propensity of the protein to aggregate under different conditions. More generally, our findings suggest a key role that specific flexible structural elements can play to ensure the existence of soluble and functional states of proteins under physiological conditions.

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 complexity of a composite amyloid fibril

The molecular structure of amyloid fibrils and the mechanism of their formation are of substantial medical and biological importance, but present an ongoing experimental and computational challenge. An early high-resolution view of amyloid-like structure was obtained on amyloid-like crystals of a small fragment of the yeast prion protein Sup35p: the peptide GNNQQNY. As GNNQQNY also forms amyloid-like fibrils under similar conditions, it has been theorized that the crystal's structural features are shared by the fibrils. Here we apply magic-angle-spinning (MAS) NMR to examine the structure and dynamics of these fibrils. Previously multiple NMR signals were observed for such samples, seemingly consistent with the presence of polymorphic fibrils. Here we demonstrate that peptides with these three distinct conformations instead assemble together into composite protofilaments. Electron microscopy (EM) of the ribbon-like fibrils indicates that these protofilaments combine in differing ways to form striations of variable widths, presenting another level of structural complexity. Structural and dynamic NMR data reveal the presence of highly restricted side-chain conformations involved in interfaces between differently structured peptides, likely comprising interdigitated steric zippers. We outline molecular interfaces that are consistent with the observed EM and NMR data. The rigid and uniform structure of the GNNQQNY crystals is found to contrast distinctly with the more complex structural and dynamic nature of these "composite" amyloid fibrils. These results provide insight into the fibril-crystal distinction and also indicate a necessary caution with respect to the extrapolation of crystal structures to the study of fibril structure and formation.

A multiscale approach to characterize the early aggregation steps of the amyloid-forming peptide GNNQQNY from the yeast prion sup-35

The self-organization of peptides into amyloidogenic oligomers is one of the key events for a wide range of molecular and degenerative diseases. Atomic-resolution characterization of the mechanisms responsible for the aggregation process and the resulting structures is thus a necessary step to improve our understanding of the determinants of these pathologies. To address this issue, we combine the accelerated sampling properties of replica exchange molecular dynamics simulations based on the OPEP coarse-grained potential with the atomic resolution description of interactions provided by all-atom MD simulations, and investigate the oligomerization process of the GNNQQNY for three system sizes: 3-mers, 12-mers and 20-mers. Results for our integrated simulations show a rich variety of structural arrangements for aggregates of all sizes. Elongated fibril-like structures can form transiently in the 20-mer case, but they are not stable and easily interconvert in more globular and disordered forms. Our extensive characterization of the intermediate structures and their physico-chemical determinants points to a high degree of polymorphism for the GNNQQNY sequence that can be reflected at the macroscopic scale. Detailed mechanisms and structures that underlie amyloid aggregation are also provided.

The formation of amyloid fibrils is associated with many neurodegenerative diseases such as Alzheimer's, Creutzfeld-Jakob, Parkinson's, the Prion disease and diabetes mellitus. In all cases, proteins misfold to form highly ordered insoluble aggregates called amyloid fibrils that deposit intra- and extracellularly and are resistant to proteases. All proteins are believed to have the instrinsic capability of forming amyloid fibrils that share common specific structural properties that have been observed by X-ray crystallography and by NMR. However, little is known about the aggregation dynamics of amyloid assemblies, and their toxicity mechanism is therefore poorly understood. It is believed that small amyloid oligomers, formed on the aggregation pathway of full amyloid fibrils, are the toxic species. A detailed atomic characterization of the oligomerization process is thus necessary to further our understanding of the amyloid oligomer's toxicity. Our approach here is to study the aggregation dynamics of a 7-residue amyloid peptide GNNQQNY through a combination of numerical techniques. Our results suggest that this amyloid sequence can form fibril-like structures and is polymorphic, which agrees with recent experimental observations. The ability to fully characterize and describe the aggregation pathway of amyloid sequences numerically is key to the development of future drugs to target amyloid oligomers.

Structural characterization of GNNQQNY amyloid fibrils by magic angle spinning NMR

Several human diseases are associated with the formation of amyloid aggregates, but experimental characterization of these amyloid fibrils and their oligomeric precursors has remained challenging. Experimental and computational analysis of simpler model systems has therefore been necessary, for instance, on the peptide fragment GNNQQNY7−13 of yeast prion protein Sup35p. Expanding on a previous publication, we report here a detailed structural characterization of GNNQQNY fibrils using magic angle spinning (MAS) NMR. On the basis of additional chemical shift assignments we confirm the coexistence of three distinct peptide conformations within the fibrillar samples, as reflected in substantial chemical shift differences. Backbone torsion angle measurements indicate that the basic structure of these coexisting conformers is an extended β-sheet. We structurally characterize a previously identified localized distortion of the β-strand backbone specific to one of the conformers. Intermolecular contacts are consistent with each of the conformers being present in its own parallel and in-register sheet. Overall the MAS NMR data indicate a substantial difference between the structure of the fibrillar and crystalline forms of these peptides, with a clearly increased complexity in the GNNQQNY fibril structure. These experimental data can provide guidance for future work, both experimental and theoretical, and provide insights into the distinction between fibril growth and crystal formation.

GNNQQNY--investigation of early steps during amyloid formation

Protein aggregation has been implicated in the pathology of several neurodegenerative diseases, and a better understanding of how it proceeds is essential for the development of therapeutic strategies. Recently, the amyloidogenic heptapeptide GNNQQNY has emerged as a molecule of choice for fundamental studies of protein aggregation. A number of experimental and computational studies have examined the structure of the GNNQQNY aggregate. Less work, however, has been aimed at understanding its aggregation pathway. In this study, we present a detailed computational analysis of such a pathway. To that end, transition path sampling Monte Carlo simulations are used to examine the dimerization process. A statistical analysis of the reaction pathways shows that the dimerization reaction proceeds via a zipping mechanism, initiated with the formation of distinct contacts at the third residue (N). Asparagine residues are found to play a key role in the early stages of aggregation. And, contrary to previous belief, it is also shown that the tyrosine terminal group is not required to stabilize the dimer. In fact, an asparagine residue leads to faster aggregation of the peptide.

HRMAS 1H NMR conformational study of the resin-bound amyloid-forming peptide GNNQQNY from the yeast prion Sup35

The conversion of soluble proteins to insoluble amyloid fibrils is associated with numerous human diseases. The peptide GNNQQNY is a short segment of the yeast prion protein Sup35 that previously has been found to form amyloid fibrils in a similar manner to the protein itself. The approach taken in this work was to attach this peptide sequence to an insoluble polymer matrix through solid phase peptide synthesis and give it the internal freedom to fold into its local conformation in an organic solvent. Observation of its monomeric structure, free from the effects of aggregation, entropic solvent effects, and neighboring molecules, was possible by two-dimensional high-resolution magic angle spinning (1)H NMR spectroscopy. Analysis of the through-bond correlations and through-space interactions observed in the spectra, combined with global energy minimization via computational studies, led to the observation that the peptide chain adopts a compact beta-like turn at the central hydrophilic residues. The technique of peptide attachment to a polymer resin and observation by NMR may allow for future study of single peptide fragments prone to aggregation.

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

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