Diverse effects on the native β-sheet of the human prion protein due to disease-associated mutations

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.

Unraveling the Histidine Tautomerism Effect on the Initial Stages of Prion Misfolding: New Insights from a Computational Perspective

The aggregation and structural conversion of normal prion peptide (PrP^C) into the pathogenic scrapie form (PrP^Sc), which can act as a seed to enhance prion amyloid fiber formation, is believed to be a crucial event in prionopathies. Previous research suggests that the prion monomer may play an important role in oligomer generation during disease pathogenesis. In the present study, extensive replica-exchange molecular dynamics (REMD) simulations were conducted to explore the conformational characteristics of the huPrP (125-160) monomer under the histidine tautomerism effect. Investigating the structural characteristics and fibrilization process is challenging because two histidine tautomers [N_ε2-H (ε) and N_δ1-H (δ)] can occur in the open neutral state. Molecular dynamics (MD) simulation outcomes have shown that the toxic εδ and δδ isomer (containing several and broader local minima) had the highest α-helix structures, with contents of 21.11% and 21.01%, respectively, and may have a strong influence on the organizational behavior of a monomeric prion. The amino acids aspartate 20 (D20)-asparagine 29 (N29) and isoleucine 15 (I15)-histidine 16 (H16), D20-arginine 27 (R27) as well as N29 formed α-helix with the highest probabilities in the δδ and εδ isomer, accordingly. On the basis of our findings, we propose the histidine tautomerization hypothesis as a new prion accumulation mechanism, which may exist to induce the formation of prion accumulates. Overall, our tautomerism hypothesis constitutes a promising perspective for enhancing understanding of prion disease pathobiology and may help in the design of a good inhibitor.

Mutation-Dependent Refolding of Prion Protein Unveils Amyloidogenic-Related Structural Ramifications: Insights from Molecular Dynamics Simulations

The main focus of prion structural biology studies is to understand the molecular basis of prion diseases caused by misfolding, and aggregation of the cellular prion protein PrP^C remains elusive. Several genetic mutations are linked with human prion diseases and driven by the conformational conversion of PrP^C to the toxic PrP^Sc. The main goal of this study is to gain a better insight into the molecular effect of disease-associated V210I mutation on this process by molecular dynamics simulations. This inherited mutation elicited copious structural changes in the β1-α1-β2 subdomain, including an unfolding of a helix α1 and the elongation of the β-sheet. These unusual structural changes likely appeared to detach the β1-α1-β2 subdomain from the α2-α3 core, an early misfolding event necessary for the conformational conversion of PrP^C to PrP^Sc. Ultimately, the unfolded α1 and its prior β1-α1 loop further engaged with unrestrained conformational dynamics and were widely considered as amyloidogenic-inducing traits. Furthermore, the resulting folding intermediate possesses a highly unstable β1-α1-β2 subdomain, thereby enhancing the aggregation of misfolded PrP^C through intermolecular interactions between frequently refolding regions. Briefly, these remarkable changes as seen in the mutant β1-α1-β2 subdomain are consistent with previous experimental results and thus provide a molecular basis of PrP^C misfolding associated with the conformational conversion of PrP^C to PrP^Sc.

Tumorigenic p53 mutants undergo common structural disruptions including conversion to α-sheet structure

The p53 protein is a commonly studied cancer target because of its role in tumor suppression. Unfortunately, it is susceptible to mutation-associated loss of function; approximately 50% of cancers are associated with mutations to p53, the majority of which are located in the central DNA-binding domain. Here, we report molecular dynamics simulations of wild-type (WT) p53 and 20 different mutants, including a stabilized pseudo-WT mutant. Our findings indicate that p53 mutants tend to exacerbate latent structural-disruption tendencies, or vulnerabilities, already present in the WT protein, suggesting that it may be possible to develop cancer therapies by targeting a relatively small set of structural-disruption motifs rather than a multitude of effects specific to each mutant. In addition, α-sheet secondary structure formed in almost all of the proteins. α-Sheet has been hypothesized and recently demonstrated to play a role in amyloidogenesis, and its presence in the reported p53 simulations coincides with the recent re-consideration of cancer as an amyloid disease.

