Molecular dynamics simulations indicate a possible role of parallel beta-helices in seeded aggregation of poly-Gln

Protein Self-Assembly at the Liquid-Surface Interface. Surface-Mediated Aggregation Catalysis

Protein self-assembly into aggregates of various morphologies is a ubiquitous phenomenon in physical chemistry and biophysics. The critical role of amyloid assemblies in the development of diseases, neurodegenerative diseases especially, highlights the importance of understanding the mechanistic picture of the self-assembly process. The translation of this knowledge to the development of efficient preventions and treatments for diseases requires designing experiments at conditions mimicking those in vivo. This Perspective reviews data satisfying two major requirements: membrane environment and physiologically low concentrations of proteins. Recent progress in experiments and computational modeling resulted in a novel model for the amyloid aggregation process at the membrane-liquid interface. The self-assembly under such conditions has a number of critical features, further understanding of which can lead to the development of efficient preventive means and treatments for Alzheimer's and other devastating neurodegenerative disorders.

Advances in the understanding of protein misfolding and aggregation through molecular dynamics simulation

Aberrant protein folding known as protein misfolding is counted as one of the striking factors of neurodegenerative diseases. The extensive range of pathologies caused by protein misfolding, aggregation and subsequent accumulation are mainly classified into either gain of function diseases or loss of function diseases. In order to seek for novel strategies for treatment and diagnosis of neurodegenerative diseases, insights into the mechanism of misfolding and aggregation is essential. A comprehensive knowledge on the factors influencing misfolding and aggregation is required as well. An extensive experimental study on protein aggregation is somewhat challenging due to the insoluble and noncrystalline nature of amyloid fibrils. Thus there has been a growing use of computational approaches including Monte Carlo simulation, docking simulation, molecular dynamics simulation in the study of protein misfolding and aggregation. The review presents a discussion on molecular dynamics simulation alone as to how it has emerged as a promising tool in the understanding of protein misfolding and aggregation in general, detailing upon three different aspects considering four misfold prone proteins in particular. It is noticeable that all four proteins considered in this review i.e prion, superoxide dismutase1, huntingtin and amyloid β are linked to chronic neurodegenerative diseases with debilitating effects. Initially the review elaborates on the factors influencing the misfolding and aggregation. Next, it addresses our current understanding of the amyloid structures and the associated aggregation mechanisms, finally, summarizing the contribution of this computational tool in the search for therapeutic strategies against the respective protein-deposition diseases.

Molecular Dynamics Simulations Applied to Structural and Dynamical Transitions of the Huntingtin Protein: A Review

Over the recent years, Huntington's disease (HD) has become widely discussed in the scientific literature especially because at the mutant level there are several contradictions regarding the aggregation mechanism. The specific role of the physiological huntingtin protein remains unknown, due to the lack of characterization of its entire crystallographic structure, making the experimental and theoretical research even harder when taking into consideration its involvement in multiple biological functions and its high affinity for different interacting partners. Different types of models, containing fewer (not more than 35 Qs) polyglutamine residues for the WT structure and above 35 Qs for the mutants, were subjected to classical or advanced MD simulations to establish the proteins' structural stability by evaluating their conformational changes. Outside the polyQ tract, there are two other regions of interest (the N17 domain and the polyP rich domain) considered to be essential for the aggregation kinetics at the mutant level. The polymerization process is considered to be dependent on the polyQ length. As the polyQ tract's dimension increases, the structures present more β-sheet conformations. Contrarily, it is also considered that the aggregation stability is not necessarily dependent on the number of Qs, while the initial stage of the aggregation seed might play the decisive role. A general assumption regarding the polyP domain is that it might preserve the polyQ structures soluble by acting as an antagonist for β-sheet formation.

Structural Biology of PrP Prions

Prion diseases are characterized by the deposition of amyloids, misfolded conformers of the prion protein. The misfolded conformation is self-replicating, by a mechanism solely enciphered in the conformation of the protein. Because of low solubility and heterogeneous aggregate sizes, the detailed atomic structure of the infectious isoform is still unknown. Progress has, however, been made, and has allowed insights into the structural and disease-related mechanisms of prions. Many structural models have been proposed, and a number of them support a consensus trimeric β-helical model, significantly more complex than simple amyloid models. There is evidence that such complexity may be a necessary property of prion structure. Knowledge of the structure of prions will provide a greater understanding of the protein isoform conversion mechanism, and could eventually lead to rationally designed intervention strategies.

Molecular dynamics analysis of the aggregation propensity of polyglutamine segments

Protein misfolding and aggregation is a pathogenic feature shared among at least ten polyglutamine (polyQ) neurodegenerative diseases. While solvent-solution interaction is a key factor driving protein folding and aggregation, the solvation properties of expanded polyQ tracts are not well understood. By using GPU-enabled all-atom molecular dynamics simulations of polyQ monomers in an explicit solvent environment, this study shows that solvent-polyQ interaction propensity decreases as the lengths of polyQ tract increases. This study finds a predominance in long-distance interactions between residues far apart in polyQ sequences with longer polyQ segments, that leads to significant conformational differences. This study also indicates that large loops, comprised of parallel β-structures, appear in long polyQ tracts and present new aggregation building blocks with aggregation driven by long-distance intra-polyQ interactions. Finally, consistent with previous observations using coarse-grain simulations, this study demonstrates that there is a gain in the aggregation propensity with increased polyQ length, and that this gain is correlated with decreasing ability of solvent-polyQ interaction. These results suggest the modulation of solvent-polyQ interactions as a possible therapeutic strategy for treating polyQ diseases.

The CAA repeat polymorphism in the ZFHX3 gene is associated with risk of coronary heart disease in a Chinese population

