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

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

Introducing a Rigid Loop Structure from Deer into Mouse Prion Protein Increases Its Propensity for Misfolding In Vitro

Prion diseases are fatal neurodegenerative disorders characterized by misfolding of the cellular prion protein (PrP^c) into the disease-associated isoform (PrP^Sc) that has increased β-sheet content and partial resistance to proteolytic digestion. Prion diseases from different mammalian species have varying propensities for transmission upon exposure of an uninfected host to the infectious agent. Chronic Wasting Disease (CWD) is a highly transmissible prion disease that affects free ranging and farmed populations of cervids including deer, elk and moose, as well as other mammals in experimental settings. The molecular mechanisms allowing CWD to maintain comparatively high transmission rates have not been determined. Previous work has identified a unique structural feature in cervid PrP, a rigid loop between β-sheet 2 and α-helix 2 on the surface of the protein. This study was designed to test the hypothesis that the rigid loop has a direct influence on the misfolding process. The rigid loop was introduced into murine PrP as the result of two amino acid substitutions: S170N and N174T. Wild-type and rigid loop murine PrP were expressed in E. coli and purified. Misfolding propensity was compared for the two proteins using biochemical techniques and cell free misfolding and conversion systems. Murine PrP with a rigid loop misfolded in cell free systems with greater propensity than wild type murine PrP. In a lipid-based conversion assay, rigid loop PrP converted to a PK resistant, aggregated isoform at lower concentrations than wild-type PrP. Using both proteins as substrates in real time quaking-induced conversion, rigid loop PrP adopted a misfolded isoform more readily than wild type PrP. Taken together, these findings may help explain the high transmission rates observed for CWD within cervids.     

Did the prion protein become vulnerable to misfolding after an evolutionary divide and conquer event?

Despite high sequence identity among mammalian prion proteins (PrPs), mammals have varying rates of susceptibility to prion disease resulting in a so-called species barrier. The species barrier follows no clear pattern, with closely related species or similar sequences being no more likely to infect each other, and remains an unresolved enigma. Variation of the conformationally flexible regions may alter the thermodynamics of the conformational change, commonly referred to as the conformational conversion, which occurs in the pathogenic process of the mammalian prion protein. A conformational ensemble scenario is supported by the species barrier in prion disease and evidence that there are strains of pathogenic prion with different conformations within species. To study how conformational flexibility has evolved in the prion protein, an investigation was undertaken on the evolutionary dynamics of structurally disordered regions in the mammalian prion protein, non-mammalian prion protein that is not vulnerable to prion disease, and remote homologs Doppel and Shadoo. Structural disorder prediction analyzed in an evolutionary context revealed that the occurrence of increased or altered conformational flexibility in mammalian PrPs coincides with key events among PrP, Doppel, and Shadoo. Comparatively rapid evolutionary dynamics of conformational flexibility in the prion protein suggest that the species barrier is not a static phenomenon. A small number of amino acid substitutions can repopulate the conformational ensemble and have a disproportionately large effect on pathogenesis.     

Relationships between protein structure and dynamics from a database of NMR-derived backbone order parameters

The amplitude of protein backbone NH group motions on a time-scale faster than molecular tumbling may be determined by analysis of (15)N NMR relaxation data according to the Lipari-Szabo model free formalism. An internet-accessible database has been compiled containing 1855 order parameters from 20 independent NMR relaxation studies on proteins whose three-dimensional structures are known. A series of statistical analyses has been performed to identify relationships between the structural features and backbone dynamics of these proteins. Comparison of average order parameters for different amino acid types indicates that amino acids with small side-chains tend to have greater backbone flexibility than those with large side-chains. In addition, the motions of a given NH group are also related to the sizes of the neighboring amino acids in the primary sequence. The secondary structural environment appears to influence backbone dynamics relatively weakly, with only subtle differences between the order parameter distributions of loop structures and regular hydrogen bonded secondary structure elements. However, NH groups near helix termini are more mobile on average than those in the central regions of helices. Tertiary structure influences are also relatively weak but in the expected direction, with more exposed residues being more flexible on average than residues that are relatively inaccessible to solvent.     

Application of the random coil index to studying protein flexibility

Protein flexibility lies at the heart of many protein–ligand binding events and enzymatic activities. However, the experimental measurement of protein motions is often difficult, tedious and error-prone. As a result, there is a considerable interest in developing simpler and faster ways of quantifying protein flexibility. Recently, we described a method, called Random Coil Index (RCI), which appears to be able to quantitatively estimate model-free order parameters and flexibility in protein structural ensembles using only backbone chemical shifts. Because of its potential utility, we have undertaken a more detailed investigation of the RCI method in an attempt to ascertain its underlying principles, its general utility, its sensitivity to chemical shift errors, its sensitivity to data completeness, its applicability to other proteins, and its general strengths and weaknesses. Overall, we find that the RCI method is very robust and that it represents a useful addition to traditional methods of studying protein flexibility. We have implemented many of the findings and refinements reported here into a web server that allows facile, automated predictions of model-free order parameters, MD RMSF and NMR RMSD values directly from backbone ¹H, ¹³C and ¹⁵N chemical shift assignments. The server is available at . Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