Biophysical characterization of oligomerization and fibrillization of the G131V pathogenic mutant of human prion protein

The pathogenesis of fatal neurodegenerative prion diseases is closely associated with the conversion of α-helix-rich cellular prion protein into β-sheet-rich scrapie form. Pathogenic point mutations of prion proteins usually promote the conformational conversion and trigger inherited prion diseases. The G131V mutation of human prion protein (HuPrP) was identified to be involved in Gerstmann-Sträussler-Scheinker syndrome. Few studies have been carried out to address the pathogenesis of the G131V mutant. Here, we addressed the effects of the G131V mutation on oligomerization and fibrillization of the full-length HuPrP(23-231) and truncated HuPrP(91-231) proteins. The G131V mutation promotes the oligomerization but alleviates the fibrillization of HuPrP, implying that the oligomerization might play a crucial role in the pathogenic mechanisms of the G131V mutant. Moreover, the flexible N-terminal fragment in either the wild-type or the G131V mutant HuPrP increases the oligomerization tendencies but decreases the fibrillization tendencies. Furthermore, this mutation significantly alters the tertiary structure of human PrPC and might distinctly change the conformational conversion tendency. Interestingly, both guanidine hydrochloride denaturation and thermal denaturation experiments showed that the G131V mutation does not significantly change the thermodynamic stabilities of the HuPrP proteins. This work may be of benefit to a mechanistic understanding of the conformational conversion of prion proteins and also provide clues for the prevention and treatment of prion diseases.

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

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

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

Structural basis for the complete resistance of the human prion protein mutant G127V to prion disease

Prion diseases are caused by the propagation of misfolded cellular prion proteins (PrPs). A completely prion disease-resistant genotype, V127M129, has been identified in Papua New Guinea and verified in transgenic mice. To disclose the structural basis of the disease-resistant effect of the G127V mutant, we determined and compared the structural and dynamic features of the G127V-mutated human PrP (residues 91-231) and the wild-type PrP in solution. HuPrP(G127V) contains α1, α2 and α3 helices and a stretch-strand (SS) pattern comprising residues Tyr128-Gly131 (SS1) and Val161-Arg164 (SS2), with extending atomic distances between the SS1 and SS2 strands, and a structural rearrangement of the Tyr128 side chain due to steric hindrance of the larger hydrophobic side chain of Val127. The extended α1 helix gets closer to the α2 and α3 helices. NMR dynamics analysis revealed that Tyr128, Gly131 and Tyr163 underwent significant conformational exchanges. Molecular dynamics simulations suggest that HuPrP(G127V) prevents the formation of stable β-sheets and dimers. Unique structural and dynamic features potentially inhibit the conformational conversion of the G127V mutant. This work is beneficial for understanding the molecular mechanisms underlying the complete resistance of the G127V mutant to prion disease and for developing new therapeutics for prion disease.

Molecular Dynamics Simulation Study on the Binding and Stabilization Mechanism of Antiprion Compounds to the "Hot Spot" Region of PrPC

