Individualized sample timing for PK studies using accessible PD data - - PowerPoint PPT Presentation

Individualized sample timing for PK studies using accessible PD data

Matthew S. Shotwell, Ph.D.

Department of Biostatistics Vanderbilt University School of Medicine Nashville, TN, USA

August 29, 2017

Introduction PK Studies

◮ serial blood sampling in individual after dose ◮ goal: estimate change in concentration in body over time ◮ timing of samples affects statistical precision ◮ efficient timing depends on uncertain parameters ◮ can use other information to reduce uncertainty

Muscle relaxants

◮ muscle relaxants (paralytics) often part of general anesthesia ◮ facilitates intubation by relaxing muscles around larynx ◮ prevents involuntary movement during surgery ◮ large initial intubating dose (0.1mg/kg) ◮ small maintenance dose every 30min (0.01mg/kg) ◮ significant PK/PD heterogeneity ◮ requires monitoring

Why study PK?

Better understanding of PK heterogeneity may:

◮ facilitate individualized dosing ◮ reduce adverse events (residual paralysis; aspiration) ◮ reduce time in OR and cost

Why optimize PK study design?

◮ maximize information in sample ◮ reduce number of subjects or samples ◮ design task is relatively cheap

Generic optimal design

• 1. identify optimal sample times a priori
• 2. sample all subjects at same times

Individualized optimal design

◮ optimize sample times for each subject ◮ need additional information about subject

Monitoring paralysis

Muscle paralysis monitored by:

◮ electrical stimulation of ulnar nerve ◮ train of four stimuli (TOF) stimuli ◮ count twitches (TOF count)

Monitoring paralysis

Monitoring paralysis

Generic vs. individualized

Generic:

• 1. identify optimal sample times a priori
• 2. sample all subjects at same times

Individualized:

• 1. monitor and record twitch counts
• 2. use twitch counts to individualize sample times
• 3. sample subject at optimal times

PK/PD model

Three compartment “biophase” model: dm1 dt = −k10m1 − k12m1 + k21m2 − k13m1 + k31m3 dm2 dt = k12m1 − k21m2 dm3 dt = k13m1 − k31m3

◮ m1 - amount of drug in blood (c1 = m1/v1) ◮ m2 - amount of drug in tissues ◮ m3 - amount of drug in “biophase” or “effect” compartment ◮ aim: estimate kinetic parameters θ = {k10, k12, k21, k13, k31}

• Stat. model: twitch count

Proportional odds logistic regression (POLR) model: logit[Pr(zil ≤ ν)] = ανi + m3(til, θi)βi

◮ ν = 0 . . . 3 ◮ i - subject index 1 . . . n ◮ l - measurement index 1 . . . q ◮ til - measurement time ◮ zil - measured twitch count ◮ ανi - intercept ◮ βi - log odds ratio

Examples

◮ TOF counts ◮ Example general vs. individualized strategies

• Stat. model: blood concentration

Additive normal error model: yij = γim1(tij, θi) + ǫij

◮ i - subject index 1 . . . n ◮ j - measurement index 1 . . . m ◮ tij - measurement time ◮ yij - measured drug concentration ◮ γi = 1/vi ◮ ǫij error

Priors π(θ) for kinetic parameters

◮ 272 surgical cases from EHR: ◮ fitted to PK/PD and POLR model described above ◮ collection of θ estimates treated as prior for θ

This type of prior:

◮ simulates a practical PK study ◮ heterogeneous patients undergoing surgery ◮ variability in dosing (amount and timing) based on clinical

need

Sample time optimization

Design sample times: ξ = {t1, · · · , tm} Information matrix M(ξ, θ) =

m

• j=1

∂m1(tj, θ) ∂θ ∂m1(tj, θ) ∂θT D-optimality (local) criterion: D(ξ, θ) = |M(ξ, θ)| ED-optimality (robust) criteria: ED(ξ) =

• D(ξ, θ)π(θ)

ED(ξ|z) =

• D(ξ, θ)π(θ|z)

◮ generic optimal design: ξgen maximizes ED(ξ) ◮ individualized optimal design: ξi maximizes ED(ξ|zi)

D-relative efficiency

D-relative efficiency: Ri = D(ξi, θi) D(ξgen, θi) 1/p “How much efficiency is gained in estimating θi by using individualized optimal design relative to generic optimal design”

Simulation study

Simulation steps:

• 1. simulate subject from prior
• 2. simulate twitch counts (POLR)
• 3. compute D-relative efficiency: generic vs. individualized
• ptimal sampling

Results (strategy 1/2)

Generic:

• 1. identify optimal sample times during first dose
• 2. sample all subjects at same times

Individualized:

• 1. monitor and record twitch counts during first dose
• 2. use twitch counts to select subject-specific sample times
• 3. sample subject at optimal times during second dose

D-relative efficiency (25th/50th/75th): 0.36/0.49/0.68

Sample times (strategy 1/2)

Results (strategy 2/2)

Generic:

• 1. identify optimal sample times during first dose
• 2. sample all subjects at same times

Individualized:

• 1. simultaneously sample and monitor during first dose
• 2. use twitch counts to continuously update/individualize times

D-relative efficiency (25th/50th/75th): 1.01/1.02/1.05

Conclusions

With regard to a fast-acting paralytic (vecuronium bromide):

◮ first dose most informative

◮ paralytics distribute and eliminate quickly ◮ uncertainty in concentration prior to second dose

◮ individualized optimal designs may be less efficient ◮ efficiency very sensitive to dosing history ◮ gains in D-efficiency for individual PK is < 10%

Future directions

◮ consider ‘local’ (vs. ‘robust’) individualization

◮ ξi maximizes D(ξ, ˆ

θi)

◮ ˆ

θi is estimate using twitch counts

◮ consider efficiency in estimating mixed-effects models

Thanks!

Thanks!