Neonatal opioid withdrawal syndrome / neonatal abstinence syndrome · NOWS / NAS

이 모델의 출발점은 다음 한 문장이다.

NOWS는 노르아드레날린성 질환인데, 우리는 이것을 µ-작용제로 치료하고, 그 둘 중 어느 것도 측정하지 않는 도구로 채점한다.

모든 것이 하나의 차이 위에 세워져 있다.

GAP = A − ITONE

  A      청반(locus coeruleus)의 신경적응 설정점(set-point).
         자궁 내에서 수 개월간 지속된 µ 수용체 점유가 써 놓은 값.
  ITONE  지금 이 순간 LC의 GIRK 통로에 걸려 있는 억제 긴장도.

GAP > 0 → LC 탈억제 → 노르아드레날린 유출 → 증상. GAP < 0 → 과진정. 이것은 "여유"가 아니다. Finnegan 채점표는 이것을 0점으로 채점하면서 수유를 억제한다.

신경적응 설정점의 사상(map)은 점유율에서 포화되어 있다:

A0 = AGAIN · FGA(GA) · EAD / (EA50A + EAD)     EA50A = 0.16

메타돈 40 mg/일에서도 태아 µ 점유 유효값 EAD는 이미 0.95를 넘는다. 따라서 4배 용량 차이가 만드는 A0 차이는 다음이 전부다:

A0는 1.6% 움직이고, 최고 점수는 0.2점 움직인다. 그러나 발현 시점은 1.5일 밀린다. 이유는 산술이다. 세척(washout)은 로그-선형이므로 용량비의 로그가 그대로 시간 이동이 된다: ln(4) × 16 h / 0.693 = 32 h.

산모 용량과 NOWS 중증도의 관계에 대해 문헌이 30년간 합의에 이르지 못한 이유가 여기서 나온다 — 관계는 중증도가 아니라 시점에 있다.

같은 아기, 같은 생물학, 규칙만 바꾼다:

교환율은 계산된다: 금단부담 1단위를 피하는 데 치료 0.39일.

두 가지가 따라 나온다. 첫째, 발표된 NOWS 치료기간의 3배 편차(4–28일)는 아기의 차이가 아니라 이유 규칙의 차이로 전부 설명된다 — 이 모델에서 생물학을 고정한 채 규칙만 바꾸면 21.5일에서 35.5일까지 나온다. 둘째, A를 안다고 해도 거의 아무것도 사지 못한다. 측정 불가능한 상태를 읽는 오라클 제어기(22.1일 / 부담 68.6)는 안정화 규칙(22.8일 / 부담 56.6)을 지배하지 못한다. 문제는 정보가 아니라, 아무도 명시적으로 정하지 않는 교환율이다.

ESC-NOW 시험은 평가 도구와 돌봄 환경을 동시에 바꿨으므로 그 효과를 귀속시킬 수 없다. 이 모델에서는 두 가지가 서로 다른 파라미터로 들어가므로 분리된다 (가상 코호트 n=200/셀):

읽는 법: 평가기준을 바꾸는 것만으로는 금단부담이 두 배가 된다(16.5 → 33.4). 돌봄을 바꾸는 것만으로는 부담이 줄어든다(16.5 → 10.8). 둘을 합치면 17.1 — 통상 돌봄과 사실상 같은 부담을, 5분의 1의 약물치료로 얻는다.

즉 ESC의 두 반쪽은 하나의 이익에 대한 독립적 기여가 아니다. 한쪽이 약물 감소를 사고, 다른 쪽이 그 값을 치른다. 돌봄 없이 채점표만 바꾸는 것은 이 모델에서 안전한 중재가 아니다. (단일 결정론적 환자에서도 같은 방향이다: ESC 12.0일 / 부담 101.0 대 Finnegan 22.8일 / 부담 56.6, 무치료 부담 192.8.)

코호트 자료에서는 항상 엉켜 있는 두 원인 — (a) 미성숙한 신경적응 능력, (b) 체중기준 용량에 비해 미성숙한 청소율 — 을 34주 아기에게 만삭 수준의 지속성 적응 풀을 강제로 부여해서 분리한다:

11.2일의 차이 중 10.0일(89%)이 낮은 적응이고 1.2일(11%)이 약동학이다. 그러나 약동학 쪽은 기간이 아니라 다른 곳에서 대가를 청구한다: 미숙아의 누적 과진정은 만삭의 48배다. 체중기준 용량은 미숙아의 치료를 짧게 만드는 것이 아니라 과진정시킨다 — 그리고 Finnegan 점수는 과진정을 채점하지 않는다.

페노바르비탈은 GIRK 합집합 밖(GABA-A, LC의 하류)에 그려져 있다. 순수 아편유사제 노출 아기에서:

서류상 4.5일 짧아지고, 실제 금단부담은 35% 늘고, 체중증가는 40%로 떨어진다. 벤조디아제핀 동시노출 아기에서는 이야기가 달라진다 (31.1일 → 23.0일, 8.1일 단축 대 순수 아편의 4.5일): 점수의 일부가 아편유사제 금단이 아닐 때만 페노바르비탈이 진짜로 유익하다.

