Simulateur pharmacocinétique du méthylphénidate

Explications du modèle

Statut du moteur : cette version utilise un moteur pédagogique unique. Pour Concerta, Ritaline LP/Ritalin LA et Quasym/Metadate CD, les paramètres de forme de libération (t_d, s, t_d1, s₁) viennent strictement de Gomeni 2017 Table 1, mais les fractions f_f sont rescalées par un facteur ×0,62 commun aux trois formes LP MPH (somme f_f1 + f_f2 = 0,62 au lieu de 1,0) et la disposition est celle du simulateur (k_e = 0,25 h⁻¹, amplitude calibrée sur Wargin 1983). Il ne s'agit donc pas d'un mode Gomeni natif strict avec kel = 0,18 h⁻¹ et V/F = 1520 L. Voir le détail de la convention de réconciliation d'échelle plus bas (paragraphe « Paramétrage par spécialité »).

Note de version : la calibration courante utilise dMphFraction = 0,97 (Modi 2000 stéréospécifique + FDA Concerta) et des seuils DAT en d-MPH = 5,8 / 9,7 ng/mL (Volkow 1998 × 0,97). Tout rapport PDF antérieur affichant des seuils 3,8 / 6,3 ng/mL ou T > 1,9 ng/mL a été généré avec un facteur historique 0,63 aujourd'hui rejeté (confusion comprimé racémique 50 % vs plasma ~97 % d-MPH après premier passage hépatique). Ces rapports sont à régénérer.

Standard métrologique : d-méthylphénidate (d-MPH)

Le simulateur rapporte par défaut les concentrations plasmatiques en d-méthylphénidate, énantiomère pharmacologiquement actif qui inhibe les transporteurs DAT et NET. Ce choix s'aligne sur les notices FDA DailyMed modernes (Ritalin LA, Concerta, Metadate CD, Focalin XR) et sur les études cliniques post-2000 qui utilisent des dosages LC-MS/MS stéréosélectifs. L'énantiomère l-MPH, présent à 50 % dans le comprimé racémique, subit un premier passage hépatique massif par CES1 (F ≈ 5–7 %) et contribue très peu aux concentrations plasmatiques utiles (Modi 2000 stéréospécifique, PMID 11028253 ; Patrick 2005, PMID 15706728).

La calibration interne du modèle (C_max des pulses, f_f Weibull) s'appuie sur l'étude historique de référence en MPH total racémique non chiral (Wargin 1983, PMID 6410043). Un facteur de conversion plasmatique d-MPH ≈ 0,97 × MPH total est appliqué à l'affichage final. Cette valeur découle des études stéréosélectives modernes sur Concerta : Modi 2000 stéréospécifique (PMID 11028253) mesure pour 18 mg une AUC l-MPH = 0,43 ng·h/mL contre d-MPH = 42,2, soit un ratio l/d ≈ 1/100 ; la notice FDA Concerta indique « l-isomer plasma concentrations are approximately 1/40 of d-isomer » sur la même formulation. Le facteur historique « 0,63 » parfois cité dans la littérature pédagogique confond la composition du comprimé racémique (50 % d-MPH) avec la composition plasmatique (~97 % d-MPH après premier passage hépatique massif du l-MPH par CES1).

Limite à connaître : ce ratio 0,97 est démontré spécifiquement pour Concerta. Pour les autres formulations (Ritalin LA SODAS, Metadate CD DIFFUCAPS, Medikinet gel) le ratio plasmatique l/d n'est pas systématiquement publié et peut différer modérément (formes IR plus anciennes : ratios l/d souvent plus élevés, jusqu'à ~1/8). De plus, les études primaires utilisées comme cibles de validation emploient des standards analytiques hétérogènes : certaines (Wargin 1983, Pentikis 2002, Markowitz 2003, Wigal 2011) rapportent du MPH total non chiral ; d'autres (Modi 2000 stéréospec, Shram 2012, Schapperer 2015, Roesch 2013) rapportent explicitement du d-MPH. Le facteur 0,97 reste une approximation pratique d'affichage, pas une équivalence physiologique exacte pour toutes les formes. Les repères DAT (Volkow 1998, Spencer 2006) sont convertis par la même règle (6 → 5,8 ng/mL ; 10 → 9,7 ng/mL). Le bouton « Standard d'affichage » en haut de la page permet de basculer entre d-MPH et MPH total racémique sans rechargement.

Cas particulier de la lisdexamfétamine (Elvanse / Xurta) : il s'agit d'une amphétamine, classe pharmacologique distincte du méthylphénidate. La molécule active dans le plasma est directement la d-amphétamine (énantiomère unique obtenu après hydrolyse enzymatique de la prodrogue par la carboxypeptidase érythrocytaire). La conversion d-MPH/total ne s'applique pas, et les seuils DAT MPH (5,8 / 9,7 ng/mL) ne sont pas pertinents pour la d-AMP (cibles thérapeutiques différentes, ~30–50 ng/mL d-AMP d'après Boellner 2010 PMID 20206783 ; mécanisme mixte DAT-inhibiteur + libérateur dopaminergique). Les courbes LDX sont donc rapportées en d-amphétamine indépendamment du standard MPH choisi. Il est déconseillé de mélanger LDX et MPH dans une même stratégie : les unités plasmatiques (d-AMP vs d-MPH) ne sont pas directement comparables sur un axe Y commun.

Équation pharmacocinétique de base

Le simulateur repose sur une superposition de fonctions de Bateman monocompartimentales (absorption d'ordre 1 + élimination d'ordre 1), calibrées sur l'adulte de référence 70 kg (Maldonado 2013, PMID 23646842) :

• demi-vie d'élimination t½ = 2,8 h (k_e = 0,25 h⁻¹) ;
• k_a du pulse IR = 1,5 h⁻¹ ;
• chaque spécialité est décomposée en un ou plusieurs pulses d'absorption, les pulses à libération prolongée (ER) recevant un k_a plus faible pour reproduire les Cmax publiés.

Paramétrage par spécialité

