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

First report of a novel polymorphism and genetic characteristics of the leporine prion protein (PRNP) gene

Transmissible spongiform encephalopathies (TSEs) have been reported in a broad spectrum of hosts. The genetic polymorphisms and characteristics of the prion protein (PRNP) gene have a vital impact on the development of TSEs. Notably, natural TSE infection cases have never been reported in rabbits, and genetic variations of the leporine PRNP gene have not been investigated to date. To identify leporine PRNP gene polymorphism, we performed amplicon sequencing in 203 rabbits. We report a novel single nucleotide polymorphism on the leporine PRNP gene. In addition, we performed a comparative analysis of amino acid sequences of prion protein (PrP) across several hosts using ClustalW2. Furthermore, we evaluated the effect of changes of unique leporine PrP amino acids with those conserved among various species using Swiss-Pdb Viewer. Interestingly, we found seven unique leporine amino acids, and the change of unique leporine amino acids with those conserved among other species, including S175N, Q221K, Q221R, A226Y, A230G, and A230S, was predicted to reduce hydrogen bonds in leporine PrP.

Rabbit PrP Is Partially Resistant to in vitro Aggregation Induced by Different Biological Cofactors

Prion diseases have been described in humans and other mammals, including sheep, goats, cattle, and deer. Since mice, hamsters, and cats are susceptible to prion infection, they are often used to study the mechanisms of prion infection and conversion. Mammals, such as horses and dogs, however, do not naturally contract the disease and are resistant to infection, while others, like rabbits, have exhibited low susceptibility. Infection involves the conversion of the cellular prion protein (PrP^C) to the scrapie form (PrP^Sc), and several cofactors have already been identified as important adjuvants in this process, such as glycosaminoglycans (GAGs), lipids, and nucleic acids. The molecular mechanisms that determine transmissibility between species remain unclear, as well as the barriers to transmission. In this study, we examine the interaction of recombinant rabbit PrP^C (RaPrP) with different biological cofactors such as GAGs (heparin and dermatan sulfate), phosphatidic acid, and DNA oligonucleotides (A1 and D67) to evaluate the importance of these cofactors in modulating the aggregation of rabbit PrP and explain the animal’s different degrees of resistance to infection. We used spectroscopic and chromatographic approaches to evaluate the interaction with cofactors and their effect on RaPrP aggregation, which we compared with murine PrP (MuPrP). Our data show that all cofactors induce RaPrP aggregation and exhibit pH dependence. However, RaPrP aggregated to a lesser extent than MuPrP in the presence of any of the cofactors tested. The binding affinity with cofactors does not correlate with these low levels of aggregation, suggesting that the latter are related to the stability of PrP at acidic pH. The absence of the N-terminus affected the interaction with cofactors, influencing the efficiency of aggregation. These findings demonstrate that the interaction with polyanionic cofactors is related to rabbit PrP being less susceptible to aggregation in vitro and that the N-terminal domain is important to the efficiency of conversion, increasing the interaction with cofactors. The decreased effect of cofactors in rabbit PrP likely explains its lower propensity to prion conversion.

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.

In Vitro Approach To Identify Key Amino Acids in Low Susceptibility of Rabbit Prion Protein to Misfolding

Prion diseases, or transmissible spongiform encephalopathies (TSEs), are a group of rare progressive neurodegenerative disorders caused by an abnormally folded prion protein (PrP^Sc). This is capable of transforming the normal cellular prion protein (PrP^C) into new infectious PrP^Sc Interspecies prion transmissibility studies performed by experimental challenge and the outbreak of bovine spongiform encephalopathy that occurred in the late 1980s and 1990s showed that while some species (sheep, mice, and cats) are readily susceptible to TSEs, others are apparently resistant (rabbits, dogs, and horses) to the same agent. To study the mechanisms of low susceptibility to TSEs of certain species, the mouse-rabbit transmission barrier was used as a model. To identify which specific amino acid residues determine high or low susceptibility to PrP^Sc propagation, protein misfolding cyclic amplification (PMCA), which mimics PrP^C-to-PrP^Sc conversion with accelerated kinetics, was used. This allowed amino acid substitutions in rabbit PrP and accurate analysis of misfolding propensities. Wild-type rabbit recombinant PrP could not be misfolded into a protease-resistant self-propagating isoform in vitro despite seeding with at least 12 different infectious prions from diverse origins. Therefore, rabbit recombinant PrP mutants were designed to contain every single amino acid substitution that distinguishes rabbit recombinant PrP from mouse recombinant PrP. Key amino acid residue substitutions were identified that make rabbit recombinant PrP susceptible to misfolding, and using these, protease-resistant misfolded recombinant rabbit PrP was generated. Additional studies characterized the mechanisms by which these critical amino acid residue substitutions increased the misfolding susceptibility of rabbit PrP.IMPORTANCE Prion disorders are invariably fatal, untreatable diseases typically associated with long incubation periods and characteristic spongiform changes associated with neuronal loss in the brain. Development of any treatment or preventative measure is dependent upon a detailed understanding of the pathogenesis of these diseases, and understanding the mechanism by which certain species appear to be resistant to TSEs is critical. Rabbits are highly resistant to naturally acquired TSEs, and even under experimental conditions, induction of clinical disease is not easy. Using recombinant rabbit PrP as a model, this study describes critical molecular determinants that confer this high resistance to transmissible spongiform encephalopathies.

