Steady-State Concentration During Constant Infusion

Why Constant Infusion Has a Capacity Limit

With a constant infusion, drug enters the body at a fixed rate. Under Michaelis-Menten elimination, however, the elimination rate increases with concentration only up to a limiting capacity. The balance between these two rates determines whether concentration can settle at a finite steady state. As the infusion rate approaches the elimination capacity, the predicted steady-state concentration rises sharply.

The earlier material on the Michaelis-Menten elimination rate law defines the saturable rate in terms of \(V_{\max}\) and \(K_m\). Here the aim is to apply that law to constant infusion, derive the steady-state solution, and identify what the model implies when the infusion rate is too high to be balanced by elimination.

Model and assumptions: Let \(C\) be the well-mixed drug concentration in a fixed-volume compartment. Drug is infused at a constant rate \(R_0\), measured in amount per time, and eliminated through one Michaelis-Menten pathway with rate \(V_{\max}C/(K_m+C)\). Here \(V_{\max}\) has units of amount per time, and \(K_m\) and \(C\) have the same concentration units. Assume the parameters and compartment volume remain constant, the infusion reaches the compartment without delay, and there are no other input or elimination processes in this model.

Deriving the Steady-State Concentration

A concentration balance is the input rate minus the elimination rate, divided by the fixed compartment volume \(V\). Thus the concentration follows

\frac{dC}{dt} = \frac{1}{V}\left(R_0-\frac{V_{\max}C}{K_m+C}\right).

A steady state is a concentration at which \(dC/dt=0\), so the infusion rate must equal the elimination rate. For a positive infusion rate, solve that balance as follows:

R_0=\frac{V_{\max}C_{\mathrm{ss}}}{K_m+C_{\mathrm{ss}}} \quad\Longrightarrow\quad R_0(K_m+C_{\mathrm{ss}})=V_{\max}C_{\mathrm{ss}} \quad\Longrightarrow\quad C_{\mathrm{ss}}=\frac{K_mR_0}{V_{\max}-R_0}.

This algebra gives a finite, nonnegative concentration only when \(0\leq R_0<V_{\max}\). At \(R_0=0\), the solution is \(C_{\mathrm{ss}}=0\). For a positive infusion, the condition is \(0<R_0<V_{\max}\). The denominator is not a mathematical detail to ignore: it expresses the finite elimination capacity that limits which infusion rates can be balanced at equilibrium.

Steady-state result: Under the stated single-compartment model, a positive constant infusion has a finite steady-state concentration if and only if \(0<R_0<V_{\max}\). That concentration is \(C_{\mathrm{ss}}=K_mR_0/(V_{\max}-R_0)\). If \(R_0\geq V_{\max}\), no finite steady-state concentration exists.

The equilibrium is stable when it exists. The net rate of change in concentration decreases as concentration increases, because

\frac{d}{dC}\left(R_0-\frac{V_{\max}C}{K_m+C}\right) = -\frac{V_{\max}K_m}{(K_m+C)^2}<0 \qquad (C\geq 0,\ K_m>0).

At concentrations below the equilibrium, elimination is less than the infusion rate, so concentration rises. Above it, elimination exceeds the infusion rate, so concentration falls. This sign change means the equilibrium is attracting in this model. The compartment volume affects how quickly concentration changes, but it does not appear in the steady-state balance and therefore does not alter \(C_{\mathrm{ss}}\).

What Happens at and Above \(V_{\max}\)?

At every finite concentration, the Michaelis-Menten elimination rate is less than \(V_{\max}\). It approaches \(V_{\max}\) as concentration becomes very large, but does not reach it at any finite concentration. Consequently, when \(R_0=V_{\max}\), infusion still exceeds elimination at every finite \(C\). The difference gets smaller as \(C\) grows, but it never becomes zero at a finite concentration.

When \(R_0>V_{\max}\), elimination cannot even approach the input rate: the net amount entering per unit time remains greater than \(R_0-V_{\max}\). Concentration therefore continues to rise in this idealized model. The distinction matters: at equality, the positive net rate diminishes toward zero as concentration increases; above capacity, the net rate remains bounded away from zero.

The formula also shows why a rate just below capacity can imply a very high concentration. Rewriting the result in terms of the fraction of maximum capacity used gives

C_{\mathrm{ss}} = K_m\frac{R_0/V_{\max}}{1-R_0/V_{\max}}.

