A novel canonical dual computational approach for prion AGAAAAGA amyloid fibril molecular modeling

Many experimental studies have shown that the prion AGAAAAGA palindrome hydrophobic region (113-120) has amyloid fibril forming properties and plays an important role in prion diseases. However, due to the unstable, noncrystalline and insoluble nature of the amyloid fibril, to date structural information on AGAAAAGA region (113-120) has been very limited. This region falls just within the N-terminal unstructured region PrP (1-123) of prion proteins. Traditional X-ray crystallography and nuclear magnetic resonance (NMR) spectroscopy experimental methods cannot be used to get its structural information. Under this background, this paper introduces a novel approach of the canonical dual theory to address the 3D atomic-resolution structure of prion AGAAAAGA amyloid fibrils. The novel and powerful canonical dual computational approach introduced in this paper is for the molecular modeling of prion AGAAAAGA amyloid fibrils, and that the optimal atomic-resolution structures of prion AGAAAAGA amyloid fibils presented in this paper are useful for the drive to find treatments for prion diseases in the field of medicinal chemistry. Overall, this paper presents an important method and provides useful information for treatments of prion diseases.     

Optimal atomic-resolution structures of prion AGAAAAGA amyloid fibrils

X-ray crystallography is a powerful tool to determine the protein 3D structure. However, it is time-consuming and expensive, and not all proteins can be successfully crystallized, particularly for membrane proteins. Although nuclear magnetic resonance (NMR) spectroscopy is indeed a very powerful tool in determining the 3D structures of membrane proteins, it is also time-consuming and costly. To the best of the authors' knowledge, there is little structural data available on the AGAAAAGA palindrome in the hydrophobic region (113-120) of prion proteins due to the noncrystalline and insoluble nature of the amyloid fibril, although many experimental studies have shown that this region has amyloid fibril forming properties and plays an important role in prion diseases. In view of this, the present study is devoted to address this problem from computational approaches such as global energy optimization, simulated annealing, and structural bioinformatics. The optimal atomic-resolution structures of prion AGAAAAGA amyloid fibils reported in this paper have a value to the scientific community in its drive to find treatments for prion diseases.     

Computer simulation study of amyloid fibril formation by palindromic sequences in prion peptides

We simulate the aggregation of large systems containing palindromic peptides from the Syrian hamster prion protein SHaPrP 113-120 (AGAAAAGA) and the mouse prion protein MoPrP 111-120 (VAGAAAAGAV) and eight sequence variations: GAAAAAAG, (AG)(4) , A8, GAAAGAAA, A10, V10, GAVAAAAVAG, and VAVAAAAVAV The first two peptides are thought to act as the Velcro that holds the parent prion proteins together in amyloid structures and can form fibrils themselves. Kinetic events along the fibrillization pathway influence the types of structures that occur and variations in the sequence affect aggregation kinetics and fibrillar structure. Discontinuous molecular dynamics simulations using the PRIME20 force field are performed on systems containing 48 peptides starting from a random coil configuration. Depending on the sequence, fibrillar structures form spontaneously over a range of temperatures, below which amorphous aggregates form and above which no aggregation occurs. AGAAAAGA forms well organized fibrillar structures whereas VAGAAAAGAV forms less well organized structures that are partially fibrillar and partially amorphous. The degree of order in the fibrillar structure stems in part from the types of kinetic events leading up to its formation, with AGAAAAGA forming less amorphous structures early in the simulation than VAGAAAAGAV. The ability to form fibrils increases as the chain length and the length of the stretch of hydrophobic residues increase. However as the hydrophobicity of the sequence increases, the ability to form well-ordered structures decreases. Thus, longer hydrophobic sequences form slightly disordered aggregates that are partially fibrillar and partially amorphous. Subtle changes in sequence result in slightly different fibril structures.     

Contribution of Specific Residues of the β-Solenoid Fold to HET-s Prion Function,Amyloid Structure and Stability

