Evidence that the 127-164 region of prion proteins has two equi-energetic conformations with beta or alpha features

Dynamics of Amyloid Formation from Simplified Representation to Atomistic Simulations

Amyloid fibril formation is an intrinsic property of short peptides, non-disease proteins, and proteins associated with neurodegenerative diseases. Aggregates of the Aβ and tau proteins, the α-synuclein protein, and the prion protein are observed in the brain of Alzheimer's, Parkinson's, and prion disease patients, respectively. Due to the transient short-range and long-range interactions of all species and their high aggregation propensities, the conformational ensemble of these devastating proteins, the exception being for the monomeric prion protein, remains elusive by standard structural biology methods in bulk solution and in lipid membranes. To overcome these limitations, an increasing number of simulations using different sampling methods and protein models have been performed. In this chapter, we first review our main contributions to the field of amyloid protein simulations aimed at understanding the early aggregation steps of short linear amyloid peptides, the conformational ensemble of the Aβ40/42 dimers in bulk solution, and the stability of Aβ aggregates in lipid membrane models. Then we focus on our studies on the interactions of amyloid peptides/inhibitors to prevent aggregation, and long amyloid sequences, including new results on a monomeric tau construct.

Detecting early stage structural changes in wild type, pathogenic and non-pathogenic prion variants using Markov state model

The conversion of prion protein from normal to scrapie followed by the aggregation and deposition of this scrapie form leads to various neurodegenerative diseases. A few studies carried out by researchers suggest that E219K prion mutant (glutamate to lysine mutation at residue position 219) is more stable than wild type protein. However a similar point mutation E200K (glutamate to lysine mutation at residue position 200) is pathogenic. In this study we have carried out detailed atomistic simulation of the wild type, pathogenic mutant E200K and E219K mutant which provides more stability. The aim of the study was to detect the early structural changes present in all the three variants which might be responsible for the stability or for their conversion from PrPC to PrPSc. MSM based analyses have been carried out to find out the differences between WT, E200K and E219K systems. Markov state model (MSM) analysis was able to predict the intermediate states which helped to understand the effect of same mutation at two different locations. The MSM analysis was able to show that the extra stability of E219K mutant may be a result of the increase in number of native contacts, strong salt bridges and less random motions. While pathogenicity of E200K mutant can be attributed to loss of some crucial salt-bridge interactions, increased random motions between helix 2 and helix 3.

Markov state model to find out the differences between WT, E200K and E219K systems.

The OPEP protein model: from single molecules, amyloid formation, crowding and hydrodynamics to DNA/RNA systems

The OPEP coarse-grained protein model has been applied to a wide range of applications since its first release 15 years ago. The model, which combines energetic and structural accuracy and chemical specificity, allows the study of single protein properties, DNA-RNA complexes, amyloid fibril formation and protein suspensions in a crowded environment. Here we first review the current state of the model and the most exciting applications using advanced conformational sampling methods. We then present the current limitations and a perspective on the ongoing developments.

Structural diversity and initial oligomerization of PrP106-126 studied by replica-exchange and conventional molecular dynamics simulations

Prion diseases are marked by cerebral accumulation of the abnormal isoform of the prion protein. A fragment of prion protein composed of residues 106–126 (PrP106–126) exhibits similar properties to full length prion and plays a key role in the conformational conversion from cellular prion to its pathogenic pattern. Soluble oligomers of PrP106–126 have been proposed to be responsible for neurotoxicity. However, the monomeric conformational space and initial oligomerization of PrP106–126 are still obscure, which are very important for understanding the conformational conversion of PrP106–126. In this study, replica exchange molecular dynamics simulations were performed to investigate monomeric and dimeric states of PrP106–126 in implicit solvent. The structural diversity of PrP106–126 was observed and this peptide did not acquire stable structure. The dimeric PrP106–126 also displayed structural diversity and hydrophobic interaction drove the dimerization. To further study initial oligomerization of PrP106–126, 1 µs conventional molecular dynamics simulations of trimer and tetramer formation were carried out in implicit solvent. We have observed the spontaneous formation of several basic oligomers and stable oligomers with high β-sheet contents were sampled in the simulations of trimer and tetramer formation. The β-hairpin formed in hydrophobic tail of PrP106–126 with residues 118–120 in turn may stabilize these oligomers and seed the formation oligomers. This study can provide insight into the detailed information about the structure of PrP106–126 and the dynamics of aggregation of monomeric PrP106–126 into oligomers in atomic level.

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.

