Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Deriving Allometric Clearance and Volume Predictions

Preclinical to Clinical Translation · Tutorial 506 of 598

Deriving Allometric Clearance and Volume Predictions

Learn to fit clearance and volume allometric relationships from animal data and calculate human predictions with units and assumptions made explicit.

Advanced 11 min read

What You'll Learn

  • Express clearance and volume as functions of body weight relative to a stated reference mass.
  • Derive the log-linear form used to estimate an allometric exponent.
  • Estimate an exponent from two species and from a multispecies regression.
  • Calculate human clearance and volume predictions from fitted equations.
  • Distinguish a fitted scaling relationship from a mechanistic explanation.
  • Recognize when data, assumptions, or extrapolation limit an allometric prediction.

Why Scale Clearance and Volume?

Preclinical PK measurements can help project human pharmacokinetic parameters, but a value measured in one animal species is not itself a human value. Body weight provides a simple way to describe how a parameter changes across species. In an allometric model, clearance and volume are represented as power functions of body weight. The exponent describes the estimated scaling relationship in the data; it is not automatically a universal biological constant.

The predictions in this tutorial are conditional on comparable parameter definitions and measurement methods across species. For example, the clearance values should refer to the same kind of clearance, and volume values should refer to the same volume parameter. A mass-based relationship is useful for translation, but it does not by itself account for every species difference or establish a mechanism.

Definition: A reference-mass allometric model for a PK parameter \(P\) is $$ P(W)=P_{\mathrm{ref}}\left(\frac{W}{W_{\mathrm{ref}}}\right)^b, $$ where \(W\) is body weight, \(W_{\mathrm{ref}}\) is a chosen reference weight, \(P_{\mathrm{ref}}\) is the parameter value at that weight, and \(b\) is the allometric exponent. The weight ratio is dimensionless, so \(P(W)\) and \(P_{\mathrm{ref}}\) have the same units. For clearance, \(P\) is \(CL\); for volume, it is \(V\).

From a Power Function to a Fitted Line

Taking logarithms of the allometric model turns its power relationship into a straight-line relationship. For any one parameter \(P\),

$$ \ln P(W)=\ln P_{\mathrm{ref}}+b\ln\left(\frac{W}{W_{\mathrm{ref}}}\right). $$

Thus, if \(x_i=\ln(W_i/W_{\mathrm{ref}})\) and \(y_i=\ln(P_i)\), the model can be written \(y_i=a+bx_i+\epsilon_i\), where \(a=\ln(P_{\mathrm{ref}})\) and \(\epsilon_i\) is the residual on the log scale. A regression of logged parameter values on logged relative weights estimates the exponent as the slope and the reference parameter as the exponential of the intercept. The same procedure is applied separately to clearance and volume; there is no requirement that their exponents be equal.

This regression interpretation relies on assumptions. The species-level values must be comparable and positive so their logarithms exist. For ordinary least squares on the log scale, the residuals are assumed to have mean zero at each \(x_i\), to have a common variance, and to be independent across observations. The assumed mean relationship is linear on the log scale. If values are species means, treating them as independent observations is an additional working assumption; biological relatedness or shared study methods can make their errors dependent. With only two species, a line can always pass through both points, so those data cannot test whether the power relationship is adequate.

The unreferenced form \(P=aW^b\) is equivalent, but its coefficient \(a\) depends on the units used for weight. The reference form makes the reference point explicit and keeps the ratio dimensionless. Changing the reference weight changes \(P_{\mathrm{ref}}\), but not the exponent if the same fitted line is expressed at a different reference weight.

Estimating an Exponent

With exactly two species, the exponent is the slope between their log-transformed values. With three or more species, a log-linear regression uses all the observations to estimate the slope and intercept. For an ordinary least-squares fit, the slope is

$$ \hat b= \frac{\sum_i(x_i-\bar{x})(y_i-\bar{y})} {\sum_i(x_i-\bar{x})^2}, \qquad \hat a=\bar{y}-\hat b\bar{x}. $$

Here \(\bar{x}\) and \(\bar{y}\) are the means of the logged relative weights and logged parameters, not the means of the original weights and parameters. The fitted reference parameter is \(\widehat{P}_{\mathrm{ref}}=\exp(\hat a)\) when natural logarithms are used. If base-10 logarithms are used instead, it is \(10^{\hat a}\). The estimated exponent is the same under either logarithm base, provided the same base is used for both axes.

