Probing the instabilities in the dynamics of helical fragments from mouse PrPC

Mutations of evolutionarily conserved aromatic residues suggest that misfolding of the mouse prion protein may commence in multiple ways

The misfolding of the mammalian prion protein from its α-helix rich cellular isoform to its β-sheet rich infectious isoform is associated with several neurodegenerative diseases. The determination of the structural mechanism by which misfolding commences, still remains an unsolved problem. In the current study, native-state hydrogen exchange coupled with mass spectrometry has revealed that the N state of the mouse prion protein (moPrP) at pH 4 is in dynamic equilibrium with multiple partially unfolded forms (PUFs) capable of initiating misfolding. Mutation of three evolutionarily conserved aromatic residues, Tyr168, Phe174, and Tyr217 present at the interface of the β2-α2 loop and the C-terminal end of α3 in the structured C-terminal domain of moPrP significantly destabilize the native state (N) of the protein. They also reduce the free energy differences between the N state and two PUFs identified as PUF1 and PUF2**. It is shown that PUF2** in which the β2-α2 loop and the C-terminal end of α3 are disordered, has the same stability as the previously identified PUF2*, but to have a very different structure. Misfolding can commence from both PUF1 and PUF2**, as it can from PUF2*. Hence, misfolding can commence and proceed in multiple ways from structurally distinct precursor conformations. The increased extents to which PUF1 and PUF2** are populated at equilibrium in the case of the mutant variants, greatly accelerate their misfolding. The results suggest that the three aromatic residues may have been evolutionarily selected to impede the misfolding of moPrP.

Intrinsic determinants of prion protein neurotoxicity in Drosophila: from sequence to (dys)function

Prion diseases are fatal brain disorders characterized by deposition of insoluble isoforms of the prion protein (PrP). The normal and pathogenic structures of PrP are relatively well known after decades of studies. Yet our current understanding of the intrinsic determinants regulating PrP misfolding are largely missing. A 3D subdomain of PrP comprising the β2-α2 loop and helix 3 contains high sequence and structural variability among animals and has been proposed as a key domain regulating PrP misfolding. We combined in vivo work in Drosophila with molecular dynamics (MD) simulations, which provide additional insight to assess the impact of candidate substitutions in PrP from conformational dynamics. MD simulations revealed that in human PrP WT the β2-α2 loop explores multiple β-turn conformations, whereas the Y225A (rabbit PrP-like) substitution strongly favors a 310-turn conformation, a short right-handed helix. This shift in conformational diversity correlates with lower neurotoxicity in flies. We have identified additional conformational features and candidate amino acids regulating the high toxicity of human PrP and propose a new strategy for testing candidate modifiers first in MD simulations followed by functional experiments in flies. In this review we expand on these new results to provide additional insight into the structural and functional biology of PrP through the prism of the conformational dynamics of a 3D domain in the C-terminus. We propose that the conformational dynamics of this domain is a sensitive measure of the propensity of PrP to misfold and cause toxicity. This provides renewed opportunities to identify the intrinsic determinants of PrP misfolding through the contribution of key amino acids to different conformational states by MD simulations followed by experimental validation in transgenic flies.

Odd-even disparity in the population of slipped hairpins in RNA repeat sequences with implications for phase separation

Low-complexity nucleotide repeat sequences, which are implicated in several neurological disorders, undergo liquid-liquid phase separation (LLPS) provided the number of repeat units, n, exceeds a critical value. Here, we establish a link between the folding landscapes of the monomers of trinucleotide repeats and their propensity to self-associate. Simulations using a coarse-grained Self-Organized Polymer (SOP) model for (CAG)n repeats in monovalent salt solutions reproduce experimentally measured melting temperatures, which are available only for small n. By extending the simulations to large n, we show that the free-energy gap, ΔGS, between the ground state (GS) and slipped hairpin (SH) states is a predictor of aggregation propensity. The GS for even n is a perfect hairpin (PH), whereas it is a SH when n is odd. The value of ΔGS (zero for odd n) is larger for even n than for odd n. As a result, the rate of dimer formation is slower in (CAG)30 relative to (CAG)31, thus linking ΔGS to RNA-RNA association. The yield of the dimer decreases dramatically, compared to the wild type, in mutant sequences in which the population of the SH decreases substantially. Association between RNA chains is preceded by a transition to the SH even if the GS is a PH. The finding that the excitation spectrum-which depends on the exact sequence, n, and ionic conditions-is a predictor of self-association should also hold for other RNAs (mRNA for example) that undergo LLPS.

Conversion of Helix 1 into a Loop in Prion Protein Misfolding

