Phenol is an industrial pollutant and its removal from industrial wastewaters is of great importance. In order to design optimised phenol removal procedures by using horseradish peroxidase-based systems, there are some points that have to be dealt with. One of the most important issues is the need for reliable kinetics as this is one of the difficulties found during process scale-up. Although simplified kinetics can be used for limited ranges of operating conditions, they are not usually reliable for the description of varying process conditions. The present work describes the implementation of a kinetic model, based on a mechanism, for the co-oxidation of phenol and 4-aminoantipyrine (Am-NH2), which is used as a chromogen agent, with hydrogen peroxide as the oxidant. The model covers not only the variation of the concentrations of all the species involved, but also the effect of temperature in the reaction. The estimation of kinetic rate constants and activation energies for the various steps in the mechanism is performed with a single optimisation procedure, and all the experimental results are described using a unique set of parameters, which, thus, is valid over an extended range of operating conditions. The mechanism allowed the determination of a reliable kinetic model which is appropriate for the range of experimental conditions used. The computational model was also tested with an independent set of experiments with different conditions from the ones for which the parameters were estimated.