The [Het-s] prion of the fungus Podospora anserina represents a good model system for studying the structure-function relationship in amyloid proteins because a high resolution solid-state NMR structure of the amyloid prion form of the HET-s prion forming domain (PFD) is available. The HET-s PFD adopts a specific β-solenoid fold with two rungs of β-strands delimiting a triangular hydrophobic core. A C-terminal loop folds back onto the rigid core region and forms a more dynamic semi-hydrophobic pocket extending the hydrophobic core. Herein, an alanine scanning mutagenesis of the HET-s PFD was conducted. Different structural elements identified in the prion fold such as the triangular hydrophobic core, the salt bridges, the asparagines ladders and the C-terminal loop were altered and the effect of these mutations on prion function, fibril structure and stability was assayed. Prion activity and structure were found to be very robust; only a few key mutations were able to corrupt structure and function. While some mutations strongly destabilize the fold, many substitutions in fact increase stability of the fold. This increase in structural stability did not influence prion formation propensity in vivo. However, if an Ala replacement did alter the structure of the core or did influence the shape of the denaturation curve, the corresponding variant showed a decreased prion efficacy. It is also the finding that in addition to the structural elements of the rigid core region, the aromatic residues in the C-terminal semi-hydrophobic pocket are critical for prion propagation. Mutations in the latter region either positively or negatively affected prion formation. We thus identify a region that modulates prion formation although it is not part of the rigid cross-β core, an observation that might be relevant to other amyloid models.     

Steric zipper of the amyloid fibrils formed by residues 109-122 of the Syrian hamster prion protein

We report the results of atomic force microscopy, Fourier-transform infrared spectroscopy, solid-state nuclear magnetic resonance, and molecular dynamics (MD) calculations for amyloid fibrils formed by residues 109-122 of the Syrian hamster prion protein (H1). Our data reveal that H1 fibrils contain no more than two β-sheet layers. The peptide strands of H1 fibrils are antiparallel with the A117 residues aligned to form a linear chain in the direction of the fibril axis. The molecular structure of the H1 fibrils, which adopts the motif of steric zipper, is highly uniform in the region of the palindrome sequence AGAAAAGA. The closest distance between the two adjacent β-sheet layers is found to be about 5 Å. The structural features of the molecular model of H1 fibrils obtained by MD simulations are consistent with the experimental results. Overall, our solid-state NMR and MD simulation data indicate that a steric zipper, which was first observed in the crystals of fibril-forming peptides, can be formed in H1 fibrils near the region of the palindrome sequence.     

Molecular dynamics simulations of alanine rich beta-sheet oligomers: Insight into amyloid formation

The aggregation observed in protein conformational diseases is the outcome of significant new beta-sheet structure not present in the native state. Peptide model systems have been useful in studies of fibril aggregate formation. Experimentally, it was found that a short peptide AGAAAAGA is one of the most highly amyloidogenic peptides. This peptide corresponds to the Syrian hamster prion protein (ShPrP) residues 113-120. The peptide was observed to be conserved in all species for which the PrP sequence has been determined. We have simulated the stabilities of oligomeric AGAAAAGA and AAAAAAAA (A8) by molecular dynamic simulations. Oligomers of both AGAAAAGA and AAAAAAAA were found to be stable when the size is 6 to 8 (hexamer to octamer). Subsequent simulation of an additional alpha-helical AAAAAAAA placed on the A8-octamer surface has revealed molecular events related to conformational change and oligomer growth. Our study addresses both the minimal oligomeric size of an aggregate seed and the mechanism of seed growth. Our simulations of the prion-derived 8-residue amyloidogenic peptide and its variant have indicated that an octamer is stable enough to be a seed and that the driving force for stabilization is the hydrophobic effect.     

The Hydrophobic Core Sequence Modulates the Neurotoxic and Secondary Structure Properties of the Prion Peptide 106-126