Cellular prion protein PrP^C consists of three α-helices, one β-sheet, and an unstructured N-terminal domain. Misfolding of this protein into the scrapie form (PrP^Sc) increases dramatically the β-sheet content. H1 is the most stable helix on PrP^C and contains an unusual number of hydrophilic amino acids. Its fate in PrP^Sc is not clear. We performed replica exchange molecular dynamics simulations on H1 alone, H1 together with an N-terminally flanking H1B1 loop and H1 in complex with other hydrophilic regions of the prion protein. In the presence of the H⁹⁹SQWNKPSKPKTNMK¹¹³ sequence, H1 is almost completely converted to a loop structure stabilized by a network of salt bridges. On the other hand, H1 retains its helical structure alone or together with the other sequences considered in this study. We carried out an additional simulation by restraining the distance between the two ends of H1, mimicking a possible geometric restriction by the rest of the protein. Even though the loop was the major conformation, a significant amount of helical structure was also observed. This suggests that the interaction with H⁹⁹SQWNKPSKPKTNMK¹¹³ is necessary for complete helix-to-loop conversion.

Molecular simulation unravels the amyloidogenic misfolding of nascent ApoA1 protein, driven by deleterious point mutations occurring in between 170-178 hotspot region

Protein ApoA1 is extensively studied for its role in lipid metabolism. Its seedy dark side of amyloid formulation remains relatively understudied yet. Due to genetic mutations, the protein pathologically misshapes into its amyloid form that gets accumulated in various organs, including the heart. To contrive effective therapeutics against this debilitating congenital disorder, it is imperative to comprehend the structural ramifications induced by mutations in APoA1's dynamic conformation. Till now, several point mutations have been implicated in ApoA1's amyloidosis, although only a handful has been examined considerably. Especially, the single nucleotide polymorphisms (SNPs) that occur in-between 170-178 mutation hotspot site of APoA1 needs to be investigated, since most of them are culpable of amyloid deposition in the heart. To that effect, in the present study, we have computationally quantified and studied the ApoA1's biomolecular modifications fostered by SNPs in the 170-178 mutation hotspot. Findings from discrete molecular dynamics simulation studies indicate that the SNPs have noticeably steered the ApoA1's behaviour from its native structural dynamics. Analysis of protein's secondary structural changes exhibits a considerable change upon mutations. Further, subjecting the protein structures to simulated thermal denaturation shows increased resistance to denaturation among mutants when compared to native. Further, normal mode analysis of protein's dynamic motion also shows discrepancy in its dynamic structural change upon SNP. These structural digressions induced by SNPs can very well be the biomolecular incendiary that drives ApoA1 into its amyloidogenesis. And, understanding these structural modifications initiates a better understanding of SNP's amyloidogenic pathology on APoA1.Communicated by Ramaswamy H. Sarma.

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.

Simulation analysis of selective alanine mutation effect on stability of human prion protein

Prion diseases are neurodegenerative disorders caused by spongiform degeneration of the brain. Understanding the fundamental mechanism of prion protein aggregation caused by mutations is very crucial to resolve the pathology of prion diseases. To help understand the roles of individual residues on the stability of the human prion protein, the computational method of free energy simulations based on atomistic molecular dynamics trajectories is applied to Phe175 → Ala, Val180 → Ala, and Val209 → Ala mutations of the human prion protein. The simulations show that all three alanine mutations destabilize the human prion protein. The calculated free energy change differences, ΔΔG, for the Phe175 → Ala, Val180 → Ala, and Val209 → Ala mutations are in good agreement with the experimental values. The significant destabilizing effects on the mutants relative to the wild-type protein arise from van der Waals terms. Furthermore, our free energy decomposition analysis shows that the major contribution to destabilizing the V180A and V209A mutants relative to the wild-type protein is originated from van der Waals interactions from residues near the mutation sites. In contrast, the contribution to destabilizing the F175A mutant is mainly caused by van der Waals interactions from residues near and far away from the mutation site. Our results show that the free energy simulation with a thermodynamic integration approach for selected alanine scanning mutations is beneficial for understanding the detailed mechanism of human prion protein destabilization, specific residues' role, and the hydrophobic effect on protein stability.Communicated by Ramaswamy H. Sarma.

Identification of Two Early Folding Stage Prion Non-Local Contacts Suggested to Serve as Key Steps in Directing the Final Fold to Be Either Native or Pathogenic

The initial steps of the folding pathway of the C-terminal domain of the murine prion protein mPrP(90–231) are predicted based on the sequential collapse model (SCM). A non-local dominant contact is found to form between the connecting region between helix 1 and β-sheet 1 and the C-terminal region of helix 3. This non-local contact nucleates the most populated molten globule-like intermediate along the folding pathway. A less stable early non-local contact between segments 120–124 and 179–183, located in the middle of helix 2, promotes the formation of a less populated molten globule-like intermediate. The formation of the dominant non-local contact constitutes an example of the postulated Nature’s Shortcut to the prion protein collapse into the native structure. The possible role of the less populated molten globule-like intermediate is explored as the potential initiation point for the folding for three pathogenic mutants (T182A, I214V, and Q211P in mouse prion numbering) of the prion protein.

Methionine oxidation within the prion protein

