Influence of pH on NMR structure and stability of the human prion protein globular domain

Experiment and molecular dynamics simulations reveal proanthocyanidin B2 and B3 can inhibit prion aggregation by different mechanisms

Prion diseases are a group of fatal neurodegenerative diseases caused by the misfolding and aggregation of prion protein (PrP), and the inhibition of PrP aggregation is one of the most effective therapeutic strategies. Proanthocyanidin B2 (PB2) and B3 (PB3), the effective natural antioxidants have been evaluated for the inhibition of amyloid-related protein aggregation. Since PrP has similar aggregation mechanism with other amyloid-related proteins, will PB2 and PB3 affect the aggregation of PrP? In this paper, experimental and molecular dynamics (MD) simulation methods were combined to investigate the influence of PB2 and PB3 on PrP aggregation. Thioflavin T assays showed PB2 and PB3 could inhibit PrP aggregation in a concentrate-dependent manner in vitro. To understand the underlying mechanism, we performed 400 ns all-atom MD simulations. The results suggested PB2 could stabilize the α2 C-terminus and the hydrophobic core of protein by stabilizing two important salt bridges R156-E196 and R156-D202, and consequently made global structure of protein more stable. Surprisingly, PB3 could not stabilize PrP, which may inhibit PrP aggregation through a different mechanism. Since dimerization is the first step of aggregation, will PB3 inhibit PrP aggregation by inhibiting the dimerization? To verify our assumption, we then explored the effect of PB3 on protein dimerization by performing 800 ns MD simulations. The results suggested PB3 could reduce the residue contacts and hydrogen bonds between two monomers, preventing dimerization process of PrP. The possible inhibition mechanism of PB2 and PB3 on PrP aggregation could provide useful information for drug development against prion diseases.Communicated by Ramaswamy H. Sarma.

Loss of residues 119 - 136, including the first β-strand of human prion protein, generates an aggregation-competent partially "open" form

In prion replication, the cellular form of prion protein (PrP^C) must undergo a full conformational transition to its disease-associated fibrillar form. Transmembrane forms of PrP have been implicated in this structural conversion. The cooperative unfolding of a structural core in PrP^C presents a substantial energy barrier to prion formation, with membrane insertion and detachment of parts of PrP presenting a plausible route to its reduction. Here, we examined the removal of residues 119 - 136 of PrP, a region which includes the first β-strand and a substantial portion of the conserved hydrophobic region of PrP, a region which associates with the ER membrane, on the structure, stability and self-association of the folded domain of PrP^C. We see an "open" native-like conformer with increased solvent exposure which fibrilises more readily than the native state. These data suggest a stepwise folding transition, which is initiated by the conformational switch to this "open" form of PrP^C.

Metadynamics simulations of ligands binding to protein surfaces: a novel tool for rational drug design

Structure-based drug design protocols may encounter difficulties to investigate poses when the biomolecular targets do not exhibit typical binding pockets. In this study, by providing two concrete examples from our labs, we suggest that the combination of metadynamics free energy methods (validated against affinity measurements), along with experimental structural information (by X-ray crystallography and NMR), can help to identify the poses of ligands on protein surfaces. The simulation workflow proposed here was implemented in a widely used code, namely GROMACS, and it could straightforwardly be applied to various drug-design campaigns targeting ligands' binding to protein surfaces.

Effects of the pathological E200K mutation on human prion protein: A computational screening and molecular dynamics approach

The human prion protein gene (PRNP) is mapped to the short arm of chromosome 20 (20pter-12). Prion disease is associated with mutations in the prion protein-encoding gene sequence. Earlier studies found that the mutation G127V in the PRNP increases protein stability. In contrast, the mutation E200K, which has the highest mutation rate in the prion protein, causes Creutzfeldt-Jakob disease (CJD) in humans and induces protein aggregation. We aimed to identify the structural mechanisms of E200k and G127V mutations causing CJD. We used a variety of bioinformatic algorithms, including SIFT, PolyPhen, I-Mutant, PhD-SNP, and SNP& GO, to predict the association of the E200K mutation with prion disease. MD simulation is performed, and graphs for root mean square deviation, root mean square fluctuation, radius of gyration, DSSP, principal component analysis, porcupine, and free energy landscape are generated to confirm and prove the stability of the wild-type and mutant protein structures. The protein is analyzed for aggregation, and the results indicate more fluctuations in the protein structure during the simulation owing to the E200K mutation; however, the G127V mutation makes the protein structure stable against aggregation during the simulation.

