In this study, the optimal seismic design methodology that can consider the structural integrity of not only the piping systems but also elasto-plastic supporting devices is developed. This methodology employs a genetic algorithm and can search the optimal conditions such as the supporting location, capacity and stiffness of the supporting devices. Here, a lead extrusion damper is treated as a typical elasto-plastic damper. Four types of evaluation functions are considered. It is found that the proposed optimal seismic design methodology is very effective and can be applied to the actual seismic design for piping systems supported by elasto-plastic dampers. The effectiveness of the evaluation functions is also clarified.