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.
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\),
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
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
The fitted reference clearance is \(0.100\ \mathrm{L/h}\) at \(0.250\ \mathrm{kg}\). Therefore,
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.100 | 0.0200 | -1 | -1.6990 |
| 1.00 | 0.110 | 0 | -0.9586 |
| 10.0 | 0.603 | 1 | -0.2197 |
| 100 | 3.31 | 2 | 0.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
The fitted intercept is
The corresponding fitted reference clearance is \(10^{-0.9591}\approx0.1099\ \mathrm{L/h}\) at \(1.00\ \mathrm{kg}\). Thus the fitted model is
At \(70.0\ \mathrm{kg}\), the weight ratio is \(70.0/1.00=70.0\), so
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
For a \(70.0\ \mathrm{kg}\) human, the prediction is
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.
Check Your Understanding
Answer each question, then compare your reasoning with the short solutions.
- 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.
- 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?
- A model predicts volume at a weight well above every observed species weight. What is the main limitation of this prediction?
- 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.