Investigation of the conformation of human prion protein in ethanol solution using molecular dynamics simulations

When the conformation of protein is changed from its natural state to a misfolded state, some diseases will happen like prion disease. Prion diseases are a set of deadly neurodegenerative diseases caused by prion protein misfolding and aggregation. Monohydric alcohols have a strong influence on the structure of protein. However, whether monohydric alcohols inhibit amyloid fibrosis remains uncertain. Here, to elucidate the effect of ethanol on the structural stability of human prion protein, molecular dynamics simulations were employed to analyze the conformational changes and dynamics characteristics of human prion proteins at different temperatures. The results show that the extension of β-sheet occurs more easily and the α-helix is more easily disrupted at high temperatures. We found that ethanol can destroy the hydrophobic interactions and make the hydrogen bonds stable, which protects the secondary structure of the protein, especially at 500 K.Communicated by Ramaswamy H. Sarma.

On the Dependence of Prion and Amyloid Structure on the Folding Environment

Currently available analyses of amyloid proteins reveal the necessity of the existence of radical structural changes in amyloid transformation processes. The analysis carried out in this paper based on the model called fuzzy oil drop (FOD) and its modified form (FOD-M) allows quantifying the role of the environment, particularly including the aquatic environment. The starting point and basis for the present presentation is the statement about the presence of two fundamentally different methods of organizing polypeptides into ordered conformations—globular proteins and amyloids. The present study shows the source of the differences between these two paths resulting from the specificity of the external force field coming from the environment, including the aquatic and hydrophobic one. The water environment expressed in the fuzzy oil drop model using the 3D Gauss function directs the folding process towards the construction of a micelle-like system with a hydrophobic core in the central part and the exposure of polarity on the surface. The hydrophobicity distribution of membrane proteins has the opposite characteristic: Exposure of hydrophobicity at the surface of the membrane protein with an often polar center (as in the case of ion channels) is expected. The structure of most proteins is influenced by a more or less modified force field generated by water through the appropriate presence of a non-polar (membrane-like) environment. The determination of the proportion of a factor different from polar water enables the assessment of the protein status by indicating factors favoring the structure it represents.

Epitope-specific anti-PrP antibody toxicity: a comparative in-silico study of human and mouse prion proteins

Despite having therapeutic potential, anti-PrP antibodies caused a major controversy due to their neurotoxic effects. For instance, treating mice with ICSM antibodies delayed prion disease onset, but both were found to be either toxic or innocuous to neurons by researchers following cross-linking PrP^C. In order to elucidate and understand the reasons that led to these contradictory outcomes, we conducted a comprehensive in silico study to assess the antibody-specific toxicity. Since most therapeutic anti-PrP antibodies were generated against human truncated recombinant PrP^91-231 or full-length mouse PrP^23-231, we reasoned that host specificity (human vs murine) of PrP^C might influence the nature of the specific epitopes recognized by these antibodies at the structural level possibly explaining the 'toxicity' discrepancies reported previously. Initially, molecular dynamics simulation and pro-motif analysis of full-length human (hu)PrP and mouse (mo)PrP 3D structure displayed conspicuous structural differences between huPrP and moPrP. We identified 10 huPrP and 6 moPrP linear B-cell epitopes from the prion protein 3D structure where 5 out of 10 huPrP and 3 out of 6 moPrP B-cell epitopes were predicted to be potentially toxic in immunoinformatics approaches. Herein, we demonstrate that some of the predicted potentially 'toxic' epitopes identified by the in silico analysis were similar to the epitopes recognized by the toxic antibodies such as ICSM18 (146-159), POM1 (138-147), D18 (133-157), ICSM35 (91-110), D13 (95-103) and POM3 (95-100). This in silico study reveals the role of host specificity of PrP^C in epitope-specific anti-PrP antibody toxicity.

Destabilization of polar interactions in the prion protein triggers misfolding and oligomerization

