Optimizing the hydrogen‐bond network in Poisson–Boltzmann equation‐based pKa calculations