Prion diseases are characterized by the self-templated misfolding of the cellular prion protein (PrPC) into infectious aggregates (PrPSc). The detailed molecular basis of the misfolding and aggregation of PrPC remains incompletely understood. It is believed that the transient misfolding of PrPC into partially structured intermediates precedes the formation of insoluble protein aggregates and is a critical component of the prion misfolding pathway. A number of environmental factors have been shown to induce the destabilization of PrPC and promote its initial misfolding. Recently, oxidative stress and reactive oxygen species (ROS) have emerged as one possible mechanism by which the destabilization of PrPC can be induced under physiological conditions. Methionine residues are uniquely vulnerable to oxidation by ROS and the formation of methionine sulfoxides leads to the misfolding and subsequent aggregation of PrPC. Here, we provide a review of the evidence for the oxidation of methionine residues in PrPC and its potential role in the formation of pathogenic prion aggregates.

Insight From Animals Resistant to Prion Diseases: Deciphering the Genotype - Morphotype - Phenotype Code for the Prion Protein

Prion diseases are a group of neurodegenerative diseases endemic in humans and several ruminants caused by the misfolding of native prion protein (PrP) into pathological conformations. Experimental work and the mad-cow epidemic of the 1980s exposed a wide spectrum of animal susceptibility to prion diseases, including a few highly resistant animals: horses, rabbits, pigs, and dogs/canids. The variable susceptibility to disease offers a unique opportunity to uncover the mechanisms governing PrP misfolding, neurotoxicity, and transmission. Previous work indicates that PrP-intrinsic differences (sequence) are the main contributors to disease susceptibility. Several residues have been cited as critical for encoding PrP conformational stability in prion-resistant animals, including D/E159 in dog, S167 in horse, and S174 in rabbit and pig PrP (all according to human numbering). These amino acids alter PrP properties in a variety of assays, but we still do not clearly understand the structural correlates of PrP toxicity. Additional insight can be extracted from comparative structural studies, followed by molecular dynamics simulations of selected mutations, and testing in manipulable animal models. Our working hypothesis is that protective amino acids generate more compact and stable structures in a C-terminal subdomain of the PrP globular domain. We will explore this idea in this review and identify subdomains within the globular domain that may hold the key to unravel how conformational stability and disease susceptibility are encoded in PrP.

A seven-residue deletion in PrP leads to generation of a spontaneous prion formed from C-terminal C1 fragment of PrP

Prions result from a drastic conformational change of the host-encoded cellular prion protein (PrP), leading to the formation of β-sheet-rich, insoluble, and protease-resistant self-replicating assemblies (PrP^Sc). The cellular and molecular mechanisms involved in spontaneous prion formation in sporadic and inherited human prion diseases or equivalent animal diseases are poorly understood, in part because cell models of spontaneously forming prions are currently lacking. Here, extending studies on the role of the H2 α-helix C terminus of PrP, we found that deletion of the highly conserved ¹⁹⁰HTVTTTT¹⁹⁶ segment of ovine PrP led to spontaneous prion formation in the RK13 rabbit kidney cell model. On long-term passage, the mutant cells stably produced proteinase K (PK)-resistant, insoluble, and aggregated assemblies that were infectious for naïve cells expressing either the mutant protein or other PrPs with slightly different deletions in the same area. The electrophoretic pattern of the PK-resistant core of the spontaneous prion (Δ^Spont) contained mainly C-terminal polypeptides akin to C1, the cell-surface anchored C-terminal moiety of PrP generated by natural cellular processing. RK13 cells expressing solely the Δ190-196 C1 PrP construct, in the absence of the full-length protein, were susceptible to Δ^Spont prions. Δ^Spont infection induced the conversion of the mutated C1 into a PK-resistant and infectious form perpetuating the biochemical characteristics of Δ^Spont prion. In conclusion, this work provides a unique cell-derived system generating spontaneous prions and provides evidence that the 113 C-terminal residues of PrP are sufficient for a self-propagating prion entity.

Unraveling the Molecular Mechanism of Prion H2 C-Terminus Misfolding by Metadynamics Simulations

Conformational transition from the normal cellular form of prion protein (PrP^C) to the pathogenic "scrapie" form (PrP^Sc) is considered to be a key event in the occurrence of prion disease. Additionally, the H2 C-terminus is widely considered to be a vital site for PrP conformational transition, which can be used as an important region to explore the potential mechanism of PrP misfolding. Therefore, to study the misfolding mechanism of PrP, 500 ns well-tempered metadynamics simulations were performed by focusing on the H2 C-terminus of PrP. For comparison, three systems were designed in total, including PrP in neutral and acidic conditions, as well as H187R mutant. The resulting free energy surfaces (FESs) obtained from metadynamics simulations reveal that acidic conditions and H187R mutation can facilitate PrP misfolding by decreasing free energy barriers for conformational transition and forming energy stable conformational states. Further analyses aimed at H2 C-terminus show that due to the increase of positive charge on residue 187 in both acidic and H187R systems, the electrostatic repulsion of residue 187 and R136/R156 increases greatly, which disrupts the electrostatic interaction network around H2 C-terminus and exposes the hydrophobic core to the solvent. Taken together, acidic conditions and H187R mutation can accelerate PrP misfolding mainly by forming more energetically stable metastable conformations with lower free energy barriers, and electrostatic network disruption involving residue 187 drives the initial misfolding of H2 C-terminus. This study provides quantitative insight into the related function of the H2 C-terminus in the PrP misfolding process, which may guide H2 C-terminus mediated drug design in the future.