The prion protein (PrP) misfolds and oligomerizes at pH 4 in the presence of physiological salt concentrations. Low pH and salt cause structural perturbations in the monomeric prion protein that lead to misfolding and oligomerization. However, the changes in stability within different regions of the PrP prior to oligomerization are poorly understood. In this study, we have characterized the local stability in PrP at high resolution using amide temperature coefficients (T_C ) measured by nuclear magnetic resonance (NMR) spectroscopy. The local stability of PrP was investigated under native as well as oligomerizing conditions. We have also studied the rapidly oligomerizing PrP variant (Q216R) and the protective PrP variant (A6). We report that at low pH, salt destabilizes PrP at several polar residues, and the hydrogen bonds in helices α2 and α3 are weakened. In addition, salt changes the curvature of the α3 helix, which likely disrupts α2-α3 contacts and leads to oligomerization. These results are corroborated by the T_C values of rapidly oligomerizing Q216R-PrP. The poly-alanine substitution in A6-PrP stabilizes α2, which prevents oligomerization. Altogether, these results highlight the importance of native polar interactions in determining the stability of PrP and reveal the structural disruptions in PrP that lead to misfolding and oligomerization.

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.

Mechanism of misfolding of the human prion protein revealed by a pathological mutation

The misfolding and aggregation of the human prion protein (PrP) is associated with transmissible spongiform encephalopathies (TSEs). Intermediate conformations forming during the conversion of the cellular form of PrP into its pathological scrapie conformation are key drivers of the misfolding process. Here, we analyzed the properties of the C-terminal domain of the human PrP (huPrP) and its T183A variant, which is associated with familial forms of TSEs. We show that the mutation significantly enhances the aggregation propensity of huPrP, such as to uniquely induce amyloid formation under physiological conditions by the sole C-terminal domain of the protein. Using NMR spectroscopy, biophysics, and metadynamics simulations, we identified the structural characteristics of the misfolded intermediate promoting the aggregation of T183A huPrP and the nature of the interactions that prevent this species to be populated in the wild-type protein. In support of these conclusions, POM antibodies targeting the regions that promote PrP misfolding were shown to potently suppress the aggregation of this amyloidogenic mutant.

Structural details of amyloid β oligomers in complex with human prion protein as revealed by solid-state MAS NMR spectroscopy

Human PrP (huPrP) is a high-affinity receptor for oligomeric amyloid β (Aβ) protein aggregates. Binding of Aβ oligomers to membrane-anchored huPrP has been suggested to trigger neurotoxic cell signaling in Alzheimer’s disease, while an N-terminal soluble fragment of huPrP can sequester Aβ oligomers and reduce their toxicity. Synthetic oligomeric Aβ species are known to be heterogeneous, dynamic, and transient, rendering their structural investigation particularly challenging. Here, using huPrP to preserve Aβ oligomers by coprecipitating them into large heteroassemblies, we investigated the conformations of Aβ(1–42) oligomers and huPrP in the complex by solid-state MAS NMR spectroscopy. The disordered N-terminal region of huPrP becomes immobilized in the complex and therefore visible in dipolar spectra without adopting chemical shifts characteristic of a regular secondary structure. Most of the well-defined C-terminal part of huPrP is part of the rigid complex, and solid-state NMR spectra suggest a loss in regular secondary structure in the two C-terminal α-helices. For Aβ(1–42) oligomers in complex with huPrP, secondary chemical shifts reveal substantial β-strand content. Importantly, not all Aβ(1–42) molecules within the complex have identical conformations. Comparison with the chemical shifts of synthetic Aβ fibrils suggests that the Aβ oligomer preparation represents a heterogeneous mixture of β-strand-rich assemblies, of which some have the potential to evolve and elongate into different fibril polymorphs, reflecting a general propensity of Aβ to adopt variable β-strand-rich conformers. Taken together, our results reveal structural changes in huPrP upon binding to Aβ oligomers that suggest a role of the C terminus of huPrP in cell signaling. Trapping Aβ(1–42) oligomers by binding to huPrP has proved to be a useful tool for studying the structure of these highly heterogeneous β-strand-rich assemblies.

The Size and Stability of Infectious Prion Aggregates Fluctuate Dynamically during Cellular Uptake and Disaggregation

Prion diseases arise when PrP^Sc, an aggregated, infectious, and insoluble conformer of the normally soluble mammalian prion protein, PrP^C, catalyzes the conversion of PrP^C into more PrP^Sc, which then accumulates in the brain leading to disease. PrP^Sc is the primary, if not sole, component of the infectious prion. Despite the stability and protease insensitivity of PrP^Sc aggregates, they can be degraded after cellular uptake. However, how cells disassemble and degrade PrP^Sc is poorly understood. In this work, we analyzed how the protease sensitivity and size distribution of PrP^Sc aggregates from two different mouse-adapted prion strains, 22L, that can persistently infect cells and 87V, that cannot, changed during cellular uptake. We show that within the first 4 h following uptake large PrP^Sc aggregates from both prion strains become less resistant to digestion by proteinase K (PK) through a mechanism that is dependent upon the acidic environment of endocytic vesicles. We further show that during disassembly, PrP^Sc aggregates from both strains become more resistant to PK digestion through the apparent removal of protease-sensitive PrP^Sc, with PrP^Sc from the 87V strain disassembled more readily than PrP^Sc from the 22L strain. Taken together, our data demonstrate that the sizes and stabilities of PrP^Sc from different prion strains change during cellular uptake and degradation, thereby potentially impacting the ability of prions to infect cells.