X-ray structural and molecular dynamical studies of the globular domains of cow, deer, elk and Syrian hamster prion proteins

Misfolded prion proteins are the cause of neurodegenerative diseases that affect many mammalian species, including humans. Transmission of the prion diseases poses a considerable public-health risk as a specific prion disease such as bovine spongiform encephalopathy can be transferred to humans and other mammalian species upon contaminant exposure. The underlying mechanism of prion propagation and the species barriers that control cross species transmission has been investigated quite extensively. So far a number of prion strains have been characterized and those have been intimately linked to species-specific infectivity and other pathophysiological manifestations. These strains are encoded by a protein-only agent, and have a high degree of sequence identity across mammalian species. The molecular events that lead to strain differentiation remain elusive. In order to contribute to the understanding of strain differentiation, we have determined the crystal structures of the globular, folded domains of four prion proteins (cow, deer, elk and Syrian hamster) bound to the POM1 antibody fragment Fab. Although the overall structural folds of the mammalian prion proteins remains extremely similar, there are several local structural variations observed in the misfolding-initiator motifs. In additional molecular dynamics simulation studies on these several prion proteins reveal differences in the local fluctuations and imply that these differences have possible roles in the unfolding of the globular domains. These local variations in the structured domains perpetuate diverse patterns of prion misfolding and possibly facilitate the strain selection and adaptation.

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

Progress on low susceptibility mechanisms of transmissible spongiform encephalopathies

Transmissible spongiform encephalopathies (TSEs), also known as prion diseases, are a group of fatal neurodegenerative diseases detected in a wide range of mammalian species. The "protein-only" hypothesis of TSE suggests that prions are transmissible particles devoid of nucleic acid and the primary pathogenic event is thought to be the conversion of cellular prion protein (PrP(C)) into the disease-associated isoform (PrP(Sc)). According to susceptibility to TSEs, animals can be classified into susceptible species and low susceptibility species. In this review we focus on several species with low susceptibility to TSEs: dogs, rabbits, horses and buffaloes. We summarize recent studies into the characteristics of low susceptibility regarding protein structure, and biochemical and genetic properties.

Molecular dynamics studies on the NMR structures of rabbit prion protein wild type and mutants: surface electrostatic charge distributions

Prion diseases are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of mammalian species such as sheep and goats, cattle, deer and elk, and humans. But for rabbits, studies have shown that they have a low susceptibility to be infected by prion diseases. This paper does molecular dynamics (MD) studies of rabbit NMR structures (of the wild type and its two mutants of two surface residues), in order to understand the specific mechanism of rabbit prion proteins (RaPrP(C)). Protein surface electrostatic charge distributions are specially focused to analyze the MD trajectories. This paper can conclude that surface electrostatic charge distributions indeed contribute to the structural stability of wild-type RaPrP(C); this may be useful for the medicinal treatment of prion diseases.

Decipher the mechanisms of rabbit's low susceptibility to prion infection

Rabbits have low susceptibility to prion infection. Studies on prion protein (PrP) from animal species of different susceptibility to prion diseases identified key amino acid residues, specific motif, and special features in rabbit prion protein (RaPrP(C)) that contribute to the stability of rabbit PrP(C) and low susceptibility to prion infection. However, there is no evidence showing that rabbits are completely resistant to prion diseases. It has been reported that the rabbit prion could be generated in vitro through protein misfolding cyclic amplification and proved to be infectious and transmissible. Here, we reviewed studies on rabbit-specific PrP structures and features in relation to rabbit's low susceptibility to prion infection.

Molecular dynamics studies on the NMR and X-ray structures of rabbit prion proteins