Case for an RNA-prion world: a hypothesis based on conformational diversity

Prions and other misfolded proteins can impart their structure and functions to normal molecules. Based upon a thorough structural assessment of RNA, prions and misfolded proteins, especially from the perspective of conformational diversity, we propose a case for co-existence of these in the pre-biotic world. Analyzing the evolution of physical aspects of biochemical structures, we put forward a case for an RNA-prion pre-biotic world, instead of, merely, the "RNA World".

Structures and thermodynamics of Alzheimer's amyloid-beta Abeta(16-35) monomer and dimer by replica exchange molecular dynamics simulations: implication for full-length Abeta fibrillation

Many proteins display a strand-loop-strand motif in their amyloid fibrillar states. For instance, the amyloid beta-protein, Abeta1-40, associated with Alzheimer's disease, displays a loop at positions 22-28 in its amyloid fibril state. It has been suggested that this loop could appear early in the aggregation process, but quantitative information regarding its presence in small oligomers remains scant. Because residues 1-15 are disordered in Abeta1-42 fibrils and Abeta10-35 forms fibrils in vitro, we select the peptide Abeta16-35, centered on residues 22-28 and determine the structures and thermodynamics of the monomer and dimer using coarse-grained implicit solvent replica exchange molecular dynamics simulations. Our simulations totalling 5 mus for the monomer and 12 micros for the dimer show no sign of strong secondary structure signals in both instances and the significant impact of dimerization on the global structure of Abeta16-35. They reveal however that the loop 22-28 acts as a quasi-independent unit in both species. The loop structure ensemble we report in Abeta16-35 monomer and dimer has high similarity to the loop formed by the Abeta21-30 peptide in solution and, to a lesser extent, to the loop found in Abeta1-40 fibrils. We discuss the implications of our findings on the assembly of full-length Abeta.

Modeling of protein misfolding in disease

A short review of the results of molecular modeling of prion disease is presented in this chapter. According to the "one-protein theory" proposed by Prusiner, prion proteins are misfolded naturally occurring proteins, which, on interaction with correctly folded proteins may induce misfolding and propagate the disease, resulting in insoluble amyloid aggregates in cells of affected specimens. Because of experimental difficulties in measurements of origin and growth of insoluble amyloid aggregations in cells, theoretical modeling is often the only one source of information regarding the molecular mechanism of the disease. Replica exchange Monte Carlo simulations presented in this chapter indicate that proteins in the native state, N, on interaction with an energetically higher structure, R, can change their conformation into R and form a dimer, R(2). The addition of another protein in the N state to R(2) may lead to spontaneous formation of a trimer, R(3). These results reveal the molecular basis for a model of prion disease propagation or conformational diseases in general.

Comparative analysis of the prion protein gene sequences in African lion

The prion protein gene of African lion (Panthera Leo) was first cloned and polymorphisms screened. The results suggest that the prion protein gene of eight African lions is highly homogenous. The amino acid sequences of the prion protein (PrP) of all samples tested were identical. Four single nucleotide polymorphisms (C42T, C81A, C420T, T600C) in the prion protein gene (Prnp) of African lion were found, but no amino acid substitutions. Sequence analysis showed that the higher homology is observed to felis catus AF003087 (96.7%) and to sheep number M31313.1 (96.2%) Genbank accessed. With respect to all the mammalian prion protein sequences compared, the African lion prion protein sequence has three amino acid substitutions. The homology might in turn affect the potential intermolecular interactions critical for cross species transmission of prion disease.

