한 줄 요약 — 경구 철 흡수는 비율이 아니라 세 인자의 곱이며, 그 세 번째 인자(장세포 배출능)는 직전 용량이 만들어낸 헵시딘이 정한다. 즉 약은 자신이 통과해야 할 문을 스스로 닫는다. 그리고 그 문 아래에서 헤모글로빈을 제한하는 것은 전혀 다른 병목(대식세포 재활용이 지배하는 골수 철 공급) 이므로, 흡수량과 Hb 반응은 고용량에서 서로 갈라진다.

철결핍성 빈혈은 "철이 부족하니 철을 주면 된다"는 질환처럼 보인다. 그런데 임상 현실은 그렇지 않다. 같은 총량을 매일 나눠 주는 것과 격일로 주는 것이 다른 결과를 내고, 정맥 철제를 두 배로 늘려도 Hb가 두 배 빨리 오르지 않으며, 염증이 있으면 철이 흡수되어도 Hb가 오르지 않는다. 이 세 가지 "이상한 일"은 모두 철 흡수가 조절되는 관문(gate)이라는 사실과 Hb 상승을 제한하는 병목이 그 관문과 다른 곳에 있다는 사실에서 나온다. 이 모델은 그 두 제약을 각각 방정식으로 쓴다.

흡수 플럭스 = A_LUM  ×  DMT1 용량  ×  FPN_ENT
               │          │            └── 장세포 배출능:
               │          │                직전 용량의 헵시딘이 결정
               │          └── 포화 (Km 11 mg)
               └── 장내 가용 철

골수는 하루 13.3 mg의 철을 쓰는데 그중 약 **85 %**가 노화 적혈구를 먹은 대식세포에서 재활용되어 온다. 장에서 오는 몫은 잘해야 하루 10 mg 안팎이다. 따라서 철 공급을 늘려도 Hb 상승 속도의 상한은 골수가 정한다.

두 기준 상태는 입력한 값이 아니라 모델의 평형해다. 정상 상태에서 출혈량 파라미터(VBLEED)만 6.375 mL/day(≈ 178 mL/주기)로 올리고 정상상태까지 적분하면 아래 IDA 상태가 나온다.

주목할 두 가지. 첫째, 망상적혈구가 빈혈에 비해 부적절하게 낮다(42 ×10⁹/L). EPO는 94 mIU/mL로 올라 있는데도 골수가 확장하지 못하기 때문이다 — 모델은 이를 F_EXP = TSAT/(TSAT + 12) 라는 게이트로 표현한다(TSAT 6.5 %에서 EPO 확장의 35 %만 실현). 둘째, 흡수량이 2.53 mg/day로 정상의 2배로 올라 있다. 결핍이 헵시딘을 억제해 관문을 연 결과다. 그런데도 병이 유지되는 것은 손실이 더 크기 때문이다. IDA 평형은 "Hb가 떨어져 출혈 철 손실이 줄어들어 마침내 흡수량과 같아지는 지점"이다.

경구 elemental iron 60 mg 1회 (IDA 기저 상태):

분획흡수율 22.0 % (13.2 mg / 60 mg) — 철결핍 여성의 관찰값 20–22 %와 일치.

핵심은 시계가 두 개라는 점이다. 헵시딘 자체는 반감기 2.5 h로 빨리 사라지지만, 장세포의 배출능은 훨씬 느린 시계로 돌아온다. 실효 회복 시간상수는 1/(KSYN_FPE + KDEG_FPE·HEP^0.6)이며, IDA 기저 헵시딘에서 약 51시간(≈2일)이다. 그래서 혈중 헵시딘이 정상으로 돌아온 뒤에도 흡수 능력은 여전히 낮다. 1943년 Hahn이 "점막 차단(mucosal block)"이라 부른 현상이 이 두 시계의 차이다.

동일한 60 mg을 조건화 용량 뒤 여러 간격에 주고, 프로브를 주지 않은 대조 시뮬레이션과의 차이로 프로브 자체의 흡수량을 측정했다.

같은 날 4시간 뒤 두 번째 용량은 4분의 1을 잃는다(Moretti 2015가 관찰한 바로 그 현상). 24시간이 지나도 완전히 회복되지 않는다.

14일간 투여, 완전 순응 가정, 용량당 elemental iron 고정:

두 지표가 정반대로 정렬된다. 간격을 늘리면 분획흡수율은 오르고(13.9 → 20.9 %) 하루 흡수량은 떨어진다(16.72 → 4.48 mg/day). 즉 격일 투여는 효율과 내성의 최적점이지, 전달량의 최적점이 아니다. 삼킨 1 mg당 얻는 철은 격일이 매일의 1.14배지만, 하루에 실제로 들어오는 철은 매일 투여가 격일의 1.75배다.

