Tutorials › Nonlinear Mixed-Effects Models (NLME) › Worked Example: Multiple Maturation Pathways

Allometric Scaling and Maturation Functions · Tutorial 399 of 797

Worked Example: Multiple Maturation Pathways

Learn to combine pathway-specific maturation with allometric scaling, calculate total clearance, and interpret how renal and hepatic contributions change with age.

Intermediate 10 min read

What You'll Learn

  • Calculate renal and hepatic clearance separately using pathway-specific maturation functions.
  • Combine pathway clearances correctly when elimination routes act in parallel.
  • Compute route contributions and fractions without rounding inconsistencies.
  • Compare age effects while distinguishing maturation from body-size scaling.
  • Identify assumptions, diagnostics, and limitations of a multiple-pathway calculation.

Why Separate Maturation Pathways?

A drug eliminated by both the kidneys and the liver need not have one common developmental profile. Renal and hepatic clearance can each increase with age, but at different rates. Applying one maturation factor to the drug’s entire clearance can therefore misrepresent both total clearance and the contribution of each route. The calculation is clearer when each pathway is scaled and matured separately, and the resulting clearances are then added.

This tutorial works through a hypothetical drug with renal and hepatic elimination. The calculations use a simple allometric size relationship and distinct maturation functions. They illustrate how to combine those components; the numerical parameter values are assumed for teaching and are not estimates from a real drug or study.

Definition: For parallel elimination pathways, total systemic clearance is the sum of the pathway clearances. In the model used here, pathway \(j\) has adult reference clearance \(CL_{j,70}\) at a body weight of 70 kg, a common allometric exponent of 0.75, and a pathway-specific maturation fraction \(M_j(a)\), where \(a\) is age in years: $$ CL_j(a,WT)=CL_{j,70}\left(\frac{WT}{70\ \mathrm{kg}}\right)^{0.75}M_j(a), \qquad CL_{\mathrm{total}}(a,WT)=\sum_j CL_j(a,WT). $$ Here \(WT\) is body weight in kg; each clearance is in L/h. \(M_j(a)\) is dimensionless and ranges from 0 to 1 in this example.

Model and Assumptions for the Example

Suppose the hypothetical drug has a mature renal reference clearance of 4 L/h and a mature hepatic reference clearance of 16 L/h, both defined at 70 kg. These reference values are already pathway-specific: they are not total clearance values to which each route is subsequently applied. The model uses the same weight exponent for both routes, but allows renal and hepatic maturation to differ.

$$ M_R(a)=\frac{a}{a+1}, \qquad M_H(a)=\frac{a}{a+6}, $$

The renal maturation function has \(TM_{50}=1\) year and the hepatic function has \(TM_{50}=6\) years. In each case, \(TM_{50}\) is the age at which that pathway’s maturation fraction is 0.5: substituting \(a=1\) into \(M_R\), or \(a=6\) into \(M_H\), gives 0.5. Both functions approach 1 as age increases. For this hypothetical calculation, age means postnatal age in years, and the functions are treated as applicable over the ages shown.

The calculation assumes that the two routes contribute additively to systemic elimination, that the reference clearances are appropriate mature values, and that body size affects both pathways through the stated common exponent. It also assumes that each pathway’s maturation depends on age through its specified function, with no additional age-dependent change in the allometric exponent. These are modeling assumptions, not consequences of the arithmetic. In an application, physiological evidence and data must support them.

Key calculation: At a given age and weight, first calculate the shared size factor, then apply the appropriate maturation fraction to each pathway’s reference clearance, and finally sum the pathway clearances. Do not average the two maturation functions and apply that average to the total unless the resulting model is justified for the intended purpose.

Why Pathway Clearances Add

The additivity follows from the definition of clearance. At a common systemic concentration \(C\), the elimination rate through pathway \(j\) is \(C\,CL_j\). If renal and hepatic elimination occur as parallel routes, their rates add:

$$ \text{total elimination rate} = C\,CL_R+C\,CL_H = C(CL_R+CL_H). $$

Total clearance is the total elimination rate divided by the same concentration \(C\). Therefore \(CL_{\mathrm{total}}=CL_R+CL_H\). This reasoning depends on the pathway clearances being defined on the same systemic basis and at the same covariate values. The fraction of total clearance contributed by route \(j\) is \(CL_j/CL_{\mathrm{total}}\), a dimensionless quantity.