Structural instability of the prion protein upon M205S/R mutations revealed by molecular dynamics simulations

The point mutations M205S and M205R have been demonstrated to severely disturb the folding and maturation process of the cellular prion protein (PrP(C)). These disturbances have been interpreted as consequences of mutation-induced structural changes in PrP, which are suggested to involve helix 1 and its attachment to helix 3, because the mutated residue M205 of helix 3 is located at the interface of these two helices. Furthermore, current models of the prion protein scrapie (PrP(Sc)), which is the pathogenic isoform of PrP(C) in prion diseases, imply that helix 1 disappears during refolding of PrP(C) into PrP(Sc). Based on molecular-dynamics simulations of wild-type and mutant PrP(C) in aqueous solution, we show here that the native PrP(C) structure becomes strongly distorted within a few nanoseconds, once the point mutations M205S and M205R have been applied. In the case of M205R, this distortion is characterized by a motion of helix 1 away from the hydrophobic core into the aqueous environment and a subsequent structural decay. Together with experimental evidence on model peptides, this decay suggests that the hydrophobic attachment of helix 1 to helix 3 at M205 is required for its correct folding into its stable native structure.

The role of the 132-160 region in prion protein conformational transitions

The native conformation of host-encoded cellular prion protein (PrP(C)) is metastable. As a result of a post-translational event, PrP(C) can convert to the scrapie form (PrP(Sc)), which emerges as the essential constituent of infectious prions. Despite thorough research, the mechanism underlying this conformational transition remains unknown. However, several studies have highlighted the importance of the N-terminal region spanning residues 90-154 in PrP folding. In order to understand why PrP folds into two different conformational states exhibiting distinct secondary and tertiary structure, and to gain insight into the involvement of this particular region in PrP transconformation, we studied the pressure-induced unfolding/ refolding of recombinant Syrian hamster PrP expanding from residues 90-231, and compared it with heat unfolding. By using two intrinsic fluorescent variants of this protein (Y150W and F141W), conformational changes confined to the 132-160 segment were monitored. Multiple conformational states of the Trp variants, characterized by their spectroscopic properties (fluorescence and UV absorbance in the fourth derivative mode), were achieved by tuning the experimental conditions of pressure and temperature. Further insight into unexplored conformational states of the prion protein, likely to mimic the in vivo structural change, was obtained from pressure-assisted cold unfolding. Furthermore, salt-induced conformational changes suggested a structural stabilizing role of Tyr150 and Phe141 residues, slowing down the conversion to a beta-sheet form.

A cylinder-shaped double ribbon structure formed by an amyloid hairpin peptide derived from the β-sheet of murine PrP: An X-ray and molecular dynamics simulation study

A structural model of the murine PrP small beta-sheet was obtained by synthesizing the RGYMLGSADPNGNQVYYRG peptide comprising the two beta-strands 127-133 and 159-164 linked by a four-residue sequence of high turn propensity. The DPNG turn sequence is a "short circuit" replacing the original protein sequence between the two strands. This 19-residue peptide spontaneously forms very long single fibrils as observed by electron microscopy. The X-ray diffraction patterns of a partially oriented sample reveals an average arrangement of the hairpin peptides into a structure which can be geometrically approximated by an empty-core cylinder. The hairpins are oriented perpendicular to the cylinder axis and a 130 A helix period is observed. Based on X-ray diffraction constraints and on more indirect general protein structure considerations, a precise and consistent fibril model was built. The structure consists of two beta-sheet ribbons wound around a cylinder and assembled into a single fibril with a hairpin orientation perpendicular to the fibril axis. Subsequent implicit and explicit solvent molecular dynamics simulations provided the final structure at atomic resolution and further insights into the stabilizing interactions. Particularly important are the zipper-like network of polar interactions between the edges of the two ribbons, including the partially buried water molecules. The hydrophobic core is not optimally compact explaining the low density of this region seen by X-ray diffraction. The present findings provide also a simple model for further investigating the sequence-stability relationship using a mutational approach with a quasi-independent consideration of the polar and apolar interactions.