Coronary heart disease (CHD) is a disease resulting from the interaction between genetic variations and environmental factors. Zinc finger homeobox 3 (ZFHX3) is a transcription factor and contains a poly-glutamine tract in a compositionally biased region that is encoded by exon 9, containing a cluster of CAG and CAA triplets followed by the polymorphic CAA repeats: (CAG)2(CAA)2(CAG)3CAACAG(CAA)nGCA. Thus, nine successive glutamine residues precede the poly-glutamine tract, encoded by the polymorphic CAA repeats. The aim of this study was to investigate the association of the CAA repeat polymorphism in exon 9 of the ZFHX3 gene with the risk of CHD in a Chinese population. The CAA repeat polymorphism was determined by polymerase chain reaction followed by DNA sequencing in 321 CHD patients. Genotype frequencies were compared using the non-parametric mood median test. Four alleles of CAG(CAA)10GCA, CAG(CAA)8GCA, CAG(CAA)9GCA, and CAG(CAA)11GCA were found in Chinese CHD patients in exon 9 of the ZFHX3 gene. The CAG(CAA)10GCA was a major allele (95.95%), and the CAG(CAA)8GCA was a minor allele (3.58%). The CAG(CAA)9GCA and CAG(CAA)11GCA were rare alleles (0.31% and 0.16%). The CAG(CAA)10GCA allele encodes a poly-glutamine tract of 19 residues. Importantly, the CHD patients homozygous for the CAG(CAA)10GCA allele had a higher risk of CHD, compared to the heterozygous patients carrying a CAG(CAA)8GCA allele. Moreover, the CAG(CAA)10GCA allele was significantly associated with hypertension, diabetes mellitus, or dyslipidemia (P < 0.05). Thus, the CAA repeat polymorphism in exon 9 of the ZFHX3 gene contributes to the CHD susceptibility in the Chinese population.

The bactofilin cytoskeleton protein BacM of Myxococcus xanthus forms an extended β-sheet structure likely mediated by hydrophobic interactions

Bactofilins are novel cytoskeleton proteins that are widespread in Gram-negative bacteria. Myxococcus xanthus, an important predatory soil bacterium, possesses four bactofilins of which one, BacM (Mxan_7475) plays an important role in cell shape maintenance. Electron and fluorescence light microscopy, as well as studies using over-expressed, purified BacM, indicate that this protein polymerizes in vivo and in vitro into ~3 nm wide filaments that further associate into higher ordered fibers of about 10 nm. Here we use a multipronged approach combining secondary structure determination, molecular modeling, biochemistry, and genetics to identify and characterize critical molecular elements that enable BacM to polymerize. Our results indicate that the bactofilin-determining domain DUF583 folds into an extended β-sheet structure, and we hypothesize a left-handed β-helix with polymerization into 3 nm filaments primarily via patches of hydrophobic amino acid residues. These patches form the interface allowing head-to-tail polymerization during filament formation. Biochemical analyses of these processes show that folding and polymerization occur across a wide variety of conditions and even in the presence of chaotropic agents such as one molar urea. Together, these data suggest that bactofilins are comprised of a structure unique to cytoskeleton proteins, which enables robust polymerization.

The structure of the infectious prion protein: experimental data and molecular models

The structures of the infectious prion protein, PrP(Sc), and that of its proteolytically truncated variant, PrP 27-30, have evaded experimental determination due to their insolubility and propensity to aggregate. Molecular modeling has been used to fill this void and to predict their structures, but various modeling approaches have produced significantly different models. The disagreement between the different modeling solutions indicates the limitations of this method. Over the years, in absence of a three-dimensional (3D) structure, a variety of experimental techniques have been used to gain insights into the structure of this biologically, medically, and agriculturally important isoform. Here, we present an overview of experimental results that were published in recent years, and which provided new insights into the molecular architecture of PrP(Sc) and PrP 27-30. Furthermore, we evaluate all published models in light of these recent, experimental data, and come to the conclusion that none of the models can accommodate all of the experimental constraints. Moreover, this conclusion constitutes an open invitation for renewed efforts to model the structure of PrP(Sc).

Engineering amyloid fibrils from β-solenoid proteins for biomaterials applications

Nature provides numerous examples of self-assembly that can potentially be implemented for materials applications. Considerable attention has been given to one-dimensional cross-β or amyloid structures that can serve as templates for wire growth or strengthen materials such as glue or cement. Here, we demonstrate controlled amyloid self-assembly based on modifications of β-solenoid proteins. They occur naturally in several contexts (e.g., antifreeze proteins, drug resistance proteins) but do not aggregate in vivo due to capping structures or distortions at their ends. Removal of these capping structures and regularization of the ends of the spruce budworm and rye grass antifreeze proteins yield micron length amyloid fibrils with predictable heights, which can be a platform for biomaterial-based self-assembly. The design process, including all-atom molecular dynamics simulations, purification, and self-assembly procedures are described. Fibril formation with the predicted characteristics is supported by evidence from thioflavin-T fluorescence, circular dichroism, dynamic light scattering, and atomic force microscopy. Additionally, we find evidence for lateral assembly of the modified spruce budworm antifreeze fibrils with sufficient incubation time. The kinetics of polymerization are consistent with those for other amyloid formation reactions and are relatively fast due to the preformed nature of the polymerization nucleus.

The diversity and relationship of prion protein self-replicating states

It has become evident that the prion protein (PrP) can form a diverse range of self-replicating structures in addition to bona fide PrP(Sc) or strain-specific PrP(Sc) variants. Some self-replicating states can be only produced in vitro, whereas others can be formed in vivo and in vitro. While transmissible, not all states that replicate in vivo are truly pathogenic. Some of them can replicate silently without causing symptoms or clinical diseases. In the current article we discuss the data on PK-digestion patterns of different self-replicating PrP states in connection with other structural data available to date and assess possible relationships between different self-replicating states. Even though different self-replicating PrP states appear to have significantly different global folding patterns, it seems that the C-terminal region exhibits a cross-β-sheet structure in all self-replicating states, as this region acquires the proteolytically most stable conformation. We also discuss the possibility of the transformation of self-replicating states and triggering of PrP(Sc) formation within the frame of the deformed templating model. The spread of silent self-replicating states is of a particular concern because they can lead to transmissible prion disease. Moreover, examples on how different replication requirements favor different states are discussed. This knowledge can help in designing conditions for selective amplification of a particular PrP state in vitro.

Transformation between α-helix and β-sheet structures of one and two polyglutamine peptides in explicit water molecules by replica-exchange molecular dynamics simulations

Aggregation of polyglutamine peptides with β-sheet structures is related to some important neurodegenerative diseases such as Huntington's disease. However, it is not clear how polyglutamine peptides form the β-sheets and aggregate. To understand this problem, we performed all-atom replica-exchange molecular dynamics simulations of one and two polyglutamine peptides with 10 glutamine residues in explicit water molecules. Our results show that two polyglutamine peptides mainly formed helix or coil structures when they are separated, as in the system with one-polyglutamine peptide. As the interpeptide distance decreases, the intrapeptide β-sheet structure sometimes appear as an intermediate state, and finally the interpeptide β-sheets are formed. We also find that the polyglutamine dimer tends to form the antiparallel β-sheet conformations rather than the parallel β-sheet, which is consistent with previous experiments and a coarse-grained molecular dynamics simulation.