For example, an infusion at 90% of \(V_{\max}\) gives \(C_{\mathrm{ss}}=9K_m\), whereas an infusion at 99% gives \(C_{\mathrm{ss}}=99K_m\). The concentration does not increase in proportion to infusion rate near the capacity limit. This is a property of the model, not evidence that \(K_m\) is a hard boundary beyond which saturation suddenly begins.

Worked Example: Calculating a Finite Steady State

Suppose a hypothetical drug is infused at \(R_0=30\) mg/h. Its elimination parameters are \(V_{\max}=50\) mg/h and \(K_m=2\) mg/L. Assume the single-compartment model described above. First check the capacity condition: \(30\) mg/h is less than \(50\) mg/h, so a finite steady state is possible.

Substitute the values into the derived equation:

C_{\mathrm{ss}} = \frac{(2\ \mathrm{mg/L})(30\ \mathrm{mg/h})} {50\ \mathrm{mg/h}-30\ \mathrm{mg/h}} = \frac{60}{20}\ \mathrm{mg/L} = 3.00\ \mathrm{mg/L}.

Check the result in the original rate balance. At \(3.00\) mg/L, the modeled elimination rate is

\frac{V_{\max}C_{\mathrm{ss}}}{K_m+C_{\mathrm{ss}}} = \frac{(50\ \mathrm{mg/h})(3.00\ \mathrm{mg/L})} {2\ \mathrm{mg/L}+3.00\ \mathrm{mg/L}} = 30.0\ \mathrm{mg/h}.

The elimination rate matches the infusion rate, as required. The calculated steady-state concentration is 3.00 mg/L under the stated hypothetical assumptions.

Worked Infusion Problems

Worked Example: Infusion Close to the Elimination Capacity

Consider a second hypothetical drug with \(V_{\max}=50\) mg/h and \(K_m=2\) mg/L, infused at \(R_0=45\) mg/h. Since \(45<50\) mg/h, a finite steady state exists, but the infusion uses 90% of the modeled maximum elimination capacity.

The steady-state concentration is

C_{\mathrm{ss}} = \frac{(2\ \mathrm{mg/L})(45\ \mathrm{mg/h})} {50\ \mathrm{mg/h}-45\ \mathrm{mg/h}} = \frac{90}{5}\ \mathrm{mg/L} = 18.0\ \mathrm{mg/L}.

This is \(9K_m\), as expected for an infusion rate at 90% of \(V_{\max}\). A direct check gives an elimination rate of \((50\ \mathrm{mg/h})(18\ \mathrm{mg/L})/(2+18\ \mathrm{mg/L})=45\) mg/h. The result illustrates the sensitivity of concentration to the remaining capacity: a modest difference between infusion rate and \(V_{\max}\) appears in the denominator and produces a much larger concentration than in the first example.

Worked Example: Finding the Infusion Rate for a Target Concentration

Suppose a hypothetical model has \(V_{\max}=60\) mg/h and \(K_m=3\) mg/L. Find the constant infusion rate that balances elimination at a specified concentration of \(10\) mg/L. This is a model calculation, not a clinical dosing recommendation.

At steady state, the infusion rate must equal the Michaelis-Menten elimination rate at the target concentration:

R_0 = \frac{V_{\max}C_{\mathrm{ss}}}{K_m+C_{\mathrm{ss}}} = \frac{(60\ \mathrm{mg/h})(10\ \mathrm{mg/L})} {3\ \mathrm{mg/L}+10\ \mathrm{mg/L}} = \frac{600}{13}\ \mathrm{mg/h} = 46.15\ \mathrm{mg/h}.

The calculated rate is less than \(V_{\max}=60\) mg/h, consistent with the condition for a finite steady state. Substituting \(R_0=600/13\) mg/h into the concentration formula also returns \(10\) mg/L. More generally, a desired finite concentration can be translated to a balancing infusion rate using the elimination rate law; the assumptions of that law still need to be appropriate for the application.

Assumptions, Checks, and Alternatives

Before interpreting a calculated concentration, verify that the model represents the situation being analyzed. In particular, \(R_0\) and \(V_{\max}\) must both be expressed as rates of amount per time, while \(K_m\) and \(C\) must use compatible concentration units. A calculation that compares an infusion rate in mg/h with \(V_{\max}\) in mg/min has not established whether the capacity condition holds. Convert units before comparing or substituting.

