Why neonatal dosing changes quickly
For a neonate, a dose calculated from weight alone can miss important changes in drug disposition. Body size affects clearance and volume, but renal filtration, transporter activity, and metabolic pathways also mature over time. Because those processes can change appreciably over weeks, a dose that was appropriate at one age may not produce the same exposure later—even if weight has changed only modestly.
A useful calculation therefore distinguishes three age measures. Gestational age at birth (GA) is the duration of gestation when the infant was born. Postnatal age (PNA) is the time elapsed since birth. Postmenstrual age (PMA) is GA at birth plus PNA. For example, an infant born at 28 weeks’ gestation who is now 4 weeks old has a PMA of 32 weeks. PMA is often a useful age scale for maturation models, but it is not a universal substitute for PNA: the appropriate age predictor depends on the drug and the evidence supporting its model.
Building a clearance estimate
The exponent 0.75 is a commonly used allometric size exponent for clearance; it should be supported by the drug model rather than assumed to fit every neonatal population. A maturation function can be specified in several ways. One standard sigmoidal form is the Hill function:
Here \(TM_{50}\) is the PMA at which the maturation factor is 0.5, and \(h\) controls the steepness of the curve. The function increases with PMA for positive \(h\), approaches 0 at very low PMA, and approaches 1 at high PMA. Its parameters and age scale must be estimated or justified for the particular drug and pathway.
A single maturation factor may be inadequate when renal and metabolic clearances mature differently. In that case, model the pathways separately and add their clearances:
This equation assumes the pathway clearances are additive, each pathway follows the stated size relationship, and the maturation functions are appropriate for the population and age range. If a pathway is not mature at the reference age, or if disease, treatment, or another covariate changes its activity, the model needs to represent that evidence rather than treating the simple equation as complete.
Worked Example: Clearance at 32 Weeks’ PMA
Hypothetical model. Suppose a drug has a mature reference clearance of 6.0 L/h at 70 kg, follows a 0.75 clearance size exponent, and has \(MF(PMA)=PMA^3/(40^3+PMA^3)\), with PMA in weeks. An infant was born at 28 weeks, is now 4 weeks old, and weighs 1.50 kg. Estimate clearance.
First calculate PMA: \(28+4=32\) weeks. The size factor is:
The maturation factor is:
Using the unrounded factors for the product gives:
The result is a model-based estimate for this hypothetical drug and infant, not a recommended clinical dose. It is much lower than the mature 70-kg reference clearance because both the size factor and the maturation factor are below 1.
From clearance to a maintenance dose
For linear pharmacokinetics at steady state, the average concentration over a dosing interval is the dose rate divided by clearance. Thus, if the intended average concentration is \(C_{\mathrm{avg,target}}\), the maintenance dose rate is \(CL\,C_{\mathrm{avg,target}}\). For an interval \(\tau\), the corresponding dose is \(CL\,C_{\mathrm{avg,target}}\tau\). This calculation assumes the target is appropriate, clearance is stable during the interval, and the drug follows linear, time-invariant pharmacokinetics over the relevant range.
Units provide a useful check: \((\mathrm{L/h})(\mathrm{mg/L})(\mathrm{h})=\mathrm{mg}\). For the clearance estimated above, a hypothetical target average concentration of 10 mg/L implies a rate of \(0.1138\times10=1.138\) mg/h. Over 24 hours, that is approximately 27.31 mg. This is a mathematical illustration only; therapeutic targets and practical regimens must be drug-specific.
Dose size and dosing interval are related but not interchangeable. Clearance determines the dose rate needed to maintain an average concentration under the assumptions above. Volume of distribution, together with clearance, affects the concentration fluctuation between doses. In a one-compartment model with first-order elimination, the elimination rate constant is \(k=CL/V\), and the half-life is \(t_{1/2}=0.693/k\). A shorter half-life generally means more rapid decline between doses, but selecting an interval also requires the intended peak, trough, route, formulation, and clinical context. An average-concentration calculation alone does not establish a safe or suitable interval.
Worked Example: Separate Renal and Metabolic Maturation
Hypothetical model. Consider the 28-week-at-birth infant from the previous example, now at PMA 32 weeks and weighing 1.50 kg. Suppose the mature 70-kg reference clearances are 4.5 L/h for renal elimination and 1.5 L/h for metabolism. Let the renal maturation function be \(PMA^3/(40^3+PMA^3)\), and the metabolic function be \(PMA^2/(50^2+PMA^2)\). Estimate total clearance.
The renal maturation factor is 0.3386243, as calculated above. The metabolic maturation factor is:
The size factor is 0.0560073. Substituting the pathway clearances and factors:
The two pathways contribute different amounts because their reference clearances and maturation functions differ. Applying one shared maturation factor to the combined 6.0 L/h reference value would impose an assumption that both pathways mature identically; this example shows why that simplification can be inappropriate.
Updating an estimate as the infant grows
When a model uses both current weight and age, a reassessment requires updating both inputs. Holding the dose fixed while clearance rises can lower average exposure; conversely, changing the dose solely in proportion to weight can miss the maturation contribution. The next example uses the same hypothetical single-pathway model as the first example to isolate those effects.
Worked Example: Recalculating Clearance and Dose Rate
Suppose the infant born at 28 weeks is now 8 weeks old, weighs 2.00 kg, and the same hypothetical model applies. PMA is \(28+8=36\) weeks. The size and maturation factors are:
Using unrounded factors, the new clearance is:
At the same hypothetical target average concentration of 10 mg/L, the dose rate is \(0.1758\times10=1.758\) mg/h, or approximately \(1.758\times24=42.19\) mg per day. This is about 1.545 times the earlier model-based daily amount of 27.31 mg. The difference reflects both increased weight and maturation; it does not imply that an actual patient’s dose should be changed by this ratio without checking the model, target, and clinical data.
For an additional half-life illustration, suppose volume is \(0.7\) L/kg and scales linearly with weight. At 1.50 kg, \(V=1.05\) L, so \(t_{1/2}=0.693(1.05/0.1138)\approx6.40\) h. At 2.00 kg, \(V=1.40\) L, so \(t_{1/2}=0.693(1.40/0.1758)\approx5.52\) h. This hypothetical model predicts a shorter half-life at the later age because clearance has increased proportionally more than volume.
Assumptions, diagnostics, and alternatives
Before using a calculation, check that the model was developed for a relevant drug, route, gestational and postnatal age range, and clinical setting. Confirm whether the age variable is PMA, PNA, or another measure; whether weight is current and measured reliably; and whether renal and metabolic pathways are represented appropriately. The calculation also assumes that the infant is not undergoing a rapid change in physiology that makes a fixed clearance estimate unsuitable over the dosing interval.
Model diagnostics should look for systematic prediction errors across PMA, PNA, weight, and clinically relevant subgroups. If observed concentrations are consistently above or below predictions at particular ages, the maturation curve, its age scale, or the pathway assumptions may be inadequate. Unexpected concentrations can also reflect timing or sampling errors, inaccurate weight, changing renal function, interacting medicines, organ dysfunction, or adherence and administration problems. A maturation function should not be extrapolated far beyond the ages represented by its supporting data simply because its equation returns a number.
A weight-only calculation is simple but assumes that size captures the relevant differences in clearance; it can be biased when maturation is substantial. A categorical neonatal dosing table can be practical, but age or weight bands create abrupt boundaries and may conceal continuous changes. A model-based calculation can represent size and maturation continuously and distinguish pathways, but it depends on suitable data, model validation, and correctly measured covariates. These are alternatives with different trade-offs, not interchangeable guarantees of accurate dosing.
Neonatal dosing is especially sensitive to model limitations because maturation can be rapid and the clinical consequences of over- or underexposure can be serious. The equations here explain how a model-based estimate is constructed; they do not replace a drug-specific validated regimen, therapeutic drug monitoring where appropriate, or clinical assessment. When organ function or treatment changes, a model that assumes stable clearance may no longer describe the infant adequately.
Check Your Understanding
Try each problem, then compare your reasoning with the solutions.
- An infant was born at 30 weeks’ gestation and is now 6 weeks old. What is the PMA, and what is the distinction between PMA and PNA?
- Using \(MF(PMA)=PMA^3/(40^3+PMA^3)\), calculate the maturation factor at 40 weeks’ PMA.
- A hypothetical model estimates \(CL=0.0800\) L/h. What dose per 12-hour interval corresponds to a target average concentration of 5.0 mg/L under linear steady-state assumptions?
- Why might combining renal and metabolic clearance under one maturation function produce a misleading neonatal estimate?
- If clearance rises while volume remains unchanged, what happens to the elimination rate constant and half-life in a one-compartment model?
Solutions. (1) PMA is \(30+6=36\) weeks. PNA is time since birth (6 weeks), whereas PMA includes gestational age at birth. (2) \(MF(40)=40^3/(40^3+40^3)=0.5000\). (3) Dose \(=0.0800\ \mathrm{L/h}\times5.0\ \mathrm{mg/L}\times12\ \mathrm{h}=4.80\) mg per 12-hour interval. (4) The pathways may have different reference clearances and maturation rates; one factor would force them to mature identically. (5) Since \(k=CL/V\), increasing clearance at fixed volume increases \(k\); because \(t_{1/2}=0.693/k\), half-life decreases.