• Ritaline LI — Bateman mono-pulse, 100 % à 0 h, k_a = 1,5 h⁻¹. T_max ≈ 1–2 h, t½ ≈ 2–3 h (Wargin 1983, PMID 6410043 ; Kimko 1999, PMID 10628897). Représentation canonique de la formulation IR.
• Ritaline LP (SODAS) — double Weibull issu de Gomeni 2017 Table 1 (Ritalin LA) : Weibull 1 (release rapide t_d = 0,92 h, s = 2,60) + Weibull 2 (release lent t_d1 = 5,55 h, s₁ = 4,18), ratio f_f1/f_f2 = 36/64 (préservé), fractions rescalées ×0,62 (f_f1 = 0,223, f_f2 = 0,397, somme = 0,62). Cibles de validation (MPH total, 20 mg adulte) : C_max2 ≈ 6,2 ng/mL (DailyMed Ritalin LA ; González 2002, PMID 11996212) ; T_max2 ≈ 5,5 h (Markowitz 2003, PMID 12648029). Profil simulé : C_max = 5,57 ng/mL, T_max = 6,55 h, AUC₂₄ ≈ 49 ng·h/mL.
• Quasym (DIFFUCAPS) — double Weibull issu de Gomeni 2017 Table 1 (Metadate CD) : Weibull 1 (release rapide t_d = 1,04 h, s = 2,87) + Weibull 2 (release lent ER t_d1 = 4,21 h, s₁ = 3,68), ratio f_f1/f_f2 = 39/61 (préservé), fractions rescalées ×0,62 (f_f1 = 0,242, f_f2 = 0,378, somme = 0,62). Cibles de validation (MPH total, 20 mg adulte) : Pentikis 2002 PMID 11931601 fasting C_max 4,6 ; SmPC EMC fed 6,4 ; T_max 4,5–5,0 h (Pentikis ; Rochdi 2004, PMID 14764313 ; González 2002). Profil simulé : C_max = 6,26 ng/mL, T_max = 5,00 h, AUC₂₄ ≈ 49 ng·h/mL.
• Medikinet (matrice gel hydrophile pH-dépendant) — Bateman bi-pulse calibré sur SmPC Medikinet XL (UK EMC 313) et Schütz 2009 (PMID 19954715) : 50 % IR avec lag = 1 h (retard de vidange gastrique post-prandial) et k_a = 1,5 → T_max1 ≈ 2,4 h ; 50 % ER libéré progressivement à partir de 2,5 h sur environ 4 h avec k_aER = 0,3 → épaulement tardif vers 6–7 h, Cmax2 < Cmax1. Modèle calibré pour la prise recommandée avec petit-déjeuner. Cibles PK publiées (20 mg avec repas) : T_max,global ≈ 2,75-2,82 h, C_max ≈ 6,4 ng/mL, second rebond ≈ 3 h après le pic IR (< C_max), t½ ≈ 3,2 h. À jeun, la matrice gel ne se structure pas correctement et le profil devient mono-phasique (Rodríguez-Cruz 2008, PMID 18793577 ; SmPC Medikinet XL EMC 313) — non représenté ici.
• Concerta (OROS) — double Weibull issu de Gomeni 2017 Table 1 (Concerta/OROS, PMID 28369788) : Weibull 1 (surcouche IR, t_d = 0,76 h, s = 3,18) + Weibull 2 (cœur osmotique, t_d1 = 6,33 h, s₁ = 3,40), ratio f_f1/f_f2 = 19/81 (strictement préservé — à distinguer du 22 %/78 % théorique de la composition du drug-overcoat rapporté dans le label FDA Concerta : les 19 %/81 % de Gomeni décrivent la fraction réellement libérée rapidement in vivo obtenue par déconvolution popPK et non la masse de MPH déposée sur la couche externe du comprimé OROS), fractions rescalées ×0,62 (f_f1 = 0,118, f_f2 = 0,502, somme = 0,62). Le modèle popPK a été développé sur Concerta puis qualifié sur plusieurs formulations LP ; les précisions paramétriques (RSE %) publiées pour la référence OROS — f_f 6,66 ; t_d 16,14 ; s 4,29 ; t_d1 4,38 ; s₁ 4,11 — concernent ce modèle de référence et ne doivent pas être lues comme une validation individuelle de chaque formulation. Cibles PK publiées (18 mg adulte, MPH total) : Modi 2000 PMID 10761165 C_max 3,7 ; FDA Concerta package insert C_max 4,4 ; Markowitz 2003 PMID 12648029 C_max 5,9 ; Wigal 2011 PMID 21663428 (pédiatrique 38 kg) C_max 6,0. Profil simulé brut : C_max = 4,80 ng/mL, T_max = 7,55 h, AUC₂₄ = 44,2 ng·h/mL — centre du faisceau littéraire multi-études.
• Elvanse / Xurta (lisdexamfétamine mésilate, LDX) — classe pharmacologique distincte du méthylphénidate : prodrogue inactive convertie en d-amphétamine (d-AMP) par hydrolyse plasmatique via la carboxypeptidase érythrocytaire (rate-limiting), sans premier passage hépatique. Bateman mono-pulse avec k_a apparent = 0,6 h⁻¹ (vitesse de conversion RBC) et k_e override = 0,063 h⁻¹ (t½ d-AMP = 11 h adulte ; ≈ 8,6 h pédiatrique). Cibles PK publiées (Elvanse 70 mg adulte, à jeun) : T_max d-AMP ≈ 3,5–4 h ; C_max d-AMP ≈ 75 ng/mL ; AUC_0-∞ d-AMP ≈ 1130 ng·h/mL ; durée d'effet clinique 12–13 h (FDA Vyvanse Prescribing Information ; Krishnan 2007, PMID 17696576 ; Pennick 2010, PMID 20520592 ; Boellner 2010, PMID 20206783). Particularités : linéarité dose–exposition stricte sur 30–70 mg, faible variabilité interindividuelle (CV ≈ 25 %), absence d'effet alimentaire significatif (un repas riche en graisse retarde T_max de 1 h sans modifier C_max/AUC), faible potentiel d'abus IV/intranasal car la prodrogue intacte est inactive (Jasinski 2009, PMID 19154100). Affichage : les courbes Elvanse/Xurta sont rapportées en d-amphétamine (énantiomère actif déjà obtenu après conversion enzymatique) ; la conversion d-MPH/total ne s'applique pas (option bypassDmphConversion). Les seuils DAT MPH (5,8 / 9,7 ng/mL) ne sont pas pertinents pour la d-AMP (cibles thérapeutiques d-AMP ≈ 30–50 ng/mL d'après Boellner 2010, mécanisme d'action mixte DAT + libérateur dopamine). Il est déconseillé de mélanger LDX et MPH dans une même stratégie : les unités plasmatiques (d-AMP vs d-MPH) ne sont pas comparables sur un axe Y commun.