The neurodegeneration seen in spongiform encephalopathies is believed to be mediated by protease-resistant forms of the prion protein (PrP). A peptide encompassing residues 106-126 of human PrP has been shown to be neurotoxic in vitro. The neurotoxicity of PrP106-126 appears to be dependent upon its adoption of an aggregated fibril structure. To examine the role of the hydrophobic core, AGAAAAGA, on PrP106-126 toxicity, we performed structure-activity analyses by substituting two or more hydrophobic residues for the hydrophilic serine residue to decrease its hydrophobicity. A peptide with a deleted alanine was also synthesized. We found all the peptides except the deletion mutant were no longer toxic on mouse cerebellar neuronal cultures. Circular dichroism analysis showed that the nontoxic PrP peptides had a marked decrease in beta-sheet structure. In addition, the mutants had alterations in aggregability as measured by turbidity, Congo red binding, and fibril staining using electron microscopy. These data show that the hydrophobic core sequence is important for PrP106-126 toxicity probably by influencing its assembly into a neurotoxic structure. The hydrophobic sequence may similarly affect aggregation and toxicity observed in prion diseases.     

Structural polymorphism in amyloids: new insights from studies with Y145Stop prion protein fibrils

The C-terminally-truncated human prion protein variant Y145Stop (or PrP23-144), associated with a familial prion disease, provides a valuable model for studying the fundamental properties of protein amyloids. In previous solid-state NMR experiments, we established that the β-sheet core of the PrP23-144 amyloid is composed of two β-strand regions encompassing residues ～113-125 and ～130-140. The former segment contains a highly conserved hydrophobic palindrome sequence, (113)AGAAAAGA(120), which has been considered essential to PrP conformational conversion. Here, we examine the role of this segment in fibrillization of PrP23-144 using a deletion variant, Δ113-120 PrP23-144, in which the palindrome sequence is missing. Surprisingly, we find that deletion of the palindrome sequence affects neither the amyloidogenicity nor the polymerization kinetics of PrP23-144, although it does alter amyloid conformation and morphology. Using two-dimensional and three-dimensional solid-state NMR methods, we find that Δ113-120 PrP23-144 fibrils contain an altered β-core extended N-terminally to residue ～106, encompassing residues not present in the core of wild-type PrP23-144 fibrils. The C-terminal β-strand of the core, however, is similar in both fibril types. Collectively, these data indicate that amyloid cores of PrP23-144 variants contain "essential" (i.e. nucleation-determining) and "nonessential" regions, with the latter being "movable" in amino acid sequence space. These findings reveal an intriguing new mechanism for structural polymorphism in amyloids and suggest a potential means for modulating the physicochemical properties of amyloid fibrils without compromising their polymerization characteristics.     

Development of the Structural Core and of Conformational Heterogeneity during the Conversion of Oligomers of the Mouse Prion Protein to Worm-like Amyloid Fibrils

Understanding how structure develops during the course of amyloid fibril formation by the prion protein is important for understanding prion diseases. Determining how conformational heterogeneity manifests itself in the fibrillar and pre-fibrillar amyloid aggregates is critical for understanding prion strain phenotypes. In this study, the formation of worm-like amyloid fibrils by the mouse prion protein has been characterized structurally by hydrogen-deuterium exchange coupled to mass spectrometry. The structural cores of these fibrils and of the oligomer on the direct pathway of amyloid fibril formation have been defined, showing how structure develops during fibril formation. The structural core of the oligomer not on the direct pathway has also been defined, allowing the delineation of the structural features that make this off-pathway oligomer incompetent to directly form fibrils. Sequence segments that exhibit multiple local conformations in the three amyloid aggregates have been identified, and the development of structural heterogeneity during fibril formation has been characterized. It is shown that conformational heterogeneity is not restricted to only the C-terminal domain region, which forms the structural core of the aggregates; it manifests itself in the N-terminal domain of the protein as well. Importantly, all three amyloid aggregates are shown to be capable of disrupting lipid membrane structure, pointing to a mechanism by which they may be toxic.     

Steric zipper formed by hydrophobic peptide fragment of Syrian hamster prion protein