Multimodal small-molecule screening for human prion protein binders

Prion disease is a rapidly progressive neurodegenerative disorder caused by misfolding and aggregation of the prion protein (PrP), and there are currently no therapeutic options. PrP ligands could theoretically antagonize prion formation by protecting the native protein from misfolding or by targeting it for degradation, but no validated small-molecule binders have been discovered to date. We deployed a variety of screening methods in an effort to discover binders of PrP, including 19F-observed and saturation transfer difference (STD) NMR spectroscopy, differential scanning fluorimetry (DSF), DNA-encoded library selection, and in silico screening. A single benzimidazole compound was confirmed in concentration-response, but affinity was very weak (Kd > 1 mm), and it could not be advanced further. The exceptionally low hit rate observed here suggests that PrP is a difficult target for small-molecule binders. Whereas orthogonal binder discovery methods could yield high-affinity compounds, non-small-molecule modalities may offer independent paths forward against prion disease.

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.

Combining molecular dynamics simulations and experimental analyses in protein misfolding

The fold of a protein determines its function and its misfolding can result in loss-of-function defects. In addition, for certain proteins their misfolding can lead to gain-of-function toxicities resulting in protein misfolding diseases such as Alzheimer's, Parkinson's, or the prion diseases. In all of these diseases one or more proteins misfold and aggregate into disease-specific assemblies, often in the form of fibrillar amyloid deposits. Most, if not all, protein misfolding diseases share a fundamental molecular mechanism that governs the misfolding and subsequent aggregation. A wide variety of experimental methods have contributed to our knowledge about misfolded protein aggregates, some of which are briefly described in this review. The misfolding mechanism itself is difficult to investigate, as the necessary timescale and resolution of the misfolding events often lie outside of the observable parameter space. Molecular dynamics simulations fill this gap by virtue of their intrinsic, molecular perspective and the step-by-step iterative process that forms the basis of the simulations. This review focuses on molecular dynamics simulations and how they combine with experimental analyses to provide detailed insights into protein misfolding and the ensuing diseases.

Structural insight into conformational change in prion protein by breakage of electrostatic network around H187 due to its protonation

A conformational change from normal prion protein(PrP^C) to abnormal prion protein(PrP^SC) induces fatal neurodegenerative diseases. Acidic pH is well-known factors involved in the conformational change. Because the protonation of H187 is strongly linked to the change in PrP stability, we examined the charged residues R156, E196, and D202 around H187. Interestingly, there have been reports on pathological mutants, such as H187R, E196A, and D202N. In this study, we focused on how an acidic pH and pathological mutants disrupt this electrostatic network and how this broken network destabilizes PrP structure. To do so, we performed a temperature-based replica-exchange molecular dynamics (T-REMD) simulation using a cumulative 252 μs simulation time. We measured the distance between amino acids comprising four salt bridges (R156–E196/D202 and H187–E196/D202). Our results showed that the spatial configuration of the electrostatic network was significantly altered by an acidic pH and mutations. The structural alteration in the electrostatic network increased the RMSF value around the first helix (H1). Thus, the structural stability of H1, which is anchored to the H2–H3 bundle, was decreased. It induces separation of R156 from the electrostatic network. Analysis of the anchoring energy also shows that two salt-bridges (R156-E196/D202) are critical for PrP stability.

Atomic insights into the effects of pathological mutants through the disruption of hydrophobic core in the prion protein

