Lp(a)는 이상하게 다루기 어려운 위험인자입니다. 인과성은 멘델 무작위화로 확립되었고, 표적은 안전하며(LPA 결손자는 건강합니다), 90% 이상 낮추는 약도 이미 있습니다. 그런데도 이 분야는 서로 모순돼 보이는 관찰들로 가득합니다.

• 스타틴은 LDL 수용체를 올리는데 Lp(a)는 오히려 올립니다.
• PCSK9 억제제는 LDL-C를 60% 낮추면서 Lp(a)는 25-30%만 낮춥니다. 같은 수용체인데 왜 다른 답이 나오나요?
• 나이아신은 Lp(a)를 20% 낮추는데 임상시험은 완전히 음성이었습니다.
• 무발라플린은 검사 방법에 따라 −85.8%로도 −70%로도 보고됩니다.
• 대동맥판막 협착에 대한 지질강하 시험(SEAS·ASTRONOMER·SALTIRE)은 전부 음성이었는데, 유전학은 강한 인과 신호를 보여줍니다.
• 멘델 무작위화가 요구하는 감소폭과 임상시험이 요구하는 감소폭이 3배 다릅니다.
• 같은 환자의 Lp(a)가 mg/dL로는 역치 미달, nmol/L로는 역치 초과로 나옵니다.

교과서는 이것들을 각각 별개의 주의사항으로 나열합니다. 이 모델은 그것들이 하나의 구조에서 나오는 산술적 귀결임을 보이려고 만들었습니다.

Lp(a)는 하나의 입자이고, 하나의 율속 단계에서 만들어지며, 세 개의 작용 팔로 병을 일으키고, 두 개의 서로 호환되지 않는 눈금으로 읽힙니다.

50개 미분방정식은 모든 시나리오에서 완전히 동일합니다. 위 일곱 가지 모순은 파라미터를 바꿔 재현한 것이 아니라, 이 구조에서 계산되어 나옵니다.

Lp(a)의 분획 이화율(FCR)은 약 0.25 pools/day 이며, 10배의 농도 범위에 걸쳐 사실상 변하지 않습니다(Rader 1993, 1994). 즉 사람 사이의 차이는 전부 생산 쪽입니다. 모델은 이것을 가정이 아니라 구조로 넣었습니다: KCAT0는 고정하고, 목표 기저 농도로부터 번역속도 KTL을 역산합니다.

생산 사슬에서 아이소폼 크기가 개입하는 지점은 단 하나, 소포체에서의 분비 전 분해입니다.

SECEFF = KSZ³ / (KSZ³ + NKIV2³)

이 항 하나가 약 5배를 만들고, 나머지는 대립유전자별 전사량 차이입니다.

이것이 청소 표적 약물 이야기의 전부이며, 모델은 이를 주장하지 않고 유도합니다.

KCAT0 = KLDLR_LPA + KOTH_LPA = 0.060 + 0.190 = 0.250 /day   (LDLR 의존 24%)
LDL   = KLDLR_LDL + KOTH_LDL = 0.300 + 0.050 = 0.350 /day   (LDLR 의존 86%)

에볼로쿠맙이 LDLR을 1.00 → 2.75로 올리면, 같은 수용체 변화가 두 지단백에 전혀 다른 답을 줍니다.

25-30%라는 관측값은 모델에 넣은 적이 없습니다. 24%라는 수용체 생물학에서 나옵니다.

같은 산술이 스타틴에서는 반대 방향으로 작동합니다.

청소 항 :  LDLR 1.00 → 1.78  →  Lp(a) −15.7%
전사 항 :  스테롤 반응요소 →  ×1.30
                              ────────────
                       순효과   +8.1%     ← 관측 +8 ~ +20%

전사 항 ESTA = 0.30은 이 엔드포인트에 맞춘 유일한 값이지만, 임의의 값이 아닙니다. 청소 산술이 주어진 이상 이보다 작은 값으로는 순증가 자체가 불가능합니다.

그리고 이것은 검증 가능한 예측을 낳습니다: LDL 수용체가 없는 환자(HoFH)에서는 상쇄해 줄 청소 이득이 없으므로 스타틴이 Lp(a)를 전사 항 전부만큼 올려야 합니다. 모델값 +25.9%.

LDLR만 올리고 LPA 전사에는 손대지 않는 약은 Lp(a)를 LDL의 약 1/4만큼 낮춰야 합니다.

에제티미브는 음성 대조가 아니라 24%라는 숫자에 대한 정량적 시험이고, 통과합니다.

