The present paper deals with shape optimization in discretized two-dimensional (2D) contact problems with Coulomb friction, where the coefficient of friction is assumed to depend on the unknown solution. Discretization of the continuous state problem leads to a system of finite-dimensional implicit variational inequalities, parametrized by the so-called design variable, that determines the shape of the underlying domain. It is shown that if the coefficient of friction is Lipschitz and sufficiently small in the C-0,C-1-norm, then the discrete state problems are uniquely solvable for all admissible values of the design variable (the admissible set is assumed to be compact), and the state variables are Lipschitzian functions of the design variable. This facilitates the numerical solution of the discretized shape optimization problem by the so-called implicit programming approach. Our main results concern sensitivity analysis, which is based on the well-developed generalized differential calculus of B. Mordukhovich and generalizes some of the results obtained in this context so far. The derived subgradient information is then combined with the bundle trust method to compute several model examples, demonstrating the applicability and efficiency of the presented approach.