Tutorials › Nonlinear Mixed-Effects Models (NLME) › Combining Allometry and Maturation for Clearance

Allometric Scaling and Maturation Functions · Tutorial 387 of 797

Combining Allometry and Maturation for Clearance

Learn to combine weight scaling and age-dependent maturation in one clearance equation, calculate pediatric predictions, and check the model’s reference values, units, and assumptions.

Intermediate 11 min read

What You'll Learn

  • Write clearance as a weight-scaling factor multiplied by an age-maturation factor.
  • Define the reference clearance as the fully mature value at a specified reference weight.
  • Calculate predicted clearance for a child using weight, age, and maturation parameters.
  • Compare predictions for children who differ in both weight and age.
  • Check relative-clearance calculations without relying on rounded intermediate predictions.
  • Recognize common unit, reference-value, and model-interpretation errors.

Why Clearance Needs Both Size and Maturation

A child’s clearance can differ from an adult’s for more than one reason. Body size affects clearance, and the physiological processes responsible for eliminating a drug may also be developmentally immature. A model with only a weight term can therefore miss age-related development; a model with only an age term can miss the effect of size. This tutorial combines the two factors in one equation.

The central modeling choice here is multiplicative: weight scaling adjusts a reference clearance for body size, and a maturation function adjusts it for developmental age. The resulting prediction is conditional on the specified reference values and maturation curve. It is not a claim that weight and age are interchangeable, or that this equation is appropriate for every drug.

Definition: Let \(CL_i\) be the predicted clearance for individual \(i\), \(WT_i\) the individual’s body weight, and \(A_i\) the age used to describe maturation. A combined allometric–maturation model is $$ CL_i = CL_{ref}\left(\frac{WT_i}{WT_{ref}}\right)^{\theta} F_{mat}(A_i), $$ where \(CL_{ref}\) is the fully mature clearance at reference weight \(WT_{ref}\), \(\theta\) is the allometric exponent, and \(F_{mat}(A_i)\) is a dimensionless maturation factor. In the examples below, \(F_{mat}(A)\) is the sigmoidal Hill function $$ F_{mat}(A)=\frac{A^h}{TM_{50}^{\,h}+A^h}, $$ with \(TM_{50}\) the age at half maturation and \(h>0\) the Hill coefficient.

The equation is dimensionally consistent when \(WT_i\) and \(WT_{ref}\) use the same weight units, and \(A_i\) and \(TM_{50}\) use the same age units. The weight ratio and maturation factor are dimensionless, so \(CL_i\) has the same units as \(CL_{ref}\). For illustration, the calculations use \(\theta=0.75\), a common allometric exponent in clearance models, but the exponent should be justified for the application rather than assumed to be universally correct.

The reference definition matters. Since \(F_{mat}(A)\) approaches 1 as age becomes large, \(CL_{ref}\) is the clearance at \(WT_{ref}\) under fully mature function. It is not necessarily the clearance at the reference weight for a child at a particular age. For example, at \(A=TM_{50}\), the maturation factor is 0.5, so predicted clearance at \(WT_{ref}\) is half of \(CL_{ref}\).

How the Two Factors Combine

The size factor is \((WT_i/WT_{ref})^\theta\). It scales clearance up or down relative to the reference weight. The maturation factor ranges between 0 and 1 for positive ages under the Hill function. It scales the fully mature value downward when function is not yet mature. Because the factors multiply, the combined prediction reflects both adjustments.

The Hill function has useful reference properties. At \(A=TM_{50}\), its numerator and denominator are equal, giving \(F_{mat}=0.5\). At ages much smaller than \(TM_{50}\), the factor is small; at ages much larger than \(TM_{50}\), it approaches 1. The age and \(TM_{50}\) must be expressed in the same units. The exponent \(h\) controls the curve’s steepness, but the form and parameter values should be supported by the modeling context.

Worked Example: Predicting Clearance at One Weight and Age