Steric zippers, where the residues of two neighboring β-sheet layers are tightly interdigitated, have been proposed as fundamental structural units of amyloid fibrils by Eisenberg and co-workers. The steric zipper formed by polypeptides containing the palindromic sequence AGAAAAGA has a distinctive feature that the distance between two interdigitated β-sheet layers is comparable to the interstrand distance of the individual β-sheet. This structural motif is of great interest in the study of prion disease because the AGAAAAGA sequence is highly conserved in prion proteins of different species. In this work, the amyloid fibrils formed by the polypeptides of PrP(113-127), viz. Ac-AGAAAAGAVVGGLGG-NH(2), are taken as the model compound to investigate the biophysical principles governing the steric zipper formation. The target fibrils adopt the structural motif of class 7 steric zipper, which is formed by stacking of antiparallel β-sheet layers with residue 117 + k forming backbone hydrogen bonds to residue 120 - k. Implication of our results in the infectivity of scrapie prion is briefly discussed.     

Formation of amyloid fibrils from β-amylase

Fibril formation has been considered a significant feature of amyloid proteins. However, it has been proposed that fibril formation is a common property of many proteins under appropriate conditions. We studied the fibril formation of β-amylase, a non-amyloid protein rich in α-helical structure, because the secondary structure of β-amylase is similar to that of prions. With the conditions for the fibril formation of prions, β-amylase proteins were converted into amyloid fibrils. The features of β-amylase proteins and fibrils are compared to prion proteins and fibrils. Furthermore, the cause of neurotoxicity in amyloid diseases is discussed.     

Solid-state NMR as a method to reveal structure and membrane-interaction of amyloidogenic proteins and peptides

It is important to understand the Amyloid fibril formation in view of numerous medical and biochemical aspects. Structural determination of amyloid fibril has been extensively studied using electron microscopy. Subsequently, solid state NMR spectroscopy has been realized to be the most important means to determine not only microscopic molecular structure but also macroscopic molecular packing. Molecular structure of amyloid fibril was first predicted to be parallel β-sheet structure, and subsequently, was further refined for Aβ(1-40) to be cross β-sheet with double layered in register parallel β-sheet structure by using solid state NMR spectroscopy. On the other hand, anti-parallel β-sheet structure has been reported to short fragments of Aβ-amyloid and other amyloid forming peptides. Kinetic study of amyloid fibril formation has been studied using a variety of methods, and two-step autocatalytic reaction mechanism used to explain fibril formation. Recently, stable intermediates or proto-fibrils have been observed by electron microscope (EM) images. Some of the intermediates have the same microscopic structure as the matured fibril and subsequently change to matured fibrils. Another important study on amyloid fibril formation is determination of the interaction with lipid membranes, since amyloid peptide are cleaved from amyloid precursor proteins in the membrane interface, and it is reported that amyloid lipid interaction is related to the cytotoxicity. Finally it is discussed how amyloid fibril formation can be inhibited. Firstly, properly designed compounds are reported to have inhibition ability of amyloid fibril formation by interacting with amyloid peptide. Secondly, it is revealed that site directed mutation can inhibit amyloid fibril formation. These inhibitors were developed by knowing the fibril structure determined by solid state NMR.     

Solid-state NMR as a method to reveal structure and membrane-interaction of amyloidogenic proteins and peptides

It is important to understand the Amyloid fibril formation in view of numerous medical and biochemical aspects. Structural determination of amyloid fibril has been extensively studied using electron microscopy. Subsequently, solid state NMR spectroscopy has been realized to be the most important means to determine not only microscopic molecular structure but also macroscopic molecular packing. Molecular structure of amyloid fibril was first predicted to be parallel beta-sheet structure, and subsequently, was further refined for Abeta(1-40) to be cross beta-sheet with double layered in register parallel beta-sheet structure by using solid state NMR spectroscopy. On the other hand, anti-parallel beta-sheet structure has been reported to short fragments of Abeta-amyloid and other amyloid forming peptides. Kinetic study of amyloid fibril formation has been studied using a variety of methods, and two-step autocatalytic reaction mechanism used to explain fibril formation. Recently, stable intermediates or proto-fibrils have been observed by electron microscope (EM) images. Some of the intermediates have the same microscopic structure as the matured fibril and subsequently change to matured fibrils. Another important study on amyloid fibril formation is determination of the interaction with lipid membranes, since amyloid peptide are cleaved from amyloid precursor proteins in the membrane interface, and it is reported that amyloid lipid interaction is related to the cytotoxicity. Finally it is discussed how amyloid fibril formation can be inhibited. Firstly, properly designed compounds are reported to have inhibition ability of amyloid fibril formation by interacting with amyloid peptide. Secondly, it is revealed that site directed mutation can inhibit amyloid fibril formation. These inhibitors were developed by knowing the fibril structure determined by solid state NMR.     