Structural transitions in the prion protein from the cellular form, PrP^C, into the pathological isoform, PrP^Sc, are regarded as the main cause of the transmissible spongiform encephalopathies, also known as prion diseases. Hence, discovering and designing effective antiprion drugs that can inhibit PrP^C to PrP^Sc conversion is regarded as a promising way to cure prion disease. Among several strategies to inhibit PrP^C to PrP^Sc conversion, stabilizing the native PrP^C via specific binding is believed to be one of the valuable approaches and many antiprion compounds have been reported based on this strategy. However, the detailed mechanism to stabilize the native PrP^C is still unknown. As such, to unravel the stabilizing mechanism of these compounds to PrP^C is valuable for the further design and discovery of antiprion compounds. In this study, by molecular dynamics simulation method, we investigated the stabilizing mechanism of several antiprion compounds on PrP^C that were previously reported to have specific binding to the "hot spot" region of PrP^C. Our simulation results reveal that the stabilization mechanism of specific binding compounds can be summarized as (I) to stabilize both the flexible C-terminal of α2 and the hydrophobic core, such as BMD42-29 and GN8; (II) to stabilize the hydrophobic core, such as J1 and GJP49; (III) to stabilize the overall structure of PrP^C by high binding affinity, as NPR-056. In addition, as indicated by the H-bond analysis and decomposition analysis of binding free energy, the residues N159 and Q160 play an important role in the specific binding of the studied compounds and all these compounds interact with PrP^C in a similar way with the key interacting residues L130 in the β1 strand, P158, N159, Q160, etc. in the α1-β2 loop, and H187, T190, T191, etc. in the α2 C-terminus although the compounds have large structural difference. As a whole, our obtained results can provide some insights into the specific binding mechanism of main antiprion compounds to the "hot spot" region of PrP^C at the molecular level and also provide guidance for effective antiprion drug design in the future.

Simulations of membrane-bound diglycosylated human prion protein reveal potential protective mechanisms against misfolding

Prion diseases are associated with the misfolding of the prion protein (PrP) from its normal cellular form (PrP^C ) to its infectious scrapie form (PrP^Sc ). Post-translational modifications in PrP in vivo can play an important role in modulating the process of misfolding. To gain more insight into the effects of post-translational modifications in PrP structure and dynamics and to test the hypothesis that such modifications can interact with the protein, we have performed molecular dynamics simulations of diglycosylated human PrP^C bound to a lipid bilayer via a glycophosphatidylinositol anchor. Multiple simulations were performed at three different pH ranges to explore pH effects on structure and dynamics. In contrast to simulations of protein-only PrP^C , no large effects were observed upon lowering the pH of the system. The protein tilted toward the membrane surface in all of the simulations and the putative PrP^Sc oligomerization sites became inaccessible, thereby offering a possible protective mechanism against PrP^Sc -induced misfolding of PrP^C .

Exploration of the Main Sites for the Transformation of Normal Prion Protein (PrPC) into Pathogenic Prion Protein (PrPsc)

Introduction

The functions and mechanisms of prion proteins (PrP^C) are currently unknown, but most experts believe that deformed or pathogenic prion proteins (PrP^Sc) originate from PrP^C, and that there may be plural main sites for the conversion of normal PrP^C into PrP^Sc. In order to better understand the mechanism of PrP^C transformation to PrP^Sc, the most important step is to determine the replacement or substitution site.

Material and Methods

BALB/c mice were challenged with prion RML strain and from 90 days post-challenge (dpc) mice were sacrificed weekly until all of them had been at 160 dpc. The ultra-structure and pathological changes of the brain of experimental mice were observed and recorded by transmission electron microscopy.

Results

There were a large number of pathogen-like particles aggregated in the myelin sheath of the brain nerves, followed by delamination, hyperplasia, swelling, disintegration, phagocytic vacuolation, and other pathological lesions in the myelin sheath. The aggregated particles did not overflow from the myelin in unstained samples. The phenomenon of particle aggregation persisted all through the disease course, and was the earliest observed pathological change.

Conclusion

It was deduced that the myelin sheath and lipid rafts in brain nerves, including axons and dendrites, were the main sites for the conversion of PrP^C to PrP^Sc, and the PrP^Sc should be formed directly by the conversion of protein conformation without the involvement of nucleic acids.

Perturbations in inter-domain associations may trigger the onset of pathogenic transformations in PrP(C): insights from atomistic simulations

