<p>This study explores an optimal control problem of the monodomain model, a coupled system of a reaction-diffusion partial differential equation and an ordinary differential equation, employing the concept of virtual control to effectively handle both control and state constraints. The state constraints are addressed through a regularization technique based on virtual control, enabling a rigorous analysis of the existence of optimal control and the derivation of the first-order optimality system for the regularized problem. We provide a detailed estimation of the regularization error in terms of the regularization parameter, facilitating the analysis of the convergence of the optimal solution of the regularized problem to that of the original problem. Finally, numerical results are presented using the semi-smooth Newton method for optimization and the finite element method for discretization.</p>