Residual Structures in the Acid-Unfolded States of Vλ6 Proteins Affect Amyloid Fibrillation

Many proteins form amyloid-like fibrils in vitro under partially or highly unfolding conditions. Recently, we showed that the residual structure in highly unfolded state is closely related to amyloid fibril formation in hen lysozyme. Thus, to better understand the role of the residual structure on amyloid fibril formation, we focused on AL amyloidosis, which results from the extracellular deposition of monoclonal immunoglobulin light-chain variable domains (V_Ls) as insoluble fibrils. We examined the relationship between the residual structure and amyloid fibril formation on three λ6 recombinant V_L (rVλ6) proteins, wild type, Jto, and Wil. Although rVλ6 proteins are highly unfolded in pH 2, ¹⁵N NMR transverse relaxation experiments revealed nonrandom structures in regions, which include some hydrophobic residues and a single disulfide bond, indicating the existence of residual structure in rVλ6 proteins. However, the residual structure of Wil was markedly disrupted compared with those of the other proteins, despite there being no significant differences in amino acid sequences. Fibrillation experiments revealed that Wil had a longer lag time for fibril formation than the others. When the single disulfide bond was reduced and alkylated, the residual structure was largely disrupted and fibril formation was delayed in all three rVλ6 proteins. It was suggested that the residual structure in highly unfolded state has a crucial role in amyloid fibril formation in many proteins, even pathogenic ones.     

A Simple Lattice Model That Captures Protein Folding,Aggregation and Amyloid Formation

The ability of many proteins to convert from their functional soluble state to amyloid fibrils can be attributed to inter-molecular beta strand formation. Such amyloid formation is associated with neurodegenerative disorders like Alzheimer''s and Parkinson''s. Molecular modelling can play a key role in providing insight into the factors that make proteins prone to fibril formation. However, fully atomistic models are computationally too expensive to capture the length and time scales associated with fibril formation. As the ability to form fibrils is the rule rather than the exception, much insight can be gained from the study of coarse-grained models that capture the key generic features associated with amyloid formation. Here we present a simple lattice model that can capture both protein folding and beta strand formation. Unlike standard lattice models, this model explicitly incorporates the formation of hydrogen bonds and the directionality of side chains. The simplicity of our model makes it computationally feasible to investigate the interplay between folding, amorphous aggregation and fibril formation, and maintains the capability of classic lattice models to simulate protein folding with high specificity. In our model, the folded proteins contain structures that resemble naturally occurring beta-sheets, with alternating polar and hydrophobic amino acids. Moreover, fibrils with intermolecular cross-beta strand conformations can be formed spontaneously out of multiple short hydrophobic peptide sequences. Both the formation of hydrogen bonds in folded structures and in fibrils is strongly dependent on the amino acid sequence, indicating that hydrogen-bonding interactions alone are not strong enough to initiate the formation of beta sheets. This result agrees with experimental observations that beta sheet and amyloid formation is strongly sequence dependent, with hydrophobic sequences being more prone to form such structures. Our model should open the way to a systematic study of the interplay between the factors that lead to amyloid formation.     

All or none fibrillogenesis of a prion peptide.

Amyloid proteins and peptides comprise a diverse group of molecules that vary both in size and amino-acid sequence, yet assemble into amyloid fibrils that have a common core structure. Kinetic studies of amyloid fibrillogenesis have revealed that certain amyloid proteins form oligomeric intermediates prior to fibril formation. We have investigated fibril formation with a peptide corresponding to residues 195-213 of the human prion protein. Through a combination of kinetic and equilibrium studies, we have found that the fibrillogenesis of this peptide proceeds as an all-or-none reaction where oligomeric intermediates are not stably populated. This variation in whether oligomeric intermediates are stably populated during fibril formation indicates that amyloid proteins assemble into a common fibrillar structure; however, they do so through different pathways.     

