As espécies reativas de oxigénio e azoto desempenham um papel importante em certas funções fisiológicas, como na inflamação, mas em situações de sobreprodução ou desregulação podem ter vários efeitos lesivos para as biomoléculas. O objetivo foi construir um modelo matemático de equações diferenciais ordinárias (ODEs) que descrevesse a rede de reações químicas entre as diferentes espécies reativas e as biomoléculas. Para tal: a) foram recolhidas da literatura e utilizadas as constantes cinéticas das várias reações para escrever as equações cinéticas de todas as reações que constituíram o sistema estudado; b) com base nas várias equações cinéticas foram construídas as equações diferenciais que descreviam a evolução temporal de cada espécie química no sistema, recorrendo à abordagem do tipo Continuously Stirred Batch Reactor (CSBR); c) otimizou-se os parâmetros do modelo de modo a obter um sistema em estado estacionário que mimetizou as condições observadas na circulação sistémica; d) foi feita a análise do modelo com recurso à análise de fluxos e de coeficientes de controlo, que identificaram quais os principais pontos de controlo do sistema. Com estes passos obteve-se um modelo matemático robusto que permitiu a simulação de situações fisiológicas como a desregulação dos sistemas de produção de espécies reativas ou o efeito de antioxidantes sobre o nível global de espécies reativas e sobre o nível de lesões induzidas às biomoléculas: ▪ As concentrações de estado estacionário das espécies tamponantes do sangue bem como da água, do hidrónio e do anião hidróxido mantiveram-se constantes em todas as variações; ▪ As variações do óxido nítrico e do superóxido ocorreram exatamente de acordo com o que se esperava; ▪ A variação da concentração do anião nitrato correspondeu à variação da concentração do peroxinitrito, indicando a composição do mesmo; ▪ Todas as espécies reativas tiveram a variação das suas concentrações como esperado em relação da variação das concentrações de estado estacionário dos antioxidantes. Contudo, não foi possível obter um modelo que descreva adequadamente o sistema de tamponamento do sangue, tendo-se identificado a origem desse erro no excessivo detalhe utilizado para descrever essas reações.