Les trois formes LP à base de MPH (Concerta, Ritaline LP, Quasym) utilisent le même moteur double Weibull avec les paramètres de forme (t_d, s, ratio f_f) issus de Gomeni 2017 Table 1. Le facteur ×0,62 appliqué aux f_f est un facteur interne de réconciliation d'échelle — pas un paramètre physiologique de bioavailability : Gomeni utilise une convention où Σf_f = 1 avec V_d ≈ 1520 L explicite, alors que ce simulateur a un V_d implicite via amp = 70/W calibré sur Wargin 1983 (Ritaline LI 10 mg → 7 ng/mL). Le facteur 0,62 est commun aux trois formes LP MPH (pas un tuning ad hoc) et préserve les ratios f_f1/f_f2 publiés. Seule Medikinet conserve un modèle Bateman bi-pulse, la matrice gel pH-dépendante n'étant pas couverte par Gomeni 2017.

La fonction de Weibull décrit la fraction cumulée de dose libérée par un processus : F(t) = 1 − exp(−((t − t_lag) / t_d)^s). Avec s = 1, elle dégénère en exponentielle et la convolution avec l'élimination donne la formule de Bateman classique ; avec s > 1, elle produit une libération sigmoïde (latence puis accélération) caractéristique des pompes osmotiques et des coatings polymériques. La convolution est évaluée numériquement (méthode du point milieu, pas 0,05 h).

Option « Ajuster aux données cliniques »

Le simulateur expose deux modes de lecture des courbes plasmatiques, sélectionnables via la case à cocher « Modèle PK → Ajuster aux données cliniques » (décochée par défaut).

Mode brut (par défaut) — référence mécanistique Gomeni 2017

Les paramètres Weibull sont strictement issus de la Table 1 de Gomeni 2017 (modèle popPK développé sur Concerta puis qualifié sur plusieurs formulations LP). Les précisions paramétriques publiées concernent le modèle de référence et ne doivent pas être interprétées comme une validation individuelle formulation par formulation : pour Ritaline LP et Quasym, l'usage des paramètres Gomeni reste une extrapolation raisonnée à partir des moyennes publiées. Le facteur de réconciliation d'échelle ×0,62 est commun aux trois formes LP MPH et ne dépend pas du C_max cible : il découle de la calibration Wargin 1983 (Ritaline LI 10 mg → 7 ng/mL adulte 70 kg) qui fixe implicitement la combinaison F × V_d du simulateur. Les C_max simulés en mode brut tombent au centre du faisceau littéraire multi-études :

• Concerta 18 mg adulte : 4,80 ng/mL — entre Modi 2000 (3,7), FDA Concerta (4,4) et Markowitz 2003 (5,9).
• Ritaline LP 20 mg adulte : 5,57 ng/mL — proche du C_max2 DailyMed (~6,2) et de González 2002 (~5).
• Quasym 20 mg adulte : 6,26 ng/mL — entre Pentikis 2002 fasting (4,6) et SmPC EMC fed (6,4).

Ce mode est recommandé pour la lecture pédagogique : il préserve la cohérence mécanistique (mêmes paramètres Gomeni pour les 3 LP) et expose la moyenne attendue d'une population publiée plutôt qu'une étude pivot particulière.

Mode ajusté — calibration sur l'étude pivot par formulation

Quand la case est cochée, l'amplitude C(t) est multipliée par un facteur empirique amp propre à chaque formulation, choisi pour aligner le C_max simulé sur la moyenne ou l'étude pivot d'une référence ciblée :

Ce que l'ajustement modifie — et ce qu'il ne modifie pas

• Amplitude C(t) → multipliée par amp : C_max et AUC sont rééchelonnés du même facteur (l'AUC ajustée est donc plus basse que l'AUC publiée quand la cible C_max est en deçà du brut, par construction).
• Forme temporelle : strictement préservée. T_max, ratios C_max1/C_max2, durée d'action et profil biphasique restent identiques entre brut et ajusté.
• Allométrie poids/sexe et variabilité interindividuelle : appliquées après la calibration ajustée. Le facteur amp ne perturbe pas les Monte Carlo ni la pharmacogénétique CES1.
• Décomposition par composantes IR/ER : préservée — l'ajustement multiplie chaque processus Weibull du même facteur, sans modifier les fractions f_f relatives.

Limites et précautions

• Calibrer sur une étude pivot peut s'éloigner d'autres références. Concerta 18 mg ajusté à 3,73 ng/mL (Modi) sous-estime FDA (-15 %) et Markowitz LC-MS/MS (-37 %). En mode brut le simulateur tombe au centre du faisceau ; en mode ajusté il calibre sur Modi seul.
• L'AUC ajustée n'est pas validée comme cible : elle suit C_max. Pour des comparaisons d'exposition cumulée, le mode brut reste préférable.
• L'ajustement est par formulation, pas par dose ni par population : un seul facteur amp par molécule. Les écarts pédiatriques (Wigal 2011 fasting) restent partiellement non corrigés.
• Aucun fitting individuel : l'ajustement ne remplace pas un dosage plasmatique réel. Le simulateur reste un outil pédagogique.