같은 논리의 반대편이 클로니딘이다. α2 수용체는 µ 수용체와 같은 GIRK 통로를 연다. 그래서 클로니딘은 GAP을 실제로 메우고, 동시에 재유도 가능한 적응 풀 AT를 다시 몰지 않는다: 6 µg/kg/일에서 22.8 → 16.7일, 부담 56.6 → 30.5. 다만 공짜는 아니다 — 12 µg/kg/일에서 과진정 적분이 0.7 → 47.2로 뛴다.

적응 설정점은 두 개의 풀이다.

A = AD + AT
  AD  지속성 풀. 자궁 내에서 설정, t½ ≈ 12일, 출생 후 재유도 불가.
  AT  재유도 가능 풀(고전적 adenylyl cyclase superactivation),
      τ ≈ 60 h, 지금 몸에 있는 어떤 µ 작용제로도 다시 몰림 — 치료약 포함.

AT를 꺼 보면(ATMAX → 0) 치료가 22.8일에서 19.0일로 줄어든다. 치료 자체가 재유도하는 적응이 그 치료 기간에 3.8일(17%)을 더한다. 이것이 "점수를 더 세게 누르는 것"이 공짜가 아닌 이유이고, 목표점수를 3점으로 낮추면 치료가 42.5일로 늘어나면서 과진정 적분이 0.7에서 87.4로 폭증하는 이유다.

경쟁적 점유    θf = xf/(1+xf+xb),  θb = xb/(1+xf+xb)
               xf = UF/EC50_eff,   xb = UB/EC50b_eff
내성 이동      EC50_eff = EC50 · (1 + KSHIFT · A)          KSHIFT = 3.0
µ 효과         EMU = RMU · (1.00·θf + 0.80·θb)
적응 구동      EAD = RMU · (1.00·θf + 0.80·0.40·θb)        ← 부프레노르핀은 두 번 등장
α2 효과        EA2 = 0.34 · Cclon / (1.6·(1+1.2·TACH) + Cclon)
GIRK 합집합    ITONE = EMU + EA2 − EMU·EA2                 ← 합이 아니라 합집합
GAP            GAP = (AD + AT) − ITONE
LC 구동        LC = LCB0·(1+0.55·(1−CARE)) + 11.9·GAP⁺·(1−0.30·CARE)

부프레노르핀이 두 번, 서로 다른 숫자로 들어가는 것이 이 모델의 구조적 약속이다. 급성 효과에는 부분작용제 천장(EMAXB = 0.80, 금단 억제라는 종말점에 대해서), 적응 구동에는 β-arrestin 동원이 약한 것(BADRV = 0.40)이 별도로 들어간다. 그래서 부프레노르핀 노출아의 낮은 A0(0.487 대 0.696)는 "중증도 상수"로 맞춘 값이 아니라 수용체 약리학에서 유도된다.

시작 역치 → 잠금(latch) → 상승(속도 제한) → 안정화 용량 잠금 → 48시간 안정 → 네 가지 이유 규칙 중 하나. $ODE 안에 들어 있으므로 프로토콜 자체가 시뮬레이션 대상이며, 위 1-2절의 교환곡선은 이 블록에서 나온다.

nows_qsp_model.dot — 188개 노드, 17개 클러스터 (산모/태반 · 신생아 PK · 성숙 · 치료 PK · µ 수용체 · 신경적응 · 청반 · 다른 약물의 작용지점 · 다물질 노출 · CNS 징후 · 자율신경 징후 · 위장관 징후 · 측정도구 · 돌봄환경 · 성장 · 프로토콜 · 임상 종말점). SVG · PNG

R 툴체인이 없는 환경에서 만들어졌으므로, mrgsolve 모델의 모든 ODE와 대수 블록을 파이썬/scipy로 독립 재구현하여 검증한다 (nows_verify_python.py). 이 스크립트는 두 가지를 한다.

1. nows_mrgsolve_model.R의 $PARAM 블록을 직접 파싱해 145개 파라미터를 전수 대조한다 → 불일치 0.
2. 61개 앵커(약동학 · 구조 불변량 · 임상)를 돌린다 → 61/61 PASS.

전체 출력: nows_verification_output.txt

• MOTHER의 누적 모르핀 9.5배 차이는 재현하지 못한다. 모델은 0.50배 (2배 차이)까지만 간다. MOTHER의 극단적 비는 치료받은 아기의 용량 차이와 치료받은 비율의 차이가 곱해진 결과인데, 이 모델의 부프레노르핀 노출 코호트는 치료 비율이 40.0%로 MOTHER의 47%와 비슷하므로, 남는 차이는 용량에서 나와야 한다. 그만큼의 용량 차이를 만들려면 BADRV를 더 낮춰야 하고 그러면 부프레노르핀 노출아의 무치료 최고점수가 임상적으로 너무 가벼워진다. 이 긴장은 해소하지 않고 남겨 둔다.
• 치료 기간의 절대값은 아기의 성질이 아니라 이유 규칙의 성질이다. 기본 시나리오의 22.8일은 문헌 범위(4–28일)의 중간이며, 어느 한 기관의 재원일수를 예측한다고 읽어서는 안 된다.
• 신경발달 결과는 예측하지 않는다. AUCGAP/AUCSED는 노출 지표이지 결과 예측자가 아니다. NOWS와 후기 발달을 잇는 관찰 문헌은 사회적 환경에 너무 교란되어 있어 보정 대상으로 쓸 수 없다.