Computational approaches to understanding protein aggregation in neurodegeneration

The generation of toxic non-native protein conformers has emerged as a unifying thread among disorders such as Alzheimer's disease, Parkinson's disease, and amyotrophic lateral sclerosis. Atomic-level detail regarding dynamical changes that facilitate protein aggregation, as well as the structural features of large-scale ordered aggregates and soluble non-native oligomers, would contribute significantly to current understanding of these complex phenomena and offer potential strategies for inhibiting formation of cytotoxic species. However, experimental limitations often preclude the acquisition of high-resolution structural and mechanistic information for aggregating systems. Computational methods, particularly those combine both all-atom and coarse-grained simulations to cover a wide range of time and length scales, have thus emerged as crucial tools for investigating protein aggregation. Here we review the current state of computational methodology for the study of protein self-assembly, with a focus on the application of these methods toward understanding of protein aggregates in human neurodegenerative disorders.

Atomistic mechanisms of huntingtin N-terminal fragment insertion on a phospholipid bilayer revealed by molecular dynamics simulations

The huntingtin protein is characterized by a segment of consecutive glutamines (Q(N)) that is responsible for its fibrillation. As with other amyloid proteins, misfolding of huntingtin is related to Huntington's disease through pathways that can involve interactions with phospholipid membranes. Experimental results suggest that the N-terminal 17-amino-acid sequence (htt(NT)) positioned just before the Q(N) region is important for the binding of huntingtin to membranes. Through all-atom explicit solvent molecular dynamics simulations, we unveil the structure and dynamics of the htt(NT)Q(N) fragment on a phospholipid membrane at the atomic level. We observe that the insertion dynamics of this peptide can be described by four main steps-approach, reorganization, anchoring, and insertion-that are very diverse at the atomic level. On the membrane, the htt(NT) peptide forms a stable α-helix essentially parallel to the membrane with its nonpolar side-chains-mainly Leu-4, Leu-7, Phe-11 and Leu-14-positioned in the hydrophobic core of the membrane. Salt-bridges involving Glu-5, Glu-12, Lys-6, and Lys-15, as well as hydrogen bonds involving Thr-3 and Ser-13 with the phospholipids also stabilize the structure and orientation of the htt(NT) peptide. These observations do not significantly change upon adding the Q(N) region whose role is rather to provide, through its hydrogen bonds with the phospholipids' head group, a stable scaffold facilitating the partitioning of the htt(NT) region in the membrane. Moreover, by staying accessible to the solvent, the amyloidogenic Q(N) region could also play a key role for the oligomerization of htt(NT)Q(N) on phospholipid membranes.

The N-terminal sequence of prion protein consists an epitope specific to the abnormal isoform of prion protein (PrP(Sc))

The conformation of abnormal prion protein (PrP^Sc) differs from that of cellular prion protein (PrP^C), but the precise characteristics of PrP^Sc remain to be elucidated. To clarify the properties of native PrP^Sc, we attempted to generate novel PrP^Sc-specific monoclonal antibodies (mAbs) by immunizing PrP-deficient mice with intact PrP^Sc purified from bovine spongiform encephalopathy (BSE)-affected mice. The generated mAbs 6A12 and 8D5 selectivity precipitated PrP^Sc from the brains of prion-affected mice, sheep, and cattle, but did not precipitate PrP^C from the brains of healthy animals. In histopathological analysis, mAbs 6A12 and 8D5 strongly reacted with prion-affected mouse brains but not with unaffected mouse brains without antigen retrieval. Epitope analysis revealed that mAbs 8D5 and 6A12 recognized the PrP subregions between amino acids 31–39 and 41–47, respectively. This indicates that a PrP^Sc-specific epitope exists in the N-terminal region of PrP^Sc, and mAbs 6A12 and 8D5 are powerful tools with which to detect native and intact PrP^Sc. We found that the ratio of proteinase K (PK)-sensitive PrP^Sc to PK-resistant PrP^Sc was constant throughout the disease time course.

Species-barrier phenomenon in prion transmissibility from a viewpoint of protein science

Transmissible spongiform encephalopathies (TSEs), or prion diseases, are fatal infectious neurodegenerative disorders. Their causative agents are prions, which are composed of disease-associated forms of prion protein (PrP(Sc)). Naturally occurring cases of TSEs are found in several mammalian species including humans, sheep, goats, minks, cattle and deer. Prions are also experimentally transmissible to other mammals such as mice, hamsters and monkeys, but interspecies transmission is often inefficient due to the 'species-barrier'. Studies have suggested that the barrier is not only simply determined by differences in amino acid sequences of cellular PrP (PrP(C)) among animal species, but also by prion strains which are closely associated with conformational properties of PrP(Sc) aggregates. Although the conformational properties of PrP(Sc) remain largely unknown, recent investigation of local structures of PrP(C) and, in particular, structural modelling of PrP(Sc) aggregates have provided molecular insight into this field. In this review, we discuss the species-barrier phenomenon in terms of the protein science.

Assessing polyglutamine conformation in the nucleating event by molecular dynamics simulations

Polyglutamine (polyQ) diseases comprise a group of dominantly inherited pathology caused by an expansion of an unstable polyQ stretch which is presumed to form β-sheets. Similar to other amyloid pathologies, polyQ amyloidogenesis occurs via a nucleated polymerization mechanism, and proceeds through energetically unfavorable nucleus whose existence and structure are difficult to detect. Here, we use atomistic molecular dynamics simulations in explicit solvent to assess the conformation of the polyQ stretch in the nucleus that initiates polyQ fibrillization. Comparison of the kinetic stability of various structures of polyQ peptide with a Q-length in the pathological range (Q40) revealed that steric zipper or nanotube-like structures (β-nanotube or β-pseudohelix) are not kinetically stable enough to serve as a template to initiate polyQ fibrillization as opposed to β-hairpin-based (β-sheet and β-sheetstack) or α-helical conformations. The selection of different structures of the polyQ stretch in the aggregation-initiating event may provide an alternative explanation for polyQ aggregate polymorphism.

Sequence-dependent stability test of a left-handed β-helix motif

