Electrostatic actuators, actuators actuated by electrostatic forces, are now widely used as sensors and switches, especially in Micro-Electro-Mechanical Systems (MEMS). Among different kinds of electrostatic actuators, the comb drive type is one of the most popular because it has a relatively large range of displacement. In design problems for electrostatic actuators, the driving force profile is of primary engineering importance. In this paper, we develop a structural optimization method for comb drive electrostatic actuators that achieves prescribed driving force profiles, based on a level set-based shape optimization method that provides optimal configurations with clear boundaries, solutions that are valid in an engineering sense. Accurate calculation of the electrostatic forces that occur on the structural boundaries during optimization is important for developing actuators that operate with prescribed driving forces. In the conventional level set-based shape optimization methods, inaccuracies in the calculation of these electrostatic forces occur because the structural boundaries are seldom aligned with the finite element method (FEM)