Our overall aim was simple: demonstrate that insights into residue coupling can be resolved by means of alchemical free energy calculations. While such insights can be directly resolved via constant-pH molecular dynamics simulations -an exceedingly powerful method -we demonstrate with both toy examples and real protein systems that similar insights can be extracted via a different approach.