Conversion of the predominantly α-helical cellular prion protein (PrP(C)) to the misfolded β-sheet enriched Scrapie form (PrP(Sc)) is a critical event in prion pathogenesis. However, the conformational triggers that lead to the isoform conversion (PrP(C) to PrP(Sc)) remain obscure, and conjectures about the role of unusually hydrophilic, short helix H1 of the C-terminal globular domain in the transition are varied. Helix H1 is anchored to helix H3 via a few stabilizing polar interactions. We have employed fully atomistic molecular dynamics simulations to study the effects triggered by a minor perturbation in the network of these non-bonded interactions in PrP(C). The elimination of just one of the key H1-H3 hydrogen bonds led to a cascade of conformational changes that are consistent with those observed in partially unfolded intermediates of PrP(C), with pathogenic mutations and in low pH environments. Our analyses reveal that the perturbation results in the enhanced conformational flexibility of the protein. The resultant enhancement in the dynamics leads to overall increased solvent exposure of the hydrophobic core residues and concomitant disruption of the H1-H3 inter-domain salt bridge network. This study lends credence to the hypothesis that perturbing the cooperativity of the stabilizing interactions in the PrP(C) globular domain can critically affect its dynamics and may lead to structural transitions of pathological relevance.

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.

Free energy and hidden barriers of the β-sheet structure of prion protein

On-the-fly free-energy parametrization is a new collective variable biasing approach akin to metadynamics with one important distinction: rather than acquiring an accelerated distribution via a history-dependent bias potential, sampling on this distribution is achieved from the beginning of the simulation using temperature-accelerated molecular dynamics. In the present work, we compare the performance of both approaches to compute the free-energy profile along a scalar collective variable measuring the H-bond registry of the β-sheet structure of the mouse Prion protein. Both methods agree on the location of the free-energy minimum, but free-energy profiles from well-tempered metadynamics are subject to a much higher degree of statistical noise due to hidden barriers. The sensitivity of metadynamics to hidden barriers is shown to be a consequence of the history dependence of the bias potential, and we detail the nature of these barriers for the prion β-sheet. In contrast, on-the-fly parametrization is much less sensitive to these barriers and thus displays improved convergence behavior relative to that of metadynamics. While hidden barriers are a frequent and central issue in free-energy methods, on-the-fly free-energy parametrization appears to be a robust and preferable method to confront this issue.

Different misfolding mechanisms converge on common conformational changes: human prion protein pathogenic mutants Y218N and E196K

Prion diseases are caused by misfolding and aggregation of the prion protein (PrP). Pathogenic mutations such as Y218N and E196K are known to cause Gerstmann-Sträussler-Scheinker syndrome and Creutzfeldt-Jakob disease, respectively. Here we describe molecular dynamics simulations of these mutant proteins to better characterize the detailed conformational effects of these sequence substitutions. Our results indicate that the mutations disrupt the wild-type native PrP(C) structure and cause misfolding. Y218N reduced hydrophobic packing around the X-loop (residues 165-171), and E196K abolished an important wild-type salt bridge. While differences in the mutation site led PrP mutants to misfold along different pathways, we observed multiple traits of misfolding that were common to both mutants. Common traits of misfolding included: 1) detachment of the short helix (HA) from the PrP core; 2) exposure of side chain F198; and 3) formation of a nonnative strand at the N-terminus. The effect of the E196K mutation directly abolished the wild-type salt bridge E196-R156, which further destabilized the F198 hydrophobic pocket and HA. The Y218N mutation propagated its effect by increasing the HB-HC interhelical angle, which in turn disrupted the packing around F198. Furthermore, a nonnative contact formed between E221 and S132 on the S1-HA loop, which offered a direct mechanism for disrupting the hydrophobic packing between the S1-HA loop and HC. While there were common misfolding features shared between Y218N and E196K, the differences in the orientation of HB and HC and the X-loop conformation might provide a structural basis for identifying different prion strains.

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).