같은 총량을 다르게 나눈 비교(각 1680 mg): 60 mg q12h는 234.0 mg를 흡수하고 120 mg q24h는 188.7 mg를 흡수한다. 여기서는 용량 포화(Km 11 mg)를 피하는 이득이 헵시딘 굴절기의 손해보다 크다. 이 둘의 균형이 용량·간격에 따라 뒤집히는 것이 이 모델의 비자명한 예측이다.

GI 증상이 순응도를 깎는 항(K_GI, EMAX_ADH)을 켜면 잦은 투여가 추가로 벌을 받는다. 완전 순응 대비 실제 흡수량:

간격이 길수록 손실이 작다. 임상 지침이 격일 투여를 권하는 이유는 흡수 효율과 내성이 같은 방향으로 움직이기 때문이며, 전달량은 그 반대 방향으로 움직인다는 사실을 이 표가 같이 보여준다.

읽을 점 몇 가지.

• 정맥 철제의 이득은 크기가 아니라 속도다. FCM 1000 mg는 12주 ΔHb에서 경구 65 mg 1일2회(+3.23)와 사실상 같은 +3.20을 낸다. 다른 것은 도달 시각이다: +2 g/dL까지 21.5일 대 37일. 망상적혈구 정점도 74 vs 60 ×10⁹/L.
• 출혈 교정 단독으로는 부족하고(+0.81), 철 보충 단독으로도 부족하다(+2.59). 둘을 합치면 +3.36으로 정맥 철제에 근접한다 — 원인 교정이 용량 증량보다 값싼 개입이라는 점을 모델이 정량화한다.
• 염증에서는 경구가 사실상 실패한다(+0.45). 정맥도 부분적 반응(+1.38)에 그친다. 시나리오 12에서 페리틴은 50.3까지 오르는데 TSAT은 12.3 %에 머문다 — 기능적 철결핍의 지표 해리가 가정이 아니라 기전에서 나온다.

500 mg → 3000 mg (6배)에서 최대 주간 Hb 기울기는 36 %만 늘고(0.72 → 0.98), 1주차 반응은 11 %만 늘지만(0.35 → 0.39), 저장철은 24배 늘어난다 (72 → 1720 mg). 시계를 쥐고 있는 것은 철 공급이 아니라 골수다. 이것이 "총 철결핍량(Ganzoni)을 채우는 것"과 "빈혈을 빨리 교정하는 것"이 서로 다른 목표인 이유다.

IL-6 40 pg/mL에서 장은 여전히 하루 4.5 mg을 흡수하고 순 철 수지도 +1.52 mg/day로 양수인데, Hb는 12주 동안 전혀 오르지 않는다(−0.07). 철이 없어서가 아니라 철이 잘못된 곳에 있기 때문이다: 대식세포 ferroportin이 닫혀 있어 저장고에 갇히고(페리틴 55), 동시에 IL-6가 EPO 생산·전구세포·적혈구 수명을 직접 억제한다. 그래서 "철을 더 주는 것"은 이 상황에서 답이 아니다 — 염증을 끄는 것이 답이다.

TMPRSS6를 0.30으로 낮춘 별도의 평형 상태(출혈 없음):

병소가 "저장철 고갈"이 아니라 "헵시딘이 부적절하게 높음"이므로 페리틴은 정상인데 경구 철은 듣지 않는다. 12주간 292 mg을 흡수했는데도 Hb는 0.32 g/dL만 올랐다. 페리틴이 정상인 미세적혈구 빈혈에서 경구 철이 듣지 않을 때 헵시딘을 재보라는 임상 권고가 여기서 유도된다.

한계: 이 모델의 IRIDA 기저 Hb(12.78)는 실제 IRIDA의 중등도 빈혈보다 경증이다. 경구 불응성이라는 정성적 결론은 유지되나 중증도는 과소평가되어 있다.

정맥 철제의 탄수화물 껍질만 바꾼 비교. 혈액학적 결과는 동일(시나리오 8 vs 10: ΔHb 둘 다 +3.20 g/dL)하지만 인(phosphate) 대사는 전혀 다르다.

Carboxymaltose는 FGF23 절단을 억제해 intact FGF23을 올리고(→ NaPi-2a/2c 억제로 인 소실, CYP27B1 억제로 1,25D 감소, 이차적 PTH 상승), derisomaltose는 그렇지 않다. 효능 지표로는 두 약을 구분할 수 없다 — 기전 모델만이 이 비대칭을 보여줄 수 있고, 이것이 QSP가 임상시험 요약표보다 더 말해주는 지점이다. FCM 750 mg를 두 번 나눠 주면 총량은 1500 mg로 늘고 저인산혈증 일수는 4.2 → 9.2일로 2배 이상 늘어난다.

