Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Sensitivity Analysis for Translation Assumptions

Preclinical to Clinical Translation · Tutorial 529 of 598

Sensitivity Analysis for Translation Assumptions

Use local and scenario-based sensitivity analysis to see how translation assumptions change predicted human exposure and response, and to prioritize assumptions for further data collection.

Advanced 12 min read

What You'll Learn

  • Define a model output as a function of translation assumptions and distinguish local from global sensitivity.
  • Quantify how an allometric exponent changes projected human clearance and exposure.
  • Examine how unbound fraction and target affinity jointly affect predicted target occupancy.
  • Trace disease-state assumptions through a predicted pharmacodynamic response.
  • Rank uncertainties without treating assumed ranges as probability distributions.
  • Diagnose limitations of one-at-a-time sensitivity analysis and select follow-up data accordingly.

Why Sensitivity Analysis Matters in Translation

A human PK/PD projection depends not only on measured animal data but also on assumptions used to transport parameters across species. An allometric exponent affects predicted clearance; an assumed unbound fraction changes the concentration available to interact with a target; target affinity influences the concentration-response relation; and disease state can alter baseline response or the size of a drug effect. If a human prediction changes substantially when one such assumption changes modestly, the assumption deserves particular attention.

Sensitivity analysis is a structured way to investigate this dependence. It does not establish which assumption is true, and it does not replace uncertainty analysis or new data. It shows how specified changes in assumptions move a specified model output, while other inputs are held fixed or varied according to an explicit plan.

Definition: Let \(Y=g(\theta_1,\ldots,\theta_k)\) be a human PK/PD prediction calculated from a model \(g\) and translation assumptions \(\theta_1,\ldots,\theta_k\). A sensitivity analysis evaluates how changes in one or more \(\theta_j\) change \(Y\). A local sensitivity describes change near a chosen reference set of assumptions; a scenario-based sensitivity evaluates predictions over specified alternative values or combinations of assumptions.

The output \(Y\) must be defined before sensitivity is assessed. It might be clearance in L/h, unbound concentration in mg/L, target occupancy as a proportion, or a biomarker response in its measurement units. Sensitivity depends on the output: an assumption that strongly affects exposure need not strongly affect a response that is already near its plateau.

Local and Scenario-Based Measures

For a differentiable model, the local absolute sensitivity of \(Y\) to an input \(\theta_j\) is its partial derivative, evaluated at a stated reference point. It gives the approximate change in \(Y\) for a small change in \(\theta_j\), with other inputs fixed. When both quantities are positive, a dimensionless local measure is the elasticity: the proportional change in output per proportional change in input.

$$ \frac{\partial Y}{\partial \theta_j} \qquad\text{and}\qquad S_j=\frac{\partial Y}{\partial \theta_j}\frac{\theta_j}{Y} =\frac{\partial \log Y}{\partial \log \theta_j} $$

An elasticity of 0.5 means that, locally, a 1% increase in the input corresponds to about a 0.5% increase in the output. This interpretation is not appropriate if \(Y\) or \(\theta_j\) is zero, and a local derivative can be misleading over a wide range or near a nonlinear transition.

A scenario-based analysis instead specifies a range or set of values for each uncertain input and recalculates the output. The resulting output range is conditional on those chosen scenarios; unless probabilities are assigned and justified, it is not a confidence interval or a probability statement. Varying one input at a time is useful for isolating its effect, but it can miss interactions. For interacting inputs, evaluate combinations as well as individual changes.

Allometric Exponent and Projected Exposure

Suppose human clearance is projected from a reference species using a power relationship. For this analysis, body weights and reference clearance are treated as fixed, while the exponent \(b\) is uncertain. This simplified calculation isolates the exponent’s effect; a full translation would also consider other PK assumptions and evidence.

$$ CL_H=CL_R\left(\frac{WT_H}{WT_R}\right)^b $$

Here \(CL_H\) and \(CL_R\) are human and reference-species clearance, respectively, and \(WT_H/WT_R\) is a dimensionless body-weight ratio. Taking the logarithm gives \(\log(CL_H)=\log(CL_R)+b\log(WT_H/WT_R)\). Therefore,

$$ \frac{\partial \log(CL_H)}{\partial b} =\log\left(\frac{WT_H}{WT_R}\right) $$

This derivative is the change in log-clearance per unit change in the exponent. The elasticity with respect to \(b\), treating \(b\) as the input, is \(b\log(WT_H/WT_R)\). This distinction matters: the derivative is not itself the percentage change for a 1% increase in \(b\).

Worked Example: Varying an Allometric Exponent

Hypothetical assumptions. A reference species weighs 20.0 kg and has clearance \(CL_R=0.400\) L/h. The projected human weight is 70.0 kg, so the weight ratio is 3.50. Examine exponents \(b=0.67,\ 0.75,\ 0.85\). Assume that, for a simplified exposure comparison, a constant input rate of 24.0 mg/h produces steady-state average concentration equal to input rate divided by clearance.