Amyloid aggregates of the prion peptide PrP106-126 are destabilised by oxidation and by the action of dendrimers

The prion protein (PrP) peptide 106-126 forms amyloid aggregates in vitro and this sequence is speculated to be involved in the formation of amyloid fibrils by the abnormally folded PrP protein (PrPSc) found in spongiform encephalopathies. It is shown here by incubation experiments in water using Thioflavin T (ThT) as a fluorescent probe for amyloid formation that changes in C-terminal charge, oxidation state and conformational stabilisation lead to large changes in amyloid forming behaviour (amyloidogenicity) of this peptide. Amyloid formation is favoured by a charged C-terminus and is strongly inhibited by oxidation. Furthermore, cationic dendrimers are shown to perturb peptide fibrillation in a process dependent on the nature of the charged groups on the dendrimer surface.

Coupled prediction of protein secondary and tertiary structure

The strong coupling between secondary and tertiary structure formation in protein folding is neglected in most structure prediction methods. In this work we investigate the extent to which nonlocal interactions in predicted tertiary structures can be used to improve secondary structure prediction. The architecture of a neural network for secondary structure prediction that utilizes multiple sequence alignments was extended to accept low-resolution nonlocal tertiary structure information as an additional input. By using this modified network, together with tertiary structure information from native structures, the Q3-prediction accuracy is increased by 7-10% on average and by up to 35% in individual cases for independent test data. By using tertiary structure information from models generated with the ROSETTA de novo tertiary structure prediction method, the Q3-prediction accuracy is improved by 4-5% on average for small and medium-sized single-domain proteins. Analysis of proteins with particularly large improvements in secondary structure prediction using tertiary structure information provides insight into the feedback from tertiary to secondary structure.

Minimal model for studying prion-like folding pathways

The Monte Carlo technique is used to simulate the energy landscape and the folding kinetics of a minimal prion-like protein model. We show that the competition between hydrogen-bonding and hydrophobic interactions yields two energetically favored secondary structures, an alpha-helix and a beta-hairpin. Folding simulations indicate that the probability of reaching the alpha-helix form from a denatured random conformation is much higher than the probability of reaching the beta-sheet form, even though the beta-sheet has a lower energy. The existence of a lower energy beta-sheet state gives the possibility for the normal alpha-helix structure to take a structural transformation into the beta-sheet structure under external influences.

Molecular evolution of the mammalian prion protein

Prion protein (PrP) sequences are until now available for only six of the 18 orders of placental mammals. A broader comparison of mammalian prions might help to understand the enigmatic functional and pathogenic properties of this protein. We therefore determined PrP coding sequences in 26 mammalian species to include all placental orders and major subordinal groups. Glycosylation sites, cysteines forming a disulfide bridge, and a hydrophobic transmembrane region are perfectly conserved. Also, the sequences responsible for secondary structure elements, for N- and C-terminal processing of the precursor protein, and for attachment of the glycosyl-phosphatidylinositol membrane anchor are well conserved. The N-terminal region of PrP generally contains five or six repeats of the sequence P(Q/H)GGG(G/-)WGQ, but alleles with two, four, and seven repeats were observed in some species. This suggests, together with the pattern of amino acid replacements in these repeats, the regular occurrence of repeat expansion and contraction. Histidines implicated in copper ion binding and a proline involved in 4-hydroxylation are lacking in some species, which questions their importance for normal functioning of cellular PrP. The finding in certain species of two or seven repeats, and of amino acid substitutions that have been related to human prion diseases, challenges the relevance of such mutations for prion pathology. The gene tree deduced from the PrP sequences largely agrees with the species tree, indicating that no major deviations occurred in the evolution of the prion gene in different placental lineages. In one species, the anteater, a prion pseudogene was present in addition to the active gene.