Molecular dynamics simulations capture the misfolding of the bovine prion protein at acidic pH

Bovine spongiform encephalopathy (BSE), or mad cow disease, is a fatal neurodegenerative disease that is transmissible to humans and that is currently incurable. BSE is caused by the prion protein (PrP), which adopts two conformers; PrPC is the native innocuous form, which is α-helix rich; and PrPSc is the β-sheet rich misfolded form, which is infectious and forms neurotoxic species. Acidic pH induces the conversion of PrPC to PrPSc. We have performed molecular dynamics simulations of bovine PrP at various pH regimes. An acidic pH environment induced conformational changes that were not observed in neutral pH simulations. Putative misfolded structures, with nonnative β-strands formed in the flexible N-terminal domain, were found in acidic pH simulations. Two distinct pathways were observed for the formation of nonnative β-strands: at low pH, hydrophobic contacts with M129 nucleated the nonnative β-strand; at mid-pH, polar contacts involving Q168 and D178 facilitated the formation of a hairpin at the flexible N-terminus. These mid- and low pH simulations capture the process of nonnative β-strand formation, thereby improving our understanding of how PrPC misfolds into the β-sheet rich PrPSc and how pH factors into the process.

Effect of charged residues in the N-domain of Sup35 protein on prion [PSI+] stability and propagation

Recent studies have shown that Sup35p prion fibrils probably have a parallel in-register β-structure. However, the part(s) of the N-domain critical for fibril formation and maintenance of the [PSI(+)] phenotype remains unclear. Here we designed a set of five SUP35 mutant alleles (sup35(KK)) with lysine substitutions in each of five N-domain repeats, and investigated their effect on infectivity and ability of corresponding proteins to aggregate and coaggregate with wild type Sup35p in the [PSI(+)] strain. Alleles sup35-M1 (Y46K/Q47K) and sup35-M2 (Q61K/Q62K) led to prion loss, whereas sup35-M3 (Q70K/Q71K), sup35-M4 (Q80K/Q81K), and sup35-M5 (Q89K/Q90K) were able to maintain the [PSI(+)] prion. This suggests that the critical part of the parallel in-register β-structure for the studied [PSI(+)] prion variant lies in the first 63-69 residues. Our study also reveals an unexpected interplay between the wild type Sup35p and proteins expressed from the sup35(KK) alleles during prionization. Both Sup35-M1p and Sup35-M2p coaggregated with Sup35p, but only sup35-M2 led to prion loss in a dominant manner. We suggest that in the fibrils, Sup35p can bind to Sup35-M1p in the same conformation, whereas Sup35-M2p only allowed the Sup35p conformation that leads to the non-heritable fold. Mutations sup35-M4 and sup35-M5 influence the structure of the prion forming region to a lesser extent, and can lead to the formation of new prion variants.

Probing early misfolding events in prion protein mutants by NMR spectroscopy

The post-translational conversion of the ubiquitously expressed cellular form of the prion protein, PrPC, into its misfolded and pathogenic isoform, known as prion or PrPSc, plays a key role in prion diseases. These maladies are denoted transmissible spongiform encephalopathies (TSEs) and affect both humans and animals. A prerequisite for understanding TSEs is unraveling the molecular mechanism leading to the conversion process whereby most α-helical motifs are replaced by β-sheet secondary structures. Importantly, most point mutations linked to inherited prion diseases are clustered in the C-terminal domain region of PrPC and cause spontaneous conversion to PrPSc. Structural studies with PrP variants promise new clues regarding the proposed conversion mechanism and may help identify "hot spots" in PrPC involved in the pathogenic conversion. These investigations may also shed light on the early structural rearrangements occurring in some PrPC epitopes thought to be involved in modulating prion susceptibility. Here we present a detailed overview of our solution-state NMR studies on human prion protein carrying different pathological point mutations and the implications that such findings may have for the future of prion research.

