A Hepatic Compartmental Model of Herbicide Bioconcentration in Nile Tilapia (Oreochromis niloticus)

Resumo Técnico: Um Modelo de Compartimentos Hepáticos da Bioconcentração de Herbicidas em Tilápia do Nilo

Definição do Problema
O fator de bioconcentração (BCF) é tradicionalmente tratado como uma quantidade agregada, fundindo a captação branquial, a distribuição sistêmica e a biotransformação hepática em uma única razão de estado estacionário. Embora útil para triagem regulatória, essa medida agregada não pode ser confrontada com dados fisiológicos independentes ou decomposta para revelar mecanismos subjacentes, como interações metabólicas em misturas químicas. Esta limitação é particularmente aguda na agricultura de cana-de-açúcar, onde a tilápia do Nilo (Oreochromis niloticus) é rotineiramente exposta conjuntamente a misturas de herbicidas (ametrina, tebutiurom, hexazinona e diuron). Os autores buscam determinar até que ponto um BCF agregado, derivado apenas da amostragem de tecido muscular, pode ser resolvido em constantes de taxa fisiologicamente interpretáveis, e avaliar se essa decomposição revela efeitos dependentes de misturas invisíveis aos modelos padrão baseados apenas em coeficientes de partição.

Metodologia
O estudo reanalisa dados de bioconcentração publicados para quatro herbicidas de cana-de-açúcar em tilápia do Nilo (Jonsson et al., 2019) utilizando um modelo cinético de três compartimentos: Água $\rightarrow$ Fígado $\rightarrow$ Músculo. Diferente da formulação original de dois compartimentos (Água $\rightarrow$ Músculo), este modelo isola explicitamente o fígado como o sítio de biotransformação.

1. Estrutura do Modelo: O sistema é governado por equações diferenciais para as concentrações no fígado e no músculo. Sob uma suposição de pseudo-estado estacionário hepático (o turnover do fígado $\ll$ a dinâmica do músculo), o sistema reduz-se a uma única EDO de segunda ordem para a concentração muscular com dois parâmetros efetivos identificáveis: captação efetiva ($K_1$) e eliminação efetiva ($K_2$).
2. Calibração: $K_1$ e $K_2$ foram ajustados via mínimos quadrados não lineares às concentrações musculares de estado estacionário (médias de platô) do estudo original através de nove combinações de composto $\times$ mistura. Os autores observam explicitamente que, para sete dessas combinações, o número de pontos de dados é igual ao número de parâmetros ajustados, tornando isso uma calibração em vez de uma validação independente.
3. Decomposição Hepática: Para decompor os parâmetros efetivos em taxas fisiológicas ($k_{up}$, $k_{fm}$, $k_{mf}$, $k_{elim}$), o modelo foi ancorado em fontes externas e independentes:
   • Eliminação ($k_{elim}$): Derivada do módulo BCFBAF do US EPA EPI Suite (usando QSARs), com um valor medido específico utilizado para o diuron.
   • Captação ($k_{up}$): Um valor de consenso foi estabelecido tomando a média geométrica de duas estimativas independentes: um QSAR molecular (Kuo & Di Toro) e uma estimativa fisiológica baseada na eficiência de extração branquial medida (McKim et al.), escalonada para a fisiologia da tilápia.
   • Distribuição: As demais constantes de taxa ($k_{fm}$, $k_{mf}$) foram resolvidas algebricamente.
4. Avaliação de Risco: A avaliação de risco determinística do estudo original foi reformulada como uma simulação de Monte Carlo (Amostragem de Hipercubo Latino) para propagar a incerteza no BCF e no peso corporal do consumidor, calculando a probabilidade de exceder a Ingestão Diária Aceitável (IDA).
5. Verificação: O modelo foi testado contra 214 determinações individuais de peixes (não apenas médias) para avaliar a identificabilidade dos parâmetros e verificar limitações estruturais no desenho experimental (ex: declínios de concentração não monotônicos).

Principais Resultados

• Ajuste do Modelo e Identificabilidade: O modelo reduzido reproduziu os BCFs empíricos com resíduos negligenciáveis. No entanto, os dados resolvidos por dia confirmaram que $K_1$ e $K_2$ são individualmente inidentificáveis em sete dos nove casos porque a amostragem ocorreu apenas na fase de platô. A capacidade do modelo de separar essas constantes dependeu inteiramente do ancoramento externo descrito acima.
• Decomposição Fisiológica: A decomposição revelou que a captação branquial ($k_{up}$) é significativamente mais rápida que o clareamento hepático ($k_{elim}$), que por sua vez é mais rápido que a distribuição fígado-músculo ($k_{fm}$). Isso apoia um mecanismo de "primeira passagem hepática forte", onde o fígado limita a fração de substância química que atinge o músculo.
• Verificações de Plausibilidade Externa:
  • As razões implícitas de partição fígado-músculo para tebutiurom, hexazinona e diuron concordaram dentro de uma ordem de magnitude com medições independentes para atrazina em tilápia.
  • O modelo falhou para a ametrina, prevendo uma razão fígado-músculo ~10 vezes maior do que os benchmarks observados. Os autores identificam a distribuição fígado-músculo da ametrina como o elemento menos restrito do modelo, provavelmente devido à sua lenta eliminação hepática.
• Efeitos de Mistura: A taxa de eliminação efetiva ($K_2$) para o tebutiurom diminuiu monotonicamente conforme o número de herbicidas co-expostos aumentava (uma redução de 3 de 34% de misturas binárias para quaternárias). Este padrão é consistente com a inibição competitiva de enzimas hepáticas, um mecanismo não detectável pelo BCF agregado original.
• Avaliação de Risco: A simulação de Monte Carlo confirmou a conclusão do estudo original: a probabilidade de exceder a IDA é negligenciável (0% em 1.000 simulações) mesmo sob cenários extremos de consumo. O BCF foi o principal driver de incerteza no cálculo da Ingestão Diária Máxima (IDM).

Significância e Alegações
O artigo afirma que reformular dados de BCF agregados em um modelo de compartimento hepático mecanístico permite:

1. Decomposição: Quebrar o BCF em constantes de taxa fisiologicamente interpretáveis ($k_{up}$, $k_{elim}$, etc.) ao ancorá-las em bases de dados independentes, em vez de deixá-las como números ajustados opacos.
2. Sensibilidade à Mistura: Detectar que o BCF não é uma propriedade fixa de um composto, mas varia com a composição da mistura, uma dependência que os QSARs baseados apenas em coeficiente de partição falham em capturar.
3. Geração de Hipóteses: Fornecer uma base estrutural para hipotetizar interações metabólicas (ex: inibição competitiva) em cenários de co-exposição, algo que a abordagem empírica original só poderia especular.
4. Transparência Metodológica: Destacar as limitações estruturais de protocolos padrão de bioconcentração de duas fases (amostragem apenas no platô e após transferência para água limpa), que tornam as taxas de captação e eliminação inidentificáveis sem restrições externas e podem mascarar interações de mistura ao remover todas as co-exposições durante a fase de depuração.

Os autores mantêm uma postura modesta, enfatizando que o modelo é uma calibração de dados existentes e não uma validação de cinética de primeira ordem, e que a decomposição depende da convergência de métodos de estimativa independentes em vez de medição direta de resíduos no fígado. O trabalho serve como uma prova de conceito para resolver dados de bioconcentração agregados em componentes mecanísticos quando apenas o tecido comestível (músculo) é amostrado.