At \(b=0.67\), \(3.50^{0.67}\approx2.3149\), so \(CL_H=0.400(2.3149)=0.9259\) L/h. At \(b=0.75\), \(3.50^{0.75}\approx2.5590\), so \(CL_H=0.400(2.5590)=1.0236\) L/h. At \(b=0.85\), \(3.50^{0.85}\approx2.9004\), so \(CL_H=0.400(2.9004)=1.1602\) L/h. Across these stated assumptions, projected clearance ranges from approximately 0.9259 to 1.1602 L/h.

The corresponding simplified concentrations are \(24.0/0.9259\approx25.92\) mg/L, \(24.0/1.0236\approx23.45\) mg/L, and \(24.0/1.1602\approx20.69\) mg/L, respectively. Thus the larger exponent produces higher clearance and lower concentration in this calculation. These values illustrate the direction and size of the exponent effect under the stated input-rate assumption; they are not a complete human PK prediction.

At the reference value \(b=0.75\), the elasticity is \(S_b=0.75\log(3.50)\approx0.940\). Locally, a 1% increase in \(b\) therefore corresponds to about a 0.940% increase in clearance. This local approximation should not be substituted for the scenario calculations across the full exponent range.

Unbound Fraction and Target Affinity

When target interaction is governed by unbound concentration, total concentration alone does not determine predicted occupancy. A simple equilibrium model defines unbound concentration as \(C_u=f_uC_{\mathrm{tot}}\), where \(f_u\) is the unbound fraction, and occupancy as \(C_u/(K_d+C_u)\), where \(K_d\) is the dissociation constant. This illustration assumes that the unbound fraction is applicable at the concentration of interest, that equilibrium is relevant, and that free concentration and \(K_d\) are expressed in the same units.

$$ O=\frac{f_u C_{\mathrm{tot}}}{K_d+f_u C_{\mathrm{tot}}} $$

Worked Example: Varying Unbound Fraction and Affinity

Hypothetical assumptions. Let total concentration be \(C_{\mathrm{tot}}=100\) ng/mL. Consider \(f_u\) values 0.02, 0.05, and 0.10, and examine target affinities represented by \(K_d=2,\ 4,\) and \(8\) ng/mL. The calculation assumes the same total concentration and the same occupancy model in every scenario.

At the central assumptions \(f_u=0.05\) and \(K_d=4\) ng/mL, \(C_u=0.05(100)=5\) ng/mL and \(O=5/(4+5)=0.5556\), or 55.56%. Holding \(K_d=4\) ng/mL fixed, the unbound-fraction scenarios give:

  • For \(f_u=0.02\), \(C_u=2\) ng/mL and \(O=2/(4+2)=0.3333\), or 33.33%.
  • For \(f_u=0.05\), \(C_u=5\) ng/mL and \(O=5/(4+5)=0.5556\), or 55.56%.
  • For \(f_u=0.10\), \(C_u=10\) ng/mL and \(O=10/(4+10)=0.7143\), or 71.43%.

Now hold \(f_u=0.05\), and thus \(C_u=5\) ng/mL, fixed while varying affinity. At \(K_d=2\) ng/mL, occupancy is \(5/(2+5)=0.7143\); at \(K_d=4\) ng/mL, it is \(5/9=0.5556\); and at \(K_d=8\) ng/mL, it is \(5/(8+5)=0.3846\). The combined corner scenarios span 20.00% occupancy at \(f_u=0.02,\ K_d=8\) ng/mL, to 83.33% at \(f_u=0.10,\ K_d=2\) ng/mL. This wide range also shows why varying each input separately around a central value may understate joint variation.

The model has a direct interaction: occupancy depends on the ratio \(f_uC_{\mathrm{tot}}/K_d\). A higher unbound fraction and a lower \(K_d\) both increase occupancy, but their effects are nonlinear and diminish as occupancy approaches one. The scenario range is conditional on the selected values and should not be read as a probability interval.

Disease-State Assumptions and Response

A translated PD prediction can depend on both baseline disease state and the magnitude of the drug effect. The following simple model illustrates how these assumptions can be separated for sensitivity analysis. It treats a biomarker as decreasing from baseline under treatment and assumes a fixed concentration and \(EC_{50}\); it does not assert that every disease or mechanism follows this form.

$$ Y=R_0-E_{\max}\frac{C}{EC_{50}+C} $$

Here \(R_0\) is the untreated disease-state biomarker value, \(E_{\max}\) is the maximum decrease in the biomarker, and \(C\) and \(EC_{50}\) use the same concentration units. With \(C\) and \(EC_{50}\) fixed, the derivatives are \(\partial Y/\partial R_0=1\) and \(\partial Y/\partial E_{\max}=-C/(EC_{50}+C)\). These provide a direct local account of how baseline and drug-effect assumptions enter the predicted endpoint.

Worked Example: Separating Baseline and Drug-Effect Assumptions

Hypothetical assumptions. Let concentration \(C=6.0\) mg/L and \(EC_{50}=4.0\) mg/L. Consider disease-state baselines \(R_0=140\) or 160 biomarker units and transported drug-effect assumptions \(E_{\max}=30\) or 45 biomarker units. Assume the same concentration-response form applies in each scenario.

