Polymer ionization differs from that for their monomeric counterparts due to intramolecular correlations. Such effects are conventionally described in terms of the site-binding model that accounts for short-range interactions between neighboring sites. With an apparent equilibrium constant for each ionizable group and the nearest-neighbor energy as adjustable parameters, the site-binding method is useful to correlate experimental titration curves when the site-site interactions are insignificant at long ranges. This work aims to describe the electrostatic behavior of weak polyelectrolytes in aqueous solutions on the basis of the intrinsic equilibrium constants of the individual ionizable groups and solution conditions underlying the thermodynamic non-ideality. A molecular thermodynamic model is proposed for the protonation of weak polyelectrolytes by incorporating classical density functional theory into the site-binding model to account for the effects of the local ionic environment on both inter-chain and intra-chain correlations. By an extensive comparison of theoretical predictions with experimental titration curves, we demonstrate that the thermodynamic model is able to quantify the ionization behavior of weak polyelectrolytes over a broad range of molecular architectures and solution conditions.