The role of MD simulations and free energy calculation methods such as those employed here may provide insight not readily accessible to conventional prediction methods. One particular insight, namely the pK a values of coupled residues, requires careful consideration of the role of nearby protonatable residues. Moreover, given that these residues are often found at the active site -frequently the target of engineered therapeutics -the relevance of this problem extends beyond basic research.