Comment un inhibiteur ou un inducteur enzymatique déplace la concentration plasmatique au fil du temps.
Trois horloges commandent une interaction CYP, et le modèle les tient séparées. L'exposition au précipitant u(t) est analytique : elle se construit vers 1 avec la demi-vie PROPRE de l'inhibiteur ou de l'inducteur et, après l'arrêt, décroît à partir de la valeur qu'elle a eu le temps d'atteindre. La dent de scie entre les prises n'est pas modélisée — une variation du résiduel au pic est sans intérêt à l'échelle de la semaine et aurait coûté un paramètre qui n'apprend rien.
Le pool enzymatique est un modèle de renouvellement : synthèse moins dégradation, plus l'inactivation irréversible lorsque cette branche est activée. L'activité est le pool multiplié par l'inhibition compétitive instantanée dans le cas d'une inhibition réversible, et le pool seul dans le cas mécanisme-dépendant. La constante d'inactivation est étalonnée de sorte que les deux branches aboutissent au même équilibre, (1 + IC)(1 − IR). L'interrupteur ne change donc que l'évolution dans le temps, jamais l'amplitude — et c'est toute la raison de sa présence, puisque c'est précisément cette confusion qui fait croire que la clarithromycine relâche son emprise une fois le dernier comprimé pris.
Le médicament objet est traité comme un compartiment avec un apport d'entretien continu, exprimé en multiples de l'équilibre de départ. La clairance varie avec l'activité selon la forme d'Ohno–Hisaka, avec un chiffre de contribution unique, CR, pour toute la clairance orale. Elle a été préférée à la branche hépatique et intestinale séparée de Fahmi et de la FDA pour trois raisons : elle est étalonnée sur 113 études d'interaction chez l'homme et les restitue, les valeurs d'IR et de CR publiées sont précisément les valeurs de cette forme, si bien que chaque chiffre par défaut a une source au lieu d'être une supposition, et une même équation couvre à la fois l'inhibition et l'induction.
La figure montre deux courbes et non une. La courbe en pointillé est le rapport d'exposition pour une dose PRISE à cet instant précis et saute immédiatement en cas d'inhibition réversible ; la courbe pleine est la concentration chez celui qui prend le médicament tous les jours et suit avec retard, à la demi-vie propre du médicament, désormais allongée. La surface teintée entre les deux est ce retard, et elle est la réponse à la question de savoir si une interaction se manifeste immédiatement ou au bout d'une semaine : une dose unique suit la courbe en pointillé, un traitement d'entretien la courbe pleine.
La bande sous la figure montre le délai jusqu'à 90 % — 3,32 demi-vies — des trois processus partiels sur le même axe en jours, de sorte qu'un trait vertical puisse être tracé entre les figures. Le repère de l'enzyme est grisé et marqué comme inactif en cas d'inhibition purement réversible sans induction, puisque le pool ne bouge alors jamais. Le repère est le temps de renouvellement du pool : en cas d'inhibition mécanisme-dépendante, le pool est démoli avec k_dég + k_inact, donc plus vite que ne le laisse croire le repère, mais il se reconstruit toujours avec k_dég seul. L'horloge de l'enzyme est donc l'horloge de la récupération, non celle de l'installation — et c'est pour la même raison que le losange, le t90 simulé, se trouve parfois à gauche du repère de l'enzyme.
L'intégration est exponentielle et non explicite à la Euler : la mise à jour est exacte pour des coefficients constants par morceaux et inconditionnellement stable, avec des sous-pas à l'intérieur de chaque pas de sortie fixés par le processus restant le plus rapide. Ce n'est pas une optimisation mais une nécessité, puisque les constantes de vitesse couvrent quatre ordres de grandeur dans l'espace des paramètres autorisé. La fenêtre temporelle est dérivée des réglages — la durée de la cure plus cinq fois la plus longue constante de temps — de sorte que l'arrêt tienne toujours dans l'image.
u(t) = 1 − e^(−kp·t) pendant la cure u(t) = u(T)·e^(−kp·(t−T)) aprèsdE/dt = k_dég·(1 + IC·u) − (k_dég + k_inact)·Ek_inact = k_dég · IR·u / (1 − IR·u) inactivation mécanisme-dépendante seulementA(t) = E·(1 − IR·u) réversible A(t) = E mécanisme-dépendanteq(t) = (1 − CR) + CR·A(t) dC/dt = k_él·(1 − q·C)Rapport d'exposition pour une dose prise en t : 1/q(t)À l'équilibre : 1 / [ (1 − CR) + CR·(1 + IC)(1 − IR) ]Plafond = 1 / (1 − CR) Nouvelle demi-vie = t½ / q_ss
Simplifications. Le modèle ne connaît aucun transporteur. La P-gp, la BCRP et l'OATP1B1 sont absentes, et ce n'est pas un détail : la rifampicine inhibe l'OATP1B1 dès la première dose et induit le CYP3A4 sur une semaine — des directions opposées sur des horloges différentes, et le pic aigu n'existe tout simplement pas dans la courbe. Près de la moitié de ce qui figure comme interaction CYP dans les RCP est pour partie une interaction de transporteur. Le modèle n'a qu'un seul pool enzymatique : la paroi intestinale et le foie sont fondus en un seul chiffre de contribution, si bien qu'un précipitant sélectif de l'intestin comme le jus de pamplemousse ne peut être représenté que comme une inhibition corps entier effective. Ce que cela coûte peut être chiffré — la rifampicine induit le CYP3A4 intestinal environ dix-neuf fois et celui du foie environ 2,4 fois, et c'est toute la raison pour laquelle le midazolam oral perd 96 % de son AUC alors que sa demi-vie est à peine plus que divisée par deux. Cette asymétrie, le modèle ne peut pas la montrer. Il calcule un médicament et un précipitant ; les patients réels en ont cinq, et les inhibiteurs se disputent la même enzyme sans que leurs effets s'additionnent. Il est linéaire : pas de saturation, pas d'auto-induction, pas de cinétique dose-dépendante. La phénytoïne est grossièrement sous-estimée, la carbamazépine est absente en tant que phénomène, et un précipitant administré avec une dose de charge — l'amiodarone en est l'exemple — s'accumule plus vite que sa demi-vie terminale ne le laisse croire, si bien que le modèle est là trop lent. Il n'a pas de métabolites, et les promédicaments retournent tout le raisonnement : un CYP2C19 inhibé donne moins d'effet du clopidogrel, un CYP2D6 inhibé moins d'effet de la codéine, du tramadol et du tamoxifène. La courbe montre toujours la molécule mère, si bien qu'une courbe qui monte signifie là un traitement moins efficace. Il n'a pas de pharmacodynamie — l'INR sous warfarine suit la concentration avec deux à trois jours de retard, car le facteur II a une demi-vie d'environ 60 heures, et la répercussion réelle arrive donc plus tard que ne le montre la courbe, par-dessus tous les autres retards. Il ne connaît pas la génétique : la CR est saisie comme un chiffre, et le modèle ignore qu'un métaboliseur lent du CYP2D6 a déjà une CR proche de zéro pour cette enzyme et ne peut donc subir aucune interaction, ou qu'un patient déjà sous inhibiteur est phénoconverti. Pas d'insuffisance rénale ni hépatique, pas d'âge, pas de poids, pas d'observance. Il donne la concentration moyenne sur l'intervalle d'administration et non les pics et les résiduels, si bien qu'une toxicité pilotée par le pic — l'allongement du QT en est l'exemple le plus net — est plus élevée que ne le montre la courbe. L'absorption est instantanée : pas de constante d'absorption, pas de délai, pas d'interaction alimentaire. Le temps de renouvellement de l'enzyme est le chiffre le plus incertain du domaine, avec des estimations publiées pour le CYP3A4 entre 10 et 140 heures, et le choix déplace l'évolution prédite de plusieurs fois ; le curseur existe précisément pour cela — l'incertitude doit se voir, non se cacher dans une constante. Les paramètres, enfin, se combinent librement et produisent alors des extrêmes arithmétiquement cohérents mais pharmacologiquement jamais observés : les plus fortes interactions réellement mesurées chez l'homme se situent autour de 20–30 fois.
La demi-vie de l'azithromycine est de près de trois jours, et cela n'a aucune importance : l'inhibition est de 11 % et le rapport reste sous 1,1×. Le plafond de la simvastatine est en revanche de 33×, et c'est ce qu'un autre antibiotique peut atteindre.