The derivation also assumes one well-mixed compartment and one elimination pathway. In a PBPK application, the concentration at an enzyme site may not equal a measured plasma concentration, and transport or distribution may prevent the one-compartment balance from representing the system. Additional inputs or elimination routes also change the balance. A predicted concentration that is implausibly high, or an infusion rate estimated at or above \(V_{\max}\), should prompt review of the model assumptions and parameter units rather than blind use of the formula.

A useful comparison is the low-concentration approximation. When \(C\) is much smaller than \(K_m\), \(K_m+C\) is approximately \(K_m\), so the elimination rate is approximately \((V_{\max}/K_m)C\), a first-order rate. Under that approximation, the steady-state concentration is approximately \(K_mR_0/V_{\max}\). The exact expression is larger for positive \(R_0\) below capacity, because its denominator is \(V_{\max}-R_0\), not \(V_{\max}\). The approximation is most reasonable when \(R_0/V_{\max}\) is small and the resulting concentration is indeed small relative to \(K_m\).

The analytical equation is convenient when the stated model applies: it directly gives the finite equilibrium and makes the capacity condition explicit. Numerical integration of the concentration-balance differential equation is an alternative when the model includes more complex inputs or rate laws, but simulation does not remove the need to check whether a finite equilibrium is possible. A more detailed PBPK model may be needed when distribution between tissues or site-specific concentrations are central to the question.

Common Mistakes and Diagnostics

Ignoring the capacity condition. Substituting \(R_0\geq V_{\max}\) into the formula can produce a zero or negative denominator and an invalid finite concentration. Check \(R_0<V_{\max}\) before interpreting the result.
Assuming a steady state exists at \(R_0=V_{\max}\). Equality of the infusion rate to the limiting capacity is not attained at finite concentration. At any finite \(C\), the elimination rate is still below \(V_{\max}\).
Treating the steady-state concentration as proportional to infusion rate. That approximation is only appropriate in the low-concentration, near-first-order range. Near \(V_{\max}\), the exact concentration rises sharply.
Mixing amount and concentration rates or incompatible units. \(V_{\max}\) and \(R_0\) must share rate units; \(K_m\) and \(C\) must share concentration units. The concentration \(C_{\mathrm{ss}}\) is not an amount in the compartment.
Applying the single-pathway result to a more complex system without checking it. If another input, elimination route, or important distribution process is present, the balance used in the derivation may not describe the system.

Key takeaway: A constant infusion can be balanced by Michaelis-Menten elimination at a finite concentration only when its rate is below \(V_{\max}\). The steady-state concentration is \(K_mR_0/(V_{\max}-R_0)\), so predictions become highly sensitive to infusion rate as it approaches the elimination capacity.

Check Your Understanding

Answer each question, then compare with the solutions below.

A hypothetical drug has \(V_{\max}=40\) mg/h, \(K_m=4\) mg/L, and \(R_0=20\) mg/h. What is its predicted steady-state concentration?
For the same \(V_{\max}\) and \(K_m\), what does the model imply if \(R_0=40\) mg/h? Explain why.
A hypothetical infusion rate is 80% of \(V_{\max}\). Express the predicted steady-state concentration in terms of \(K_m\).
For \(V_{\max}=72\) mg/h, \(K_m=6\) mg/L, and \(C_{\mathrm{ss}}=6\) mg/L, calculate the balancing infusion rate.
Why does compartment volume affect the time course toward steady state but not the steady-state concentration in this model?

Solutions. (1) The infusion is below capacity, so \(C_{\mathrm{ss}}=(4\ \mathrm{mg/L})(20\ \mathrm{mg/h})/(40-20\ \mathrm{mg/h})=4.00\) mg/L. (2) There is no finite steady state: at every finite concentration, elimination is below \(V_{\max}\), so it cannot exactly balance an infusion of \(40\) mg/h. (3) With \(R_0/V_{\max}=0.8\), \(C_{\mathrm{ss}}=K_m(0.8)/(1-0.8)=4K_m\). (4) \(R_0=(72\ \mathrm{mg/h})(6\ \mathrm{mg/L})/(6+6\ \mathrm{mg/L})=36.0\) mg/h. (5) Volume divides the net rate of amount input into a rate of concentration change in the differential equation. At steady state, the net rate of amount input is zero, so volume does not enter the balance that determines concentration.