Worked Example: Estimating Clearance Scaling From Two Species

Suppose a hypothetical study reports clearance of \(0.100\ \mathrm{L/h}\) at \(0.250\ \mathrm{kg}\) and \(1.00\ \mathrm{L/h}\) at \(10.0\ \mathrm{kg}\). Estimate the exponent and predict clearance at \(70.0\ \mathrm{kg}\). Use the \(0.250\ \mathrm{kg}\) animal as the reference mass.

Under the two-point power model, the ratio of clearances equals the ratio of weights raised to the exponent. Taking logarithms gives

$$ \hat b =\frac{\ln(1.00/0.100)}{\ln(10.0/0.250)} =\frac{\ln(10)}{\ln(40)} =\frac{2.3026}{3.6889} \approx 0.6242. $$

The fitted reference clearance is \(0.100\ \mathrm{L/h}\) at \(0.250\ \mathrm{kg}\). Therefore,

$$ \widehat{CL}_{70} =0.100\left(\frac{70.0}{0.250}\right)^{0.6242} =0.100(280)^{0.6242} \approx 3.37\ \mathrm{L/h}. $$

The prediction is \(3.37\ \mathrm{L/h}\) under this fitted two-point relationship. It is an extrapolation beyond the heavier observed weight of \(10.0\ \mathrm{kg}\), so the calculation does not show that the same relationship continues to human body weight. With only two species, there are no additional observations to assess the fit.

Worked Example: Fitting Clearance With Multispecies Data

Consider these hypothetical species-level clearance values. Use \(W_{\mathrm{ref}}=1.00\ \mathrm{kg}\), fit an unweighted least-squares line using base-10 logarithms, and predict clearance at \(70.0\ \mathrm{kg}\).

Body weight \(W\) (kg)Clearance \(CL\) (L/h)\(x=\log_{10}(W/1.00)\)\(y=\log_{10}(CL)\)
0.1000.0200-1-1.6990
1.000.1100-0.9586
10.00.6031-0.2197
1003.3120.5198

The mean logged weight is \(\bar{x}=0.5000\), and the mean logged clearance is approximately \(\bar{y}=-0.5894\). Substitution in the least-squares slope formula gives

$$ \hat b =\frac{\sum_i(x_i-0.5000)(y_i+0.5894)} {\sum_i(x_i-0.5000)^2} \approx\frac{3.6977}{5.000} =0.7395. $$

The fitted intercept is

$$ \hat a=\bar{y}-\hat b\bar{x} =-0.5894-(0.7395)(0.5000) \approx -0.9591. $$

The corresponding fitted reference clearance is \(10^{-0.9591}\approx0.1099\ \mathrm{L/h}\) at \(1.00\ \mathrm{kg}\). Thus the fitted model is

$$ \widehat{CL}(W)=0.1099\left(\frac{W}{1.00\ \mathrm{kg}}\right)^{0.7395}\ \mathrm{L/h}. $$

At \(70.0\ \mathrm{kg}\), the weight ratio is \(70.0/1.00=70.0\), so

$$ \widehat{CL}_{70} =0.1099(70.0)^{0.7395} \approx 2.54\ \mathrm{L/h}. $$

This is an estimate from the stated regression and units, not evidence that the fitted exponent applies to every species or compound. In this example, \(70.0\ \mathrm{kg}\) lies within the observed weight range, although it is still a prediction at a weight where no species was observed.

Predicting Volume

The same derivation and fitting procedure applies to volume. The volume parameter must be defined consistently across species: mixing, for example, different volume concepts would make the fitted relationship difficult to interpret. The following example uses a deliberately simple dataset to show how the regression and human projection work when the values follow a power law exactly.

Worked Example: Estimating Volume Scaling and Predicting a Human Value

Suppose hypothetical volume measurements are \(0.0600\), \(0.600\), \(6.00\), and \(60.0\ \mathrm{L}\) at body weights \(0.100\), \(1.00\), \(10.0\), and \(100\ \mathrm{kg}\), respectively. Fit the reference-mass model at \(1.00\ \mathrm{kg}\) and predict volume at \(70.0\ \mathrm{kg}\).

The base-10 logged weights are \(-1,0,1,2\), while the logged volumes are approximately \(-1.2218,-0.2218,0.7782,1.7782\). Each increase of one unit in logged weight corresponds to an increase of one unit in logged volume, so the least-squares slope is \(\hat b=1.000\). The intercept is \(-0.2218\), giving a reference volume of \(10^{-0.2218}\approx0.6000\ \mathrm{L}\) at \(1.00\ \mathrm{kg}\). The fitted relationship is therefore