The left-handed β-helix (LHBH) is an intriguing, rare structural pattern in polypeptides that has been implicated in the formation of amyloid aggregates. We used accurate all-atom replica-exchange molecular dynamics (REMD) simulations to study the relative stability of diverse sequences in the LHBH conformation. Ensemble-average coordinates from REMD served as a scoring criterion to identify sequences and threadings optimally suited to the LHBH, as in a fold recognition paradigm. We examined the repeatability of our REMD simulations, finding that single simulations can be reliable to a quantifiable extent. We find expected behavior for the positive and negative control cases of a native LHBH and intrinsically disordered sequences, respectively. Polyglutamine and a designed hexapeptide repeat show remarkable affinity for the LHBH motif. A structural model for misfolded murine prion protein was also considered, and showed intermediate stability under the given conditions. Our technique is found to be an effective probe of LHBH stability, and promises to be scalable to broader studies of this and potentially other novel or rare motifs. The superstable character of the designed hexapeptide repeat suggests theoretical and experimental follow-ups.

Mouse prion protein (PrP) segment 100 to 104 regulates conversion of PrP(C) to PrP(Sc) in prion-infected neuroblastoma cells

Prion diseases are characterized by the replicative propagation of disease-associated forms of prion protein (PrP(Sc); PrP refers to prion protein). The propagation is believed to proceed via two steps; the initial binding of the normal form of PrP (PrP(C)) to PrP(Sc) and the subsequent conversion of PrP(C) to PrP(Sc). We have explored the two-step model in prion-infected mouse neuroblastoma (ScN2a) cells by focusing on the mouse PrP (MoPrP) segment 92-GGTHNQWNKPSKPKTN-107, which is within a region previously suggested to be part of the binding interface or shown to differ in its accessibility to anti-PrP antibodies between PrP(C) and PrP(Sc). Exchanging the MoPrP segment with the corresponding chicken PrP segment (106-GGSYHNQKPWKPPKTN-121) revealed the necessity of MoPrP residues 99 to 104 for the chimeras to achieve the PrP(Sc) state, while segment 95 to 98 was replaceable with the chicken sequence. An alanine substitution at position 100, 102, 103, or 104 of MoPrP gave rise to nonconvertible mutants that associated with MoPrP(Sc) and interfered with the conversion of endogenous MoPrP(C). The interference was not evoked by a chimera (designated MCM2) in which MoPrP segment 95 to 104 was changed to the chicken sequence, though MCM2 associated with MoPrP(Sc). Incubation of the cells with a synthetic peptide composed of MoPrP residues 93 to 107 or alanine-substituted cognates did not inhibit the conversion, whereas an anti-P8 antibody recognizing the above sequence in PrP(C) reduced the accumulation of PrP(Sc) after 10 days of incubation of the cells. These results suggest the segment 100 to 104 of MoPrP(C) plays a key role in conversion after binding to MoPrP(Sc).

Fibrillogenesis of huntingtin and other glutamine containing proteins

This chapter focuses on the aggregation of glutamine containing peptides and proteins with an emphasis on huntingtin protein, whose aggregation leads to the development of Huntington's disease. The kinetics that leads to the formation of amyloids, the structure of aggregates of various types and the morphological mechanical properties of amyloid fibrils are described. The kinetics of amyloid fibril formation has been proposed to follow a nucleation dependent polymerization model, dependent upon the size of the nucleus. This model and the effect of the polyglutamine length on the nucleus size are reviewed. Aggregate structure is characterized at two different levels. The atomic-scale resolution structure of fibrillar and crystalline aggregates of polyglutamine containing proteins and peptides was determined by X-ray crystallography and solid-state nuclear magnetic resonance (NMR). The chapter outlines the results obtained by both these techniques. Atomic force microscopy (AFM) was instrumental in elucidating the morphology of fibrils, their organization and assembly. The chapter also discusses the high stability of amyloid fibrils, including their mechanical properties as revealed by AFM.

Origins and evolution of the HET-s prion-forming protein: searching for other amyloid-forming solenoids

The HET-s prion-forming domain from the filamentous fungus Podospora anserina is gaining considerable interest since it yielded the first well-defined atomic structure of a functional amyloid fibril. This structure has been identified as a left-handed beta solenoid with a triangular hydrophobic core. To delineate the origins of the HET-s prion-forming protein and to discover other amyloid-forming proteins, we searched for all homologs of the HET-s protein in a database of protein domains and fungal genomes, using a combined application of HMM, psi-blast and pGenThreader techniques, and performed a comparative evolutionary analysis of the N-terminal alpha-helical domain and the C-terminal prion-forming domain of HET-s. By assessing the tandem evolution of both domains, we observed that the prion-forming domain is restricted to Sordariomycetes, with a marginal additional sequence homolog in Arthroderma otae as a likely case of horizontal transfer. This suggests innovation and rapid evolution of the solenoid fold in the Sordariomycetes clade. In contrast, the N-terminal domain evolves at a slower rate (in Sordariomycetes) and spans many diverse clades of fungi. We performed a full three-dimensional protein threading analysis on all identified HET-s homologs against the HET-s solenoid fold, and present detailed structural annotations for identified structural homologs to the prion-forming domain. An analysis of the physicochemical characteristics in our set of structural models indicates that the HET-s solenoid shape can be readily adopted in these homologs, but that they are all less optimized for fibril formation than the P. anserina HET-s sequence itself, due chiefly to the presence of fewer asparagine ladders and salt bridges. Our combined structural and evolutionary analysis suggests that the HET-s shape has "limited scope" for amyloidosis across the wider protein universe, compared to the 'generic' left-handed beta helix. We discuss the implications of our findings on future identification of amyloid-forming proteins sharing the solenoid fold.

Structural characterization of polyglutamine fibrils by solid-state NMR spectroscopy

Protein aggregation via polyglutamine stretches occurs in a number of severe neurodegenerative diseases such as Huntington's disease. We have investigated fibrillar aggregates of polyglutamine peptides below, at, and above the toxicity limit of around 37 glutamine residues using solid-state NMR and electron microscopy. Experimental data are consistent with a dry fibril core of at least 70-80 Å in width for all constructs. Solid-state NMR dipolar correlation experiments reveal a largely β-strand character of all samples and point to tight interdigitation of hydrogen-bonded glutamine side chains from different sheets. Two approximately equally frequent populations of glutamine residues with distinct sets of chemical shifts are found, consistent with local backbone dihedral angles compensating for β-strand twist or with two distinct sets of side-chain conformations. Peptides comprising 15 glutamine residues are present as single extended β-strands. Data obtained for longer constructs are most compatible with a superpleated arrangement with individual molecules contributing β-strands to more than one sheet and an antiparallel assembly of strands within β-sheets.

Highly polar environments catalyze the unfolding of PrP C helix 1