Prion diseases, traditionally referred to as transmissible spongiform encephalopathies (TSEs), are invariably fatal and highly infectious neurodegenerative diseases that affect a wide variety of mammalian species, 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. 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)), and the conversion of PrP(C) to PrP(Sc) is believed to involve conformational change from a predominantly α-helical protein to one rich in β-sheet structure. Such a conformational change may be amenable to study by molecular dynamics (MD) techniques. For rabbits, classical studies show that they have a low susceptibility to be infected by PrP(Sc), but recently it was reported that rabbit prions can be generated through saPMCA (serial automated Protein Misfolding Cyclic Amplification) in vitro and the rabbit prion is infectious and transmissible. In this paper, we first do a detailed survey on the research advances of rabbit prion protein (RaPrP) and then we perform MD simulations on the NMR and X-ray molecular structures of rabbit prion protein wild-type and mutants. The survey shows to us that rabbits were not challenged directly in vivo with other known prion strains and the saPMCA result did not pass the test of the known BSE strain of cattle. Thus, we might still look rabbits as a prion resistant species. MD results indicate that the three α-helices of the wild-type are stable under the neutral pH environment (but under low pH environment the three α-helices have been unfolded into β-sheets), and the three α-helices of the mutants (I214V and S173N) are unfolded into rich β-sheet structures under the same pH environment. In addition, we found an interesting result that the salt bridges such as ASP201-ARG155, ASP177-ARG163 contribute greatly to the structural stability of RaPrP.

Prion protein oligomer and its neurotoxicity

The prion diseases, also known as transmissible spongiform encephalopathies, are fatal neurodegenerative disorders. According to the 'protein only' hypothesis, the key molecular event in the pathogenesis of prion disease is the conformational conversion of the host-derived cellular prion protein (PrP(C)) into a misfolded form (scrapie PrP, PrP(Sc)). Increasing evidence has shown that the most infectious factor is the smaller subfibrillar oligomers formed by prion proteins. Both the prion oligomer and PrP(Sc) are rich in β-sheet structure and resistant to the proteolysis of proteinase K. The prion oligomer is soluble in physiologic environments whereas PrP(Sc) is insoluble. Various prion oligomers are formed in different conditions. Prion oligomers exhibited more neurotoxicity both in vitro and in vivo than the fibrillar forms of PrP(Sc), implying that prion oligomers could be potential drug targets for attacking prion diseases. In this article, we describe recent experimental evidence regarding prion oligomers, with a special focus on prion oligomer formation and its neurotoxicity.

N-terminal helix-cap in α-helix 2 modulates β-state misfolding in rabbit and hamster prion proteins

Susceptibility of a particular species to prion disease is affected by small differences in the sequence of PrP and correlates with the propensity of its PrP to assume the β-state. A helix-cap motif in the β2-α2-loop of native α-helical rabbit PrP, a resistant species, contains sequence differences that influence intra- and interspecies transmission. To determine the effect the helix-cap motif on β-state refolding propensity, we mutated S170N, S174N, and S170N/S174N of the rabbit PrP helix-cap to resemble that of hamster PrP and conversely, N170S, N174S, and N170S/N174S of hamster PrP to resemble the helix-cap of rabbit PrP. High-resolution crystal structures (1.45-1.6 Å) revealed that these mutations ablate hydrogen-bonding interactions within the helix-cap motif in rabbit PrP(C). They also alter the β-state-misfolding propensity of PrP; the serine mutations in hamster PrP decrease the propensity up to 35%, whereas the asparagine mutations in rabbit PrP increase it up to 42%. Rapid dilution of rabbit and hamster into β-state buffer conditions causes quick conversion to β-state monomers. Kinetic monitoring using size-exclusion chromatography showed that the monomer population decreases exponentially mirrored by an increase in an octameric species. The monomer-octamer transition rates are faster for hamster than for rabbit PrP. The N170S/N174S mutant of hamster PrP has a smaller octamer component at the endpoint compared to the wild-type, whereas the kinetics of octamer formation in mutant and wild-type rabbit PrP are comparable. These findings demonstrate that the sequence of the β2-α2 helix-cap affects refolding to the β-state and subsequently, may influence susceptibility to prion disease.

Naturally prion resistant mammals: a utopia?

Each known abnormal prion protein (PrP (Sc) ) is considered to have a specific range and therefore the ability to infect some species and not others. Consequently, some species have been assumed to be prion disease resistant as no successful natural or experimental challenge infections have been reported. This assumption suggested that, independent of the virulence of the PrP (Sc) strain, normal prion protein (PrP (C) ) from these 'resistant' species could not be induced to misfold. Numerous in vitro and in vivo studies trying to corroborate the unique properties of PrP (Sc) have been undertaken. The results presented in the article "Rabbits are not resistant to prion infection" demonstrated that normal rabbit PrP (C) , which was considered to be resistant to prion disease, can be misfolded to PrP (Sc) and subsequently used to infect and transmit a standard prion disease to leporids. Using the concept of species resistance to prion disease, we will discuss the mistake of attributing species specific prion disease resistance based purely on the absence of natural cases and incomplete in vivo challenges. The BSE epidemic was partially due to an underestimation of species barriers. To repeat this error would be unacceptable, especially if present knowledge and techniques can show a theoretical risk. Now that the myth of prion disease resistance has been refuted it is time to re-evaluate, using the new powerful tools available in modern prion laboratories, whether any other species could be at risk.