Molecular dynamics studies on 3D structures of the hydrophobic region PrP(109-136)

Prion diseases, traditionally referred to as transmissible spongiform encephalopathies, are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of mammalian species, manifesting as scrapie in sheep, bovine spongiform encephalopathy (or 'mad-cow' disease) in cattle, and Creutzfeldt-Jakob disease, Gerstmann-Strussler-Scheinker syndrome, fatal familial insomnia (FFI), and Kulu in humans, etc. These neurodegenerative diseases are caused by the conversion from a soluble normal cellular prion protein (PrP(C)) into insoluble abnormally folded infectious prions (PrP(Sc)). The hydrophobic region PrP(109-136) controls the formation of diseased prions: the normal PrP(113-120) AGAAAAGA palindrome is an inhibitor/blocker of prion diseases and the highly conserved glycine-xxx-glycine motif PrP(119-131) can inhibit the formation of infectious prion proteins in cells. This article gives detailed reviews on the PrP(109-136) region and presents the studies of its three-dimensional structures and structural dynamics.

Systematic investigation of predicted effect of nonsynonymous SNPs in human prion protein gene: a molecular modeling and molecular dynamics study

Nonsynonymous mutations in the human prion protein (HuPrP) gene contribute to the conversion of HuPrP(C) to HuPrP(Sc) and amyloid formation which in turn leads to prion diseases such as familial Creutzfeldt-Jakob disease and Gerstmann-Straussler-Scheinker disease. In order to better understand and predict the role of HuPrP mutations, we developed the following procedure: first, we consulted the Human Genome Variation database and dbSNP databases, and we reviewed literature for the retrieval of aggregation-related nsSNPs of the HuPrP gene. Next, we used three different methods - Polymorphism Phenotyping (PolyPhen), PANTHER, and Auto-Mute - to predict the effect of nsSNPs on the phenotype. We compared the predictions against experimentally reported effects of these nsSNPs to evaluate the accuracy of the three methods: PolyPhen predicted 17 out of 22 nsSNPs as "probably damaging" or "possibly damaging"; PANTHER predicted 8 out of 22 nsSNPs as "Deleterious"; and Auto-Mute predicted 9 out of 20 nsSNPs as "Disease". Finally, structural analyses of the native protein against mutated models were investigated using molecular modeling and molecular dynamics (MD) simulation methods. In addition to comparing predictor methods, our results show the applicability of our procedure for the prediction of damaging nsSNPs. Our study also elucidates the obvious relationship between predicted values of aggregation-related nsSNPs in HuPrP gene and molecular modeling and MD simulations results. In conclusion, this procedure would enable researchers to select outstanding candidates for extensive MD simulations in order to decipher more details of HuPrP aggregation. An animated interactive 3D complement (I3DC) is available in Proteopedia at http://proteopedia.org/w/Journal:JBSD:34.

Reversibility of prion misfolding: insights from constant-pH molecular dynamics simulations

The prion protein (PrP) is the cause of a group of diseases known as transmissible spongiform encephalopathies (TSEs). Creutzfeldt-Jakob disease and bovine spongiform encephalopathy are examples of TSEs. Although the normal form of PrP (PrP(C)) is monomeric and soluble, it can misfold into a pathogenic form (PrP(Sc)) that has a high content of β-structure and can aggregate forming amyloid fibrils. The mechanism of conversion of PrP(C) into PrP(Sc) is not known but different triggers have been proposed. It can be catalyzed by a PrP(Sc) sample, or it can be induced by an external factor, such as low pH. The pH effect on the structure of PrP was recently studied by computational methods [Campos et al. J. Phys. Chem. B 2010, 114, 12692-12700], and an evident trend of loss of helical structure was observed with pH decrease, together with a gain of β-structures. In particular, one simulation at pH 2 showed an evident misfolding transition. The main goal of the present work was to study the effects of a change in pH to 7 in several transient conformations of this simulation, in order to draw some conclusions about the reversibility of PrP misfolding. Although the most significant effect caused by the change of pH to 7 was a global stabilization of the protein structure, we could also observe that some conformational transitions induced by pH 2 were reversible in many of our simulations, namely those started from the early moments of the misfolding transition. This observation is in good agreement with experiments showing that, even at pH as low as 1.7, it is possible to revert the misfolding process [Bjorndahl et al. Biochemistry 2011, 50, 1162-1173].