A Structural Overview of the Vertebrate Prion Proteins

Among the diseases caused by protein misfolding is the family associated with the prion protein (PrP). This is a small extracellular membrane-anchored molecule of yet unknown function. Understanding how PrP folds both into its cellular and pathological forms is thought to be crucial for explaining protein misfolding in general and the specific role of PrP in disease. Since the first structure determination, an increasing number of structural studies of PrP have become available, showing that the protein is formed by a flexible N-terminal region and a highly conserved globular C-terminal domain. We review here the current knowledge on PrP structure. We focus on vertebrate PrPs and analyse in detail the similarities and the differences among the coordinates of the C-terminal domain of PrP from different species, in search for understanding the mechanism of disease-causing mutations and the molecular bases of species barrier.Key Words: amyloid, NMR, prion, scrapie, structure, X-ray     

Single-molecule approaches to prion protein misfolding

The structural conversion of the prion protein PrP into a transmissible, misfolded form is the central element of prion disease, yet there is little consensus as to how it occurs. Key aspects of conversion into the diseased state remain unsettled, from details about the earliest stages of misfolding such as the involvement of partially- or fully-unfolded intermediates to the structure of the infectious state. Part of the difficulty in understanding the structural conversion arises from the complexity of the underlying energy landscapes. Single molecule methods provide a powerful tool for probing complex folding pathways as in prion misfolding, because they allow rare and transient events to be observed directly. We discuss recent work applying single-molecule probes to study misfolding in prion proteins, and what it has revealed about the folding dynamics of PrP that may underlie its unique behavior. We also discuss single-molecule studies probing the interactions that stabilize non-native structures within aggregates, pointing the way to future work that may help identify the microscopic events triggering pathogenic conversion. Although single-molecule approaches to misfolding are relatively young, they have a promising future in prion science.     

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

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

A functional role for correlated motion in the N-terminal RNA-binding domain of human U1A protein

The N-terminal RNA-binding domain of the human U1A protein (RBD1) undergoes local conformational changes upon binding to its target RNA. Here, the wild-type RBD1 and two mutants are examined with molecular dynamics simulations that are analyzed using the reorientational eigenmode dynamics (RED) formalism. The results reveal changes in the magnitude and extent of coupled intra-domain motions resulting from single amino acid substitutions. Interpretation of the novel RED results and corresponding NMR relaxation data suggests that the loss of collective motions in the mutants could account for their weak RNA-binding.     

Studies on the structural stability of rabbit prion probed by molecular dynamics simulations of its wild-type and mutants

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases that affect humans and animals. Rabbits are the only mammalian species reported to be resistant to infection from prion diseases isolated from other species (Vorberg et al., 2003). Fortunately, the NMR structure of rabbit prion (124-228) (PDB entry 2FJ3), the NMR structure of rabbit prion protein mutation S173N (PDB entry 2JOH) and the NMR structure of rabbit prion protein mutation I214V (PDB entry 2JOM) were released recently. This paper studies these NMR structures by molecular dynamics simulations. Simulation results confirm the structural stability of wild-type rabbit prion, and show that the salt bridge between D177 and R163 greatly contributes to the structural stability of rabbit prion protein.     

Nascent β Structure in the Elongated Hydrophobic Region of a Gerstmann–Sträussler–Scheinker PrP Allele

Prion diseases are neurodegenerative disorders caused by the misfolding of the cellular prion protein (PrP^C). Gerstmann–Sträussler–Scheinker syndrome is an inherited prion disease with one early-onset allele (HRdup) containing an eight-amino-acid insertion; this LGGLGGYV insert is positioned after valine 129 (human PrP^C sequence) in a hydrophobic tract in the natively disordered region. Here we have characterized the structure and explored the molecular motions and dynamics of HRdup PrP and a control allele. High-resolution NMR data suggest that the core of HRdup has a canonical PrP^C structure, yet a nascent β-structure is observed in the flexible elongated hydrophobic region of HRdup. In addition, using mouse PrP^C sequence, we observed that a methionine/valine polymorphism at codon 128 (equivalent of methionine/valine 129 in human sequence) and oligomerization caused by high protein concentration affects conformational exchange dynamics at residue G130. We hypothesize that with the β-structure at the N-terminus, the hydrophobic region of HRdup can adopt a fully extended configuration and fold back to form an extended β-sheet with the existing β-sheet. We propose that these structures are early chemical events in disease pathogenesis.     