이 데이터는 동시에 조립 단계가 LDL 기질에 포화되어 있어야 함을 강제합니다. 조립이 LDL에 1차라면 에제티미브가 Lp(a)를 20% 가까이 더 낮췄어야 하기 때문입니다. (이것은 실제로 초판 모델의 결함이었고, 아래 §검증에서 다룹니다.)

유리 apo(a)를 별도 상태변수로 둔 것이 이 모델의 두 번째 구조적 결정입니다. 조립 억제제(무발라플린)는 mRNA를 건드리지 않으므로 apo(a)는 계속 만들어지고, 입자가 되지 못한 apo(a)가 혈장에 쌓입니다.

intact Lp(a)   250 → 39.0 nmol/L   (−84.4%)     ← KRAKEN −85.8%
유리 apo(a)    2.60 → 10.66 nmol/L (×4.09)

그래서 유리 apo(a)까지 잡아내는 검사는 효과를 과소보고할 수밖에 없습니다. 방향은 정확히 맞습니다. 그런데 크기는 맞지 않습니다.

즉 모델은 이렇게 말합니다: KRAKEN의 불일치가 닫히려면 전통적 검사가 유리 apo(a)에 대해 입자 내 apo(a) 대비 약 3.8배의 몰당 반응을 가져야 한다. 이것은 측정 가능한 양이고, 찾아본 범위에서 보고된 적이 없습니다. RFREE는 3.8이 아니라 1.0으로 두었습니다 — 표를 예쁘게 만들려고 튜닝하는 대신, 모델이 KRAKEN을 덜 설명한다는 사실을 그대로 남겼습니다.

Lp(a)의 병리학적 화물은 전부 입자당 1개입니다.

아이소폼 크기에 따라 변하는 것은 apo(a) 크링글 단백 질량뿐이고, 그것은 병리학적으로 불활성입니다. 그러므로 nmol/L은 단위 취향이 아니라 기전적으로 옳은 눈금입니다.

그런데 모델을 돌려 보니, 흔히 지목되는 범인이 범인이 아니었습니다.

단위 환산은 거의 문제가 없습니다. 항체가 문제입니다.

TRUE 질량이 똑같이 50 mg/dL인 환자 다섯 명:

질량 검사는 정확히 작은 아이소폼 환자 — 유전적으로 가장 위험이 높은 환자 — 를 놓치고, 큰 아이소폼 환자를 가장 세게 표시합니다.

덤으로: 지침이 짝지어 제시하는 50 mg/dL / 125 nmol/L 은 서로 등가가 아닙니다. TRUE 질량 50 mg/dL은 아이소폼에 따라 126-139 nmol/L 이므로, 몰 역치가 언제나 더 포괄적입니다.

Lp(a) 250 nmol/L, 작은 아이소폼 환자에서:

보고되는 LDL-C  127.3 mg/dL
              = 진짜 LDL-C 100.0  +  Lp(a)-콜레스테롤 27.3
                                     (보고값의 21%는 LDL이 아님)

Dahlén 보정을 (편향된) 질량 검사값으로 수행하면 → 108.5 mg/dL
                                    진짜 LDL-C 100.0 대비 +8.5 mg/dL 과소보정

역치 오류와 같은 환자에서 같은 방향으로 어긋납니다.

위험을 시간상수가 다른 두 성분이 나릅니다. 그게 전부입니다.

62년 시뮬레이션 결과 (평생 250 vs 25 nmol/L):

65세 시점 HR                              2.29   (PAV 44.0 vs 29.8%, VULN 1.448 vs 1.000)
60세부터 80% 감소, 5년 적분 위험 감소     20.0%
평생 저Lp(a) 대비 같은 5년 구간           53.1%
                                    ─────────
                              비율      2.65배

이 2.65배는 보정계수로 넣은 값이 아닙니다. 두 개의 시간상수가 만든 값입니다.

퍼센트가 아니라 절대 감소량이 위험을 움직입니다.

AIM-HIGH·HPS2-THRIVE 참가자의 Lp(a) 중앙값(약 30 mg/dL)에서 20% 감소는 6 mg/dL 남짓입니다. 나이아신은 효력이 부족한 Lp(a) 약이 아니라, 어떤 퍼센트로도 성공할 수 없는 집단에 투여된 충분히 강한 약이었습니다.