Insights into the Bidirectional Properties of the Sheep-Deer Prion Transmission Barrier

The large chronic wasting disease (CWD)-affected cervid population in the USA and Canada, and the risk of the disease being transmitted to humans through intermediate species, is a highly worrying issue that is still poorly understood. In this case, recombinant protein misfolding cyclic amplification was used to determine, in vitro, the relevance of each individual amino acid on cross-species prion transmission. Others and we have found that the β2-α2 loop is a key modulator of transmission barriers between species and markedly influences infection by sheep scrapie, bovine spongiform encephalopathy (BSE), or elk CWD. Amino acids that differentiate ovine and deer normal host prion protein (PrP^C) and associated with structural rigidity of the loop β2-α2 (S173N, N177T) appear to confer resistance to some prion diseases. However, addition of methionine at codon 208 together with the previously described rigid loop substitutions seems to hide a key in this species barrier, as it makes sheep recombinant prion protein highly susceptible to CWD-induced misfolding. These studies indicate that interspecies prion transmission is not only governed just by the β2-α2 loop amino acid sequence but also by its interactions with the α3-helix as shown by substitution I208M. Transmissible spongiform encephalopathies, characterized by long incubation periods and spongiform changes associated with neuronal loss in the brain, have been described in several mammalian species appearing either naturally (scrapie in sheep and goats, bovine spongiform encephalopathy in cattle, chronic wasting disease in cervids, Creutzfeldt-Jakob disease in humans) or by experimental transmission studies (scrapie in mice and hamsters). Much of the pathogenesis of the prion diseases has been determined in the last 40 years, such as the etiological agent or the fact that prions occur as different strains that show distinct biological and physicochemical properties. However, there are many unanswered questions regarding the strain phenomenon and interspecies transmissibility. To assess the risk of interspecies transmission between scrapie and chronic wasting disease, an in vitro prion propagation method has been used. This technique allows to predict the amino acids preventing the transmission between sheep and deer prion diseases.

Structural Determinants of the Prion Protein N-Terminus and Its Adducts with Copper Ions

The N-terminus of the prion protein is a large intrinsically disordered region encompassing approximately 125 amino acids. In this paper, we review its structural and functional properties, with a particular emphasis on its binding to copper ions. The latter is exploited by the region’s conformational flexibility to yield a variety of biological functions. Disease-linked mutations and proteolytic processing of the protein can impact its copper-binding properties, with important structural and functional implications, both in health and disease progression.

Decoupling of size and shape fluctuations in heteropolymeric sequences reconciles discrepancies in SAXS vs. FRET measurements

Unfolded states of proteins and native states of intrinsically disordered proteins (IDPs) populate heterogeneous conformational ensembles in solution. The average sizes of these heterogeneous systems, quantified by the radius of gyration (R_G ), can be measured by small-angle X-ray scattering (SAXS). Another parameter, the mean dye-to-dye distance (R_E ) for proteins with fluorescently labeled termini, can be estimated using single-molecule Förster resonance energy transfer (smFRET). A number of studies have reported inconsistencies in inferences drawn from the two sets of measurements for the dimensions of unfolded proteins and IDPs in the absence of chemical denaturants. These differences are typically attributed to the influence of fluorescent labels used in smFRET and to the impact of high concentrations and averaging features of SAXS. By measuring the dimensions of a collection of labeled and unlabeled polypeptides using smFRET and SAXS, we directly assessed the contributions of dyes to the experimental values R_G and R_E For chemically denatured proteins we obtain mutual consistency in our inferences based on R_G and R_E , whereas for IDPs under native conditions, we find substantial deviations. Using computations, we show that discrepant inferences are neither due to methodological shortcomings of specific measurements nor due to artifacts of dyes. Instead, our analysis suggests that chemical heterogeneity in heteropolymeric systems leads to a decoupling between R_E and R_G that is amplified in the absence of denaturants. Therefore, joint assessments of R_G and R_E combined with measurements of polymer shapes should provide a consistent and complete picture of the underlying ensembles.

Mechanism of Unfolding of Human Prion Protein

Misfolding and aggregation of prion proteins are associated with several neurodegenerative diseases. Therefore, understanding the mechanism of the misfolding process is of enormous interest in the scientific community. It has been speculated and widely discussed that the native cellular prion protein (PrP^C) form needs to undergo substantial unfolding to a more stable PrP^C* state, which may further oligomerize into the toxic scrapie (PrP^Sc) form. Here, we have studied the mechanism of the unfolding of the human prion protein (huPrP) using a set of extensive well-tempered metadynamics simulations. Through multiple microsecond-long metadynamics simulations, we find several possible unfolding pathways. We show that each pathway leads to an unfolded state of lower free energy than the native state. Thus, our study may point to the signature of a PrP^C* form that corresponds to a global minimum on the conformational free-energy landscape. Moreover, we find that these global minima states do not involve an increased β-sheet content, as was assumed to be a signature of PrP^Sc formation in previous simulation studies. We have further analyzed the origin of metastability of the PrP^C form through free-energy surfaces of the chopped helical segments to show that the helices, particularly H2 and H3 of the prion protein, have the tendency to form either a random coil or a β-structure. Therefore, the secondary structural elements of the prion protein are only weakly stabilized by tertiary contacts and solvation forces so that relatively weak perturbations induced by temperature, pressure, pH, and so forth can lead to substantial unfolding with characteristics of intrinsically disordered proteins.

