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

기술 요약: 틸라피아의 제초제 생물농축에 관한 간 구획 모델

문제 제id (Problem Statement)
생물농축계수(BCF)는 전통적으로 아가미 흡수, 전신 분포, 간 대사를 하나의 정상 상태 비율로 혼합하여 다루는 집계된 양으로 취급된다. 이러한 집계된 척도는 규제 심사에는 유용하지만, 독립적인 생리학적 데이터와 교차 검증하거나 화학 물질 혼합물에서의 대사 상호작용과 같은 기저 메커니즘을 분해하여 밝혀낼 수 없다. 이러한 한계는 사탕수수 농업 지역에서 나일 틸라피아(Oreochromis niloticus)가 네 가지 제초제(아메트린, 테부티우론, 헥사지논, 디우론)에 공동 노출되는 상황에서 특히 두드러진다. 저자들은 근육 조직 샘플링만으로 도출된 집계 BCF가 얼마나 생리학적으로 해석 가능한 속도 상수(rate constants)로 분해될 수 있는지 결정하고, 이 분해가 표준 분배 계수 전용 모델에서는 보이지 않는 혼합물 의존적 효과를 드러내는지 평가하고자 한다.

방법론 (Methodology)
본 연구는 네 가지 사탕수수 제초제에 대한 기존의 생물농축 데이터(Jonsson et al., 2019)를 3구획 운동론적 모델(물 $\rightarrow$ 간 $\rightarrow$ 근육)을 사용하여 재분석한다. 기존의 2구획(물 $\rightarrow$ 근육) 공식과 달리, 이 모델은 간을 대사의 장소로 명시적으로 격리한다.

1. 모델 구조: 시스템은 간과 근육 농도에 대한 미분 방정식에 의해 지배된다. 간의 유사 정상 상태 가정(간의 회전율 $\ll$ 근육의 역학) 하에, 시스템은 두 개의 식별 가능한 유효 파라미터(유효 흡수 $K_1$, 유효 제거 $K_2$)를 가진 단일 2차 미분 방정식으로 축소된다.
2. 보정 (Calibration): $K_1$과 $K_2$는 9가지 화합물 $\times$ 혼합물 조합에 걸친 원래 연구의 정상 상태 근육 농도(고원 상태 평균값)를 통해 비선형 최소 제곱법으로 피팅되었다. 저자들은 이 중 7개 조합의 경우 데이터 포인트의 수가 피팅된 파라미터의 수와 동일하므로, 이는 검증이 아닌 보정임을 명시하였다.
3. 간 분해 (Hepatic Decomposition): 유효 파라미터를 생리학적 속도($k_{up}$, $k_{fm}$, $k_{mf}$, $k_{elim}$)로 분해하기 위해, 본 모델은 다음과 같은 외부의 독립적인 출처에 기반을 두었다:
   • 제거 ($k_{elim}$): QSAR를 사용한 US EPA EPI Suite(BCFBAF 모듈)로부터 도출되었으며, 디우론의 경우 특정 측정값이 사용되었다.
   • 흡수 ($k_{up}$): 두 가지 독립적인 추정치의 기하 평균을 취하여 합의된 값을 설정하였다: 분자 QSAR(Kuo & Di Toro)와 틸라피아 생리에 맞춰 스케일링된 측정된 아가미 추출 효율(McKim et al.)에 기반한 생리학적 추정치.
   • 분포: 나머지 속도 상수($k_{fm}$, $k_{mf}$)는 대수적으로 계산되었다.
4. 위험 평가: 원래 연구의 결정론적 위험 평가를 몬테카를로 시뮬레이션(라틴 하이퍼큐브 샘플링)으로 재구성하여, BCF와 소비자 체중의 불확실성을 전파하고 허용 일일 섭취량(ADI)을 초과할 확률을 계산하였다.
5. 검증: 파라미터 식별 가능성을 평가하고 실험 설계의 구조적 한계(예: 비단조적 농도 감소)를 확인하기 위해 214개의 개별 물고기 결정치(평균값이 아닌)를 대상으로 모델을 테스트하였다.

주요 결과 (Key Results)

