The level set method is used for shape optimization of the energy functional for the Signorini problem. The boundary variations technique is used in order to derive the shape gradients of the energy functional. The conical differentiability of solutions with respect to the boundary variations is exploited. The topology modifications during the optimization process are identified by means of an asymptotic analysis. The topological derivatives of the energy shape functional are employed for the topology variations in the form of small holes. The derivation of topological derivatives is performed within the framework proposed in (Sokołowski and Żochowski, 2003). Numerical results confirm that the method is efficient and gives better results compared with the classical shape optimization techniques.