Suppose a hypothetical drug has \(CL_{ref}=5.00\) L/h at \(WT_{ref}=70.0\) kg. Use \(\theta=0.75\), \(TM_{50}=1.00\) year, and \(h=2\). Predict clearance for a child weighing 10.0 kg at age 0.500 years. Assume the age supplied is the appropriate developmental age for this hypothetical calculation.

First calculate the weight factor: $$ \left(\frac{10.0\text{ kg}}{70.0\text{ kg}}\right)^{0.75} =\left(\frac{1}{7}\right)^{0.75} \approx 0.2324. $$ The units cancel in the ratio.

Next calculate maturation: $$ F_{mat}(0.500) =\frac{(0.500\text{ y})^2}{(1.00\text{ y})^2+(0.500\text{ y})^2} =\frac{0.250}{1.250} =0.2000. $$

Multiply the reference clearance by both factors: $$ CL_i=5.00\text{ L/h}\times 0.2324\times 0.2000 \approx 0.2324\text{ L/h}. $$ The prediction is much lower than 5.00 L/h because both the weight factor and the maturation factor are below 1. This is a model-based prediction, not a measured value.

Comparing Children Who Differ in Weight and Age

A useful way to compare two predictions is to divide one model equation by the other. The common reference clearance cancels, leaving the relative contributions of weight and maturation. This calculation helps separate the arithmetic of the two factors without treating their effects as additive.

Derivation: For individuals 1 and 2 under the same model parameters, $$ \frac{CL_2}{CL_1} = \frac{CL_{ref}(WT_2/WT_{ref})^\theta F_{mat}(A_2)} {CL_{ref}(WT_1/WT_{ref})^\theta F_{mat}(A_1)} = \left(\frac{WT_2}{WT_1}\right)^\theta \frac{F_{mat}(A_2)}{F_{mat}(A_1)}. $$ The reference clearance and reference weight cancel. This result assumes the same \(\theta\) and maturation function apply to both individuals.

Worked Example: Combining Different Weights and Ages

Use the same hypothetical model: \(CL_{ref}=5.00\) L/h at 70.0 kg, \(\theta=0.75\), \(TM_{50}=1.00\) year, and \(h=2\). Compare a child weighing 10.0 kg at 0.500 years with a child weighing 20.0 kg at 2.00 years.

For the 10.0-kg child, the weight factor is approximately 0.232368, and the maturation factor is \(0.250/(1.000+0.250)=0.2000\). Thus $$ CL_1=5.00\times 0.232368\times 0.2000 \approx 0.232368\text{ L/h}. $$

For the 20.0-kg child, the weight factor is $$ \left(\frac{20.0}{70.0}\right)^{0.75} =\left(\frac{2}{7}\right)^{0.75} \approx 0.390795, $$ and the maturation factor is $$ F_{mat}(2.00)=\frac{(2.00)^2}{(1.00)^2+(2.00)^2} =\frac{4.00}{5.00}=0.8000. $$ Therefore, $$ CL_2=5.00\times 0.390795\times 0.8000 \approx 1.56318\text{ L/h}. $$

Using unrounded factors, the relative prediction is $$ \frac{CL_2}{CL_1} =\left(\frac{20.0}{10.0}\right)^{0.75}\frac{0.8000}{0.2000} =1.68179\times 4.000 \approx 6.727. $$ The second child’s predicted clearance is about 6.727 times the first child’s under this model. If the displayed predictions are divided after rounding, the result can differ slightly; comparisons should use unrounded values when available.

Checking the Model Before Interpreting Predictions

The calculation is only meaningful if its inputs and structural assumptions fit the intended use. Before evaluating predictions, state the assumed clearance reference, weight reference, exponent, age definition, and maturation parameters. In particular, the age measure must be appropriate to the application and must use the same units as \(TM_{50}\). The equation as written also assumes that the same maturation function multiplies the allometric term across the individuals being compared.

The multiplication has a specific interpretation: each factor modifies clearance relative to the fully mature reference. It does not say that the effects are independent in every physiological sense, nor does a prediction by itself establish a causal effect of weight or age. It is a chosen model structure. If the observed process is not represented adequately by a single maturation function, or if the reference clearance and parameters are poorly supported, numerical precision in the calculation does not make the prediction reliable.

