Correct protonation states and relevant waters = better computational simulations?