The present paper proposes a new mixed-fidelity method to optimize the shape of ships using genetic algorithms (GA) and potential flow codes to evaluate the hydrodynamics of variant hull forms, enhanced by a surrogate model based on an Artificial Neural Network (ANN) to account for viscous effects. The performance of the variant hull forms generated by the GA is evaluated for calm water resistance using potential flow methods which are quite fast when they run on modern computers. However, these methods do not take into account the viscous effects which are dominant in the stern region of the ship. Solvers of the Reynolds-Averaged Navier-Stokes Equations (RANS) should be used in this respect, which, however, are too time-consuming to be used for the evaluation of some hundreds of variants within the GA search. In this study, a RANS solver is used prior to the execution of the GA to train an ANN in modeling the effect of stern design geometrical parameters only. Potential flow results, accounting for the geometrical design parameters of the rest of the hull, are combined with the aforementioned trained meta-model for the final hull form evaluation. This work concentrates on the provision of a more reliable framework for the evaluation of hull form performance in calm water without a significant increase of the computing time.