Water distribution networks are high energy and low efficiency systems, where water pressure is frequently reduced by dissipation valves to limit leakage. The dissipation produced by the valves can be converted to energy production to increase the efficiency and reduce the energy impact of networks. If valves are replaced by turbines or pumps as turbines (PATs), they can both reduce pressure and produce energy. This study focuses on the optimal location of PATs within a water distribution network in order to both produce energy and reduce leakage. A new optimization model is developed consisting of several linear and non-linear constraints and a newly proposed objective function, where the turbine installation costs as well as the energy production and the economic saving due to the reduction of leakage can be accounted all together. The case study shows that the application of the mathematical model to a synthetic network ensures better results, in terms of both energy production and water saving, in comparison to other procedures.