Modulation of prion polymerization and toxicity by rationally designed peptidomimetics

Misfolding and aggregation of cellular prion protein is associated with a large array of neurological disorders commonly called the transmissible spongiform encephalopathies. Designing inhibitors against prions has remained a daunting task owing to limited information about mechanism(s) of their pathogenic self-assembly. Here, we explore the anti-prion properties of a combinatorial library of bispidine-based peptidomimetics (BPMs) that conjugate amino acids with hydrophobic and aromatic side chains. Keeping the bispidine unit unaltered, a series of structurally diverse BPMs were synthesized and tested for their prion-modulating properties. Administration of Leu- and Trp-BPMs delayed and completely inhibited the amyloidogenic conversion of human prion protein (HuPrP), respectively. We found that each BPM induced the HuPrP to form unique oligomeric nanostructures differing in their biophysical properties, cellular toxicities and response to conformation-specific antibodies. While Leu-BPMs were found to stabilize the oligomers, Trp-BPMs effected transient oligomerization, resulting in the formation of non-toxic, non-fibrillar aggregates. Yet another aromatic residue, Phe, however, accelerated the aggregation process in HuPrP. Molecular insights obtained through MD (molecular dynamics) simulations suggested that each BPM differently engages a conserved Tyr 169 residue at the α2–β2 loop of HuPrP and affects the stability of α2 and α3 helices. Our results demonstrate that this new class of molecules having chemical scaffolds conjugating hydrophobic/aromatic residues could effectively modulate prion aggregation and toxicity.

Detailed computational analysis revealed mutation V210I on PrP induced conformational conversion on β2–α2 loop and α2–α3

The development of fatal transmissible spongiform encephalopathies (TSE) is associated with the conformational conversion of the normal cellular prion protein, PrP^C, into its pathogenic isoform, PrP^Sc.

Partially Unfolded Forms of the Prion Protein Populated under Misfolding-promoting Conditions: CHARACTERIZATION BY HYDROGEN EXCHANGE MASS SPECTROMETRY AND NMR

The susceptibility of the cellular prion protein (PrP(C)) to convert to an alternative misfolded conformation (PrP(Sc)), which is the key event in the pathogenesis of prion diseases, is indicative of a conformationally flexible native (N) state. In the present study, hydrogen-deuterium exchange (HDX) in conjunction with mass spectrometry and nuclear magnetic resonance spectroscopy were used for the structural and energetic characterization of the N state of the full-length mouse prion protein, moPrP(23-231), under conditions that favor misfolding. The kinetics of HDX of 34 backbone amide hydrogens in the N state were determined at pH 4. In contrast to the results of previous HDX studies on the human and Syrian hamster prion proteins at a higher pH, various segments of moPrP were found to undergo different extents of subglobal unfolding events at pH 4, a pH at which the protein is known to be primed to misfold to a β-rich conformation. No residual structure around the disulfide bond was observed for the unfolded state at pH 4. The N state of the prion protein was observed to be at equilibrium with at least two partially unfolded forms (PUFs). These PUFs, which are accessed by stochastic fluctuations of the N state, have altered surface area exposure relative to the N state. One of these PUFs resembles a conformation previously implicated to be an initial intermediate in the conversion of monomeric protein into misfolded oligomer at pH 4.

Molecular Mechanism of the Misfolding and Oligomerization of the Prion Protein: Current Understanding and Its Implications

Prion diseases, also known as transmissible spongiform encephalopathies, make up a group of fatal neurodegenerative disorders linked with the misfolding and aggregation of the prion protein (PrP). Although it is not yet understood how the misfolding of PrP induces neurodegeneration, it is widely accepted that the formation of misfolded prion protein (termed PrP(Sc)) is both the triggering event in the disease and the main component of the infectious agent responsible for disease transmission. Despite the clear involvement of PrP(Sc) in prion diseases, the exact composition of PrP(Sc) is not yet well-known. Recent studies show that misfolded oligomers of PrP could, however, be responsible for neurotoxicity and/or infectivity in the prion diseases. Hence, understanding the molecular mechanism of formation of the misfolded oligomers of PrP is critical for developing an understanding about the prion diseases and for developing anti-prion therapeutics. This review discusses recent advances in understanding the molecular mechanism of misfolded oligomer formation by PrP and its implications for the development of anti-prion therapeutics.

Structural and dynamic properties of the human prion protein

