Computation of the pK of liquid water using coordination constraints