Why Renal Clearance Needs More Than a Size Adjustment
A child's kidneys are smaller than an adult's, but kidney size alone does not determine how quickly they filter or secrete a drug. Renal function also changes during development. Glomerular filtration rate (GFR) matures with age, while the capacity for active tubular secretion depends on the growth and functional development of the relevant transport pathways. These processes can mature on different schedules.
A useful PBPK representation therefore separates three contributions: body-size scaling, organ growth, and functional maturation. For filtration, body surface area (BSA) is a practical size term, paired with a GFR maturation factor. For secretion, kidney growth can scale the amount of functional tissue, while a separate maturation factor represents transporter capacity per unit tissue. These are modelling choices, not claims that BSA or kidney mass alone directly determines function.
A Separable Model for Filtration and Secretion
Let \(\mathrm{GFR}_{\mathrm{ref,index}}\) be a mature reference GFR indexed to 1.73 m2. For a child with BSA \(BSA_{\mathrm{child}}\), one simple model is:
This applies BSA scaling and functional maturation as separate factors. The BSA ratio is dimensionless, so the predicted GFR has the same units as the indexed reference, such as mL/min. The maturation factor should be supported by an appropriate developmental model or data; it should not be inserted a second time if the chosen reference or scaling already incorporates age-related maturation.
For secretion, define \(CL_{\mathrm{int,sec,ref}}\) as the mature reference intrinsic secretory clearance, in mL/min, relative to unbound plasma concentration. Let \(M_{\mathrm{kidney,child}}\) and \(M_{\mathrm{kidney,ref}}\) be the relevant total kidney masses. Under the simplifying assumption that capacity is proportional to kidney mass, a separable model is:
This uses organ growth to scale the amount of tissue and \(F_{\mathrm{sec}}(a)\) to scale its functional capacity per unit mass. Proportionality to kidney mass is a modelling assumption: it may not hold when nephron structure, transporter localization, or the distribution of functional tissue changes with age. If that assumption is unsuitable, replace the mass ratio with a better-supported measure of secretory tissue or capacity.
For a low drug concentration, assume filtration and secretion are linear and their contributions add. Let \(f_u\) be the unbound fraction in plasma, and assume it is constant over the concentration range considered. If tubular reabsorption is negligible, or has already been addressed separately, renal plasma clearance is:
The filtration term follows because filtration acts on unbound drug: the filtration rate is \(f_u C_p\,\mathrm{GFR}\), where \(C_p\) is total plasma concentration. The linear secretory rate is \(f_u C_p\,CL_{\mathrm{int,sec}}\). Dividing each rate by \(C_p\) gives the corresponding contribution to clearance. This convention defines \(CL_{\mathrm{int,sec}}\) relative to unbound plasma concentration; the factor \(f_u\) converts it to clearance relative to total plasma concentration.
Worked Example: Scaling GFR in a Hypothetical Infant
Suppose a hypothetical infant has BSA \(0.40\ \mathrm{m}^2\), and an assumed \(F_{\mathrm{GFR}}(a)=0.45\). Use a mature reference GFR of \(120\ \mathrm{mL/min/1.73\ m^2}\). Assume this indexed reference is mature, BSA is the selected size-scaling term, and the maturation factor is not already included in the reference.
The model gives:
The unadjusted BSA-scaled value is \(120(0.40/1.73)=27.75\ \mathrm{mL/min}\). Multiplication by the assumed maturation factor reduces it to \(12.49\ \mathrm{mL/min}\). This result is conditional on the chosen reference, BSA, and maturation factor; it is not a measured value for an actual infant.
Worked Example: Scaling Secretory Capacity Separately
For a hypothetical child, suppose total kidney mass is \(40\ \mathrm{g}\), compared with a mature reference of \(300\ \mathrm{g}\), and \(F_{\mathrm{sec}}(a)=0.60\). Let mature reference intrinsic secretory clearance be \(150\ \mathrm{mL/min}\). Assume capacity is proportional to total kidney mass and the age-specific maturation factor applies per unit mass.
First calculate the kidney-mass ratio: \(40/300=0.1333\). Then:
The result is an intrinsic clearance relative to unbound plasma concentration. It is not yet the secretion contribution to clearance relative to total plasma concentration; that contribution also depends on \(f_u\). Keeping these quantities distinct avoids multiplying by binding twice or omitting it entirely.
Combining the Renal Contributions
The combined model follows from adding rates, not from adding maturation factors. Filtration and active secretion are parallel routes of removal under the stated linear assumptions. Their rates can be added, and dividing the combined rate by total plasma concentration yields the sum of their clearances.
A short derivation makes the assumptions explicit. Write total plasma concentration as \(C_p\). The filtration rate is \(f_u C_p\,\mathrm{GFR}\), and, in the linear range, the secretory rate is \(f_u C_p\,CL_{\mathrm{int,sec}}\). Their sum is the total renal removal rate when there is no additional reabsorption term. Dividing by \(C_p\) gives:
If net reabsorption is important, it cannot simply be ignored: represent its contribution with a suitable mechanism or data-supported net-clearance term. Likewise, secretion may not remain linear when transporter-mediated flux approaches saturation. The additive equation is not a universal renal-clearance identity; it is the result for this specified model.
Worked Example: Calculating Combined Renal Clearance
Use the hypothetical infant values from the first example, \(\mathrm{GFR}=12.49\ \mathrm{mL/min}\), and the hypothetical secretory clearance from the second example, \(CL_{\mathrm{int,sec}}=12.00\ \mathrm{mL/min}\). Suppose \(f_u=0.20\), constant over the modelled range, and assume negligible reabsorption and linear secretion.
The filtration contribution relative to total plasma concentration is:
The secretory contribution is:
Thus:
Under these assumptions, the predicted renal plasma clearance is \(4.898\ \mathrm{mL/min}\). Filtration contributes slightly more than secretion in this calculation. That comparison is specific to the hypothetical values and should not be generalized to other drugs or ages.
Checking the Model and Choosing an Alternative
Before using the combined result in a PBPK simulation, check whether the input quantities and assumptions match the drug, age range, and intended output. In particular, confirm that the reference GFR is indexed in the units used by the BSA ratio, that the selected maturation factor corresponds to the same developmental-age definition, and that the kidney-mass ratio and secretory reference describe compatible tissue and capacity. An implausible discontinuity across ages or a predicted clearance inconsistent with available renal-function or excretion data signals a need to re-examine these choices.
The low-concentration secretion assumption also needs a check. For a single transporter described by Michaelis–Menten kinetics, secretory rate is \(V_{\max}C_u/(K_m+C_u)\), where \(C_u=f_u C_p\). At concentrations much smaller than \(K_m\), the rate is approximately \((V_{\max}/K_m)C_u\), so \(CL_{\mathrm{int,sec}}=V_{\max}/K_m\) is approximately constant. At concentrations approaching or exceeding \(K_m\), this linear approximation fails: the apparent secretion clearance depends on concentration. A PBPK model must then represent the saturable process rather than treating the linear clearance as constant.
A single \(F_{\mathrm{sec}}(a)\) is also a simplification when several transporters contribute. Tubular secretion requires coordinated steps, including movement into and out of tubular cells. If those steps have different age profiles, a single factor may misrepresent the effective pathway. More detailed transporter-specific or pathway-specific representations may be preferable when suitable data are available. A measured age-specific net renal clearance can be a useful empirical alternative, but it may not separate filtration from secretion and may not transfer reliably to different binding conditions or drug concentrations.
| Approach | What it represents | Key limitation |
|---|---|---|
| Separated filtration and secretion model | Size and maturation factors applied to distinct renal pathways | Needs pathway-specific inputs and explicit assumptions about binding, reabsorption, and linearity |
| Age-specific empirical renal clearance | Observed net clearance at specified ages and conditions | May combine processes that respond differently to age, concentration, or binding |
| Saturable secretion model | Concentration-dependent transporter-mediated flux | Requires parameters such as \(V_{\max}\) and \(K_m\), and appropriate developmental scaling |
These approaches are not interchangeable in every application. The separated model is useful when the aim is to represent filtration and secretion as distinct mechanisms and the required parameters are defensible. Empirical clearance may be more appropriate when the available evidence supports the net process but not its components. A saturable representation is needed when concentrations make the linear approximation untenable.
Common Mistake: Treating Organ Growth as Maturation
A common mistake is to multiply adult secretory clearance by a kidney-mass ratio and assume that this fully captures development. This assumes that every unit of kidney tissue has mature transporter capacity at every age. The opposite mistake is to apply a maturation factor without any size term, implicitly giving a small child's kidneys the adult amount of functional tissue. The correction is to state and justify separate size and functional factors, then combine them only as specified by the model. A second frequent error is to apply \(f_u\) to a clearance that is already defined relative to total plasma concentration; check the clearance definition before applying binding.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- A child's BSA is one-half of the reference value, and the assumed GFR maturation factor is 0.40. What fraction of indexed mature reference GFR does the model predict?
- Kidney mass is one-fifth of the reference mass and \(F_{\mathrm{sec}}(a)=0.50\). What fraction of mature reference intrinsic secretory clearance is predicted under the proportional-mass model?
- Why is \(F_{\mathrm{GFR}}(a)\) not automatically interchangeable with \(F_{\mathrm{sec}}(a)\)?
- A drug has \(f_u=0.10\), GFR \(=20\ \mathrm{mL/min}\), and \(CL_{\mathrm{int,sec}}=30\ \mathrm{mL/min}\). Under the linear assumptions and with negligible reabsorption, calculate \(CL_R\).
- What changes in the model if transporter-mediated secretion is near saturation?
Solutions. (1) The predicted fraction is \((1/2)(0.40)=0.20\), or 20% of indexed mature reference GFR. (2) It is \((1/5)(0.50)=0.10\), or 10% of mature reference intrinsic secretory clearance. (3) Filtration and transporter function can have different developmental determinants and trajectories, so their maturation factors must be supported separately. (4) \(CL_R=0.10(20)+0.10(30)=2+3=5.00\ \mathrm{mL/min}\). (5) A constant linear \(CL_{\mathrm{int,sec}}\) no longer describes secretion; represent the concentration-dependent flux, for example with a Michaelis–Menten model, and scale its capacity developmentally as appropriate.