Temperature-Induced Misfolding in Prion Protein: Evidence of Multiple Partially Disordered States Stabilized by Non-Native Hydrogen Bonds.

The structural basis of pathways of misfolding of a cellular prion (PrPC) into the toxic scrapie form (PrPSC) and identification of possible intermediates (e.g., PrP*) still eludes us. In this work, we have used a cumulative ∼65 µs of replica exchange molecular dynamics simulation data to construct the conformational free energy landscapes and capture the structural and thermodynamic characteristics associated with various stages of the thermal denaturation process in human prion protein. The temperature-dependent free energy surfaces consist of multiple metastable states stabilized by non-native contacts and hydrogen bonds, thus rendering the protein prone to misfolding. We have been able to identify metastable conformational states with high ß-content (∼30-40%) and low α-content (∼10-20%) that might be precursors of PrPSC oligomer formation. These conformations also involve participation of the unstructured N-terminal domain, and its role in misfolding has been investigated. All the misfolded or partially unfolded states are quite compact in nature despite having large deviations from the native structure. Although the number of native contacts decreases dramatically at higher temperatures, the radius of gyration and number of intraprotein hydrogen bonds and contacts remain relatively unchanged, leading to stabilization of the misfolded conformations by non-native interactions. Our results are in good agreement with the established view that the C-terminal regions of the second and third helices (H2 and H3, respectively) of mammal prions might be the Achilles heels of their stability, while separation of B1-H1-B2 and H2-H3 domains seems to play a key role, as well.

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 (PrPC) form needs to undergo substantial unfolding to a more stable PrPC* state, which may further oligomerize into the toxic scrapie (PrPSc) 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 PrPC* 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 PrPSc formation in previous simulation studies. We have further analyzed the origin of metastability of the PrPC 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.

Replica Exchange Molecular Dynamics Study of Dimerization in Prion Protein: Multiple Modes of Interaction and Stabilization.

The pathological forms of prions are known to be a result of misfolding, oligomerization, and aggregation of the cellular prion. While the mechanism of misfolding and aggregation in prions has been widely studied using both experimental and computational tools, the structural and energetic characterization of the dimer form have not garnered as much attention. On one hand dimerization can be the first step toward a nucleation-like pathway to aggregation, whereas on the other hand it may also increase the conformational stability preventing self-aggregation. In this work, we have used extensive all-atom replica exchange molecular dynamics simulations of both monomer and dimer forms of a mouse prion protein to understand the structural, dynamic, and thermodynamic stability of dimeric prion as compared to the monomeric form. We show that prion proteins can dimerize spontaneously being stabilized by hydrophobic interactions as well as intermolecular hydrogen bonding and salt bridge formation. We have computed the conformational free energy landscapes for both monomer and dimer forms to compare the thermodynamic stability and misfolding pathways. We observe large conformational heterogeneity among the various modes of interactions between the monomers and the strong intermolecular interactions may lead to as high as 20% ß-content. The hydrophobic regions in helix-2, surrounding coil regions, terminal regions along with the natively present ß-sheet region appear to actively participate in prion-prion intermolecular interactions. Dimerization seems to considerably suppress the inherent dynamic instability observed in monomeric prions, particularly because the regions of structural frustration constitute the dimer interface. Further, we demonstrate an interesting reversible coupling between the Q160-G131 interaction (which leads to inhibition of ß-sheet extension) and the G131-V161 H-bond formation.