PK/PD Translation
The arithmetic between a number on a plate and a number on a label. It decides whether a compound
can reach its target concentration at a dose people can tolerate — and it is where a surprising
number of programmes discover, late, that it cannot.
No installation, no network, no key. Every calculation here is closed-form or a fixed-step
numerical integration in the standard library.
Reference: FDA, Estimating the Maximum Safe Starting Dose in Initial Clinical Trials for
Therapeutics in Adult Healthy Volunteers (2005) — the Km factors are its Table 1.
Read references/nca-and-parameters.md before quoting a PK
parameter, references/interspecies-scaling.md before
converting a dose between species, and
references/translating-potency-to-dose.md before
turning an IC50 into a dose — that one is judgement, not syntax.
The four scripts
| Script | Answers |
|---|
nca.py | What are this profile's clearance, volume, and half-life? |
pk_compartmental.py | What does this regimen look like at steady state? |
allometry.py | What is the human equivalent of this animal dose? |
exposure_margin.py | Does the projected exposure cover the target, and is it safe? |
Free drug, or the answer is wrong by 1/fu
This is the one to get right. Only unbound drug engages the target. Comparing a total plasma
concentration against a free-drug IC50 is the most common translation error, and at 99% protein
binding it is a hundredfold error in the dangerous direction:
C_free = C_total x fu 1000 nM total at fu = 0.01 -> 10 nM free
Against a 100 nM IC50 that is tenfold under-coverage, not the tenfold coverage the total number
suggests. exposure_margin.py coverage takes --fu and reports both.
A related trap: 50% inhibition is rarely enough. Occupancy is C/(C+Ki), so 90% needs nine times
Ki and 95% needs nineteen. Ask what fraction the biology requires, and for how long.
Non-compartmental analysis
python skills/pkpd-translation/scripts/nca.py --times 0.25,0.5,1,2,4,6,8,12,18,24,36,48 \
--conc 9.75,9.51,9.05,8.19,6.70,5.49,4.49,3.01,1.65,0.91,0.27,0.08 --dose 100 --route iv
Validated against an analytic one-compartment case (D = 100, V = 10, k = 0.1), this returns
CL = 1.000000, Vz = 10.000000, t½ = 6.931472, AUC(0–∞) = 100.000000 — exact.
Three defaults that matter. Linear-up / log-down trapezoid, because a straight line between
two descending points sits above an exponential and inflates AUC. Back-extrapolation to C0 for
IV bolus, because sampling never starts at t = 0 and the missing first segment was 2.5% of AUC in
the case above. And CL/F and Vz/F labelling after oral dosing, because calling them CL and V
assumes complete bioavailability.
nca.py warns when the terminal fit is poor or when more than 20% of AUC(0–∞) is extrapolated
past the last sample — above that, CL and Vz are guesses.
Steady state
python skills/pkpd-translation/scripts/pk_compartmental.py steady --dose 100 --cl 5 --v 50 --tau 24
half_life_h accumulation_ratio time_to_95pct_ss_h css_max css_min css_avg auc_tau
6.93147 1.09977 29.9573 2.19954 0.199538 0.833333 20
Accumulation depends on half-life and interval, never on dose — R = 1/(1 − e^(−kτ)). Doubling
the dose moves every concentration and leaves the ratio alone. Time to steady state depends only
on half-life, about 4.3 half-lives to 95%, which is why loading doses exist.
Note Css,min of 0.2 against Css,avg of 0.83: the average hides a fourfold swing. Report all three.
Species conversion
python skills/pkpd-translation/scripts/allometry.py hed --species rat --species dog --dose 50
species km divide_by hed_mg_kg hed_total_mg_60kg
rat 6 6.167 8.10811 486.486
dog 20 1.85 27.027 1621.62
Never transfer mg/kg directly. Metabolic rate scales with surface area, so a rat dose divides
by about 6.2 and a dog dose by 1.8. Skipping that step overdoses the first human cohort by exactly
that factor.
fih takes NOAELs across species, converts each, and applies a safety factor to the lowest —
the most sensitive species sets the dose, not the average. The default factor of 10 should rise
for steep dose–response, irreversible toxicity, or a novel target. For biologics and agonists,
MABEL usually governs instead of the NOAEL route; TGN1412 started 500-fold below the NOAEL HED and
still nearly killed six volunteers.
Coverage and margin
python skills/pkpd-translation/scripts/exposure_margin.py coverage \
--dose 100 --cl 5 --v 50 --tau 24 --target-conc 0.05 --fu 0.02
python skills/pkpd-translation/scripts/exposure_margin.py margin \
--tox-cmax 12000 --eff-cmax 400
Both sides of a margin must be the same metric and the same binding basis. A total-drug
toxicology Cmax against a free-drug efficacy target inflates the margin by 1/fu.
Four ways this arithmetic misleads
- The efficacy driver and the tolerability driver are often different metrics. Time above
threshold versus Cmax argues for smaller, more frequent doses; the reverse argues the opposite.
- Single-dose PK does not give the steady-state trough, and the trough is what determines
whether coverage holds through the interval.
- Everything assumes linear PK. Saturable metabolism, absorption, or binding breaks
superposition, and then none of the accumulation arithmetic applies.
- This is a typical subject, not a population. Real between-subject variability is severalfold
and is not represented anywhere here.
When to stop using these scripts
When the site of action is not plasma (brain, tumour, intracellular), when transporters or
first-pass metabolism dominate, when you need a population rather than a typical subject, or when
the PD is delayed or tolerance-developing. Use PK-Sim/MoBi (Open Systems Pharmacology, open
source) for PBPK, and NONMEM, Monolix, or nlmixr2 for population PK/PD.
Composing with the rest of the bundle
admet-prediction → here: predicted clearance, half-life, and plasma protein binding as inputs.chembl → here: a measured IC50 or Ki to translate into a target concentration.openfda → alongside: the approved label's dose and interval as a real-world anchor.clinicaltrials → alongside: what doses have actually been taken into humans.free-energy-perturbation → before: a potency number worth translating.
Reporting results honestly
Give the free concentration and the fu behind it. Name the exposure metric you believe drives
efficacy and say why. Give Cmax, Cmin, and Cavg, not one. Give margins on a stated matched basis.
Say that this is a one- or two-compartment plasma approximation for a typical subject with no
population variability — it is planning arithmetic, not a regulatory submission.
