Modèle à deux compartiments avec dose de charge, dose d'entretien et accumulation en cas d'insuffisance rénale.
Le médicament entre dans un compartiment central (V1 — celui que l'on mesure), s'échange avec un compartiment périphérique (V2 — le tissu) par la clairance intercompartimentale Q, et n'est éliminé qu'à partir du compartiment central, avec la clairance CL. La concentration plasmatique est A1/V1 et c'est la seule que l'on puisse mesurer ; la concentration tissulaire A2/V2 est tracée sur la figure du bas mais ne figure jamais sur un résultat de laboratoire.
Les deux microconstantes k10 = CL/V1 et k12 = Q/V1, k21 = Q/V2 donnent les constantes hybrides α et β comme racines d'une équation du second degré. La demi-vie de distribution est ln2/α et la demi-vie terminale ln2/β — elles diffèrent d'un facteur 20 à 60 dans les préréglages, et c'est de cette différence que traite toute la simulation.
La fonction rénale est mise à l'échelle linéairement selon l'approche de Dettli : CL = poids · cl · [(1 − fe) + fe · ClCr/100], où fe est la fraction rénale et cl la clairance par kilo à la valeur de référence, une ClCr de 100. La référence est 100 parce que Cockcroft–Gault, pour un homme de 40 ans pesant 70 kg avec une créatinine sérique de 90 µmol/L, donne 95 mL/min. Cette approche reproduit exactement le modèle de digoxine de Jelliffe — la CL à une ClCr de 25 vaut 4,41 L/h par les deux voies — et la régression de Matzke pour la vancomycine au dixième près.
La prescription est une suite de segments à vitesse de perfusion constante, et sur chaque segment le système est résolu en forme close : la réponse à entrée nulle plus la réponse à état initial nul. Aucune intégration numérique n'est utilisée. La raison en est la raideur — à l'extrémité du curseur de volume, α vaut 120 000 fois β, et un RK4 explicite au pas dont la figure a besoin y est instable. La solution analytique n'a aucune limite de stabilité : chaque segment est exact et le pas ne sert qu'à la résolution du tracé.
La dose de charge remplace la première dose d'entretien, elle ne s'y ajoute pas. La courbe pleine et la courbe fantôme ont donc exactement les mêmes moments d'administration et ne diffèrent que sur un seul point — la taille de la première dose. C'est ce qui rend la comparaison honnête.
Les lectures à l'état d'équilibre viennent du point fixe de la solution périodique, x = (I − Φ(τ))⁻¹ · b(τ), et non du dernier intervalle de la figure. L'horizon est de cinq demi-vies, c'est-à-dire 97 % de l'état d'équilibre, et une concentration résiduelle lue à cet endroit aurait dépendu de la longueur qu'aurait eue l'axe — et aurait été contaminée par une dose de charge qui, par définition, ne modifie pas l'état d'équilibre.
L'horizon est dérivé : cinq demi-vies terminales ou quatre intervalles d'administration, arrondis vers le haut à un nombre entier d'intervalles, avec un plafond à 21 jours. Que l'axe change de longueur quand la demi-vie change fait partie de la leçon — la gentamicine reçoit quatre jours et la digoxine en insuffisance rénale dix-huit. Qu'il s'agisse d'un nombre entier d'intervalles est une exigence de calcul et non un détail esthétique : lue exactement à cinq demi-vies, la « concentration résiduelle » de la digoxine s'établit à 1,42 au lieu de 0,50 µg/L, soit 184 % d'erreur sur le seul chiffre d'après lequel un médecin agit.
L'heure de prélèvement interprétable la plus précoce est dérivée de α − β : l'exigence que le terme de distribution soit inférieur à 5 % du terme d'élimination donne une seule formule, et de celle-ci découle la règle de prélèvement apprise pour chaque médicament — 6,0 h pour la digoxine, 2,7 h pour la vancomycine, 1,9 h pour la gentamicine, 0,9 h pour la lidocaïne.
Chaque médicament porte avec lui ce qu'est sa fenêtre — un pic, une concentration résiduelle ou un taux — et ce contre quoi la dose d'entretien doit être calculée, et ce sont là deux questions et non une. MD = CL · C_cible · τ maintient la concentration MOYENNE à l'objectif et ne vaut donc que pour les fenêtres de taux, la digoxine et la lidocaïne. Les 16–24 mg/L de la gentamicine sont un pic relevé une demi-heure après la fin de la perfusion ; les 10–20 mg/L de la vancomycine sont une concentration résiduelle, tandis que l'objectif de dosage est, depuis Rybak 2020, une AUC24 de 400–600. La cinétique est linéaire, donc chaque lecture à l'état d'équilibre est proportionnelle à la dose : un calcul de référence avec une unité de dose de 1 donne la mesure par unité de dose, et la dose est l'objectif divisé par celle-ci — la même ligne calcule le pic, la résiduelle, la moyenne et l'AUC. Calculée comme une moyenne, la proposition pour la gentamicine atteignait 2 883 mg contre les 490 du préréglage lui-même ; calculée contre le pic à l'heure du prélèvement, elle devient 470 mg, et l'AUC24 se place à 78, à l'intérieur de la fenêtre propre du médicament, 70–120. Trois lectures qui supposent que la concentration doit se trouver dans la fenêtre — la part de l'intervalle passée dans la fenêtre, le délai avant que la résiduelle atteigne la fenêtre, l'écart du prélèvement par rapport à la résiduelle — sont laissées vides devant une fenêtre de pic au lieu d'être affichées : en administration unique quotidienne, la concentration doit quitter la fenêtre, et l'écart du prélèvement par rapport à une résiduelle de 0,03 mg/L atteignait +39 822 %.
dA1/dt = R(t) − (CL/V1)·A1 − (Q/V1)·A1 + (Q/V2)·A2dA2/dt = (Q/V1)·A1 − (Q/V2)·A2 C = A1/V1, Cp = A2/V2α, β = ½·[ s ± √(s² − 4·k10·k21) ], s = k10 + k12 + k21t½,dist = ln2/α t½,élim = ln2/βCL(ClCr) = poids · cl · [ (1 − fe) + fe · ClCr/100 ]LD = V · C_cible MD = CL · C_cible · τ (objectif de moyenne)Css,moy = (MD/τ)/CL AUC24 = dose quotidienne/CLRac = 1/(1 − e^(−β·τ))t_interprétable = ln( (A/B)/0,05 ) / (α − β), A = (α−k21)/(α−β), B = (k21−β)/(α−β)horizon = ceil( max(5·t½,élim, 4·τ) / τ ) · τ, plafond 21 jourspic/résiduelle : MD = C_cible / C(unité de dose 1) AUC : MD = AUC24,cible · CL · τ/24
Simplifications. Voie intraveineuse exclusivement : F = 1, pas d'absorption, pas d'effet de premier passage. La digoxine per os ne peut pas être montrée, et la vancomycine per os contre Clostridioides difficile dans l'intestin est ici conceptuellement impossible. La cinétique est linéaire de bout en bout — la clairance ne se sature pas, si bien que la phénytoïne, l'éthanol et les salicylés à forte dose ne peuvent pas être modélisés, et c'est justement la phénytoïne qui est le médicament le plus important pour lequel le raisonnement sur l'intervalle d'administration s'effondre. Chaque paramètre est constant sur toute la durée : le patient du modèle à une ClCr de 20 a une ClCr de 20 depuis avant la première dose et l'aura pour toujours, ce qui est la seule chose qui ne soit jamais vraie du patient dont les concentrations montent. Il n'y a pas de compartiment d'effet, et le modèle ne peut donc pas dire quelle dose de charge est la bonne — il calcule V1 · C_cible et Vss · C_cible côte à côte, et laquelle s'applique dépend de l'endroit où siège le récepteur : le canal sodique s'équilibre avec le plasma en quelques minutes, la Na/K-ATPase siège dans le compartiment périphérique. La concentration cible de la dose de charge est le milieu de la fenêtre, quelle que soit cette fenêtre ; la dose d'entretien, en revanche, se calcule contre ce que la fenêtre EST — un pic pour la gentamicine, une AUC pour la vancomycine, un taux pour la digoxine et la lidocaïne — et un pic se lit à l'heure de prélèvement que la fenêtre suppose, une demi-heure après la fin de la perfusion, et non au Css,max du modèle, qui pour la gentamicine s'établit à 31,9 mg/L là où la même dose donne 20,9 à une heure. Cette heure de prélèvement est une convention inscrite dans le profil du médicament et non quelque chose que le modèle puisse déduire : d'une fenêtre citée sans sa règle de prélèvement on ne peut pas tirer une dose. La gentamicine porte deux objectifs publiés à la fois, un pic de 16–24 mg/L et une AUC24 de 70–120, et le modèle dose contre le pic ; la vancomycine porte une fenêtre de résiduelle et un objectif d'AUC, et dose contre l'AUC parce qu'une résiduelle de 15 mg/L donne une AUC24 de 651 — au-dessus du plafond, c'est-à-dire précisément le surdosage contre lequel le document de consensus a été écrit. Trois lectures sont donc laissées vides devant une fenêtre de pic — la part de l'intervalle passée dans la fenêtre, le délai avant que la résiduelle atteigne la fenêtre, l'écart du prélèvement par rapport à la résiduelle : toutes trois mesurent que la concentration SE TROUVE dans la fenêtre, et en administration unique quotidienne elle ne le doit pas. Le volume est maintenu constant lorsque la fonction rénale change, ce qui est faux pour la digoxine, dont le volume à l'état d'équilibre diminue d'environ un tiers en insuffisance rénale marquée — la dose de charge de la digoxine doit donc être abaissée en insuffisance rénale bien que la dose de charge soit indépendante de la clairance. La clairance diminue linéairement avec la clairance de la créatinine et avec rien d'autre : pas de sécrétion tubulaire saturable, pas d'interaction par la P-glycoprotéine, pas d'insuffisance hépatique, pas d'effet du débit cardiaque sur les médicaments dont la clairance est limitée par le débit, pas de modification de la liaison aux protéines. Pas de dialyse — l'hémodialyse intermittente est une clairance allumée quatre heures et éteinte quarante-quatre, et le modèle n'en a qu'une seule, continue. La concentration est totale et non libre ; en cas d'hypoalbuminémie, la concentration totale de phénytoïne se lit basse alors que la fraction libre est normale, ce qui est la plus fréquente erreur d'interprétation d'un dosage médicamenteux à l'hôpital et une chose que le modèle ne reproduit pas du tout. La dose de charge fractionnée n'existe pas : une charge, administrée une fois, alors que la digoxine se charge en réalité en trois étapes et l'amiodarone sur des semaines, précisément pour éviter le pic. Les chiffres sont des moyennes de population et non ce patient-ci — les sept patients de Koup s'étalaient sur un facteur 2,5 en volume et 3,5 en demi-vie, et le modèle n'a aucune variabilité interindividuelle ; c'est exactement pour cela que le suivi thérapeutique pharmacologique existe. Pas de pharmacodynamie : la concentration n'est pas l'effet, la fenêtre thérapeutique est tracée parce que les sources en rapportent une et non parce qu'une limite franchie déclencherait quoi que ce soit. Un pic au-dessus du seuil toxique qui se situe dans la phase de distribution n'est pas une toxicité — 500 µg de digoxine en 20 minutes portent la concentration centrale au-dessus de 4 µg/L et restent une prescription normale, puisque l'organe cible ne s'est pas encore équilibré. L'amiodarone est délibérément impossible à montrer : un volume de distribution de 66 L/kg et une demi-vie terminale de 40–55 jours se situent hors du curseur de volume comme du plafond d'horizon, et un axe des x d'un semestre rendrait illisible tout autre cas. Les unités sont celles du médicament et le modèle ne les contrôle pas — le calcul se fait en unité de dose par litre, si bien qu'une dose de vancomycine saisie alors que les étiquettes d'un préréglage de digoxine sont encore en place donne un chiffre faux d'un facteur mille sans le moindre avertissement. La digoxine est en outre rapportée en nmol/L dans les pays nordiques et en µg/L en Allemagne et en France ; le modèle calcule en µg/L, et le facteur est 1,28. La clairance de la créatinine est ici celle de Cockcroft–Gault et non un DFG estimé par CKD-EPI : chaque règle de posologie citée a été établie contre la première, et les deux diffèrent le plus chez précisément les patients où cela compte. Le poids n'est pas un seul chiffre — la vancomycine se dose sur le poids corporel total, la gentamicine sur le poids ajusté, la digoxine sur le poids idéal, et le modèle a un seul curseur de poids qu'il applique à tout. Rien n'est animé : les courbes sont retracées quand un curseur est déplacé et ne bougent jamais d'elles-mêmes, si bien que prefers-reduced-motion n'a aucun effet ici.
Le cas ordinaire : 1500 mg de charge en 90 minutes, puis 1000 mg toutes les douze heures. La concentration résiduelle tombe dans la fenêtre et l'AUC24 au milieu de l'intervalle cible de 400–600.
Le prélèvement a été fait pendant la phase de distribution. Il reflète la redistribution et non l'élimination, et n'est donc pas une concentration résiduelle.