Although intrinsically approximated, classical three-site rigid water models are widely used in molecular dynamics (MD) simulations. One key requirement is that such simplified models can capture fairly well the properties of water even relying on a reduced number of parameters. In such a representation, the interaction potential is centered on three sites, each of which corresponds to one of the atoms in the water molecule (O, H, H). Early versions of these models, including, e.g., TIP3P and SPC, were originally parameterized to accurately reproduce basic thermodynamic properties, e.g., density and enthalpy of vaporization under standard conditions. Despite their age, these models continue to be extensively utilized in classical MD simulations, and most general-purpose forcefields are parametrized on them. With the increase in computing power, it has become possible to perform high-throughput parameterization, often in an automatic fashion, by considering a large set of experimental observables under different conditions as the reference data to fit.