We consider optimal control problems for the bidomain equations of cardiacelectrophysiology together with two-variable ionic models, e.g. theRogers–McCulloch model. After ensuring the existence of global minimizers, we provide arigorous proof for the system of first-order necessary optimality conditions. The proof isbased on a stability estimate for the primal equations and an existence theorem for weaksolutions of the adjoint system.