Prion diseases are a family of infectious amyloid diseases affecting human and animals. Prion propagation in transmissible spongiform encephalopathies is associated with the unfolding and conversion of normal cellular prion protein into its pathogenic scrapie form. Understanding the fundamentals of prion protein aggregation caused by mutations is crucial to unravel the pathology of prion diseases. To help understand the contributions of individual residues to the stability of the human prion protein, we have carried out free energy simulations based on atomistic molecular dynamics trajectories. We focus on Met → Ala mutations at positions 205, 206 and 213, which are mostly buried residues located on helix 3 of the protein. The simulations predicted that all three mutations destabilize the prion protein. Changes in unfolding free energy upon mutation, ∆∆G, are 3.10 ± 0.79, 2.00 ± 0.26 and 3.06 ± 0.66 kcal/mol for M205A, M206A and M213A, respectively, in excellent agreement with the corresponding experimental values of 3.09 ± 0.28, 1.50 ± 0.34 and 3.12 ± 0.27 kcal/mol [T. Hart et al. (2009) PNAS 106, 5651–5656]. Component analysis indicates that the major contributions to the loss of protein stability arise from van der Waals interactions for the M205A and M206A mutations, and from van der Waals and covalent energy terms for M213A. Interestingly, while free energy contributions from a majority of residues neighboring the mutation sites tend to stabilize the wild type, there are a few residues stabilizing the mutant side chains. Our results show that this approach to free energy calculation can be very useful for understanding the detailed mechanism of human prion protein stability.