inst/modeldb/pharmacodynamics/indirect_1cpt_inhi_kin.R

In nlmixr2lib: A Model Library for 'nlmixr2'

indirect_1cpt_inhi_kin  <- function() {
  description <- "One compartment indirect response model with inhibition of kin."
  ini({
    lka  <- 0.45 ; label("Absorption rate (Ka)")
    lvc  <- log(90) ; label("Central volume of distribution (Vc)")
    lkel <- 0.534; label("elimination rate (1/d)")
    lIC50 <- 0.67; label("Drug concentration producing 50% of maximum inhibition at effect site (IC50)")
    limax <- 0.56; label("Maximum inhibitory factor attributed to drug (Imax)")
    lkin <- 0.48; label("Zero-order rate constant for production of drug response(1/d)")
    lkout <- 0.34; label("First-order rate constant for loss of drug response")
    lfdepot <- 0.4; label("Bioavailability (F)")
    propSd <- 0.5 ; label("Proportional residual error (fraction)")
  })
  model({
    ka  <- exp(lka)
    vc  <- exp(lvc)
    kel <- exp(lkel)
    IC50 <- exp(lIC50)
    imax <- exp(limax)
    kin <- exp(lkin)
    kout <- exp(lkout)
    fdepot   <- exp(lfdepot)
    
    d/dt(depot)      <- -ka*depot
    f(depot)         <- fdepot
    d/dt(central)    <- ka*depot -kel*central
    Cc <-  central/vc
    
    d/dt(effect) <- kin*(1-imax*Cc/(Cc + IC50)) - kout*effect
    Cc ~ prop(propSd)
  })
}