With a common allometric exponent, the size factor is shared across routes. It changes the magnitude of both clearances at a given weight, but cancels from their ratio. Consequently, in this particular model, the renal and hepatic fractions depend on age through the matured reference clearances \(CL_{R,70}M_R(a)\) and \(CL_{H,70}M_H(a)\). If the pathways had different size exponents, that cancellation would not hold.

Worked Example: Total Clearance at Four Years

Calculate the two pathway clearances and total clearance for a hypothetical 4-year-old child weighing 12 kg. Use the parameters specified above.

Step 1: Calculate the size factor. The weight ratio is \(12/70\), so

$$ \left(\frac{12}{70}\right)^{0.75}=0.266417289\ldots $$

Step 2: Evaluate maturation. At age \(a=4\) years,

$$ M_R(4)=\frac{4}{4+1}=0.8, \qquad M_H(4)=\frac{4}{4+6}=0.4. $$

Step 3: Calculate each pathway and add. Keeping the size factor through the calculation gives

$$ CL_R=4(0.266417289\ldots)(0.8) =0.852535325\ldots\ \mathrm{L/h}, $$
$$ CL_H=16(0.266417289\ldots)(0.4) =1.705070650\ldots\ \mathrm{L/h}. $$

Thus

$$ CL_{\mathrm{total}} =0.852535325\ldots+1.705070650\ldots =2.557605975\ldots\ \mathrm{L/h}. $$

The total is 2.55761 L/h at the precision shown, or 2.558 L/h to four significant figures. The renal share, calculated from the unrounded pathway clearances, is \(0.852535325\ldots/2.557605975\ldots=1/3\), or 33.3%. The hepatic share is \(2/3\), or 66.7%. The pathway fractions sum to 100%, as they should.

Comparing Pathway Contributions Across Ages

A useful way to see the effect of the distinct maturation rates is to compare the two modeled pathways at the same reference weight, 70 kg. This is a standardized calculation to isolate age-related maturation; it is not a claim that one individual maintains that weight at all ages. At 70 kg, the size factor is exactly 1, so each pathway clearance is its reference value multiplied by its maturation fraction.

AgeRenal maturationRenal clearance (L/h)Hepatic maturationHepatic clearance (L/h)Total (L/h)
0.5 years0.33331.3330.076921.2312.564
4 years0.80003.2000.40006.4009.600
12 years0.92313.6920.666710.6714.36

For example, at 0.5 years the renal fraction is \(1.3333/2.5641=0.5200\), or 52.0%, while the hepatic fraction is 48.0%. Although renal maturation is more advanced at this age, the assumed mature hepatic reference clearance is four times the renal reference clearance. The pathway with the larger maturation fraction therefore need not have the larger clearance. At 4 years the renal and hepatic clearances are 3.2 and 6.4 L/h, respectively, so their shares are one-third and two-thirds. These proportions depend on both the reference clearances and the maturation functions.

Worked Example: A Younger Child and the Route Fractions

For a second hypothetical calculation, consider age 0.5 years and weight 17.5 kg. Find each clearance, total clearance, and the percentage contributed by each route.

The size factor is

$$ \left(\frac{17.5}{70}\right)^{0.75} =(0.25)^{0.75} =0.353553391. $$

At 0.5 years, the maturation fractions are

$$ M_R(0.5)=\frac{0.5}{1.5}=\frac{1}{3}, \qquad M_H(0.5)=\frac{0.5}{6.5}=\frac{1}{13}. $$

Substituting these values gives

$$ CL_R=4(0.353553391)\left(\frac{1}{3}\right) =0.471404521\ \mathrm{L/h}, $$
$$ CL_H=16(0.353553391)\left(\frac{1}{13}\right) =0.435911865\ \mathrm{L/h}. $$

Therefore \(CL_{\mathrm{total}}=0.471404521+0.435911865=0.907316386\) L/h, or 0.9073 L/h to four significant figures. For the renal share, the common size factor cancels:

$$ \frac{CL_R}{CL_{\mathrm{total}}} =\frac{4/3}{4/3+16/13} =\frac{52}{100} =0.5200. $$

The renal contribution is 52.0%, and the hepatic contribution is 48.0%. The hepatic route has the lower maturation fraction, but its larger mature reference clearance makes its absolute clearance comparable to the renal clearance in this example.

Worked Example: Nearing the Mature Reference Values

At age 12 years and a standardized weight of 70 kg, calculate total clearance and pathway shares. This isolates maturation from the effect of a particular child’s weight.