Prion diseases involve the conformational conversion of the cellular prion protein (PrP(C)) to its misfolded pathogenic form (PrP(Sc)). To better understand the structural mechanism of this conversion, we performed extensive all-atom, explicit-solvent molecular-dynamics simulations for three structures of the wild-type human PrP (huPrP) at different pH values and temperatures. Residue 129 is polymorphic, being either Met or Val. Two of the three structures have Met in position 129 and the other has Val. Lowering the pH or raising the temperature induced large conformational changes of the C-terminal globular domain and increased exposure of its hydrophobic core. In some simulations, HA and its preceding S1-HA loop underwent large displacements. The C-terminus of HB was unstable and sometimes partially unfolded. Two hydrophobic residues, Phe-198 and Met-134, frequently became exposed to solvent. These conformational changes became more dramatic at lower pH or higher temperature. Furthermore, Tyr-169 and the S2-HB loop, or the X-loop, were different in the starting structures but converged to common conformations in the simulations for the Met-129, but not the Val-129, protein. α-Strands and β-strands formed in the initially unstructured N-terminus. α-Strand propensity in the N-terminus was different between the Met-129 and Val129 proteins, but β-strand propensity was similar. This study reveals detailed structural and dynamic properties of huPrP, providing insight into the mechanism of the conversion of PrP(C) to PrP(Sc).

Hydrogen-exchange mass spectrometry for the study of intrinsic disorder in proteins

Amide hydrogen/deuterium exchange detected by mass spectrometry (HXMS) is seeing wider use for the identification of intrinsically disordered parts of proteins. In this review, we discuss examples of how discovery of intrinsically disordered regions and their removal can aid in structure determination, biopharmaceutical quality control, the characterization of how post-translational modifications affect weak structuring of disordered regions, the study of coupled folding and binding, and the characterization of amyloid formation. This article is part of a Special Issue entitled: Mass spectrometry in structural biology.

Mouse prion protein polymorphism Phe-108/Val-189 affects the kinetics of fibril formation and the response to seeding: evidence for a two-step nucleation polymerization mechanism

Prion diseases are fatal neurodegenerative disorders associated with the polymerization of the cellular form of prion protein (PrP(C)) into an amyloidogenic β-sheet infectious form (PrP(Sc)). The sequence of host PrP is the major determinant of host prion disease susceptibility. In mice, the presence of allele a (Prnp(a), encoding the polymorphism Leu-108/Thr-189) or b (Prnp(b), Phe-108/Val-189) is associated with short or long incubation times, respectively, following infection with PrP(Sc). The molecular bases linking PrP sequence, infection susceptibility, and convertibility of PrP(C) into PrP(Sc) remain unclear. Here we show that recombinant PrP(a) and PrP(b) aggregate and respond to seeding differently in vitro. Our kinetic studies reveal differences during the nucleation phase of the aggregation process, where PrP(b) exhibits a longer lag phase that cannot be completely eliminated by seeding the reaction with preformed fibrils. Additionally, PrP(b) is more prone to propagate features of the seeds, as demonstrated by conformational stability and electron microscopy studies of the formed fibrils. We propose a model of polymerization to explain how the polymorphisms at positions 108 and 189 produce the phenotypes seen in vivo. This model also provides insight into phenomena such as species barrier and prion strain generation, two phenomena also influenced by the primary structure of PrP.

Helices 2 and 3 are the initiation sites in the PrP(C) → PrP(SC) transition

It is established that prion protein is the sole causative agent in a number of diseases in humans and animals. However, the nature of conformational changes that the normal cellular form, PrP(C), undergoes in its conversion to a self-replicating state is still not fully understood. The ordered C-terminus of PrP(C) proteins has three helices (H1-H3). Here, we use statistical coupling analysis (SCA) to infer covariations at various locations using a family of evolutionarily related sequences and the response of mouse and human PrP(C)s to mechanical force to decipher the initiation sites for the transition from PrP(C) to an aggregation-prone PrP* state. Sequence-based SCA predicts that the clustered residues in nonmammals are localized in the stable core (near H1) of PrP(C), whereas in mammalian PrP(C), they are localized in frustrated helices H2 and H3 where most of the pathogenic mutations are found. Force-extension curves and free energy profiles as a function of extension of mouse and human PrP(C) in the absence of a disulfide (SS) bond between residues Cys179 and Cys214, generated by applying mechanical force to the ends of the molecule, show a sequence of unfolding events starting first with rupture of H2 and H3. This is followed by disruption of structure in two strands. Helix H1, stabilized by three salt bridges, resists substantial force before unfolding. Force extension profiles and the dynamics of rupture of tertiary contacts also show that even in the presence of an SS bond the instabilities in most of H3 and parts of H2 still determine the propensity to form the PrP* state. In mouse PrP(C) with an SS bond, there are ∼10 residues that retain their order even at high forces. Both SCA and single-molecule force simulations show that in the conversion from PrP(C) to PrP(SC) major conformational changes occur (at least initially) in H2 and H3, which because of their sequence compositions are frustrated in the helical state. Implications of our findings for the structural model for the scrapie form of PrP(C) are discussed.

Massive conformation change in the prion protein: Using dual-basin structure-based models to find misfolding pathways

We employ all-atom structure-based models with a force field with multiple energetic basins for the C-terminal (residues 166-226) of the mammalian prion protein. One basin represents the known alpha-helical (αH) structure while the other represents the same residues in a left-handed beta-helical (LHBH) conformation. The LHBH structure has been proposed to help describe one class of in vitro grown fibrils, as well as possibly self-templating the conversion of normal cellular prion protein to the infectious form. Yet, it is unclear how the protein may make this global rearrangement. Our results demonstrate that the conformation changes are not strongly limited by large-scale geometry modification and that there may exist an overall preference for the LHBH conformation. Furthermore, our model presents novel intermediate trapping conformations with twisted LHBH structure.

