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 us 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 beta-content (similar to 30-40%) and low alpha-content (similar to 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.