Quand utiliser quoi ?

• Brut : explorations pédagogiques, comparaison entre formulations, lecture de la moyenne attendue dans une population, démonstration de la mécanistique galénique.
• Ajusté : reproduction de l'étude pivot d'une formulation donnée (typiquement Modi 2000 pour Concerta, lecture de notice ou de littérature focalisée), validation visuelle contre une référence unique citée.

Validation externe contre la littérature

Cette section compare les valeurs publiées au modèle clinique brut et à son option d'ajustement clinique. Les comparaisons MPH sont exprimées en MPH total racémique afin de rester cohérent avec la majorité des études pivots ; la lisdexamfétamine est rapportée en d-amphétamine. Dose unique à 8 h, adulte 70 kg sauf mention contraire.

Lecture importante : le moteur brut (case « Ajuster aux données cliniques » décochée — état par défaut) applique strictement les paramètres Weibull Gomeni 2017 Table 1 pour les 3 LP MPH, sans tuning. L'option ajustée multiplie l'amplitude C(t) par un facteur empirique par formulation pour caler le C_max sur une étude pivot citée (Modi 2000 pour Concerta, DailyMed C_max2 pour Ritaline LP, SmPC pour Medikinet). La forme temporelle (T_max, profil biphasique, ratios IR/ER) est strictement préservée dans les deux modes.

Légende : vert = cohérent avec la plage publiée ; orange = acceptable mais à interpréter avec prudence ; rouge = limite connue ou comparaison non directement interprétable.

Conclusion de validation : le moteur brut Gomeni 2017 (3 LP en double Weibull, Medikinet en Bateman bi-pulse) tombe au centre du faisceau littéraire multi-études pour Concerta (4,80 entre Modi 3,7 / FDA 4,4 / Markowitz 5,9), Ritaline LP (5,57 vs DailyMed C_max2 6,2) et Quasym (6,26 entre Pentikis 4,6 et SmPC 6,4). L'option « Ajuster aux données cliniques » est utile pour caler une simulation sur une étude pivot précise (Modi 2000 pour Concerta, DailyMed C_max2 pour Ritaline LP, etc.) : elle rapproche le C_max de la cible citée mais réduit l'AUC du même facteur. Elle reste donc optionnelle : la lecture mécanistique brute est plus fidèle aux moyennes inter-études, l'ajustement plus fidèle à une référence ponctuelle.

Synthèse par moyenne clinique multi-études

Pour les formulations où plusieurs études rapportent le même paramètre, cette table résume la moyenne clinique et la plage observée. Elle ne remplace pas les comparaisons détaillées ci-dessus : elle sert de lecture rapide quand les études sont hétérogènes.

Code couleur des valeurs moteur : fond vert = proche de la moyenne ou dans la plage clinique ; fond jaune = écart modéré ; fond rouge / fond rouge = hors centre de la plage ou biais visible.

Références conservées pour cette validation

• Wargin W, Patrick K, Kilts C, et al. Pharmacokinetics of methylphenidate in man, rat and monkey. J Pharmacol Exp Ther 1983;226:382-386. PMID 6410043.
• Kimko HC, Cross JT, Abernethy DR. Pharmacokinetics and clinical effectiveness of methylphenidate. Clin Pharmacokinet 1999;37:457-470. PMID 10628897.
• Modi NB, Lindemulder B, Gupta SK. Single- and multiple-dose pharmacokinetics of an oral once-a-day osmotic controlled-release OROS methylphenidate formulation. J Clin Pharmacol 2000;40:379-388. PMID 10761165.
• Markowitz JS, Straughn AB, Patrick KS. Pharmacokinetics of methylphenidate after oral administration of two modified-release formulations in healthy adults. Clin Pharmacokinet 2003;42(4):393-401. PMID 12648029. DOI 10.2165/00003088-200342040-00007.
• González MA, Pentikis HS, Anderl N, et al. Methylphenidate bioavailability from two extended-release formulations. Int J Clin Pharmacol Ther 2002;40:175-184. PMID 11996212.
• Pentikis HS, Simmons RD, Benedict MF, Hatch SJ. Methylphenidate bioavailability in adults when an extended-release multiparticulate formulation is administered sprinkled on food or as an intact capsule. J Am Acad Child Adolesc Psychiatry 2002. PMID 11931601.
• Rochdi M, González MA, Hirshey Dirksen SJ. Dose-proportional pharmacokinetics of a methylphenidate extended-release formulation. Int J Clin Pharmacol Ther 2004;42:285-292.
• Wang Y, Lee L, Somma R, et al. In vitro dissolution and in vivo oral absorption of methylphenidate from extended-release formulations. Biopharm Drug Dispos 2004. PMID 14872557.
• Schütz H, Hiemke C, Doepfmer S, et al. Lack of bioequivalence between two methylphenidate extended-release formulations. Int J Clin Pharmacol Ther 2009;47(12):761-769. PMID 19954715. DOI 10.5414/CPP47761.
• Wigal SB, Kollins SH, Childress AC, Squires L. A 13-hour laboratory school study of lisdexamfetamine dimesylate in school-aged children with ADHD. Child Adolesc Psychiatry Ment Health 2009;3:17. PMID 19508726 ; and Wigal SB et al. OROS methylphenidate food-effect / PK study. J Child Adolesc Psychopharmacol 2011;21:127-133. PMID 21663428. DOI 10.1089/cap.2010.0083.
• Gomeni R, Bressolle-Gomeni F, Spencer TJ, Faraone SV, Fang L, Babiskin A. Model-based approach for optimizing extended-release methylphenidate formulations. Clin Pharmacol Ther 2017;102:951-960. PMID 28369788. DOI 10.1002/cpt.684.
• Krishnan SM, Pennick M, Stark JG. Metabolism, distribution and elimination of lisdexamfetamine dimesylate. Clin Drug Investig 2008;28:745-755. PMID 17696576.
• Boellner SW, Stark JG, Krishnan S, Zhang Y. Pharmacokinetics of lisdexamfetamine dimesylate and its active metabolite d-amphetamine in children with ADHD. Clin Ther 2010;32:252-264. PMID 20206783.
• Ermer JC, Dennis K, Haffey MB, Doll WJ, Sandefer EP, Buckwalter M. Intranasal and oral lisdexamfetamine dimesylate: pharmacokinetics and pharmacodynamics. Clin Drug Investig 2011-2012.
• Ritalin LA, Metadate CD, Concerta, Vyvanse and Medikinet / Equasym XL official prescribing information or SmPC documents.