Destabilization of prion protein induces a conformational change from normal prion protein (PrPC) to abnormal prion protein (PrPSC). Hydrophobic interaction is the main driving force for protein folding, and critically affects the stability and solvability. To examine the importance of the hydrophobic core in the PrP, we chose six amino acids (V176, V180, T183, V210, I215, and Y218) that make up the hydrophobic core at the middle of the H2-H3 bundle. A few pathological mutants of these amino acids have been reported, such as V176G, V180I, T183A, V210I, I215V, and Y218N. We focused on how these pathologic mutations affect the hydrophobic core and thermostability of PrP. For this, we ran a temperature-based replica-exchange molecular dynamics (T-REMD) simulation, with a cumulative simulation time of 28 μs, for extensive ensemble sampling. From the T-REMD ensemble, we calculated the protein folding free energy difference between wild-type and mutant PrP using the thermodynamic integration (TI) method. Our results showed that pathological mutants V176G, T183A, I215V, and Y218N decrease the PrP stability. At the atomic level, we examined the change in pair-wise hydrophobic interactions from valine-valine to valine-isoleucine (and vice versa), which is induced by mutation V180I, V210I (I215V) at the 180th-210th (176th-215th) pair. Finally, we investigated the importance of the π-stacking between Y218 and F175.

Atomic insight into prion disorder: An intricate detail gained by 0.5 μs molecular dynamics simulation of preventive G127V and deleterious D178V mutation in prion protein

In this study we are looking into two contradicting mutations found in prion protein (PrP) viz G127V and D178V, that are reportedly protective and pathogenic, respectively. Despite significant advances in comprehension of the role of pathogenic mutations, the role of protective mutation in amyloid fold inhibition still lacks a substantial basis. To understand the structural basis of protective mutation, molecular dynamics simulation coupled with protein-protein docking and molecular mechanics/Poisson-Boltzmann surface area analysis was used to understand the instant structural variability brought about by these mutations alone and in combination on PrP and prion-prion complex. Atomic-scale investigations successfully revealed that the binding pattern of prion-prion varies differentially in protective and pathogenic mutations with secondary structure showing distinct contrasting patterns, which could supposedly be a critical factor for differential prion behavior in protective and pathogenic mutations. Considering the reported role of an amyloid fold in prion-prion binding, the contrasting pattern has given us a lead in comprehending the role of these mutations and has been used in this study to look for small molecules that can inhibit amyloid fold for prion-prion interaction in pathogenic mutant carrying PrP.

How graphene affects the misfolding of human prion protein: A combined experimental and molecular dynamics simulation study

As the broad application of graphene in the biomedical field, it is urgent and important to evaluate how the graphene affects the structure and function of the proteins in our body, especially the amyloid-related proteins. Prion protein, as a typical amyloid protein, it misfolding and aggregation will lead to serious prion diseases. To explore if graphene promotes or inhibits the formation of amyloid, here, we combined the experimental and molecular dynamics (MD) simulation methods to study the influence of graphene on the globular domain of prion protein (PrP_117-231). The results from fluorescence quenching and circular dichroism spectrum showed that the addition of graphene changed the secondary structure of prion protein largely, mainly reflecting in the reduced α-helix structure and the increased coil structure, indicating graphene may strengthen the misfolding inclination of prion. To further uncover the mechanism of conformational change of prion under the induction of graphene, the all-atoms MD simulations in explicit solvent were performed. Our simulations suggest that prion protein can be quickly and tightly adsorbed onto graphene together with the weak conformational rearrangement and may reorient when approaching the surface. The Van der Waals' force drive the adsorption process. In the induction of graphene, H1 and S2-H2 loop regions of prion become unstable and prion begins to misfold partially. Our work shows that graphene can induce the misfolding of prion protein and may cause the potential risk to biosystems.

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

Liquid and Hydrogel Phases of PrPC Linked to Conformation Shifts and Triggered by Alzheimer's Amyloid-β Oligomers

Protein phase separation by low-complexity, intrinsically disordered domains generates membraneless organelles and links to neurodegeneration. Cellular prion protein (PrP^C) contains such domains, causes spongiform degeneration, and is a receptor for Alzheimer's amyloid-β oligomers (Aβo). Here, we show that PrP^C separates as a liquid phase, in which α-helical Thr become unfolded. At the cell surface, PrP^C Lys residues interact with Aβo to create a hydrogel containing immobile Aβo and relatively mobile PrP^C. The Aβo/PrP hydrogel has a well-defined stoichiometry and dissociates with excess Aβo. NMR studies of hydrogel PrP^C reveal a distinct α-helical conformation for natively unfolded amino-terminal Gly and Ala residues. Aβo/PrP hydrogel traps signal-transducing mGluR5 on the plasma membrane. Recombinant PrP^C extracts endogenous Aβo from human Alzheimer's soluble brain lysates into hydrogel, and a PrP^C antagonist releases Aβo from endogenous brain hydrogel. Thus, coupled phase and conformational transitions of PrP^C are driven by Aβ species from Alzheimer's disease.

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.