Amyloid Properties of the Mouse Egg Zona Pellucida

The zona pellucida (ZP) surrounding the oocyte is an extracellular fibrillar matrix that plays critical roles during fertilization including species-specific gamete recognition and protection from polyspermy. The mouse ZP is composed of three proteins, ZP1, ZP2, and ZP3, all of which have a ZP polymerization domain that directs protein fibril formation and assembly into the three-dimensional ZP matrix. Egg coats surrounding oocytes in nonmammalian vertebrates and in invertebrates are also fibrillar matrices and are composed of ZP domain-containing proteins suggesting the basic structure and function of the ZP/egg coat is highly conserved. However, sequence similarity between ZP domains is low across species and thus the mechanism for the conservation of ZP/egg coat structure and its function is not known. Using approaches classically used to identify amyloid including conformation-dependent antibodies and dyes, X-ray diffraction, and negative stain electron microscopy, our studies suggest the mouse ZP is a functional amyloid. Amyloids are cross-β sheet fibrillar structures that, while typically associated with neurodegenerative and prion diseases in mammals, can also carry out functional roles in normal cells without resulting pathology. An analysis of the ZP domain from mouse ZP3 and ZP3 homologs from five additional taxa using the algorithm AmylPred 2 to identify amyloidogenic sites, revealed in all taxa a remarkable conservation of regions that were predicted to form amyloid. This included a conserved amyloidogenic region that localized to a stretch of hydrophobic amino acids previously shown in mouse ZP3 to be essential for fibril assembly. Similarly, a domain in the yeast protein α-agglutinin/Sag 1p, that possesses ZP domain-like features and which is essential for mating, also had sites that were predicted to be amyloidogenic including a hydrophobic stretch that appeared analogous to the critical site in mouse ZP3. Together, these studies suggest that amyloidogenesis may be a conserved mechanism for ZP structure and function across billions of years of evolution.     

Structures for amyloid fibrils

Alzheimer's disease and Creutzfeldt-Jakob disease are the best-known examples of a group of diseases known as the amyloidoses. They are characterized by the extracellular deposition of toxic, insoluble amyloid fibrils. Knowledge of the structure of these fibrils is essential for understanding the process of pathology of the amyloidoses and for the rational design of drugs to inhibit or reverse amyloid formation. Structural models have been built using information from a wide variety of techniques, including X-ray diffraction, electron microscopy, solid state NMR and EPR. Recent advances have been made in understanding the architecture of the amyloid fibril. Here, we describe and compare postulated structural models for the mature amyloid fibril and discuss how the ordered structure of amyloid contributes to its stability.     

Hydration effects on the HET-s prion and amyloid-beta fibrillous aggregates, studied with three-dimensional molecular theory of solvation

We study the thermodynamic properties of the experimental fragments of the amyloid fibril made of the HET-s prion proteins (the infectious element of the filamentous fungus Podospora anserina) and of amyloid-β proteins (the major component of Alzheimer's disease-associated plaques) by using the three-dimensional molecular theory of solvation. The full quantitative picture of hydration effects, including the hydration thermodynamics and hydration structure around the fragments, is presented. For both the complexes, the hydration entropic effects dominate, which results in the entropic part offsetting the unfavorable energetic part of the free energy change upon the association. This is in accord with the fact that the hydrophobic cooperativity plays an essential role in the formation of amyloid fibrils. By calculating the partial molar volume of the proteins, we found that the volume change upon the association in both the systems is large and positive, with the implication that high pressure causes destabilization of the fibril. This observation is in good agreement with the recent experimental results. We also found that both the HET-s and amyloid-β pentamers have loose intermolecular packing with voids. The three-dimensional molecular theory of solvation predicts that water molecules can be locked in the interior cavities along the fibril axis for both the HET-s and amyloid-β proteins. We provide a detailed molecular picture of the structural water localized in the interior of the fibrils. Our results suggest that the interior hydration plays an important role in the structural stability of fibrils.