An ever-growing population together with globally depleting water resources pose immense stresses for water supply systems. Desalination technologies can reduce these stresses by generating fresh water from saline water sources. Reverse osmosis (RO), as the industry leading desalination technology, typically involves a complex network of membrane modules that separate unwanted particles from water. The optimal design and operation of these complex RO systems can be computationally expensive. In this work, we present a modeling and optimization strategy for addressing the optimal operation of an industrial-scale RO plant. We employ a feed-forward artificial neural network (ANN) surrogate modeling representation with rectified linear units as activation functions to capture the membrane behavior accurately. Several ANN set-ups and surrogate models are presented and evaluated, based on collected data from the H2Oaks RO desalination plant in South-Central Texas. The developed ANN is then transformed into a mixed-integer linear programming formulation for the purpose of minimizing energy consumption while maximizing water utilization. Trade-offs between the two competing objectives are visualized in a Pareto front, where indirect savings can be uncovered by comparing energy consumption for an array of water recoveries and feed flows.