파이썬 재구현과 앵커 통과 과정에서 다섯 가지 실질적 결함이 드러났고, 전부 구조를 고쳐서 해결했다(파라미터로 덮지 않았다).

1. 치료가 신경적응을 영구히 유지하는 덫. 단일 풀 모형에서 dA/dt = KAON·EMU·(1−A) − KAOFF·A는 EMU ≈ 0.47에서 A* = 0.62에 고정점을 가지므로, 아기가 어떤 용량에서도 이유(離乳)를 시작하지 못하고 40일 넘게 같은 점수(6.0)에 머물렀다. 지속성/재유도 두 풀로 분리해서 해결.
2. µ 점유가 너무 늦게 풀림. 내성 이동(KSHIFT)이 없을 때 잔류 메타돈 5 ng/mL이 여전히 66%를 점유해서 발현이 4일 뒤로 밀렸다. 신경적응된 신생아의 농도-효과 곡선이 오른쪽으로 이동한다는 사실을 넣어 해결.
3. 클로니딘이 모르핀을 완전히 대체해 버림. EC50 0.45 ng/mL·Emax 0.75에서 6 µg/kg/일이 치료 자체를 불필요하게 만들었다. EC50 1.6 ng/mL·Emax 0.34로 고쳐 보조제 수준으로 되돌림.
4. 모유 용량이 20배 과대. 산모 일일용량에 직접 비례시킨 초기 식이 0.63 mg/kg/일을 신생아에게 주고 있었다. 상대영아용량(RID, 체중보정 산모용량의 2.5%)으로 다시 쓰자 4.5 ng/mL로 떨어졌고, 이는 Begg 2001이 보고한 영아 혈장 2–9 ng/mL와 일치한다.
5. 미숙아 실험이 실험이 아니었음. A0가 제대로 임신주수에 의존하지 않아 34주와 39주가 같은 A0를 갖고 있었다. 적응 능력에 임신주수 시그모이드 (GA50A 33.5주, Hill 8)를 넣어 해결.

또한 적분기가 초기 스텝 추정에서 멈추는 문제가 있어 first_step을 명시했다. 이것은 모델이 아니라 수치 문제였다.

신경적응 블록은 인간 신생아에서 측정된 적이 없다. 구조(자궁 내에서 설정되는 지속성 풀 + 현재 노출이 다시 모는 재유도 가능 풀)는 설치류 청반 문헌과 임상 경과에서 추론한 것이고, 두 시간상수 — 지속성 풀 12일, 재유도 풀 60시간 — 이 이 모델에서 가장 노출된 숫자다.

원리적으로는 식별 가능하다. 서로 다른 이유 규칙에서 얻은 연속 점수가 이 두 상수를 구분한다: 지수 이유(현재 용량의 −10%/일)와 절대 계단 이유 (표준 용량의 −10%/일)는 지속성 풀의 반감기에 대해 서로 다른 민감도를 갖는다. 같은 기관에서 두 규칙을 무작위 배정하고 4시간마다의 점수 궤적을 기록하면 12일이라는 값이 맞는지 몇 주 안에 알 수 있다. 그 자료가 없다면 위의 치료일수 절대값은 모두 조건부다.

부수적으로, AUCSED 지표에는 알려진 한계가 있다. A가 0으로 내려간 뒤에도 남아 있는 소량의 아편 긴장도(예: 모유 수유 팔)는 GAP을 음수로 만들고 따라서 "과진정"으로 적분된다. 그 경우 아기는 진정된 것이 아니라 단지 아주 낮은 수준으로 아편에 노출되어 있는 것이다. 모유 수유 시나리오의 AUCSED 79.1은 이렇게 읽어야 한다.

# 지도
dot -Tsvg nows_qsp_model.dot -o nows_qsp_model.svg
dot -Tpng -Gdpi=150 nows_qsp_model.dot -o nows_qsp_model.png

# 검증 (numpy + scipy)
python3 nows_verify_python.py
python3 nows_population_sim.py

# 모델과 대시보드 (R: mrgsolve, shiny, ggplot2, dplyr, tidyr, DT)
Rscript -e 'source("nows_mrgsolve_model.R"); print(run_all())'
Rscript -e 'shiny::runApp("nows_shiny_app.R")'

교육 및 연구 목적의 QSP 모델이다. 공개 문헌으로 구성했으나 독립적으로 검증·인증되지 않았으며 실제 임상 의사결정, 처방, 규제 제출에 사용해서는 안 된다. 특히 위의 어떤 용량, 역치, 이유(離乳) 일정도 임상 프로토콜이 아니다.