Un sistema de ecuaciones cinéticas que describe los cambios en la concentración de especies reactivas de oxígeno (ROS) en soluciones acuosas de proteínas se obtuvo a partir del análisis de reacciones químicas que involucran oxígeno singlete. Aplicando la condición de estacionariedad de los productos intermedios al sistema, determinamos la dependencia funcional de la concentración de peróxido de hidrógeno en la concentración de proteína bajo la acción de radiación térmica y láser. Se encontró una solución analítica aproximada al sistema no lineal de ecuaciones diferenciales que definen la dinámica de la concentración de ROS. Para las soluciones acuosas de albúmina sérica bovina (BSA) y globulina gamma bovina (BGG), se determinaron los órdenes y constantes de velocidad de las reacciones que describen las conversiones de ROS minimizando la suma de desviaciones al cuadrado de las funciones encontradas al resolver tanto los problemas estáticos como dinámicos de las dependencias medidas experimentalmente. Al resolver el problema de optimización, se utilizó el algoritmo de Levenberg-Marquardt.