Strand swapping regulates the iron-sulfur cluster in the diabetes drug target mitoNEET

MitoNEET is a recently identified diabetes drug target that coordinates a transferable 2Fe-2S cluster, and additionally contains an unusual strand swap. In this manuscript, we use a dual basin structure-based model to predict and characterize the folding and functionality of strand swapping in mitoNEET. We demonstrate that a strand unswapped conformation is kinetically accessible and that multiple levels of control are employed to regulate the conformational dynamics of the system. Environmental factors such as temperature can shift route preference toward the unswapped pathway. Additionally we see that a region recently identified as contributing to frustration in folding acts as a regulatory hinge loop that modulates conformational balance. Interestingly, strand unswapping transfers strain specifically to cluster-coordinating residues, opening the cluster-coordinating pocket. Strengthening contacts within the cluster-coordinating pocket opens a new pathway between the swapped and unswapped conformation that utilizes cracking to bypass the unfolded basin. These results suggest that local control within distinct regions affect motions important in regulating mitoNEET's 2Fe-2S clusters.

Role of water in protein aggregation and amyloid polymorphism

A variety of neurodegenerative diseases are associated with amyloid plaques, which begin as soluble protein oligomers but develop into amyloid fibrils. Our incomplete understanding of this process underscores the need to decipher the principles governing protein aggregation. Mechanisms of in vivo amyloid formation involve a number of coconspirators and complex interactions with membranes. Nevertheless, understanding the biophysical basis of simpler in vitro amyloid formation is considered important for discovering ligands that preferentially bind regions harboring amyloidogenic tendencies. The determination of the fibril structure of many peptides has set the stage for probing the dynamics of oligomer formation and amyloid growth through computer simulations. Most experimental and simulation studies, however, have been interpreted largely from the perspective of proteins: the role of solvent has been relatively overlooked in oligomer formation and assembly to protofilaments and amyloid fibrils. In this Account, we provide a perspective on how interactions with water affect folding landscapes of amyloid beta (Aβ) monomers, oligomer formation in the Aβ16-22 fragment, and protofilament formation in a peptide from yeast prion Sup35. Explicit molecular dynamics simulations illustrate how water controls the self-assembly of higher order structures, providing a structural basis for understanding the kinetics of oligomer and fibril growth. Simulations show that monomers of Aβ peptides sample a number of compact conformations. The formation of aggregation-prone structures (N*) with a salt bridge, strikingly similar to the structure in the fibril, requires overcoming a high desolvation barrier. In general, sequences for which N* structures are not significantly populated are unlikely to aggregate. Oligomers and fibrils generally form in two steps. First, water is expelled from the region between peptides rich in hydrophobic residues (for example, Aβ16-22), resulting in disordered oligomers. Then the peptides align along a preferred axis to form ordered structures with anti-parallel β-strand arrangement. The rate-limiting step in the ordered assembly is the rearrangement of the peptides within a confining volume. The mechanism of protofilament formation in a polar peptide fragment from the yeast prion, in which the two sheets are packed against each other and create a dry interface, illustrates that water dramatically slows self-assembly. As the sheets approach each other, two perfectly ordered one-dimensional water wires form. They are stabilized by hydrogen bonds to the amide groups of the polar side chains, resulting in the formation of long-lived metastable structures. Release of trapped water from the pore creates a helically twisted protofilament with a dry interface. Similarly, the driving force for addition of a solvated monomer to a preformed fibril is water release; the entropy gain and favorable interpeptide hydrogen bond formation compensate for entropy loss in the peptides. We conclude by offering evidence that a two-step model, similar to that postulated for protein crystallization, must also hold for higher order amyloid structure formation starting from N*. Distinct water-laden polymorphic structures result from multiple N* structures. Water plays multifarious roles in all of these protein aggregations. In predominantly hydrophobic sequences, water accelerates fibril formation. In contrast, water-stabilized metastable intermediates dramatically slow fibril growth rates in hydrophilic sequences.

Energy landscape of the prion protein helix 1 probed by metadynamics and NMR

The characterization of the structural dynamics of proteins, including those that present a substantial degree of disorder, is currently a major scientific challenge. These dynamics are biologically relevant and govern the majority of functional and pathological processes. We exploited a combination of enhanced molecular simulations of metadynamics and NMR measurements to study heterogeneous states of proteins and peptides. In this way, we determined the structural ensemble and free-energy landscape of the highly dynamic helix 1 of the prion protein (PrP-H1), whose misfolding and aggregation are intimately connected to a group of neurodegenerative disorders known as transmissible spongiform encephalopathies. Our combined approach allowed us to dissect the factors that govern the conformational states of PrP-H1 in solution, and the implications of these factors for prion protein misfolding and aggregation. The results underline the importance of adopting novel integrated approaches that take advantage of experiments and theory to achieve a comprehensive characterization of the structure and dynamics of biological macromolecules.

