There is a need for an efficient and accurate computational method to identify the effects of single- and multiple-residue mutations on the stability and reactivity of proteins. Such a method should ideally be consistent and yet applicable in a widespread manner, i.e. it should be applied to various proteins under the same parameter settings, and have good predictive power for all of them.