Effects of pH and aggregation in the human prion conversion into scrapie form: a study using molecular dynamics with excited normal modes

There are two different prion conformations: (1) the cellular natural (PrP^C) and (2) the scrapie (PrP^Sc), an infectious form that tends to aggregate under specific conditions. PrP^C and PrP^Sc are widely different regarding secondary and tertiary structures. PrP^Sc contains more and longer β-strands compared to PrP^C. The lack of solved PrP^Sc structures precludes a proper understanding of the mechanisms related to the transition between cellular and scrapie forms, as well as the aggregation process. In order to investigate the conformational transition between PrP^C and PrP^Sc, we applied MDeNM (molecular dynamics with excited normal modes), an enhanced sampling simulation technique that has been recently developed to probe large structural changes. These simulations yielded new structural rearrangements of the cellular prion that would have been difficult to obtain with standard MD simulations. We observed an increase in β-sheet formation under low pH (≤ 4) and upon oligomerization, whose relevance was discussed on the basis of the energy landscape theory for protein folding. The characterization of intermediate structures corresponding to transition states allowed us to propose a conversion model from the cellular to the scrapie prion, which possibly ignites the fibril formation. This model can assist the design of new drugs to prevent neurological disorders related to the prion aggregation mechanism.

Prion Proteins Without the Glycophosphatidylinositol Anchor: Potential Biomarkers in Neurodegenerative Diseases

Prion protein (PrP) is a biomolecule that is involved in neuronal signaling, myelinization, and the development of neurodegenerative diseases. In the cell, PrP is shed by the ADAM10 protease. This process generates PrP molecules that lack glycophosphatidylinositol anchor, and these molecules incorporate into toxic aggregates and neutralize toxic oligomers. Due to this dual role, these molecules are important biomarkers for neurodegenerative diseases. In this review, we present shed PrP as a potential biomarker, with a focus on PrP226*, which may be the main biomarker for predicting neurodegenerative diseases in humans.

The function of the cellular prion protein in health and disease

The essential role of the cellular prion protein (PrP^C) in prion disorders such as Creutzfeldt-Jakob disease is well documented. Moreover, evidence is accumulating that PrP^C may act as a receptor for protein aggregates and transduce neurotoxic signals in more common neurodegenerative disorders, such as Alzheimer's disease. Although the pathological roles of PrP^C have been thoroughly characterized, a general consensus on its physiological function within the brain has not yet been established. Knockout studies in various organisms, ranging from zebrafish to mice, have implicated PrP^C in a diverse range of nervous system-related activities that include a key role in the maintenance of peripheral nerve myelination as well as a general ability to protect against neurotoxic stimuli. Thus, the function of PrP^C may be multifaceted, with different cell types taking advantage of unique aspects of its biology. Deciphering the cellular function(s) of PrP^C and the consequences of its absence is not simply an academic curiosity, since lowering PrP^C levels in the brain is predicted to be a powerful therapeutic strategy for the treatment of prion disease. In this review, we outline the various approaches that have been employed in an effort to uncover the physiological and pathological functions of PrP^C. While these studies have revealed important clues about the biology of the prion protein, the precise reason for PrP^C's existence remains enigmatic.

A 13C-detected 15N double-quantum NMR experiment to probe arginine side-chain guanidinium 15Nη chemical shifts

Arginine side-chains are often key for enzyme catalysis, protein-ligand and protein-protein interactions. The importance of arginine stems from the ability of the terminal guanidinium group to form many key interactions, such as hydrogen bonds and salt bridges, as well as its perpetual positive charge. We present here an arginine ¹³C^ζ-detected NMR experiment in which a double-quantum coherence involving the two ¹⁵N^η nuclei is evolved during the indirect chemical shift evolution period. As the precession frequency of the double-quantum coherence is insensitive to exchange of the two ¹⁵N^η; this new approach is shown to eliminate the previously deleterious line broadenings of ¹⁵N^η resonances caused by the partially restricted rotation about the C^ζ-N^ε bond. Consequently, sharp and well-resolved ¹⁵N^η resonances can be observed. The utility of the presented method is demonstrated on the L99A mutant of the 19 kDa protein T4 lysozyme, where the measurement of small chemical shift perturbations, such as one-bond deuterium isotope shifts, of the arginine amine ¹⁵N^η nuclei becomes possible using the double-quantum experiment.

Implications of peptide assemblies in amyloid diseases