Prions: Generation and Spread Versus Neurotoxicity

Neurodegenerative diseases are characterized by the aggregation of misfolded proteins in the brain. Among these disorders are the prion diseases, which are transmissible, and in which the misfolded proteins (“prions”) are also the infectious agent. Increasingly, it appears that misfolded proteins in Alzheimer and Parkinson diseases and the tauopathies also propagate in a “prion-like” manner. However, the association between prion formation, spread, and neurotoxicity is not clear. Recently, we showed that in prion disease, protein misfolding leads to neurodegeneration through dysregulation of generic proteostatic mechanisms, specifically, the unfolded protein response. Genetic and pharmacological manipulation of the unfolded protein response was neuroprotective despite continuing prion replication, hence dissociating this from neurotoxicity. The data have clear implications for treatment across the spectrum of these disorders, targeting pathogenic processes downstream of protein misfolding.     

Molecular dynamics studies on the buffalo prion protein

It was reported that buffalo is a low susceptibility species resisting to transmissible spongiform encephalopathies (TSEs) (same as rabbits, horses, and dogs). TSEs, also called prion diseases, are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of species (except for rabbits, dogs, horses, and buffalo), manifesting as scrapie in sheep and goats; bovine spongiform encephalopathy (BSE or “mad–cow” disease) in cattle; chronic wasting disease in deer and elk; and Creutzfeldt–Jakob diseases, Gerstmann–Sträussler–Scheinker syndrome, fatal familial insomnia, and Kulu in humans etc. In molecular structures, these neurodegenerative diseases are caused by the conversion from a soluble normal cellular prion protein (PrP^C), predominantly with α-helices, into insoluble abnormally folded infectious prions (PrP^Sc), rich in β-sheets. In this article, we studied the molecular structure and structural dynamics of buffalo PrP^C (BufPrP^C), in order to understand the reason why buffalo is resistant to prion diseases. We first did molecular modeling of a homology structure constructed by one mutation at residue 143 from the NMR structure of bovine and cattle PrP(124–227); immediately we found that for BufPrP^C(124–227), there are five hydrogen bonds (HBs) at Asn143, but at this position, bovine/cattle do not have such HBs. Same as that of rabbits, dogs, or horses, our molecular dynamics studies also revealed there is a strong salt bridge (SB) ASP178–ARG164 (O–N) keeping the β2–α2 loop linked in buffalo. We also found there is a very strong HB SER170–TYR218 linking this loop with the C-terminal end of α-helix H3. Other information, such as (i) there is a very strong SB HIS187–ARG156 (N–O) linking α-helices H2 and H1 (if mutation H187R is made at position 187, then the hydrophobic core of PrP^C will be exposed (L.H. Zhong (2010). Exposure of hydrophobic core in human prion protein pathogenic mutant H187R. Journal of Biomolecular Structure and Dynamics 28(3), 355–361)), (ii) at D178, there is a HB Y169–D178 and a polar contact R164–D178 for BufPrP^C instead of a polar contact Q168–D178 for bovine PrP^C (C.J. Cheng, & V. Daggett. (2014). Molecular dynamics simulations capture the misfolding of the bovine prion protein at acidic pH. Biomolecules 4(1), 181–201), (iii) BufPrP^C owns three 3₁₀ helices at 125–127, 152–156, and in the β2–α2 loop, respectively, and (iv) in the β2–α2 loop, there is a strong π–π stacking and a strong π–cation F175–Y169–R164.(N)NH2, has been discovered.     

Production of cattle lacking prion protein

Prion diseases are caused by propagation of misfolded forms of the normal cellular prion protein PrP(C), such as PrP(BSE) in bovine spongiform encephalopathy (BSE) in cattle and PrP(CJD) in Creutzfeldt-Jakob disease (CJD) in humans. Disruption of PrP(C) expression in mice, a species that does not naturally contract prion diseases, results in no apparent developmental abnormalities. However, the impact of ablating PrP(C) function in natural host species of prion diseases is unknown. Here we report the generation and characterization of PrP(C)-deficient cattle produced by a sequential gene-targeting system. At over 20 months of age, the cattle are clinically, physiologically, histopathologically, immunologically and reproductively normal. Brain tissue homogenates are resistant to prion propagation in vitro as assessed by protein misfolding cyclic amplification. PrP(C)-deficient cattle may be a useful model for prion research and could provide industrial bovine products free of prion proteins.     

Fibril formation of the rabbit/human/bovine prion proteins

