Recently, Model Predictive Control (MPC) for buildings has undergone an intensive research. Usually, according to the international standards, a static range for the air temperature represents the thermal comfort which is being kept making use of MPC while minimizing the energy consumption. On contrary, this paper deals with the optimization of the trade-off between energy consumption and Predicted Mean Vote (PMV) index which, opposed to the static temperature range, describes user comfort directly. PMV index is a nonlinear function of various quantities, which makes the problem more difficult to solve. The paper will show the main differences in MPC problem formulation, propose a tractable approximation strategy and compare the control performance both to the conventional and typical predictive control strategies. The approximation of PMV computation will be shown to be sufficiently precise and moreover, such a formulation keeps the MPC optimization problem convex. Finally, it will be shown that the proposed PMV based optimal control problem formulation shifts the savings potential of typical MPC by additional 10% while keeping the comfort at a desired level.