This paper is concerned with the state-constrained optimal control of the threedimensional thermistor problem, a fully quasilinear coupled system of a parabolic and elliptic PDE with mixed boundary conditions. This system models the heating of a conducting material by means of direct current. Local existence, uniqueness and continuity for the state system are derived by employing maximal parabolic regularity in the fundamental theorem of Prüss. Global solutions are addressed, which includes analysis of the linearized state system via maximal parabolic regularity, and existence of optimal controls is shown if the temperature gradient is under control. The adjoint system involving measures is investigated using a duality argument. These results allow to derive first-order necessary conditions for the optimal control problem in form of a qualified optimality system. The theoretical findings are illustrated by numerical results.for all ξ ∈ X. Moreover, in the proof of Theorem 4.5, ∂ θ B(θ, u) was found to be a topological isomorphism between the spaces X and L s (J; W −1,q ∅ ) × X s and consequently ∂ θ B(θ, u) * is also a topological isomorphism between L s (J; W 1,q ) × X s and X . In particular, for every f ∈ X there exists a unique p = p f ∈ L s (J; W 1,q ) × X s such that ∂ θ B(θ, u) * p = f. Hence, settinḡ