Worked Example: Isolating the Age Adjustment at a Fixed Weight

For a second hypothetical illustration, let \(CL_{ref}=6.00\) L/h at 70.0 kg, with \(\theta=0.75\), \(TM_{50}=2.00\) years, and \(h=2\). Compare two children who each weigh 35.0 kg, one aged 1.00 year and the other aged 2.00 years.

The common weight factor is $$ \left(\frac{35.0}{70.0}\right)^{0.75} =(0.500)^{0.75} \approx 0.594604. $$ At age 1.00 year, maturation is $$ F_{mat}(1.00)=\frac{1.00^2}{2.00^2+1.00^2} =\frac{1}{5}=0.2000. $$ The prediction is \(6.00\times0.594604\times0.2000\approx0.7135\) L/h.

At age 2.00 years, maturation is \(4.00/(4.00+4.00)=0.5000\). The prediction is \(6.00\times0.594604\times0.5000\approx1.784\) L/h. The ratio is \(0.5000/0.2000=2.500\), consistent with \(1.784/0.7135\approx2.500\) using unrounded predictions. Because weight is held fixed, the difference comes entirely from the maturation factor in this model.

Diagnostics, Alternatives, and Limitations

A basic computational check is to verify the two factors separately before multiplying them. The weight factor should equal 1 at \(WT_i=WT_{ref}\), and the maturation factor should equal 0.5 at \(A_i=TM_{50}\). At positive ages, the Hill maturation factor should lie between 0 and 1. A result outside that range usually signals a coding, unit, or parameter error for this particular function. Also check that the final clearance has the same units as \(CL_{ref}\), and that calculations at the reference weight and a very mature age approach \(CL_{ref}\).

A common mistake is to use a reference clearance that already represents a child’s partially mature clearance, then multiply by a maturation factor a second time. That applies the maturation adjustment twice. The correction is to define \(CL_{ref}\) explicitly as the fully mature clearance at \(WT_{ref}\), or to reformulate the reference parameter consistently if a different reference is intended.

Another pitfall is to add the weight and maturation factors, for example by using \(CL_{ref}\{(WT_i/WT_{ref})^\theta+F_{mat}(A_i)\}\). That expression does not preserve the intended reference interpretation and does not apply each factor as a proportional adjustment. Under the model taught here, use their product. A size-only model omits developmental change; a maturation-only model omits weight scaling. Each can be considered as a simpler alternative when scientifically justified, but neither includes both adjustments.

The combined equation is not a guarantee that predictions are accurate outside the range supported by data. Its parameters may be uncertain, the chosen age function may not represent the relevant biology, and a single maturation factor may be inadequate when several elimination processes contribute differently. The present calculations show how to apply a specified model; they do not select among maturation functions or establish that one is empirically preferable.

Key takeaway: A combined clearance prediction is the fully mature reference clearance multiplied by a dimensionless weight-scaling factor and a dimensionless maturation factor. Define the references and units first, compute each factor separately, and compare predictions using unrounded values.

Check Your Understanding

Answer each question, then check the solutions below.

  1. In the combined equation, what does \(CL_{ref}\) represent?
  2. If \(A=TM_{50}\), what is the value of the Hill maturation factor?
  3. Using \(\theta=0.75\), how does weight scaling change when weight increases from 10 kg to 20 kg?
  4. Why must age and \(TM_{50}\) use the same units?
  5. A calculation gives a maturation factor of 1.20 for a positive age using the stated Hill function. Name a likely problem.

Solutions. (1) \(CL_{ref}\) is clearance at \(WT_{ref}\) under fully mature function. (2) The numerator and denominator are equal, so \(F_{mat}=0.5000\). (3) The weight-factor ratio is \((20/10)^{0.75}=2^{0.75}\approx1.682\), so the weight contribution increases by a factor of about 1.682. (4) The Hill function compares age quantities in a ratio; mismatched units change the numerical ratio and therefore the predicted maturation. (5) For positive age and positive parameters the stated function is below 1, so check the formula, units, or inputs.