The first alpha-helix (H1) likely plays an important role in the conversion of the cellular prion protein (PrP(C)) into its pathogenic isoform (PrP(Sc)). In this conversion, H1 may either have to unfold or may represent a site of intermolecular contact. A recent molecular dynamics simulation suggested that H1 can unfold if it is detached from the protein core (Hirschberger et al. in Biophys J 90:3908, 2006). It has been hypothesized that the high dielectric constant epsilon (S) of the bulk water environment facilitates the unfolding of H1. To check this hypothesis, we performed a number of replica exchange molecular dynamics simulations of an H1 peptide in solvents of different epsilon (S). We found that the equilibrium helix fraction in water is less than 40%, in agreement with previous experimental findings, and that the helix unfolds much faster in water than in less polar solvents. The kinetically stabilizing effect of the organic solvents is largely unspecific and correlates well with their dielectric constant epsilon (S).

LpxA: a natural nanotube

UDP-N-acetylglucosamine 3-O-acyltransferase is a protein with a left-handed parallel beta-helix, which is a natural nanotube. They are associated with unusual high stability. To identify the reason behind the structural stability of beta-helical nanotubular structure, we have performed a total of 4 mus molecular dynamics simulations of the protein in implicit solvent at four different temperatures and monitored the unfolding pathway. The correlation in movement between different regions of the nanotubular structure has been identified from the dynamical cross-correlation map and contribution of some specific residues towards unfolding transition has been identified by principal component analysis. Difference in stability of the three loop regions has also been characterized. Construction of the unfolding conformational energy landscape identifies the probable intermediates that can appear in the unfolding pathway of the protein.

Aggregation of fibrils and plaques in amyloid molecular systems

Amyloidlike proteins form highly organized aggregates, such as fibrils and plaques, preceded by the assembly of a wide range of unstructured oligomers and protofibrils. Despite their importance in a number of human neurodegenerative diseases, a comprehensive understanding of their kinetics and thermodynamics is still missing. We investigate, by computer simulations, a realistic model of amyloid molecules interacting via the experimentally determined Derjaguin-Landau-Verwey-Overbeek potential and derive its phase diagram. We show that fibrils and plaques, along with their precursors, correspond to different equilibrium and metastable thermodynamics phases and discuss the dynamical mechanisms leading to the nucleation and self-assembly of large scale structures.

Heterologous stacking of prion protein peptides reveals structural details of fibrils and facilitates complete inhibition of fibril growth

Fibrils play an important role in the pathogenesis of amyloidosis; however, the underlying mechanisms of the growth process and the structural details of fibrils are poorly understood. Crucial in the fibril formation of prion proteins is the stacking of PrP monomers. We previously proposed that the structure of the prion protein fibril may be similar as a parallel left-handed beta-helix. The beta-helix is composed of spiraling rungs of parallel beta-strands, and in the PrP model residues 105-143 of each PrP monomer can contribute two beta-helical rungs to the growing fibril. Here we report data to support this model. We show that two cyclized human PrP peptides corresponding to residues 105-124 and 125-143, based on two single rungs of the left-handed beta-helical core of the human PrP(Sc) fibril, show spontaneous cooperative fibril growth in vitro by heterologous stacking. Because the structural model must have predictive value, peptides were designed based on the structure rules of the left-handed beta-helical fold that could stack with prion protein peptides to stimulate or to block fibril growth. The stimulator peptide was designed as an optimal left-handed beta-helical fold that can serve as a template for fibril growth initiation. The inhibiting peptide was designed to bind to the exposed rung but frustrate the propagation of the fibril growth. The single inhibitory peptide hardly shows inhibition, but the combination of the inhibitory with the stimulatory peptide showed complete inhibition of the fibril growth of peptide huPrP-(106-126). Moreover, the unique strategy based on stimulatory and inhibitory peptides seems a powerful new approach to study amyloidogenic fibril structures in general and could prove useful for the development of therapeutics.

Analysis of the sequence and structural features of the left-handed beta-helical fold

The left-handed parallel beta-helix (LbetaH) is a structurally repetitive, highly regular, and symmetrical fold formed by coiling of elongated beta-sheets into helical "rungs." This canonical fold has recently received interest as a possible solution to the fibril structure of amyloid and as a building block of self-assembled nanotubular structures. In light of this interest, we aimed to understand the structural requirements of the LbetaH fold. We first sought to determine the sequence characteristics of the repeats by analyzing known structures to identify positional preferences of specific residues types. We then used molecular dynamics simulations to demonstrate the stabilizing effect of successive rungs and the hydrophobic core of the LbetaH. We show that a two-rung structure is the minimally stable LbetaH structure. In addition, we defined the structure-based sequence preference of the LbetaH and undertook a genome-wide sequence search to determine the prevalence of this unique protein fold. This profile-based LbetaH search algorithm predicted a large fraction of LbetaH proteins from microbial origins. However, the relative number of predicted LbetaH proteins per specie was approximately equal across the genomes from prokaryotes to eukaryotes.

Folding of polyglutamine chains

Long polyglutamine chains have been associated with a number of neurodegenerative diseases. These include Huntington's disease, where expanded polyglutamine (PolyQ) sequences longer than 36 residues are correlated with the onset of symptoms. In this paper we study the folding pathway of a 54-residue PolyQ chain into a beta-helical structure. Transition path sampling Monte Carlo simulations are used to generate unbiased reactive pathways between unfolded configurations and the folded beta-helical structure of the polyglutamine chain. The folding process is examined in both explicit water and an implicit solvent. Both models reveal that the formation of a few critical contacts is necessary and sufficient for the molecule to fold. Once the primary contacts are formed, the fate of the protein is sealed and it is largely committed to fold. We find that, consistent with emerging hypotheses about PolyQ aggregation, a stable beta-helical structure could serve as the nucleus for subsequent polymerization of amyloid fibrils. Our results indicate that PolyQ sequences shorter than 36 residues cannot form that nucleus, and it is also shown that specific mutations inferred from an analysis of the simulated folding pathway exacerbate its stability.

Left handed beta helix models for mammalian prion fibrils

