Comparison of the unfolding and oligomerization of human prion protein under acidic and neutral environments by molecular dynamics simulations

Cellular prion protein gene polymorphisms linked to differential scrapie susceptibility correlate with distinct residue connectivity between secondary structure elements

The conformational conversion of the cellular prion protein (PrP^C) to the misfolded and aggregated isoform, termed scrapie prion protein (PrP^Sc), is key to the development of a group of neurodegenerative diseases known as transmissible spongiform encephalopathies (TSEs). Although the conversion mechanism is not fully understood, the role of gene polymorphisms in varying susceptibilities to prion diseases is well established. In ovine, specific gene polymorphisms in PrP^C alter prion disease susceptibility: the Valine136-Glutamine171 variant (Susceptible structure) displays high susceptibility to classical scrapie while the Alanine136-Arginine171 variant (Resistant structure) displays reduced susceptibility. The opposite trend has been reported in atypical scrapie. Despite the differentiation between classical and atypical scrapie, a complete understanding of the effect of polymorphisms on the structural dynamics of PrP^C is lacking. From our structural bioinformatics study, we propose that polymorphisms locally modulate the network of residue interactions in the globular C-terminus of the ovine recombinant prion protein while maintaining the overall fold. Although the two variants we examined exhibit a densely connected group of residues that includes both β-sheets, the β2-α2 loop and the N-terminus of α-helix 2, only in the Resistant structure do most residues of α-helix 2 belong to this group. We identify the structural role of Valine136Alanine and Glutamine171Arginine: modulation of residue interaction networks that affect the connectivity between α-helix 2 and α-helix 3. We propose blocking interactions of residue 171 as a potential target for the design of therapeutics to prevent efficient PrP^C misfolding. We discuss our results in the context of initial PrP^C conversion and extrapolate to recently proposed PrP^Sc structures.Communicated by Ramaswamy H. Sarma.

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.