A hydro-kinetic scheme for water-like fluids, based on a lattice version of the Boltzmann equation, is presented and applied to the popular TIP3P model of liquid water. By proceeding in much larger steps than molecular dynamics, the scheme proves to be very effective in attaining global minima of classical pair potentials with directional and radial interactions, as confirmed by further simulations using the three-dimensional Ben-Naim water potential. The scheme is shown to reproduce the propensity of water to form nearly four hydrogen bonds per molecule, as well as their statistical distribution in the presence of thermal fluctuations, at a linear cost of computational time with the system size.