Warfarin, DOAK og heparin over tid — opstart, seponering og reversering.
Fire lag, ren fremad-Euler med dt = 0,05 h over et vindue på mindst fjorten døgn. Lag et er lægemiddeleksponeringen: ét kompartment med førsteordensabsorption for alt, der gives oralt eller subkutant, hvor nyrefunktionen skalerer clearance lineært i den renalt udskilte andel — nul for warfarin, 0,27 for apixaban, 0,33 for rivaroxaban, 0,50 for edoxaban, 0,80 for dabigatran og 0,56 for enoxaparin. Absorptionshastigheden løses ud af den publicerede Tmax i stedet for at angives frit, fordi Tmax følger af absorptionen og eliminationen tilsammen.
Ufraktioneret heparin får i stedet mætbar elimination af Michaelis–Menten-type, fordi den dosisafhængige halveringstid er hele pointen med UFH. Modellen giver effektiv halveringstid 65, 81 og 162 minutter ved 0,2, 0,5 og 2,0 IE/mL, inden for de 30–150 minutter, Hirsh rapporterer for 25–400 E/kg, og infusionen 18 E/kg/h lander på 0,50 IE/mL — midt i det terapeutiske 0,3–0,7.
Koncentrationsskalaen er kalibreret, ikke udledt. Publicerede V/F og CL/F stammer fra tokompartmenttilpasninger og kan ikke samtidig reproducere top, dal og terminal halveringstid i en enkompartmentmodel. Derfor køres en referencesimulering pr. middel ved standarddosis og CrCl 100, og hele sporet skaleres, så dets steady state-top bliver den publicerede Cmax,ss. Det, der står på et laboratoriesvar, er ng/mL — så skal figurens akse være de samme ng/mL. Dalværdien falder ud af modellen og vises over for den publicerede som selvkontrol: apixaban 103 mod 103, dabigatran 103 mod 91, rivaroxaban 53 mod 44, edoxaban 42 mod 36 ng/mL.
Lag to gælder kun warfarin og er en Emax-kurve med Hills eksponent 1,15 på S-warfarinkoncentrationen. IC50 er sat til 0,10 mg/L, hvilket giver INR omkring 2,5 ved 5 mg/døgn. Warfarinfølsomhedsskyderen skalerer IC50 omvendt og repræsenterer VKORC1-genotype, CYP2C9-genotype og interagerende lægemidler i ét eneste tal.
Lag tre er selve simuleringen: en indirekte responsmodel med hæmmet produktion, én ligning pr. vitamin K-afhængigt protein, hvor den samme syntesblokade møder fem forskellige ure. FVII har halveringstid 5 timer, protein C 8, FIX 24, FX 40 og FII 67 — det sidste ikke et lærebogstal, men Shapiro og Martinez' måling med jodmærket protrombin. At FVII ligger på 25 procent, mens FII stadig står på 80, er ikke en detalje, men hele lektionen.
Lag fire er aflæsningerne. INR beregnes ud af faktorniveauerne som en vægtet geometrisk protrombinkompleksaktivitet med vægtene 0,45 for FVII, 0,35 for FX og 0,20 for FII: FII er altså den mindst betydningsfulde af de tre for INR, selv om den er den vigtigste for antitrombotisk effekt. INR på DOAK modelleres som den artefakt, den er, lineært i koncentrationen, med hældninger, der gør, at rivaroxaban når INR 2,1 ved et helt normalt topniveau, mens apixaban standser på 1,09. APTT-ratioen er lineær for UFH og krumlineær for dabigatran, hvor den flader ud og bliver ufølsom ved dalværdi.
Reverseringerne er mekanismer, ikke specialtilfælde. K-vitamin ophæver blokaden i stedet for at tilføre faktorer. 4F-PCC er et spring i faktorniveauerne, hvorefter den tilførte aktivitet nedbrydes med hver faktors egen halveringstid — og derfor vender INR tilbage inden for et halvt døgn uden en eneste ekstra linje kode. Specifikke antidoter neutraliserer en andel af den cirkulerende koncentration. En kombination, der ikke kan virke, beregnes ikke, men afvises: protamin på apixaban gør ingenting, og modellen siger hvorfor.
dAg/dt = −ka·Ag + Dosis·δ(t − t_n) dAc/dt = ka·Ag − (CL/V)·Ac C = Ac/VCL(CrCl) = CL_ref · [ (1 − f_renal) + f_renal · CrCl/100 ]UFH: dC/dt = R(t)/V − Vmax·C/(Km + C) V 0,07 L/kg · Vmax 0,77 IE/mL/h · Km 1,0 IE/mLE(t) = C^1,15 / (IC50^1,15 + C^1,15) IC50 = 0,10 mg/L S-warfarin / (følsomhed/100)dF_i/dt = k_i · [ 100·(1 − E_eff(t)) − F_i ] k_i = ln2 / t½_it½: FVII 5 h · protein C 8 h · FIX 24 h · FX 40 h · FII 67 hPCA = (FVII/100)^0,45 · (FX/100)^0,35 · (FII/100)^0,20 INR = PCA^(−0,66)DOAK: INR = 1 + s·C/100 s: rivaroxaban 0,45 · edoxaban 0,30 · dabigatran 0,20 · apixaban 0,05APTT-ratio: UFH 1 + 2,1·C dabigatran 1 + 0,75·√(C/100)K-vitamin: E_eff(t) = E(t)·(1 − K(t)) K(t) = K_max·(1 − e^(−(t−t_rev)/τ_on))·e^(−(t−t_rev)/τ_off)4F-PCC: F_i ← min(100, F_i + a_i·dosis) a: FII 1,6 · FVII 1,2 · FIX 1,0 · FX 1,8 · protein C 1,2 %-point pr. E/kgAntidot: C_fri(t) = C_tot(t)·(1 − R(t)) R = R_max under infusionen, derefter R_max·e^(−ln2·(t−t_slut)/t½)
Forenklinger. Modellen beregner lægemiddelniveauer og faktorniveauer over tid. Den beregner ikke blødningsrisiko, tromboserisiko eller udfald, og den kan ikke sige, at en patient kommer til at bløde eller få en blodprop. Der beregnes ingen trombingenerering: der findes ingen ETP, ingen cellebaseret koagulationsmodel, ingen trombocytter, intet fibrinogen, ingen von Willebrand-faktor og ingen fibrinolyse. Faktorniveauerne er et surrogat for hæmostatisk evne, og sammenhængen mellem dem er ikke lineær. Antitrombin mangler helt, hvilket betyder, at heparinresistens er umulig i modellen — en patient med antitrombinmangel svarer her lige så godt som alle andre, hvilket vedkommende ikke gør i virkeligheden. INR er en kalibreret empirisk funktion af tre faktorniveauer, ikke en simuleret koagulationsprøve: ISI, reagens og instrument modelleres ikke, samme plasma giver i virkeligheden INR, der adskiller sig 0,3–0,5 mellem laboratorier, og PT-artefakten ved DOAK er stærkt reagensafhængig. Alle vitamin K-afhængige proteiner antages at blive hæmmet med samme andel af deres syntese, hvilket målte niveauer ikke viser: sammenlignet med Harrison 1997 giver modellen en FVII, der ligger omkring tyve procentpoint for LAVT ved 36 timer (32 mod 54 procent på 5 mg), mens FII rammer inden for to enheder. Restitutionen efter seponering er for langsom af samme grund: modellen skal bruge 145 timer til INR under 1,2, hvor White i 1995 målte 96–115, fordi FX og FII bærer 0,55 af INR-vægten med halveringstider på 40 og 67 timer, samtidig med at S-warfarinets egen udskillelse holder syntesen blokeret endnu et døgn. Efter 4F-PCC lander modellen på INR omkring 1,6, hvor Sarode i 2013 så 62 procent nå 1,3 eller lavere — og indsætter man Sarodes egne målte faktorniveauer i INR-afbildningen, giver den 1,50, så afvigelsen sidder i afbildningen og ikke i PCC-koefficienterne. Farmakokinetikken er populationsgennemsnit: variationskoefficienten for DOAK-niveauer mellem patienter er 30–40 procent, og warfarinets dosisbehov varierer tredivefold, så en kurve er en median og ikke en patient. Enoxaparinets fordelingsvolumen er fast 5,24 L ifølge TIMI 11A, hvilket gør topniveauet proportionalt med vægten på en måde, der overdriver forskellen hos den meget tunge. Interaktioner modelleres kun som en følsomhedsknap for warfarin; P-gp- og CYP3A4-hæmmere hæver DOAK-niveauerne 1,5–2,5 gange, og rifampicin halverer dem, og intet af det findes pr. lægemiddel. Intet leversvigt, intet absorptionssvigt, ingen manglende adhærens, ingen dialyse — dabigatran er dialyserbart, de øvrige er det ikke, og den forskel ses ikke her. Ingen HIT, intet antifosfolipidsyndrom, ingen mekaniske klapper, ingen graviditet og ingen børn. Reverseringsmidlerne er beskrevet som en neutraliseret andel med empirisk aftagende varighed, ikke som bindingskinetik: andexanets egen protrombotiske virkning modelleres slet ikke, selv om ANNEXA-4 rapporterede tromboembolisk hændelse hos 10 procent og død hos 14 procent inden for 30 dage, og dabigatranets tilbagevenden efter idarucizumab sker her hos alle, hvor RE-VERSE AD så den hos omkring en femtedel. Sikkerhedstærsklerne 30 og 50 ng/mL er konsensusgrænser, ikke udfaldsvaliderede; spinale hæmatomer er forekommet ved meget lave niveauer, og at søjlen slutter, betyder, at koncentrationen har passeret et aftalt tal, intet andet. Modellen træffer ingen beslutninger. Den viser niveauer over tid.
Den sædvanlige opstart hos en patient med atrieflimren, normal nyrefunktion og uden interagerende lægemidler. INR bliver terapeutisk længe før protrombinet er det.