We propose models for in vitro grown mammalian prion protein fibrils based upon left handed beta helices formed both from the N-terminal and C-terminal regions of the proteinase resistant infectious prion core. The C-terminal threading onto a beta-helical structure is almost uniquely determined by fixing the cysteine disulfide bond on a helix corner. In comparison to known left handed helical peptides, the resulting model structures have similar stability attributes including relatively low root mean square deviations in all atom molecular dynamics, substantial side-chain-to-side-chain hydrogen bonding, good volume packing fraction, and low hydrophilic/hydrophobic frustration. For the N-terminus, we propose a new threading of slightly more than two turns, which improves upon the above characteristics relative to existing three turn beta-helical models. The N-terminal and C-terminal beta helices can be assembled into eight candidate models for the fibril repeat units, held together by large hinge (order 30 residues) domain swapping, with three amenable to fibril promoting domain swapping via a small (five residue) hinge on the N-terminal side. Small concentrations of the metastable C-terminal beta helix in vivo might play a significant role in templating the infectious conformation and in enhancing conversion kinetics for inherited forms of the disease and explain resistance (for canines) involving hypothesized coupling to the methionine 129 sulfur known to play a role in human disease.

Insights into structure, stability, and toxicity of monomeric and aggregated polyglutamine models from molecular dynamics simulations

Nine genetically inherited neurodegenerative diseases are linked to abnormal expansions of a polyglutamine (polyQ) encoding region. Over the years, several structural models for polyQ regions have been proposed and confuted. The cross-beta-spine steric zipper motif, identified recently for the GNNQQNY peptide, represents an attractive model for amyloid fibers formed by polyQ fragments. Here we report a detailed molecular dynamics investigation of polyQ models assembled by cross-beta-spine steric zipper motifs. Our simulations indicate clearly that these assemblies are very stable. Glutamine side chains contribute strongly to the overall stability of the models by fitting perfectly within the zipper. In contrast to GNNQQNY zipper motifs, hydrogen bonding interactions provide a significant contribution to the overall stability of polyQ models. Molecular dynamics simulations carried out on monomeric polyQ forms (composed by 40-60 residues) show clearly that they can also assume structures stabilized by steric zipper motifs. Based on these findings, we build monomeric polyQ models that can explain recent data on the toxicity exerted by these species. In a more general context, our data suggests that polyQ models with interdigitated side chains can provide a structural rationale to several literature experiments on polyQ formation, stability, and toxicity.

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.

Molecular dynamics simulation study on the structural stabilities of polyglutamine peptides

It is known that Huntington's disease patients commonly have glutamine (Q) repeat sequences longer than a critical length in the coding area of Huntingtin protein in their genes. As the polyglutamine (polyQ) region becomes longer than the critical length, the disease occurs and Huntingtin protein aggregates, both in vitro and in vivo, as suggested by experimental and clinical data. The determination of polyglutamine structure is thus very important for elucidation of the aggregation and disease mechanisms. Here, we perform molecular dynamics calculations on the stability of the structure based on the beta-helix structure suggested by Perutz et al. (2002) [Perutz, M.F., Finch, J.T., Berriman, J., Lesk, A., 2002. Amyloid fibers are water-filled nanotubes. Proc. Natl. Acad. Sci. USA 99, 5591]. We ensure that perfect hydrogen bonds are present between main chains of the beta-helix based on the previous studies, and perform simulations of stretches with 20, 25, 30, 37 and 40 glutamine residues (20Q, 25Q, 30Q, 37Q and 40Q) for the Perutz models with 18.5 and 20 residues per turn (one coil). Our results indicate that the structure becomes more stable with the increase of repeated number of Q, and there is a critical Q number of around 30, above which the structure of the Perutz model is kept stable. In contrast to previous studies, we started molecular dynamics simulations from conformations in which the hydrogen bonds are firmly formed between stacked main chains. This has rendered the initial beta-helix structures of polyQ much more stable for longer time, as compared to those proposed previously. Model calculations for the initial structures of polyQ dimer and tetramer have also been carried out to study a possible mechanism for aggregation.

Flanking polyproline sequences inhibit beta-sheet structure in polyglutamine segments by inducing PPII-like helix structure

Polyglutamine (poly(Q)) expansion is associated with protein aggregation into beta-sheet amyloid fibrils and neuronal cytotoxicity. In the mutant poly(Q) protein huntingtin, associated with Huntington's disease, both aggregation and cytotoxicity may be abrogated by a polyproline (poly(P)) domain flanking the C terminus of the poly(Q) region. To understand structural changes that may occur with the addition of the poly(P) sequence, we synthesized poly(Q) peptides with 3-15 glutamine residues and a corresponding set of poly(Q) peptides flanked on the C terminus by 11 proline residues (poly(Q)-poly(P)), as occurs in the huntingtin sequence. The shorter soluble poly(Q) peptides (three or six glutamine residues) showed polyproline type II-like (PPII)-like helix conformation when examined by circular dichroism spectroscopy and were monomers as judged by size-exclusion chromatography (SEC), while the longer poly(Q) peptides (nine or 15 glutamine residues) showed a beta-sheet conformation by CD and defined oligomers by SEC. Soluble poly(Q)-poly(P) peptides showed PPII-like content but SEC showed poorly defined, overlapping oligomeric peaks, and as judged by CD these peptides retained significant PPII-like structure with increasing poly(Q) length. More importantly, addition of the poly(P) domain increased the threshold for fibril formation to approximately 15 glutamine residues. X-ray diffraction, electron microscopy, and film CD showed that, while poly(Q) peptides with >or=6 glutamine residues formed beta-sheet-rich fibrils, only the longest poly(Q)-poly(P) peptide (15 glutamine residues) did so. From these and other observations, we propose that poly(Q) domains exist in a "tug-of-war" between two conformations, a PPII-like helix and a beta-sheet, while the poly(P) domain is conformationally constrained into a proline type II helix (PPII). Addition of poly(P) to the C terminus of a poly(Q) domain induces a PPII-like structure, which opposes the aggregation-prone beta-sheet. These structural observations may shed light on the threshold phenomenon of poly(Q) aggregation, and support the hypothesized evolution of "protective" poly(P) tracts adjacent to poly(Q) aggregation domains.

Studying the structural properties of polyalanine and polyglutamine peptides

Poly-(Ala) and poly-(Gln) peptides have important biological effects, and can cause various human illnesses and neurodegenerative diseases. Conformational analysis of these homo-oligopeptides (HOPs) was carried out by simulated annealing in order to identify their structural properties regarding secondary structures and intramolecular H-bonding patterns. Poly-(Ala) and poly-(Gln) peptides composed of 7, 10, 14 or 20 amino acids were modelled in both charged and terminally blocked forms. In the case of conformers derived from simulated annealing calculations, the presence of various secondary structural elements (different types of beta-turns, alpha-helix, 3(10)-helix, poly-proline II helix, parallel and antiparallel beta-strands) was investigated. Moreover, the intramolecular H-bonding patterns formed either between the backbone atoms for both HOPs or between the backbone and side-chain atoms for the poly-(Gln) peptides were examined. Our results showed that different secondary structural elements (type I and type III beta-turns, alpha-helix, 3(10)-helix, antiparallel beta-strand) could be observed in both poly-(Ala) and poly-(Gln) peptides and, according to their presence, characteristic H-bonding patterns formed mainly by i<--i+3 and i<--i+4 H-bonds could be found.

