Destabilization of prion protein induces a conformational change from normal prion protein (PrPC) to abnormal prion protein (PrPSC). Hydrophobic interaction is the main driving force for protein folding, and critically affects the stability and solvability. To examine the importance of the hydrophobic core in the PrP, we chose six amino acids (V176, V180, T183, V210, I215, and Y218) that make up the hydrophobic core at the middle of the H2-H3 bundle. A few pathological mutants of these amino acids have been reported, such as V176G, V180I, T183A, V210I, I215V, and Y218N. We focused on how these pathologic mutations affect the hydrophobic core and thermostability of PrP. For this, we ran a temperature-based replica-exchange molecular dynamics (T-REMD) simulation, with a cumulative simulation time of 28 μs, for extensive ensemble sampling. From the T-REMD ensemble, we calculated the protein folding free energy difference between wild-type and mutant PrP using the thermodynamic integration (TI) method. Our results showed that pathological mutants V176G, T183A, I215V, and Y218N decrease the PrP stability. At the atomic level, we examined the change in pair-wise hydrophobic interactions from valine-valine to valine-isoleucine (and vice versa), which is induced by mutation V180I, V210I (I215V) at the 180th–210th (176th–215th) pair. Finally, we investigated the importance of the π-stacking between Y218 and F175.