• 모델 적합도 및 식별 가능성: 축약된 모델은 경험적 BCF를 잔차 없이 재현하였다. 그러나 일별 데이터는 샘플링이 고원 단계에서만 이루어졌기 때문에 9개 사례 중 7개에서 $K_1$과 $K_2$가 개별적으로 식별 불가능함을 확인시켜 주었다. 이 상수들을 분리하는 모델의 능력은 전술한 외부 앵커링(anchoring)에 전적으로 의존하였다.
• 생리학적 분해: 분해 결과, 아가미 흡수($k_{up}$)는 간 제거($k_{elim}$)보다 현저히 빠르고, 간 제거는 간-근육 분포($k_{fm}$)보다 빠르다. 이는 간이 근육에 도달하는 화학 물질의 비율을 제한하는 "강한 간 초회 통과(strong hepatic first-pass)" 메커니즘을 뒷받침한다.
• 외부 타당성 검토:
  • 테부티우론, 헥사지논, 디우론의 암시된 간-근육 분배 비율은 틸라피아 내 아트라진에 대한 독립적인 측정치와 10배 이내의 차이를 보였다.
  • 모델은 아메트린에 대해 실패하였는데, 관찰된 벤치마크보다 약 10배 높은 간-근육 비율을 예측하였다. 저자들은 아메트린의 간-근육 분포가 모델에서 가장 제약이 적은 요소이며, 이는 아마도 느린 간 제거 때문인 것으로 보고 있다.
• 혼합물 효과: 테부티우론의 유효 제거율($K_2$)은 공동 노출된 제초제의 수가 증가함에 따라 단조 감소하였다(이진 혼합물에서 4구획 혼합물로 갈 때 34% 감소). 이 패턴은 간 효소의 경쟁적 저해를 나타내며, 이는 원래의 집계 BCF로는 탐지할 수 없는 메커니즘이다.
• 위험 평가: 몬테카를로 시뮬레이션은 원래 연구의 결론을 확인해주었다. 극단적인 소비 시나리오에서도 ADI 초과 확률은 무시할 수 있는 수준(1,000회 시뮬레이션 중 0%)이었다. BCF는 최대 일일 섭취량(MDI) 계산에서 불확실성의 주요 동인이었다.

의의 및 주장 (Significance and Claims)
본 논문은 집계된 BCF를 메커니즘적인 간 구획 모델로 재구성함으로써 다음과 같은 것이 가능하다고 주장한다:

1. 분해: BCF를 독립적인 데이터베이스에 고정함으로써, 이를 불투명한 피팅 숫자로 남겨두는 대신 생리학적으로 해석 가능한 속도 상수($k_{up}$, $k_{elim}$ 등)로 분해할 수 있다.
2. 혼합물 민감도: BCF가 화합물의 고정된 특성이 아니라 혼합물 조성에 따라 변하는 특성임을 감지할 수 있으며, 이는 분배 계수 전용 QSAR가 포착하지 못하는 의존성이다.
3. 가설 생성: 공동 노출 시나리오에서 대사 상호작용(예: 경쟁적 저해)을 가설화할 수 있는 구조적 근거를 제공하며, 이는 원래의 경험적 접근 방식으로는 추측에 그쳤던 부분이다.
4. 방법론적 투명성: (샘플링이 고원 단계와 깨끗한 물로의 전이 후에만 이루어지는) 표준 2단계 생물농축 프로토콜의 구조적 한계를 강조한다. 이러한 프로토콜은 외부 제약 없이는 흡수 및 제거율을 식별 불가능하게 만들며, 탈출(depuration) 단계 동안 모든 공동 노출을 제거함으로써 혼합물 상호작용을 은폐할 수 있다.

저자들은 이 모델이 1차 운동론의 검증이라기보다는 기존 데이터에 대한 보정임을 강조하며, 분해가 직접적인 간 잔류물 측정보다는 독립적인 추정 방법들의 수렴에 의존한다는 점을 들어 신중한 입장을 유지한다. 이 연구는 오직 식용 조직(근육)만을 샘플링했을 때 집계된 생물농축 데이터를 메커니즘적 구성 요소로 분해할 수 있음을 보여주는 개념 증명 역할을 한다.