Relative and regional stabilities of the hamster, mouse, rabbit, and bovine prion proteins toward urea unfolding assessed by nuclear magnetic resonance and circular dichroism spectroscopies

The residue-specific urea-induced unfolding patterns of recombinant prion proteins from different species (bovine, rabbit, mouse, and Syrian hamster) were monitored using high-resolution (1)H nuclear magnetic resonance (NMR) spectroscopy. Protein constructs of different lengths, and with and without a His tag attached at the N-terminus, were studied. The various species showed different overall sensitivities toward urea denaturation with stabilities in the following order: hamster ≤ mouse < rabbit < bovine protein. This order is in agreement with recent circular dichroism (CD) spectroscopic measurements for several species [Khan, M. Q. (2010) Proc. Natl. Acad. Sci. U.S.A.107, 19808-19813] and for the bovine protein presented herein. The [urea](1/2) values determined by CD spectroscopy parallel those of the most stable residues observed by NMR spectroscopy. Neither the longer constructs containing an additional hydrophobic region nor the His tag influenced the stability of the structured domain of the constructs studied. The effect of the S174N mutation in rabbit PrP(C) was also investigated. The rank order of the regional stabilities within each protein remained the same for all species. In particular, the residues in the β-sheet region in all four species were more sensitive to urea-induced unfolding than residues in the α2 and α3 helical regions. These observations indicate that the regional specific unfolding pattern is the same for the four mammalian prion proteins studied but militate against the idea that PrP(Sc) formation is linked with the global stability of PrP(C).

Solution structure and dynamics of the I214V mutant of the rabbit prion protein

BACKGROUND

The conformational conversion of the host-derived cellular prion protein (PrP(C)) into the disease-associated scrapie isoform (PrP(Sc)) is responsible for the pathogenesis of transmissible spongiform encephalopathies (TSEs). Various single-point mutations in PrP(C)s could cause structural changes and thereby distinctly influence the conformational conversion. Elucidation of the differences between the wild-type rabbit PrP(C) (RaPrP(C)) and various mutants would be of great help to understand the ability of RaPrP(C) to be resistant to TSE agents.

METHODOLOGY/PRINCIPAL FINDINGS

We determined the solution structure of the I214V mutant of RaPrP(C)(91-228) and detected the backbone dynamics of its structured C-terminal domain (121-228). The I214V mutant displays a visible shift of surface charge distribution that may have a potential effect on the binding specificity and affinity with other chaperones. The number of hydrogen bonds declines dramatically. Urea-induced transition experiments reveal an obvious decrease in the conformational stability. Furthermore, the NMR dynamics analysis discloses a significant increase in the backbone flexibility on the pico- to nanosecond time scale, indicative of lower energy barrier for structural rearrangement.

CONCLUSIONS/SIGNIFICANCE

Our results suggest that both the surface charge distribution and the intrinsic backbone flexibility greatly contribute to species barriers for the transmission of TSEs, and thereby provide valuable hints for understanding the inability of the conformational conversion for RaPrP(C).

Unique structural characteristics of the rabbit prion protein

Rabbits are one of the few mammalian species that appear to be resistant to transmissible spongiform encephalopathies due to the structural characteristics of the rabbit prion protein (RaPrP(C)) itself. Here, we determined the solution structures of the recombinant protein RaPrP(C)-(91-228) and its S173N variant and detected the backbone dynamics of their structured C-terminal domains-(121-228). In contrast to many other mammalian PrP(C)s, loop 165-172, which connects β-sheet-2 and α-helix-2, is well-defined in RaPrP(C). For the first time, order parameters S(2) are obtained for residues in this loop region, indicating that loop 165-172 of RaPrP(C) is highly ordered. Compared with the wild-type RaPrP(C), less hydrogen bonds form in the S173N variant. The NMR dynamics analysis reveals a distinct increase in the structural flexibility of loop 165-172 and helix-3 after the S173N substitution, implying that the S173N substitution disturbs the long range interaction of loop 165-172 with helix-3, which further leads to a marked decrease in the global conformational stability. Significantly, RaPrP(C) possesses a unique charge distribution, carrying a continuous area of positive charges on the surface, which is distinguished from other PrP(C)s. The S173N substitution causes visible changes of the charge distribution around the recognition sites for the hypothetical protein X. Our results suggest that the ordered loop 165-172 and its interaction with helix-3, together with the unique distribution of surface electrostatic potential, significantly contribute to the unique structural characteristics of RaPrP(C).