$$ \widehat{V}(W)=0.6000\left(\frac{W}{1.00\ \mathrm{kg}}\right)^{1.000}\ \mathrm{L}. $$

For a \(70.0\ \mathrm{kg}\) human, the prediction is

$$ \widehat{V}_{70} =0.6000\left(\frac{70.0}{1.00}\right)^{1.000} =42.0\ \mathrm{L}. $$

The exact fit here is a feature of the constructed numbers, not a general property of volume data. In real multispecies data, observations need not lie exactly on a straight line after logging.

Assumptions, Checks, and Alternatives

A fitted line should be treated as a useful approximation only when the data support its use. Before calculating a prediction, check that weights and parameters are positive, units and parameter definitions are consistent, and the fitted relationship is being used over a defensible range. A plot of logged parameter against logged weight can reveal curvature or an observation far from the general pattern. Residuals on the log scale can reveal systematic departures from the assumed line or a spread that changes with weight. Such patterns warn that the stated log-linear model may not adequately describe the data; the equation alone does not resolve the cause.

Common alternatives answer somewhat different questions. A conventional fixed exponent can be used when a project has a justified prespecified scaling rule; it avoids estimating a slope from a small dataset, but it may not represent the compound and species data at hand. A two-species calculation is transparent and requires little data, but provides no independent information about curvature, scatter, or influential observations. A multispecies regression estimates the exponent from more observations, but depends on the regression assumptions and on whether those observations provide independent, comparable information. More mechanistic approaches can represent physiological processes directly, but require additional data and assumptions; an empirical allometric fit is not a substitute for those mechanisms.

If the model is fitted to logarithms, exponentiating the fitted line gives the fitted geometric-scale relationship. When log-scale residuals are random with mean zero, this corresponds to the conditional median on the original scale under a lognormal error model. It is not automatically the conditional arithmetic mean. If the intended prediction is an arithmetic mean, retransformation may require a correction based on the residual distribution. Do not silently label an exponentiated log-scale fit as an arithmetic mean prediction.

Allometric projections also have limitations beyond the regression calculation. A body-weight relationship does not account explicitly for species differences in enzymes, transporters, binding, or physiology. The fitted exponent and reference parameter summarize the data used; neither proves a causal mechanism. Predictions outside the observed weight range are particularly dependent on the assumption that the same power relationship continues there. These considerations are reasons to state the data range, fitted parameter definition, units, and prediction target alongside a projected human value.

Common pitfall: Do not calculate a human parameter by multiplying an animal value by a weight ratio without raising that ratio to the fitted exponent. In the reference model, the correct prediction is \(P_{\mathrm{ref}}(W/W_{\mathrm{ref}})^{\hat b}\). Also retain enough precision in \(\hat b\) and the intercept while calculating; round the final prediction to a sensible number of significant figures.
Key takeaway: Allometric translation estimates a power relationship between body weight and a consistently defined PK parameter. Fit the slope and intercept on the log scale, convert the intercept back to the parameter scale at a stated reference weight, and apply the resulting equation with units and extrapolation limits made explicit.

Check Your Understanding

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

  1. A parameter is \(0.200\ \mathrm{L/h}\) at \(0.500\ \mathrm{kg}\) and \(0.800\ \mathrm{L/h}\) at \(8.00\ \mathrm{kg}\). Write the two-point expression for the exponent.
  2. In a log-linear fit using \(x=\log_{10}(W/W_{\mathrm{ref}})\) and \(y=\log_{10}(CL)\), what do the slope and intercept represent?
  3. A model predicts volume at a weight well above every observed species weight. What is the main limitation of this prediction?
  4. Why can two species determine an exponent but not provide a check of whether a power relationship is adequate?

Solutions. (1) \(\hat b=\ln(0.800/0.200)/\ln(8.00/0.500)=\ln(4)/\ln(16)=0.500\). (2) The slope estimates the allometric exponent; the intercept is the logged clearance at \(W_{\mathrm{ref}}\), so exponentiating it gives the reference clearance. (3) It is an extrapolation that assumes the fitted power relationship continues beyond the observed range. (4) Two points determine a unique straight line on the log scale, leaving no additional observations to assess departures from that line.

Clinical Biostats