Prion diseases are infectious fatal neurodegenerative diseases including Creutzfeldt-Jakob disease in humans and bovine spongiform encephalopathy in cattle. The misfolding and conversion of cellular PrP in such mammals into pathogenic PrP is believed to be the key procedure. Rabbits are among the few mammalian species that exhibit resistance to prion diseases, but little is known about the molecular mechanism underlying such resistance. Here, we report that the crowding agents Ficoll 70 and dextran 70 have different effects on fibrillization of the recombinant full-length PrPs from different species: although these agents dramatically promote fibril formation of the proteins from human and cow, they significantly inhibit fibrillization of the rabbit protein by stabilizing its native state. We also find that fibrils formed by the rabbit protein contain less β-sheet structure and more α-helix structure than those formed by the proteins from human and cow. In addition, amyloid fibrils formed by the rabbit protein do not generate a proteinase K-resistant fragment of 15–16-kDa, but those formed by the proteins from human and cow generate such proteinase K-resistant fragments. Together, these results suggest that the strong inhibition of fibrillization of the rabbit PrP by the crowded physiological environment and the absence of such a protease-resistant fragment for the rabbit protein could be two of the reasons why rabbits are resistant to prion diseases.     

The intrinsic stability of the human prion β-sheet region investigated by molecular dynamics

Human prion diseases are neurodegenerative disorders associated to the misfolding of the prion protein (PrP). Common features of prion disorders are the fibrillar amyloid deposits and the formation of prefibrillar oligomeric species also suggested as the origin of cytotoxicity associated with diseases. Although the process of PrP misfolding has been extensively investigated, many crucial aspects of this process remain unclear. We have here carried out a molecular dynamics study to evaluate the intrinsic dynamics of PrP β-sheet, a region that is believed to play a crucial role in prion aggregation. Moreover, as this region mediates protein association in dimeric assemblies frequently observed in prion crystallographic investigations, we also analyzed the dynamics of these intermolecular interactions. The extensive sampling of replica exchange shows that the native antiparallel β-structure of the prion is endowed with a remarkable stability. Therefore, upon unfolding, the persistence of a structured β-region may seed molecular association and influence the subsequent phases of the aggregation process. The analysis of the four-stranded β-sheet detected in the dimeric assemblies of PrP shows a tendency of this region to form dynamical structured states. The impact on the β-sheet structure and dynamics of disease associated point mutations has also been evaluated.     

Different misfolding mechanisms converge on common conformational changes

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.     

Sequence-dependent Prion Protein Misfolding and Neurotoxicity

Prion diseases are neurodegenerative disorders caused by misfolding of the normal prion protein (PrP) into a pathogenic “scrapie” conformation. To better understand the cellular and molecular mechanisms that govern the conformational changes (conversion) of PrP, we compared the dynamics of PrP from mammals susceptible (hamster and mouse) and resistant (rabbit) to prion diseases in transgenic flies. We recently showed that hamster PrP induces spongiform degeneration and accumulates into highly aggregated, scrapie-like conformers in transgenic flies. We show now that rabbit PrP does not induce spongiform degeneration and does not convert into scrapie-like conformers. Surprisingly, mouse PrP induces weak neurodegeneration and accumulates small amounts of scrapie-like conformers. Thus, the expression of three highly conserved mammalian prion proteins in transgenic flies uncovered prominent differences in their conformational dynamics. How these properties are encoded in the amino acid sequence remains to be elucidated.     

Proteome analysis of cultivated vascular smooth muscle cells from a CADASIL patient

Cerebral autosomal dominant arteriopathy with subcortical infarcts and leukoencephalopathy (CADASIL) is a vascular dementing disease caused by mutations in the NOTCH3 gene, most which are missense mutations leading to an uneven number of cysteine residues in epidermal growth factor-like repeats in the extracellular domain of Notch3 receptor (N3ECD). CADASIL is characterized by degeneration of vascular smooth muscle cells (VSMC) and accumulation of N3ECD on the VSMCs of small and middle-sized arteries. Recent studies have demonstrated that impairment of Notch3 signaling is not the primary cause of the disease. In the present study we used proteomic analysis to characterize the protein expression pattern of a unique material of genetically genuine cultured human CADASIL VSMCs. We identified 11 differentially expressed proteins, which are involved in protein degradation and folding, contraction of VSMCs, and cellular stress. Our findings indicate that misfolding of Notch3 may cause endoplasmic reticulum stress and activation of unfolded protein response, leading to increased reactive oxygen species and inhibition of cell proliferation. In addition, upregulation of contractile proteins suggests an alteration in the signaling system of VSMC contraction. The accumulation of N3ECD on the cell surface possibly upregulates the angiotensin II regulatory feedback loop and thereby enhances the readiness of the cells to respond to angiotensin II stimulation.