The concentration fraction is \(6.0/(4.0+6.0)=0.600\). For \(R_0=140\) and \(E_{\max}=30\), the prediction is \(Y=140-30(0.600)=122\) units. For \(R_0=140\) and \(E_{\max}=45\), it is \(Y=140-45(0.600)=113\) units. At \(R_0=160\), the predictions are \(160-30(0.600)=142\) units and \(160-45(0.600)=133\) units, respectively. The four scenarios therefore span 113 to 142 units.

Holding \(E_{\max}\) fixed, changing baseline from 140 to 160 units raises the predicted endpoint by 20 units. Holding \(R_0\) fixed, changing \(E_{\max}\) from 30 to 45 units lowers the endpoint by \(15(0.600)=9\) units. Under this model the two effects add, but in more complex models disease state may also alter potency, turnover, or other parameters. Such dependencies require joint scenarios rather than an assumption that each input acts independently.

How to Prioritize Assumptions and Detect Problems

A useful analysis starts from a defined reference prediction, varies defensible assumptions, and reports the resulting change in an interpretable output with units. A large output range identifies an assumption with substantial influence over the chosen range; it does not establish that the assumption is likely to be wrong. Prioritization should also consider whether the uncertainty is reducible, whether relevant data can be collected, and whether resolving it could change a scientific or development decision.

1
Define the output and reference assumptions.
State whether the target is clearance, exposure, occupancy, or a response endpoint. Record the values held fixed and their units.
2
Specify plausible alternatives.
Use evidence-based ranges or discrete scenarios for exponents, binding, affinity, and disease-state parameters. Do not label a range probabilistic unless a probability model has been supplied.
3
Recalculate and compare.
First vary inputs individually to see their isolated effects; then evaluate plausible combinations where interactions are possible. Report both input assumptions and output changes.
4
Link influence to a data opportunity.
Ask which influential assumption can be better measured and whether doing so could materially change the prediction or its use.

Several diagnostics help identify an unreliable sensitivity analysis. Recompute a few scenarios independently or with a second implementation, check that units are consistent, and verify expected directions where the model implies them. For example, in the allometric equation with a weight ratio greater than one, increasing \(b\) must increase projected clearance. If a response changes abruptly or reverses direction, inspect for coding errors, parameter boundaries, and nonlinear transitions before interpreting the result. If the output is non-monotone, report enough scenarios to show that pattern rather than only its minimum and maximum.

One-at-a-time analysis is transparent but can miss interactions such as the joint effect of \(f_u\) and \(K_d\). A local derivative is efficient near a reference point but can misrepresent a broad range, especially for saturating response models. Scenario analysis can expose nonlinearity and interactions but depends on selected values. When defensible probability distributions are available and the goal is probabilistic uncertainty quantification, simulation-based methods can propagate joint uncertainty; variance-based methods can further allocate output variance under their assumptions. These approaches answer different questions from a sensitivity analysis using unweighted scenarios.

Common Pitfall and Limitations

Common mistake: “The exponent generated the largest output range, so it is the most likely assumption to be wrong.”

Correction: A sensitivity range measures impact across the values chosen, not the probability that those values occur. An assumption can be influential but already well established, or uncertain but have little effect on the selected output. Prioritize by considering both influence and how much relevant evidence could reduce the uncertainty.

Sensitivity analysis is conditional on the model structure. It cannot expose a missing mechanism merely by varying parameters already included in the model, and it cannot validate a cross-species relationship. The selected output and range also matter: sensitivity of clearance may not rank assumptions in the same order as sensitivity of response. For these reasons, preserve the model context, show assumptions transparently, and avoid presenting one sensitivity ranking as universal.

Check Your Understanding

Answer each question, then compare your reasoning with the brief solutions.

  1. If \(CL_H=CL_Rr^b\) with \(r=2\), does increasing \(b\) raise or lower \(CL_H\)? Why?
  2. At \(C_u=6\) ng/mL and \(K_d=3\) ng/mL, calculate occupancy. What is the occupancy if \(K_d\) doubles while \(C_u\) stays fixed?
  3. In the disease-state model, how much does \(Y\) change when \(R_0\) rises by 8 units while all other inputs are fixed?
  4. Why is a scenario output range not automatically a 95% interval?
  5. Name one reason to supplement one-at-a-time sensitivity analysis with joint scenarios.

Solutions. (1) It rises because \(r>1\), so \(r^b\) increases with \(b\). (2) Occupancy is \(6/(3+6)=0.6667\), or 66.67%; if \(K_d=6\) ng/mL, it is \(6/(6+6)=0.5000\), or 50.00%. (3) The derivative with respect to \(R_0\) is 1, so \(Y\) increases by 8 units. (4) The range uses chosen scenarios and has no probability coverage unless a probability model and interval procedure are defined. (5) Joint scenarios can reveal interactions, such as the combined influence of unbound fraction and affinity on occupancy.

Clinical Biostats