같은 표가 Lp(a)HORIZON(≥70 mg/dL)과 OCEAN(a)-Outcomes(≥200 nmol/L)의 등록 기준이 왜 그렇게 정해졌는지도 설명합니다. 기저 선택이 퍼센트 효력을 압도합니다.

판막 석회화는 문턱(Hill n=4)을 넘으면 Lp(a)와 무관하게 자기영속합니다.

dVCALC/dt = KCA · ( VIC_OST  +  KSELFR · VCALC⁴/(KSELFH⁴ + VCALC⁴) )
                   ↑ Lp(a) 의존          ↑ Lp(a) 무관, 문턱형

62년(18세 → 80세) 시뮬레이션:

30세에 시작하면 칼슘의 87%를 막습니다. 60세에 시작하면 — 입자의 80%를 제거하는 약을 20년 써도 — 15%입니다.

이것이 SEAS·ASTRONOMER·SALTIRE가 모두 음성이었던 이유이자, 동시에 판막 적응증은 아무리 강력한 약이라도 이차예방 설계로는 성공할 수 없다는 예측입니다.

Lp(a) → OxPL → 단핵구 NF-κB → IL-6 → LPA 프로모터의 IL-6 반응요소 → Lp(a)

모델이 계산한 개방 루프 이득 g = 0.019, 증폭 1.02. 임상적으로 무시할 수준입니다. 루프를 강조하는 서술은 과장입니다.

무시할 수 없는 것은 기저 IL-6 톤입니다. IL-6가 2 pg/mL에 불과해도 프로모터 배수는 이미 1.33이므로, 기저 LPA 전사의 약 25%가 IL-6 구동입니다. IL-6 차단이 제거하는 것은 루프가 아니라 이 톤입니다.

하나의 방정식, IL-6 수치만 바꿔서:

50개 ODE 전체를 파이썬에서 고정 스텝 RK4로 독립 재구현하여 대조했습니다. 아래 다섯 가지는 그 과정에서 발견되어 수정된 실제 결함입니다.

이 밖에 위험 지수 항을 유한하게 제한하고, 평균 압력차 출력을 판막이 거의 막힌 극한에서 보호했습니다.

# 모델 + 20개 시나리오 + 4개 분석 런
Rscript lpa_mrgsolve_model.R

# 대시보드
shiny::runApp("lpa_shiny_app.R")

# 지도 재생성
dot -Tsvg lpa_qsp_model.dot -o lpa_qsp_model.svg
dot -Tpng -Gdpi=150 lpa_qsp_model.dot -o lpa_qsp_model.png

정직하게 남겨 둔 네 가지입니다. 전부 파라미터로 노출되어 있어 새 데이터가 나오면 바로 시험할 수 있습니다.

1. KRAKEN 검사 불일치의 크기. 방향은 재현하지만 크기는 재현하지 못합니다. 닫으려면 RFREE ≈ 3.8이어야 하며, 이는 측정 가능하지만 미보고 값입니다.
2. 조립 위치. 세포외(White 1994) 대 세포내(Frischmann 2012) 논쟁은 미해결입니다. 모델은 세포외를 채택했고, 세포내가 맞다면 무발라플린의 유리 apo(a) 상승 예측이 달라집니다.
3. Lp(a) 질량 중 콜레스테롤 분율. 고전적 0.30 대 최근 0.17-0.25. FCHOL이 이 불확실성을 그대로 노출하며, Dahlén 보정의 과대/과소 여부가 여기에 달려 있습니다.
4. 차단된 조립이 아껴진 apoB를 혈장 LDL로 돌려보내는가 (FRECY). 펠라카르센 시험의 apoB 감소폭은 FRECY가 0에 가까움을 시사합니다.

그리고 하지 않은 것들: 성별 구분 없음, 중성지방/LPL 축 없음, 대립유전자 2개가 아닌 단일 발현 아이소폼 가정, 플라크 기하 없음, rs3798220-아스피린 상호작용 없음. 항섬유소용해 팔은 시험관 내 생물학이 아무리 강력해도 멘델 무작위화가 정맥혈전증과의 연관을 보이지 않으므로 의도적으로 낮은 가중치(WARM3 = 0.25)를 두었습니다. 모델은 생물학이 아니라 근거를 따릅니다.

교육 및 연구 목적의 정성적·반정량적 QSP 모델입니다. 공개 문헌과 임상시험 데이터를 바탕으로 구성했으나 독립적으로 검증·인증되지 않았으며, 실제 임상 의사결정, 처방, 규제 제출에 사용해서는 안 됩니다.