Minimizing NVH and friction-induced power losses is becoming paramount in the design of geared transmissions. The aim of this paper is to present an automatic methodology to explore Pareto-optimal designs of bevel gears when minimization of noise and frictional losses is essential. In the first part, a semi-empirical model to estimate frictional power losses under elasto-hydrodynamic lubrication is described. The model has been validated against experimental data available in the literature in previous works by the authors. The efficiency calculation is coupled with a state-of-the-art loaded tooth contact analysis (LTCA) tool to obtain accurate predictions of the instantaneous load shared by the mating tooth pairs during the meshing cycle. In the second part, an automatic framework based on multi-objective optimization (MOO) is presented where the tooth micro-geometry is systematically designed. The design variables are represented by few coefficients of a polynomial basis that embodies the tooth flank ease-off topography. To ensure manufacturability, the polynomial modifications are projected onto the feasible set of the machine-tool envelopes. This step is achieved through a state-of-the-art identification algorithm that the authors have developed in previous work. Frictional losses are estimated with the aforementioned model, whereas the NVH level is measured by the loaded transmission error (LTE), directly available from the simulation tool. The maximum contact pressures are limited by the material properties, thus proper nonlinear constraints are prescribed. Application to a test case involving the design of a spiral bevel gearset reveals that the methodology presented allows the designer to obtain Pareto-optimal solutions in a systematic and automatic manner.