PK voie orale, équations différentielles, métabolite (1er ordre, MM)
ONGOING
Rédaction en cours…
parent (1 cpt) + métabolite (1 cpt) via 1er ordre
$PROB
…/…
$DATA …CMT…
;1 pour amount
;2 pour parent et 3 pour métabolite
;parent et metabolite avec même unité de concentration (mole)
$SUBROUTINE ADVAN13 TRANS1 TOL=3 ; previously ADVAN6
$MODEL NCOMP=3
COMP=(DIG, DEFDOSE)
COMP=(CENTRAL, DEFOBS)
COMP=(METAB)
$PK
V2 = popV2*EXP(ETA(1))
…/…
S2=V2
S3=V2
; La quantité de métabolite et son volume ne sont pas identifiables simultanément.
; => FIXER V3 à une valeur raisonnable
; V3=V2 OU
; S3=1 ; equiv. V3=1L
; CL=V2*(K20+K23)
; Fraction de metabolisation apparente
; FMET=K20/(K20+K23)
; On peut préférer estimer FMET ou K23
; mais => V3 soit fixé
; K23 et V3 ne sont pas dissociables
; OU fixer la fraction de dose convertie
$DES
C2 = A(2)/S2
C3 = A(3)/S3
DADT(1) = -KA*A(1)
DADT(2) = KA*A(1) – K23*A(2) – K20*A(2)
DADT(3) = K23*A(2) – K30*A(3)
$THETA
(50.0) ; popV2 ; L
( 0.1) ; popKA ; 1/h
( 0.1) ; popK23 ; 1/h
( 0.1) ; popK20 ; 1/h
( 0.1) ; popK30 ; 1/h
(50.0) ; popV3 ; L
(0, .1) ; PROP
(0, 50.0) ; ADD
$ERROR
CP = A(2)/S2
CM = A(3)/S3
IF (CMT.EQ.2) THEN
Y = F*EXP(ERR(1))
ELSE
Y = F*EXP(ERR(2))
ENDIF
IPRED=F
IRES=DV-IPRED
IWRES=IRES/W
…/…
[b]parent (1 cpt) + métabolite (1 cpt) via Michaelis Menten[/b]
$PROB parent + metabolite MM
…/…
$DATA …CMT…
;1 pour amount
;2 pour parent et 3 pour métabolite
;parent et metabolite avec même unité de concentration (mole)
$SUBROUTINE ADVAN6 TRANS1 TOL=3
$MODEL NCOMP=3
COMP=(DIG, DEFDOSE)
COMP=(CENTRAL, DEFOBS)
COMP=(METAB)
$PK
V2 = popV2*EXP(ETA(1))
…/…
S2=V2
S3=V2
; V2=V3 (V3 non estimable)
; OR S3=1 with V3=1L
$DES
C2 = A(2)/S2
C3 = A(3)/S3
DADT(1) = -KA*A(1)
DADT(2) = +KA*A(1) – VMAX*C2/(KM+C2) – K20*A(2)
DADT(3) = +VMAX*C2/(KM+C2) – K30*A(3)
$THETA
(50.0) ; popV2 ; L
( 0.1) ; popKA ; 1/h
( 0.1) ; popVMAX; A/h
( 0.1) ; popKM ; A/L
( 0.1) ; popK20 ; 1/h
( 0.1) ; popK30 ; 1/h
(50.0) ; popV3 ; L
(0, 0.1) ; PROP
(0,50.0) ; ADD
$ERROR
CP = A(2)/S2
CM = A(3)/S3
IF (CMT.EQ.2) THEN
Y = F*EXP(ERR(1))
ELSE
Y = F*EXP(ERR(2))
ENDIF
IPRED=F
IRES=DV-IPRED
IWRES=IRES/W
…/…