Neurodegenerative disorders and type 2 diabetes are global epidemics compromising the quality of life of millions worldwide, with profound social and economic implications. Despite the significant differences in pathology - much of which are poorly understood - these diseases are commonly characterized by the presence of cross-β amyloid fibrils as well as the loss of neuronal or pancreatic β-cells. In this review, we document research progress on the molecular and mesoscopic self-assembly of amyloid-beta, alpha synuclein, human islet amyloid polypeptide and prions, the peptides and proteins associated with Alzheimer's, Parkinson's, type 2 diabetes and prion diseases. In addition, we discuss the toxicities of these amyloid proteins based on their self-assembly as well as their interactions with membranes, metal ions, small molecules and engineered nanoparticles. Through this presentation we show the remarkable similarities and differences in the structural transitions of the amyloid proteins through primary and secondary nucleation, the common evolution from disordered monomers to alpha-helices and then to β-sheets when the proteins encounter the cell membrane, and, the consensus (with a few exceptions) that off-pathway oligomers, rather than amyloid fibrils, are the toxic species regardless of the pathogenic protein sequence or physicochemical properties. In addition, we highlight the crucial role of molecular self-assembly in eliciting the biological and pathological consequences of the amyloid proteins within the context of their cellular environments and their spreading between cells and organs. Exploiting such structure-function-toxicity relationship may prove pivotal for the detection and mitigation of amyloid diseases.

Sequence-dependent dynamical instability of the human prion protein: a comparative simulation study

The present study aimed to explore the most probable regions of the human prion protein backbone for which the initial steps of conformational transitions as a result of intrinsic and extrinsic perturbing factors on the protein structure can be assigned. A total of 0.3-μs molecular dynamics simulations on several analog structures of the protein have been performed. To mimic the impact of the extrinsic and intrinsic destructive parameters on the dynamical characteristics of the protein, mild acidic conditions and R208H mutation have been simulated. The findings indicated that distribution of conformational flexibilities along the protein chain was almost independent of the induced perturbing factors, and was mostly centralized on certain distinct parts of the structure comprising residues 132-145 and 187-203. Analyses also revealed that the segment comprising residues 187-203 may be considered as a peptide sequence, possessing high potential to start the initial steps of conformational rearrangements due to the induced physicochemical alterations. Sequence alignment and molecular dynamics data also revealed that segment 178-203 prefers to accommodate in extended structures rather than α-helices. Region 178-203 may be considered as a peptide switch capable of initiating the conformational transitions due to the introduced modifications and perturbing parameters.

Solution NMR Spectroscopy in Target-Based Drug Discovery

Solution NMR spectroscopy is a powerful tool to study protein structures and dynamics under physiological conditions. This technique is particularly useful in target-based drug discovery projects as it provides protein-ligand binding information in solution. Accumulated studies have shown that NMR will play more and more important roles in multiple steps of the drug discovery process. In a fragment-based drug discovery process, ligand-observed and protein-observed NMR spectroscopy can be applied to screen fragments with low binding affinities. The screened fragments can be further optimized into drug-like molecules. In combination with other biophysical techniques, NMR will guide structure-based drug discovery. In this review, we describe the possible roles of NMR spectroscopy in drug discovery. We also illustrate the challenges encountered in the drug discovery process. We include several examples demonstrating the roles of NMR in target-based drug discoveries such as hit identification, ranking ligand binding affinities, and mapping the ligand binding site. We also speculate the possible roles of NMR in target engagement based on recent processes in in-cell NMR spectroscopy.

Polymorphism at 129 dictates metastable conformations of the human prion protein N-terminal β-sheet

We study the thermodynamic stability of the native state of the human prion protein using a new free-energy method, replica-exchange on-the-fly parameterization. This method is designed to overcome hidden-variable sampling limitations to yield nearly error-free free-energy profiles along a conformational coordinate. We confirm that all four (M129V, D178N) polymorphs have a ground-state conformation with three intact β-sheet hydrogen bonds. Additionally, they are observed to have distinct metastabilities determined by the side-chain at position 129. We rationalize these findings with reference to the prion "strain" hypothesis, which links the variety of transmissible spongiform encephalopathy phenotypes to conformationally distinct infectious prion forms and classifies distinct phenotypes of sporadic Creutzfeldt-Jakob disease based solely on the 129 polymorphism. Because such metastable structures are not easily observed in structural experiments, our approach could potentially provide new insights into the conformational origins of prion diseases and other pathologies arising from protein misfolding and aggregation.

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.

Insight into Early-Stage Unfolding of GPI-Anchored Human Prion Protein