Mouse-adapted ovine scrapie prion strains are characterized by different conformers of PrPSc

The agent responsible for prion disease may exist in different forms, commonly referred to as strains, with each carrying the specific information that determines its own distinct biological properties, such as incubation period and lesion profile. Biological strain typing of ovine scrapie isolates by serial passage in conventional mice has shown some diversity in ovine prion strains. However, this biological diversity remains poorly supported by biochemical prion strain typing. The protein-only hypothesis predicts that variation between different prion strains in the same host is manifest in different conformations adopted by PrPSc. Here we have investigated the molecular properties of PrPSc associated with two principal Prnp(a) mouse-adapted ovine scrapie strains, namely, RML and ME7, in order to establish biochemical prion strain typing strategies that may subsequently be used to discriminate field cases of mouse-passaged ovine scrapie isolates. We used a conformation-dependent immunoassay and a conformational stability assay, together with Western blot analysis, to demonstrate that RML and ME7 PrPSc proteins show distinct biochemical and physicochemical properties. Although RML and ME7 PrPSc proteins showed similar resistance to proteolytic digestion, they differed in their glycoform profiles and levels of proteinase K (PK)-sensitive and PK-resistant isoforms. In addition, the PK-resistant core (PrP27-30) of ME7 was conformationally more stable following exposure to guanidine hydrochloride or Sarkosyl than was RML PrP27-30. Our data show that mouse-adapted ovine scrapie strains can be discriminated by their distinct conformers of PrPSc, which provides a basis to investigate their diversity at the molecular level.

A Lumry-Eyring nucleated polymerization model of protein aggregation kinetics: 1. Aggregation with pre-equilibrated unfolding

A mathematical model is presented of the kinetics of non-native protein aggregation that combines Lumry-Eyring and nucleated polymerization (LENP) descriptions. The LENP model is solved for cases in which aggregation rates are slow compared to folding-unfolding equilibration and is shown to be a generalization of a number of previously proposed nucleation-and-growth models for non-native and native protein aggregation. The model solutions exhibit a number of qualitative kinetic regimes. Each regime has a characteristic set of experimental signatures that are related to the relative rates of growth and nucleation as well as to the threshold size at which aggregates condense to form higher-order structures or other phases. Approximate model solutions provide practical rate equations that can be regressed against typical experimental kinetic data to obtain mechanistic parameters characterizing the aggregation pathway. In all kinetic regimes, it is found that observed rate coefficients (kobs) or half-lives (t50) obtained from extent-of-reaction measurements are convolutions of more than one stage in the pathway unless purely seeded growth occurs. Despite this convolution, the combination of apparent reaction order (time domain) and the scaling of kobs or t50 with initial protein concentration provides a means to determine a value for the dominant nucleus size in each case. Additional information, such as equilibrium unfolding thermodynamics and the limiting aggregate size distribution, are required to further deconvolute kobs into intrinsic contributions from nucleation, growth, and conformational changes. The model and analysis are expected to be generally applicable to a wide range of proteins and polypeptides that form non-native aggregates.

Effects of chain length on the aggregation of model polyglutamine peptides: molecular dynamics simulations

Aggregation in the brain of polyglutamine-containing proteins is either a cause or an associated symptom of nine hereditary neurodegenerative disorders including Huntington's disease. The molecular level mechanisms by which these proteins aggregate are still unclear. In an effort to shed light on this important phenomenon, we are investigating the aggregation of model polyglutamine peptides using molecular-level computer simulation with a simplified model of polyglutamine that we have developed. This model accounts for the most important types of intra- and inter-molecular interactions-hydrogen bonding and hydrophobic interactions-while allowing the folding process to be simulated in a reasonable time frame. The model is used to examine the folding of isolated polyglutamine peptides 16, 32, and 48 residues long and the folding and aggregation of systems of 24 model polyglutamine peptides 16, 24, 32, 36, 40, and 48 residues long. Although the isolated polyglutamine peptides did form some alpha and beta backbone-backbone hydrogen bonds they did not have as many of these bonds as they would have if they had folded into a complete alpha helix or beta sheet. In one of the simulations on the isolated polyglutamine peptide 48 residues long, we observed a structure that resembles a beta helix. In the multi-chain simulations we observed amorphous aggregates at low temperatures, ordered aggregates with significant beta sheet character at intermediate temperatures, and random coils at high temperatures. We have found that the temperature at which the model peptides undergo the transition from amorphous aggregates to ordered aggregates and the temperature at which the model peptides undergo the transition from ordered aggregates to random coils increase with increasing chain length. Our finding that the stability of the ordered aggregates increases as the peptide chain length increases may help to explain the experimentally observed relation between polyglutamine tract length and aggregation in vitro and disease progression in vivo. We have also observed in our simulations that the optimal temperature for the formation of beta sheets increases with chain length up to 36 glutamine residues but not beyond. Equivalently, at fixed temperature we find a transition from a region dominated by random coils at chain lengths less than 36 to a region dominated by relatively ordered beta sheet structures at chain lengths greater than 36. Our finding of this critical chain length of 36 glutamine residues is interesting because a critical chain length of 37 glutamine residues has been observed experimentally.

Identification of anti-prion compounds as efficient inhibitors of polyglutamine protein aggregation in a zebrafish model

Several neurodegenerative diseases, including Huntington disease (HD), are associated with aberrant folding and aggregation of polyglutamine (polyQ) expansion proteins. Here we established the zebrafish, Danio rerio, as a vertebrate HD model permitting the screening for chemical suppressors of polyQ aggregation and toxicity. Upon expression in zebrafish embryos, polyQ-expanded fragments of huntingtin (htt) accumulated in large SDS-insoluble inclusions, reproducing a key feature of HD pathology. Real time monitoring of inclusion formation in the living zebrafish indicated that inclusions grow by rapid incorporation of soluble htt species. Expression of mutant htt increased the frequency of embryos with abnormal morphology and the occurrence of apoptosis. Strikingly, apoptotic cells were largely devoid of visible aggregates, suggesting that soluble oligomeric precursors may instead be responsible for toxicity. As in nonvertebrate polyQ disease models, the molecular chaperones, Hsp40 and Hsp70, suppressed both polyQ aggregation and toxicity. Using the newly established zebrafish model, two compounds of the N'-benzylidene-benzohydrazide class directed against mammalian prion proved to be potent inhibitors of polyQ aggregation, consistent with a common structural mechanism of aggregation for prion and polyQ disease proteins.