Exploring the essential collective dynamics of interacting proteins: application to prion protein dimers

Essential collective dynamics is a promising and robust approach for analysing the slow motions of macromolecules from short molecular dynamics trajectories. In this study, an extension of the method to treat a collection of interacting protein molecules is presented. The extension is applied to investigate the effects of dimerization on the collective dynamics of ovine prion protein molecules in two different arrangements. Examination of the structural plasticity shows that aggregation has a restricting effect on the local mobility of the prion protein molecules in the interfacial regions. Domain motions of the two dimeric ovine prion protein conformations are distinctly different and can be related to interatomic correlations at the interface. Correlated motions are among the slow collective modes extensively analysed by considering both main-chain and side-chain atoms. Correlation maps reveal the existence of a vast network of dynamically correlated side groups, which extends beyond individual subunits via interfacial interconnections. The network is formed by a core of hydrophobic side chains surrounded by hydrophilic groups at the periphery. The relevance of these findings are discussed in the context of mutations associated with prion diseases. The binding free energy of the dimeric conformations is evaluated to probe their thermodynamic stability. The descriptions afforded by the analysis of the essential collective dynamics of the prion dimers are consistent with their binding free energies. The agreement validates the extension of the methodology and provides a means of interpreting the collective dynamics in terms of the thermodynamic stability of ovine prion proteins.

Dynamic diagnosis of familial prion diseases supports the β2-α2 loop as a universal interference target

BACKGROUND

Mutations in the cellular prion protein associated to familial prion disorders severely increase the likelihood of its misfolding into pathogenic conformers. Despite their postulation as incompatible elements with the native fold, these mutations rarely modify the native state structure. However they variably have impact on the thermodynamic stability and metabolism of PrP(C) and on the properties of PrP(Sc) aggregates. To investigate whether the pathogenic mutations affect the dynamic properties of the HuPrP(125-229) α-fold and find possible common patterns of effects that could help in prophylaxis we performed a dynamic diagnosis of ten point substitutions.

METHODOLOGY/PRINCIPAL FINDINGS

Using all-atom molecular dynamics simulations and novel analytical tools we have explored the effect of D178N, V180I, T183A, T188K, E196K, F198S, E200K, R208H, V210I and E211Q mutations on the dynamics of HuPrP(125-228) α-fold. We have found that while preserving the native state, all mutations produce dynamic changes which perturb the coordination of the α2-α3 hairpin to the rest of the molecule and cause the reorganization of the patches for intermolecular recognition, as the disappearance of those for conversion inhibitors and the emergence of an interaction site at the β2-α2 loop region.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that pathogenic mutations share a common pattern of dynamical alterations that converge to the conversion of the β2-α2 loop into an interacting region that can be used as target for interference treatments in genetic diseases.

Molecular dynamics as an approach to study prion protein misfolding and the effect of pathogenic mutations

Computer simulation of protein dynamics offers unique high-resolution information that complements experiment. Using experimentally derived structures of the natively folded prion protein (PrP), physically realistic dynamics and conformational changes can be simulated, including the initial steps of misfolding. By introducing mutations in silico, the effect of pathogenic mutations on PrP conformation and dynamics can be assessed. Here, we briefly introduce molecular dynamics methods and review the application of molecular dynamics simulations to obtain insight into various aspects of the PrP, including the mechanism of misfolding, the response to changes in the environment, and the influence of disease-related mutations.