Prion diseases are fatal neurodegenerative disorders, which are characterized by the accumulation of misfolded prion protein (PrPSc) converted from a normal host cellular prion protein (PrPC). Experimental studies suggest that PrPC is enriched with α-helical structure, whereas PrPSc contains a high proportion of β-sheet. In this study, we report the impact of N-glycosylation and the membrane on the secondary structure stability utilizing extensive microsecond molecular dynamics simulations. Our results reveal that the HB (residues 173 to 194) C-terminal fragment undergoes conformational changes and helix unfolding in the absence of membrane environments because of the competition between protein backbone intramolecular and protein-water intermolecular hydrogen bonds as well as its intrinsic instability originated from the amino acid sequence. This initiation of the unfolding process of PrPC leads to a subsequent increase in the length of the HB-HC loop (residues 195 to 199) that may trigger larger rigid body motions or further unfolding around this region. Continuous interactions between prion protein and the membrane not only constrain the protein conformation but also decrease the solvent accessibility of the backbone atoms, thereby stabilizing the secondary structure, which is enhanced by N-glycosylation via additional interactions between the N-glycans and the membrane surface.

Conformational dynamics of a hydrophobic prion fragment (113-127) in different pH and osmolyte solutions

Prion diseases are characterized by a conformational change in prion protein from its native state into beta-sheet rich aggregates that are neurotoxic. The central domain that contain a highly conserved hydrophobic region of the protein play an important role in the toxicity. The conformation of the proteins is largely influenced by various solvent environments. Here we report results of study of hydrophobic prion fragment peptide PrP(113-127) under different pH and osmolytes solution conditions. The secondary structure and the folding of PrP(113-127) was determined using circular dichroism and fluorescence spectroscopic methods. The results indicate that PrP(113-127) adopts a random coil conformation in aqueous buffer at neutral pH and that converted into beta sheet on aging. Even though the initial random coil conformation was similar in different pH conditions, the acidic as well as basic pH conditions delays the conformational transition to beta sheet. FRET results indicate that the distance between N and C-terminal regions increased on aging due to unfolding by self-assembly of the peptide into an organized beta sheet structure. Presence of osmolytes, prevented or decelerated the aggregation process of PrP(113-127) peptide.

Protective V127 prion variant prevents prion disease by interrupting the formation of dimer and fibril from molecular dynamics simulations

Recent studies uncovered a novel protective prion protein variant: V127 variant, which was reported intrinsically resistant to prion conversion and propagation. However, the structural basis of its protective effect is still unknown. To uncover the origin of the protective role of V127 variant, molecular dynamics simulations were performed to explore the influence of G127V mutation on two key processes of prion propagation: dimerization and fibril formation. The simulation results indicate V127 variant is unfavorable to form dimer by reducing the main-chain H-bond interactions. The simulations of formed fibrils consisting of β1 strand prove V127 variant will make the formed fibril become unstable and disorder. The weaker interaction energies between layers and reduced H-bonds number for V127 variant reveal this mutation is unfavorable to the formation of stable fibril. Consequently, we find V127 variant is not only unfavorable to the formation of dimer but also unfavorable to the formation of stable core and fibril, which can explain the mechanism on the protective role of V127 variant from the molecular level. Our findings can deepen the understanding of prion disease and may guide the design of peptide mimetics or small molecule to mimic the protective effect of V127 variant.

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.

PrP charge structure encodes interdomain interactions

Almost all proteins contain charged residues, and their chain distribution is tailored to fulfill essential ionic interactions for folding, binding and catalysis. Among proteins, the hinged two-domain chain of the cellular prion protein (PrP^C) exhibits a peculiar charge structure with unclear consequences in its structural malleability. To decipher the charge design role, we generated charge-reverted mutants for each domain and analyzed their effect on conformational and metabolic features. We found that charges contain the information for interdomain interactions. Use of dynamic light scattering and thermal denaturation experiments delineates the compaction of the α-fold by an electrostatic compensation between the polybasic 23–30 region and the α3 electronegative surface. This interaction increases stability and disfavors fibrillation. Independently of this structural effect, the N-terminal electropositive clusters regulate the α-cleavage efficiency. In the fibrillar state, use of circular dichroism, atomic-force and fluorescence microscopies reveal that the N-terminal positive clusters and the α3 electronegative surface dictate the secondary structure, the assembly hierarchy and the growth length of the fibril state. These findings show that the PrP charge structure functions as a code set up to ensure function and reduce pathogenic routes.

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.

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.