Isolation and characterization of a proteinase K-sensitive PrPSc fraction

Recent studies have shown that a sizable fraction of PrPSc present in prion-infected tissues is, contrary to previous conceptions, sensitive to digestion by proteinase K (PK). This finding has important implications in the context of diagnosis of prion disease, as PK has been extensively used in attempts to distinguish between PrPSc and PrPC. Even more importantly, PK-sensitive PrPSc (sPrPSc) might be essential to understand the process of conversion and aggregation of PrPC leading to infectivity. We have isolated a fraction of sPrPSc. This material was obtained by differential centrifugation at an intermediate speed of Syrian hamster PrPSc obtained through a conventional procedure based on ultracentrifugation in the presence of detergents. PK-sensitive PrPSc is completely degraded under standard conditions (50 mug/mL of proteinase K at 37 degrees C for 1 h) and can also be digested with trypsin. Centrifugation in a sucrose gradient showed sPrPSc to correspond to the lower molecular weight fractions of the continuous range of oligomers that constitute PrPSc. PK-sensitive PrPSc has the ability to convert PrPC into protease-resistant PrPSc, as assessed by the protein misfolding cyclic amplification assay (PMCA). Limited proteolysis of sPrPSc using trypsin allows for identification of regions that are particularly susceptible to digestion, i.e., are partially exposed and flexible; we have identified as such the regions around residues K110, R136, R151, K220, and R229. PK-sensitive PrPSc isolates should prove useful for structural studies to help understand fundamental issues of the molecular biology of PrPSc and in the quest to design tests to detect preclinical prion disease.

Structures of soluble amyloid oligomers from computer simulations

Alzheimer's, Parkinson's, and Creutzfeldt-Jakob's neurodegenerative diseases are all linked with the assembly of normally soluble proteins into amyloid fibrils. Because of experimental limitations, structural characterization of the soluble oligomers, which form early in the process of fibrillogenesis and are cytotoxic, remains to be determined. In this article, we study the aggregation paths of seven chains of the shortest amyloid-forming peptide, using an activitated method and a reduced atomic representation. Our simulations show that disordered KFFE monomers ultimately form three distinct topologies of similar energy: amorphous oligomers, incomplete rings with beta-barrel character, and cross-beta-sheet structures with the meridional but not the equatorial X-ray fiber reflections. The simulations also shed light on the pathways from misfolded aggregates to fibrillar-like structures. They also underline the multiplicity of building blocks that can lead to the formation of the critical nucleus from which rapid growth of the fibril occurs.

Polyglutamine repeats and beta-helix structure: molecular dynamics study

Neurodegenerative diseases are often associated with the formation of highly insoluble aggregates. Despite the efforts devoted to the characterization of these aggregates, their structure remains elusive. Several neurodegenerative diseases are characterized by the expansion of CAG repeats, which code for Gln. Among the structural models proposed for the aggregates observed in polyQ-linked diseases, the nanotube beta-helix model proposed by Perutz and colleagues Proc Natl Acad Sci U S A 2002;99:5591-5595 has been influential. In the present study, the stability of this beta-helix model has been investigated by performing molecular dynamics simulations on polyQ fragments of different lengths. The results indicate that models shorter than two full beta-helix turns are unstable and collapse toward irregular structures. On the other hand, longer beta-helix models, containing more than 40 residues, achieve a dynamic regular structure. This finding is in line with the observed threshold of Gln repeats (approximately 40) correlated with the insurgence of the disease. Notably, the structure of the final state of the models longer than 40 residues strictly depends on their size. A compact stable ellipsoidal structure is formed by the model made of two full helical turns (41 residues), whereas water filled tubular structures emerge from simulation on longer polypeptides. These results have been interpreted taking into account the experimental data on polyQ aggregates. A structural interpretation of the literature data has been proposed by assuming that different beta-helical models are involved in the different stages of the aggregation process.

Photo-induced crosslinking of prion protein oligomers and prions

Prion diseases are caused by a unique type of infectious agent, which is thought to consist of a misfolded beta-sheeted form of the alpha-helical cellular prion protein (PrPC). This misfolded isoform (PrPSc) tends to form insoluble amyloid-like aggregates, impeding classical structural analysis by X-ray crystallography or NMR. Intermolecular crosslinking may provide a means of stabilizing notoriously elusive oligomers for further analysis and may be used for analyzing aggregate architecture by characterising intermolecular contact sites. Using a photo-induced crosslinking method (PICUP), aggregates of recombinant PrP (rPrP) and PrPSc were linked at interacting surfaces via amino acid side chains. The degree of crosslinking within PrP aggregates was adjustable using varying light intensities and could efficiently be monitored by fluorescence correlation spectroscopy. Specific intermolecular crosslinking of PrPSc molecules was achieved even in crude brain homogenate. Functional studies showed that stabilized aggregates of rPrP did not loose their capacity to induce further protein aggregation and crosslinking of PrPSc did not alter significantly the level of infectivity, indicating that photo-induced covalent linkage of PrPSc does not destruct surfaces important for prion propagation.

Rare Large Scale Subdomain Motions in Prion Protein can Initiate Aggregation

The prion protein is thought to induce prion diseases by changing its conformation from the cellular form, PrP(C), into the infectious Scrapie-form, PrP(Sc). Little is known about the structural and dynamical features of this conformational change. We here introduce a novel concept that involves rare large scale motions between the subdomains beta1-alpha1-beta2 and alpha2-alpha3 in the carboxy-terminal, globular part of PrP. The interface between these two subdomains carries most pathogenic mutations known to be associated with prion diseases. Based on computational simulations as well as experimental results we propose that such a large scale motion subsequently destabilizes large parts of the cellular conformer PrP(C), thus, rendering it prone to structural rearrangements, including aggregation of now partially unfolded parts of the PrP sequence. We hypothesize that such large scale motions occur as a rare event even under equilibrium conditions and that the interaction of such partially destabilized PrP(C)-conformers, which we named PrP(C*), contributes to the formation of pathogenic oligomeric species of the prion protein.