Mise à l'échelle allométrique (poids du patient)

Si un poids est saisi, le modèle applique un ajustement allométrique aux paramètres pharmacocinétiques. Cette section explique pourquoi et comment.

Pourquoi un V_d plus petit donne une C_max plus élevée

L'équation de base du modèle monocompartimental après dose orale D s'écrit :

C(t) = (D · F / V_d) × (k_a / (k_a − k_e)) × (e^−k_e·t − e^−k_a·t)

Le facteur D · F / V_d est l'amplitude : il fixe la hauteur globale de la courbe. Le terme entre parenthèses ne dépend que du temps et des constantes de vitesse, pas du sujet. À dose constante, C_max est donc inversement proportionnelle à V_d. Si V_d est divisé par deux, C_max est doublée — c'est de la dilution physique simple : la même quantité de molécule active est répartie dans un volume plus petit chez un patient plus léger.

Comment V_d dépend du poids — base allométrique

D'après la théorie allométrique (Anderson & Holford 2008, PMID 17914927) appliquée à la pharmacocinétique de population (Yang 2014, PMID 24515100 ; Mahmood 2013, PMID 23588537) :

V_d(W) = V_d,70 · (W / 70)^α_V      CL(W) = CL₇₀ · (W / 70)^α_CL

Avec :

• α_V ≈ 1,0 — V_d proportionnel au poids. L'eau corporelle totale et les tissus de distribution suivent le poids quasi linéairement.
• α_CL ≈ 0,75 — clairance selon la loi 3/4 de Kleiber. La clairance hépatique et rénale (par métabolisme et filtration) augmente moins vite que le poids, suivant la loi métabolique générale des mammifères.
• Conséquence : t½ ∝ V_d/CL ∝ W^0,25 — la demi-vie augmente très modestement avec le poids (≈ 2,3 h chez un enfant de 30 kg, ≈ 3,0 h chez un adulte de 100 kg, cohérent avec Wargin 1983, PMID 6410043).

Paramètres allométriques retenus pour le méthylphénidate

Cas particulier de la lisdexamfétamine

L'allométrie standard α_V = 1,0 surestime l'amplification de la C_max chez l'enfant pour la LDX. La cause est mécanistique : l'hydrolyse de la prodrogue par la carboxypeptidase érythrocytaire est saturable (cinétique michaelienne, K_m dépassé aux doses cliniques usuelles 30–70 mg). La masse globulaire totale ne limite donc pas linéairement la conversion avec le poids ; elle plafonne. Données mesurées :

• Boellner 2010 (PMID 20206783) — enfant 6–12 ans (~30 kg), LDX 30 mg : C_max d-AMP ≈ 53,2 ng/mL.
• Krishnan 2007 (PMID 17696576) — adulte 70 kg, LDX 30 mg : C_max d-AMP ≈ 33 ng/mL.
• Ratio enfant/adulte ≈ 1,6 (mesuré) vs ≈ 2,3 (prédit par α_V = 1,0).
• Ermer 2010 (PMID 21090838) — confirme : LDX présente la plus faible variabilité PK liée au poids de tous les stimulants longue durée.

Le simulateur applique donc à la classe LDX un scaling allométrique atténué via le paramètre weightScalingFactor = 0,15, qui transforme l'amplitude selon la formule :

amp_LDX = 1 + 0,15 × ((70/W) − 1)

Cela correspond à un exposant α_V apparent ≈ 0,2 sur la plage 25–70 kg. Le ratio simulé enfant/adulte (~1,2 à 30 kg) est conservateur par rapport à Boellner 2010 (PMID 20206783, ratio observé ~1,6), ce qui convient à un outil pédagogique évitant de sur-estimer les concentrations pédiatriques. Ermer 2010 La k_e override (0,063 h⁻¹) reste soumise au k_Ratio standard (α_CL = 0,75), reproduisant l'évolution de la t½ d-AMP avec l'âge (10,4 h adulte → 8,6 h pédiatrique ; FDA Vyvanse PI).

Synthèse des paramètres et leurs bornes

• Amplitude de la courbe MPH ∝ 70 / poids (α_V = 1,0), borné [0,4 ; 3,0].
• Amplitude de la courbe LDX = 1 + 0,15 × ((70/W) − 1), borné par la même contrainte sur (70/W).
• Constante d'élimination k_e × (poids/70)^−0,25 (α_CL = 0,75), borné [0,7 ; 1,3]. Appliqué uniformément à toutes les molécules.
• Les bornes protègent des extrapolations extrêmes : nourrissons < 10 kg (immaturité enzymatique CES1 non modélisée) et obèses > 200 kg (V_d ne suit plus le poids total mais plutôt la masse maigre LBW).

Limites connues du modèle allométrique simple

1. Maturation enzymatique CES1 : chez le jeune enfant (< 5 ans), la carboxylestérase 1 (qui hydrolyse le MPH en acide ritalinique) est immature. La clairance réelle est inférieure à la prédiction α_CL = 0,75. Non modélisé (Yang 2014, PMID 24515100).
2. Composition corporelle (obésité) : chez les sujets très obèses, V_d ne suit plus le poids total mais plutôt la masse maigre. Le modèle surestime alors V_d et sous-estime C_max.
3. Polymorphismes CES1 : variants comme G143E peuvent multiplier C_max par 2–4× indépendamment du poids (Stage 2017, PMID 28087982 ; van der Schans 2023, PMID 37851403). Modélisés via la mixture pharmacogénétique activable dans la barre des options (cf. section « Variabilité interindividuelle » ci-dessous).
4. Co-administration éthanol : inhibition CES1 + transestérification → augmentation C_max +40 % et AUC +25 % (Markowitz 2009, PMID 20533759). Non modélisé.