The maturation fractions are \(M_R(12)=12/13=0.9230769\) and \(M_H(12)=12/18=2/3\). Since the size factor at 70 kg is 1,

$$ CL_R=4\left(\frac{12}{13}\right)=3.6923077\ \mathrm{L/h}, \qquad CL_H=16\left(\frac{2}{3}\right)=10.6666667\ \mathrm{L/h}. $$

The total is \(3.6923077+10.6666667=14.3589744\) L/h, or 14.36 L/h. The renal share is \(3.6923077/14.3589744=0.2571\), or 25.71%; the hepatic share is 74.29%. Both pathways are closer to their mature reference values than at younger ages, but neither has reached its asymptote in this model. The mature total at 70 kg is \(4+16=20\) L/h.

Diagnostics, Alternatives, and Limitations

A calculation can be arithmetically correct yet represent an unsuitable model. Check each pathway separately before interpreting the total:

  • Check units and reference values. Both reference clearances must be in the same clearance units and defined at the same reference weight. The maturation functions are dimensionless, and the weight ratio must compare weights in the same units.
  • Check the maturation values. For the functions used here and nonnegative age, each fraction should be between 0 and 1 and should increase with age. A value above 1 or a decreasing value indicates an input or implementation error for this model.
  • Check the sum and shares. Total clearance must equal the sum of pathway clearances. Nonnegative route shares should sum to 1 (or 100%), allowing only small discrepancies from rounding.
  • Check the age scale and domain. The example defines age in years from birth. Using a different age scale without redefining the maturation parameters changes the calculation. A simple curve should not be extrapolated beyond a physiologically supported domain merely because it gives a numerical answer.
  • Compare predictions with evidence. Implausible pathway contributions or a poor fit across ages may indicate unsuitable reference values, maturation functions, or the assumption of a common size exponent. The equations alone cannot establish physiological validity.

One alternative is a single maturation function applied to total clearance. That can be a reasonable simpler model if evidence supports a common developmental pattern for the aggregate clearance. It does not, by itself, describe distinct renal and hepatic contributions. Another alternative is to use different functional forms or size exponents for different routes, but that adds parameters and requires adequate data or physiological justification. The present example fixes all parameters, so it demonstrates calculation rather than parameter estimation, uncertainty quantification, or proof that one model is preferable.

Common mistake: “Renal maturation is 0.8 and hepatic maturation is 0.4, so total maturation is 0.6; I can multiply total mature clearance by 0.6.” This takes an unweighted average and generally changes the answer. At age 4, the correct calculation is \(4(0.8)+16(0.4)=9.6\) L/h at 70 kg. Applying 0.6 to the total mature clearance gives \(0.6(4+16)=12\) L/h, which is not the modeled total. The correction is to mature each pathway using its own function and reference clearance, then add the resulting clearances.

Check Your Understanding

Use the hypothetical pathway model in this tutorial. Give units where appropriate; solutions follow the questions.

  1. At 70 kg and age 1 year, calculate renal, hepatic, and total clearance.
  2. At age 6 years, what are the two maturation fractions? Explain why they are not equal.
  3. At age 4 years and 12 kg, what is the renal percentage contribution using unrounded pathway clearances?
  4. Why does the shared size factor cancel when calculating pathway shares in this example?
  5. A calculation gives negative hepatic clearance but positive renal clearance. Name one model check that should be made before interpreting the total.

Solutions. 1. At 1 year, \(M_R=1/2\) and \(M_H=1/7\). At 70 kg, \(CL_R=4(1/2)=2\) L/h and \(CL_H=16(1/7)=2.286\) L/h, so total clearance is 4.286 L/h. 2. \(M_R(6)=6/7=0.8571\), while \(M_H(6)=6/12=0.5000\). The functions have different \(TM_{50}\) values, 1 and 6 years. 3. The pathway clearances are proportional to \(4(0.8)=3.2\) and \(16(0.4)=6.4\); the renal fraction is \(3.2/(3.2+6.4)=1/3\), or 33.3%. 4. Both clearances include the same positive factor \((WT/70)^{0.75}\); it appears in numerator and denominator of the fraction and cancels. 5. Check the inputs, units, and implementation of the maturation function. In this model, for nonnegative age, \(M_H(a)\) is between 0 and 1, and the reference clearance, size factor, and maturation fraction are all nonnegative. A negative result is inconsistent with those assumptions and should not be interpreted as physiological clearance.