경구 치료는 12주에 Hb 결손의 53 %를 메우면서 조직 철 결손은 18 %만 메운다. 정맥 철제는 반대로 조직 철을 먼저(105 %) 채운다. 골수가 혈중 철을 우선적으로 가져가기 때문이다. 빈혈이 교정되었는데도 피로·하지불안이 남는 환자와 빈혈이 없는데도 정맥 철제로 피로가 개선되는 환자(Krayenbuehl 2011, FAIR-HF)가 같은 구조에서 나온다.

37개 상태변수:

세 가지 구조적 선택이 이 모델의 성격을 결정한다.

1. 세포수와 Hb 질량을 별도의 두 사슬로 추적한다. 그래야 MCH·CHr가 계산에서 나오고, "적은 수의 정상 세포"와 "많은 수의 저색소 세포"를 구분할 수 있다.
2. Ferroportin에 두 개의 시계를 준다. 장세포(KSYN_FPE 0.004/h, 느림)는 점막 차단의 24–48시간 기억을, 대식세포(KSYN_FPR 0.035/h, 빠름)는 염증성 저철혈증의 수시간 반응을 만든다. 하나의 시계로는 둘 다 재현할 수 없다.
3. EPO의 골수 확장 효과에 철 게이트를 씌운다(F_EXP). 이것이 없으면 모델은 빈혈에 EPO로 반응해 세포를 무한히 만들어 내며 망상적혈구가 실제보다 높아진다.

BODY_FE(t) = BODY_FE(0) + CUM_IV + CUM_ABS − CUM_LOSS를 12주 시뮬레이션에서 검증했다. 오차 ~1×10⁻¹² % (배정도 정밀도 한계). 장세포 풀(A_ENT, A_EFT)은 의도적으로 체내 철에서 제외했다 — 정점막으로 들어왔다가 세포와 함께 탈락한 철은 체내로 들어온 적이 없으므로, 이를 세면 양변이 함께 부풀려진다.

모든 결과는 두 개의 독립 구현으로 계산해 대조했다: mrgsolve(C++/LSODA)와 ida_reference_model.py의 순수 파이썬 고정스텝 RK4 구현이다. 두 구현은 코드를 공유하지 않으므로, 일치한다는 사실은 결과가 특정 솔버의 거동이 아니라 방정식 자체에서 나온다는 증거가 된다. 예시(단일 60 mg): 혈청철 정점 224.2 vs 224.2 µg/dL, 헵시딘 정점 1.243 vs 1.243 ng/mL, 분획흡수율 21.97 vs 21.97 %, 프로브 +24 h 91.5 vs 91.5 %, 14일 q24h 흡수량 147.70 vs 147.70 mg.

# ODE 모델 + 15개 시나리오 + 검증 노트 전체 실행
Rscript -e 'options(ida.report=TRUE); source("ida_mrgsolve_model.R")'

# 인터랙티브 대시보드 (8탭)
Rscript -e 'shiny::runApp("ida_shiny_app.R")'

# R 없이 재현 (표준 라이브러리만 사용)
python3 ida_reference_model.py

# 지도 재렌더링 (newrank 필요)
dot -Tsvg ida_qsp_model.dot -o ida_qsp_model.svg
dot -Tpng -Gdpi=150 ida_qsp_model.dot -o ida_qsp_model.png

다른 환자로 옮길 때는 BV_L, PV_L, VBLEED, DIET, IL6_IN을 바꾼 뒤 반드시 재평형시켜야 한다(무투여로 2–3년 적분). 기본 초기값은 위 표의 IDA 평형 상태이며, 다른 파라미터에서는 그 값이 평형이 아니다. Shiny 앱의 "기저상태 재평형" 체크박스가 이 작업을 수행한다.

문헌 대조와 명시적 한계는 ida_references.md의 마지막 두 절에 정리했다. 요약하면,

• 재현된 것: 단회 분획흡수율(22.0 % vs 관찰 20–22 %), 24시간 지속되는 헵시딘 상승(1.36배), 같은 날 재투여 흡수 저하(75.4 %), 망상적혈구 정점 시기(12–14일), FCM 저인산혈증의 크기·시기와 derisomaltose에서의 부재, 정상/IDA 기준치 전반.
• 과소평가된 것: 격일 대 매일의 분획흡수율 비가 모델 0.87 vs 관찰 0.75 — 본 모델은 굴절기의 벌칙을 보수적으로 추정한다. 헵시딘 순환 반감기를 2.5 h로 고정한 채 장세포 회복 시계만으로 24시간 기억을 만들면 깊이와 지속성이 상충하기 때문이다.
• 그 외: 혈청철 정점이 다소 늦고(5.8 h vs 2–5 h), IRIDA 표현형이 경증이며, GI–순응도 항은 반정량적이고, 60 kg 여성 단일 체격으로 보정되어 있다.

⚠️ 교육 및 연구 목적의 모델입니다. 실제 진료·처방·규제 제출에 사용할 수 없습니다.