Variabilité interindividuelle (Monte Carlo)

Le simulateur peut afficher une bande de variabilité interindividuelle conforme aux conventions de pharmacométrie de population (popPK). Trois modes sont disponibles dans la barre d'options :

Mode « Aucune » (défaut)

Affichage de la courbe typique de la population, ajustée par les covariables connues : poids (allométrie), sexe (multiplicateur ×1,12 chez la femme pour le MPH ; pas d'effet pour LDX). Elle sert de repère pédagogique pour les données saisies, sans prédire le profil individuel réel.

Mode « Standard » (log-normale homogène)

Simulation Monte Carlo de N = 500 patients virtuels, paramètres pharmacocinétiques tirés selon une loi log-normale corrélée. Les coefficients de variation (CV) utilisés sont des ordres de grandeur popPK (incluant implicitement une part attribuable à CES1 mais sans la séparer) :

Une matrice de corrélation est appliquée par décomposition de Cholesky pour respecter les corrélations physiologiques entre paramètres (ρ ≈ 0,4 entre V_d et CL via la perfusion hépatique commune ; Sun 2015).

Mode « Pharmacogénétique CES1 » (mixture pure)

Ce mode affiche uniquement une approximation de la contribution d'un phénotype CES1 réduit, sans superposer la variance log-normale. Il est conçu de façon orthogonale au mode standard pour permettre une lecture pédagogique isolée de l'effet pharmacogénétique. Mixture trimodale illustrative :

Sources : Stage 2017 (PMID 28087982) ; Shi 2022 (PMC9391289) ; Patrick 2007 (PMC3188424). Fréquences phénotypiques agrégées, pas uniquement G143E.

Aspect visuel attendu du mode CES1 : la bande affichée est discrète à 3 niveaux plutôt que continue. Les percentiles P5 à P75 sont tous superposés sur la courbe typique (correspondant aux métaboliseurs normaux), le P95 est décalé à × 1,8 (zone phénotype CES1 intermédiaire illustratif), et le maximum atteint × 3,2 (zone phénotype CES1 lent illustratif). Cette représentation « en marches » est volontaire : elle isole un mécanisme pharmacogénétique plausible sans le confondre avec une prédiction génotypique individuelle.

Comparaison entre les deux modes :

• Standard : bande continue (variabilité agrégée, log-normale), montre la dispersion réelle attendue dans une population réelle
• CES1 : bande discrète (mixture pure), montre uniquement ce qu'apporte spécifiquement le polymorphisme CES1 G143E, isolé de toute autre variance

Pour avoir la simulation populationnelle complète (variance log-normale + mixture CES1 superposées), il faudrait un mode hybride non actuellement implémenté. Les deux modes orthogonaux doivent donc être lus comme deux éclairages pédagogiques séparés : le standard pour l'ordre de grandeur de dispersion, le CES1 pour la contribution pharmacogénétique potentielle.

Mode CES1 pour la lisdexamfétamine : non applicable. La conversion de la prodrogue par carboxypeptidase érythrocytaire ne dépend pas de CES1 et n'a pas de polymorphisme cliniquement significatif documenté (Krishnan 2007, Ermer 2010). Si vous activez le mode CES1 avec une stratégie LDX/Xurta, aucune bande ne sera affichée.

Pourquoi le simulateur ne demande pas votre génotype

Le génotypage CES1 n'est pas pratiqué en routine en France (non remboursé, non recommandé en première intention). Le simulateur illustre donc une distribution pharmacogénétique populationnelle : chaque patient virtuel se voit affecter un phénotype au hasard selon les fréquences ci-dessus. La bande affichée est une aide de lecture populationnelle, pas une assignation de génotype au patient réel.

Si vous suspectez cliniquement un phénotype atypique (sur-réponse à faible dose, effets indésirables marqués, exposition inattendue), ce mode peut aider à discuter une hypothèse CES1 parmi d'autres. Cette information ne remplace pas l'évaluation clinique ; elle peut orienter la réflexion vers une molécule contournant CES1 (lisdexamfétamine) ou vers un dosage plasmatique de confirmation lorsque celui-ci est disponible.

Stabilité des bandes — seed déterministe

Le générateur de nombres pseudo-aléatoires (Mulberry32) est seed depuis un hash des inputs (poids, sexe, mode, doses, horaires). Tant que vos entrées sont identiques, les bandes sont parfaitement reproductibles entre rafraîchissements. Modifier une dose ou un horaire change le seed et donc la simulation.

Lecture des stats enrichies

Quand un mode de variabilité est actif, le panneau de stats affiche les valeurs sous le format médiane (P25–P75) :

• P50 (médiane) — valeur centrale de la distribution simulée
• P25–P75 (intervalle interquartile) — 50 % des patients de la population sont dans cette plage
• Les bandes graphiques (fines au centre, plus larges aux extrêmes) correspondent à la P5–P95 (90 % de la population)

Limites du modèle Monte Carlo

1. Indépendance log-normale par paramètre (corrélations modérées via Cholesky, mais pas de corrélations pleines paramètre × covariable).
2. Pas de variabilité inter-occasion (IOV) : le modèle suppose une stabilité jour-à-jour pour le même patient.
3. Pas d'erreur de mesure plasmatique (ε résiduelle) : le simulateur affiche l'exposition vraie, pas la concentration mesurée par dosage.
4. Mode « pharmacogénétique CES1 » non applicable à la lisdexamfétamine (carboxypeptidase érythrocytaire indépendante de CES1).
5. Mélange MPH + LDX dans une même stratégie : la variabilité n'est pas affichée (unités hétérogènes d-MPH vs d-AMP).

Bibliographie spécifique au modèle de variabilité

• Sun L, Mills R, Sadler BM, Rege B. Population pharmacokinetics of methylphenidate ER beads in pediatric ADHD. Drugs R D 2015;15:169-76. PMID 26060393. — CV V_d 30 %, CV CL 37 %, matrice de corrélation popPK pédiatrique.
• Nemec K, Schubert-Zsilavecz M, Abdel-Tawab M. Population pharmacokinetics of methylphenidate in adults emphasizing CES1 variants. Clin Transl Sci 2017;10:276-84. PMC5351003. — CV CL adulte ~45 %, dont ~50 % attribuable à CES1.
• Stage C, Jürgens G, Guski LS, et al. The impact of CES1 genotypes on methylphenidate pharmacokinetics. Br J Clin Pharmacol 2017;83:1506-14. PMID 28087982. — effet majeur chez les porteurs de variants CES1 réducteurs ; fréquence allélique G143E ~1,6 % en population européenne ; les multiplicateurs ×1,8 / ×3,2 utilisés ici sont des approximations illustratives.
• Shi J, Wang X, Eyler RF, et al. PBPK modeling to predict methylphenidate exposure affected by CES1 pharmacogenetics, drug-drug interactions, and sex. J Pharm Sci 2022;111:3109-20. PMC9391289. — synthèse intégrée variance × covariables ; analyse Monte Carlo PBPK validée.
• Patrick KS, Straughn AB, Minhinnett RR, et al. Influence of ethanol and gender on methylphenidate pharmacokinetics. Clin Pharmacol Ther 2007;81:346-53. PMID 17698715. — effet sexe ×1,12 sur Cmax/AUC chez la femme, MPH.
• van der Schans J, Çiçek R, Vardar S, et al. High-dose methylphenidate and CES1 genetic variability: case series. J Clin Psychopharmacol 2023;43:520-4. PMID 37851403. — illustrations cliniques de la variabilité G143E.
• Anderson BJ, Holford NH. Mechanism-based concepts of size and maturity in pharmacokinetics. Annu Rev Pharmacol Toxicol 2008;48:303-32. PMID 17914927. — référence théorique allométrie + matrices de variance-covariance.
• Beal SL, Sheiner LB. NONMEM User's Guides. Icon Development Solutions, 1992. — convention popPK standard pour log-normale + Cholesky.

Repères plasmatiques d'efficacité (occupation DAT)

Les deux lignes horizontales dérivent des études PET d'occupation du transporteur de la dopamine (DAT) striatal. Volkow 1998 (PMID 9766762) rapporte les concentrations en MPH total racémique (méthode antérieure à la généralisation des dosages chiraux) ; Spencer 2006 (PMID 16513858) rapporte les concentrations en d-MPH. Pour homogénéiser l'affichage, le simulateur stocke les seuils en MPH total et applique la conversion plasmatique d-MPH ≈ 0,97 × total quand l'unité d-MPH est affichée. L'écart pratique entre les deux conventions est faible (~3 %), ce qui rend les seuils visuellement comparables, mais la coexistence de deux standards analytiques dans la littérature DAT-PET reste à garder à l'esprit. Ce sont des repères pharmacodynamiques, pas des seuils universels d'efficacité clinique :

• 5,8 ng/mL d-MPH (équivalent 6 ng/mL MPH total) ≈ 50 % d'occupation DAT — repère d'effet pharmacologique apparent (Volkow 1998, PMID 9766762).
• 9,7 ng/mL d-MPH (équivalent 10 ng/mL MPH total) ≈ 80 % d'occupation DAT — repère de forte occupation DAT, au-delà duquel le gain clinique peut plafonner et les effets indésirables devenir plus probables (Spencer 2006, PMID 16513858).

Pourquoi les seuils sont indépendants du poids

Les seuils plasmatiques DAT (Volkow 1998, Spencer 2006) sont des propriétés intrinsèques du couple molécule–récepteur : 6 ng/mL de MPH total dans le plasma produisent ~50 % d'occupation DAT, que le sujet pèse 30 kg ou 100 kg. Le simulateur conserve donc des seuils fixes indépendamment du poids saisi.

Cette convention conserve le comportement attendu du modèle : pour une même dose, un patient plus léger atteint des concentrations plasmatiques plus élevées (V_d ∝ poids), sa courbe dépasse plus largement le repère fixe, et la durée au-dessus du repère augmente. Exemple Ritaline LP 18 mg + 10 mg :

• Adulte 60 kg : C_max ≈ 12,0 ng/mL d-MPH ; T > 5,8 ng/mL ≈ 10,8 h.
• Enfant 38 kg : C_max ≈ 18,0 ng/mL d-MPH (+50 %) ; T > 5,8 ng/mL ≈ 12,4 h (+1,6 h).

Cette lecture évite l'option historique (repères également scalés par amp) qui aurait masqué l'effet de la dilution dans V_d.

Seuils spécifiques à la lisdexamfétamine (d-amphétamine)

Pour la classe amphétamine (Elvanse / Xurta), les seuils MPH-DAT ne sont pas pertinents : la d-amphétamine agit par mécanisme mixte (inhibiteur DAT/NET + libérateur dopaminergique vésiculaire VMAT-2), sans relation simple PET-DAT/effet. Les seuils retenus sont empiriques, dérivés des études exposition-réponse cliniques en ADHD :

• 30 ng/mL d-AMP — seuil thérapeutique. C_max pédiatrique ≈ 53,2 ng/mL avec Elvanse 30 mg associée à un effet clinique significatif sur l'échelle SKAMP-Combined (Boellner 2010, PMID 20206783). La durée d'efficacité corrèle avec le temps passé au-dessus de ce seuil (Ermer 2016, PMID 27021968).
• 80 ng/mL d-AMP — saturation clinique pratique. Correspond approximativement à la C_max de la dose maximale approuvée chez l'adulte (Elvanse 70 mg, C_max ≈ 75 ng/mL ; FDA Vyvanse PI). L'analyse exposition-réponse en population pédiatrique japonaise (Tsuda 2020, PMID 33082099) montre un plateau d'efficacité dès l'exposition minimale, sans gain proportionnel au-delà.

Le basculement automatique entre seuils MPH et d-AMP se fait selon les molécules actives dans la stratégie : si une stratégie mélange MPH et LDX, les seuils ne sont pas affichés (mélange d'unités d-MPH et d-AMP non comparables).

Adaptation au poids spécifique à Xurta / Elvanse (scaling atténué)

Contrairement au méthylphénidate (Vd ∝ poids^1, scaling allométrique standard), la pharmacocinétique clinique de la lisdexamfétamine est peu dépendante du poids. Deux mécanismes l'expliquent :

1. L'hydrolyse de la prodrogue par la carboxypeptidase érythrocytaire est saturable (cinétique de Michaelis–Menten). À doses cliniques (30–70 mg), le système est partiellement saturé chez l'adulte ; chez l'enfant, la masse globulaire totale réduite ne diminue pas proportionnellement la vitesse de conversion. Cela aplatit la dépendance au poids attendue par scaling allométrique pur.
2. La liaison aux protéines plasmatiques de la d-amphétamine est très faible (~16–20 %), réduisant la sensibilité de la fraction libre aux variations de composition corporelle (Pennick 2010, PMID 20520592).

Données cliniques de référence :

• Boellner 2010 (PMID 20206783) — enfants 6–12 ans (~30 kg) avec Elvanse 30 mg : C_max d-AMP ≈ 53,2 ng/mL.
• Krishnan 2007 (PMID 17696576) — adultes (~70 kg) avec LDX 30 mg : C_max d-AMP ≈ 33 ng/mL.
• Ratio enfant/adulte observé ≈ 1,6, alors que le scaling allométrique pur (70/W) prédirait ≈ 2,3.

Implémentation dans le simulateur : la classe LDX utilise weightScalingFactor = 0,15, qui atténue le facteur d'amplitude selon la formule : amp = 1 + 0,15 × ((70/W) − 1). Concrètement :

Cette correction est cohérente avec la posologie officielle Xurta / Elvanse (FDA Vyvanse PI ; RCP ANSM) qui n'est pas exprimée en mg/kg : la dose initiale est de 30 mg/jour à tout âge (à partir de 6 ans), titrée jusqu'à 70 mg maximum selon la réponse clinique et la tolérance, indépendamment du poids. Pour le méthylphénidate au contraire, la posologie pédiatrique est usuellement calculée en mg/kg (typiquement 0,5–1,5 mg/kg/jour MPH total), ce qui justifie le scaling allométrique standard pour cette classe.

k_e et demi-vie d-AMP — la t½ d-AMP varie modérément avec l'âge (10,4 h adulte vs 8,6 h pédiatrique selon FDA Vyvanse PI ; rapport 1,21). Le simulateur applique le facteur kRatio = (W/70)^(−0,25) à la k_e override de 0,063 h⁻¹, ce qui reproduit cette variation : à 30 kg, k_e = 0,075 h⁻¹ → t½ = 9,2 h, conforme à Boellner 2010 (PMID 20206783).

Facteurs de variabilité non modélisés

La variabilité interindividuelle réelle peut atteindre ±50 %. Les principaux déterminants identifiés dans la littérature mais non intégrés, ou seulement illustrés de façon simplifiée, sont les suivants.

Sur l'absorption

• Alimentation. Un repas augmente l'AUC du d-MPH (IR et SR) sans modifier significativement la vitesse d'absorption (Chan 1983, PMID 11587491). Pour l'OROS Concerta, l'effet est cliniquement négligeable (Auiler 2002, PMID 12240794). Cas particulier du Medikinet : la notice impose une prise avec petit-déjeuner car la matrice gel a besoin du bol alimentaire pour ralentir la libération ; à jeun, le profil ER s'estompe et se rapproche d'un profil IR.
• pH gastro-intestinal et microbiote. Le MPH est stable de pH 2 à 8, mais sa biodisponibilité luminale dépend du pH intestinal, lui-même modulé par le régime et la flore (Khaled 2021, PMID 34451830). Effet individuel non prédictible.
• Altération de la galénique. Broyer ou ouvrir une matrice OROS détruit le mécanisme de libération prolongée et produit un profil type IR avec risque de Cmax très élevé.
• Antiacides / IPP. En modifiant le pH gastrique, ils peuvent accélérer la libération des formes ER à matrice pH-dépendante (principalement Medikinet).

Sur l'élimination

• Polymorphisme CES1 (facteur majeur). La carboxylestérase 1 métabolise ~80 % du MPH en acide ritalinique inactif. Le variant rs71647871 (G143E) — fréquence 3–4 % chez les Caucasiens — réduit l'activité enzymatique : les porteurs hétérozygotes présentent une AUC du d-MPH augmentée de +149 % (Stage 2017, PMID 28087982). D'autres variants et diplotypes CES1 ont été identifiés comme covariables significatives (Nemec 2017, PMC5351003 ; Shi 2022, PMC9391289 ; van der Schans 2023, PMID 37851403).
• Éthanol. L'éthanol inhibe CES1 et induit une transestérification produisant de l'éthylphénidate (métabolite actif). La co-ingestion augmente la C_max du d-MPH d'environ +40 % et l'AUC de +25 % (Markowitz 2009, PMID 20533759 ; Patrick 2007, PMC3188424), avec renforcement de l'effet euphorisant.
• Sexe. AUC ~20 % plus élevée chez la femme (Patrick 2007, PMC3188424), attribuée à une activité CES1 modulée par les hormones sexuelles.
• Interactions médicamenteuses. Contre-indication absolue avec les IMAO (risque de crise hypertensive). Effets modestes avec les ISRS/IRSN via compétition pour CES1 (fluoxétine surtout).
• Fonction hépatique. Peu documentée ; pas d'ajustement posologique standardisé. Le métabolite (acide ritalinique) est excrété à ~78–97 % dans les urines, donc la fonction rénale influence peu le d-MPH actif.