Discordant and chameleon sequences: their distribution and implications for amyloidogenicity

Identification of ambiguous encoding in protein secondary structure is paramount to develop an understanding of key protein segments underlying amyloid diseases. We investigate two types of structurally ambivalent peptides, which were hypothesized in the literature as indicators of amyloidogenic proteins: discordant α-helices and chameleon sequences. Chameleon sequences are peptides discovered experimentally in different secondary-structure types. Discordant α-helices are α-helical stretches with strong β-strand propensity or prediction. To assess the distribution of these features in known protein structures, and their potential role in amyloidogenesis, we analyzed the occurrence of discordant α-helices and chameleon sequences in nonredundant sets of protein domains (n = 4263) and amyloidogenic proteins extracted from the literature (n = 77). Discordant α-helices were identified if discordance was observed between known secondary structures and secondary-structure predictions from the GOR-IV and PSIPRED algorithms. Chameleon sequences were extracted by searching for identical sequence words in α-helices and β-strands. We defined frustrated chameleons and very frustrated chameleons based on varying degrees of total β propensity ≥α propensity. To our knowledge, this is the first study to discern statistical relationships between discordance, chameleons, and amyloidogenicity. We observed varying enrichment levels for some categories of discordant and chameleon sequences in amyloidogenic sequences. Chameleon sequences are also significantly enriched in proteins that have discordant helices, indicating a clear link between both phenomena. We identified the first set of discordant-chameleonic protein segments we predict may be involved in amyloidosis. We present a detailed analysis of discordant and chameleons segments in the family of one of the amyloidogenic proteins, the Prion Protein.

NMR structure and CD titration with metal cations of human prion alpha2-helix-related peptides

The 173–195 segment corresponding to the helix 2 of the C-globular prion protein domain could be one of several “spots” of intrinsic conformational flexibility. In fact, it possesses chameleon conformational behaviour and gathers several disease-associated point mutations. We have performed spectroscopic studies on the wild-type fragment 173–195 and on its D178N mutant dissolved in trifluoroethanol to mimic the in vivo system, both in the presence and in the absence of metal cations. NMR data showed that the structure of the D178N mutant is characterized by two short helices separated by a kink, whereas the wild-type peptide is fully helical. Both peptides retained these structural organizations, as monitored by CD, in the presence of metal cations. NMR spectra were however not in favour of the formation of definite ion-peptide complexes. This agrees with previous evidence that other regions of the prion protein are likely the natural target of metal cation binding.

Comparative analysis of essential collective dynamics and NMR-derived flexibility profiles in evolutionarily diverse prion proteins

Collective motions on ns-μs time scales are known to have a major impact on protein folding, stability, binding and enzymatic efficiency. It is also believed that these motions may have an important role in the early stages of prion protein misfolding and prion disease. In an effort to accurately characterize these motions and their potential influence on the misfolding and prion disease transmissibility we have conducted a combined analysis of molecular dynamic simulations and NMR-derived flexibility measurements over a diverse range of prion proteins. Using a recently developed numerical formalism, we have analyzed the essential collective dynamics (ECD) for prion proteins from 8 different species including human, cow, elk, cat, hamster, chicken, turtle and frog. We also compared the numerical results with flexibility profiles generated by the random coil index (RCI) from NMR chemical shifts. Prion protein backbone flexibility derived from experimental NMR data and from theoretical computations show strong agreement with each other, demonstrating that it is possible to predict the observed RCI profiles employing the numerical ECD formalism. Interestingly, flexibility differences in the loop between second beta strand (S2) and the second alpha helix (HB) appear to distinguish prion proteins from species that are susceptible to prion disease and those that are resistant. Our results show that the different levels of flexibility in the S2-HB loop in various species are predictable via the ECD method, indicating that ECD may be used to identify disease resistant variants of prion proteins, as well as the influence of prion proteins mutations on disease susceptibility or misfolding propensity.

Influence of pH on the human prion protein: insights into the early steps of misfolding

Transmissible spongiform encephalopathies, or prion diseases, are caused by misfolding and aggregation of the prion protein PrP. Conversion from the normal cellular form (PrP(C)) or recombinant PrP (recPrP) to a misfolded form is pH-sensitive, in that misfolding and aggregation occur more readily at lower pH. To gain more insight into the influence of pH on the dynamics of PrP and its potential to misfold, we performed extensive molecular-dynamics simulations of the recombinant PrP protein (residues 90-230) in water at three different pH regimes: neutral (or cytoplasmic) pH (∼7.4), middle (or endosomal) pH (∼5), and low pH (<4). We present five different simulations of 50 ns each for each pH regime, amounting to a total of 750 ns of simulation time. A detailed analysis and comparison with experiment validate the simulations and lead to new insights into the mechanism of pH-induced misfolding. The mobility of the globular domain increases with decreasing pH, through displacement of the first helix and instability of the hydrophobic core. At middle pH, conversion to a misfolded (PrP(Sc)-like) conformation is observed. The observed changes in conformation and stability are consistent with experimental data and thus provide a molecular basis for the initial steps in the misfolding process.