Tutorials › Nonlinear Mixed-Effects Models (NLME) › Two-Compartment Microconstants and Their Directions

One-, Two- and Three-Compartment Structural Models · Tutorial 25 of 797

Two-Compartment Microconstants and Their Directions

Learn how each two-compartment microconstant describes a directed first-order flow of drug amount, and why none is itself a clearance.

Beginner 9 min read

What You'll Learn

  • Identify the source and destination represented by k10, k12, and k21.
  • Write the two-compartment amount-balance equations from the directed flows.
  • Calculate an instantaneous amount flow from a rate constant and compartment amount.
  • Check units, signs, and total amount balance in a numerical example.
  • Distinguish a first-order rate constant from clearance and relate them when volume is specified.

Why Direction Matters

A two-compartment model describes drug amount in a central compartment and a peripheral compartment. The central compartment is the site of direct input and elimination in the standard model; drug can also move between the central and peripheral compartments in either direction. The three microconstants \(k_{10}\), \(k_{12}\), and \(k_{21}\) specify the rates of these directed processes. Their subscripts are not interchangeable labels: they tell us which compartment supplies the flow and where it goes.

A rate constant is not an amount flow by itself. It is a proportionality factor: multiply it by the amount currently in the source compartment to obtain the instantaneous amount transferred or eliminated per unit time. This distinction is essential when interpreting parameter values and writing amount-balance equations.

Definition: Label the central compartment 1 and the peripheral compartment 2. In the standard linear two-compartment model, \(k_{10}\) is the first-order elimination rate constant from compartment 1 to the outside of the modeled system; \(k_{12}\) is the first-order transfer rate constant from compartment 1 to compartment 2; and \(k_{21}\) is the first-order transfer rate constant from compartment 2 to compartment 1. The subscript 0 denotes an external elimination sink, not a third drug compartment.

From Directed Flows to Amount Balance

Let \(A_c(t)\) be the drug amount in the central compartment and \(A_p(t)\) the amount in the peripheral compartment at time \(t\). Assume each compartment is well mixed, the three rate constants are nonnegative and constant over the interval being modeled, and each process is first order in the amount in its source compartment. In this model, elimination occurs from the central compartment only. Any input is represented by \(R_{\mathrm{in}}(t)\), an amount entering the central compartment per unit time.

Under these assumptions, the instantaneous flows are \(k_{10}A_c\) out of the modeled system, \(k_{12}A_c\) from central to peripheral, and \(k_{21}A_p\) from peripheral to central. An amount balance adds flows entering a compartment and subtracts flows leaving it. Thus the central amount gains input and return from the peripheral compartment, while losing drug to elimination and peripheral distribution.

$$ \frac{dA_c}{dt} = R_{\mathrm{in}}(t)-k_{10}A_c-k_{12}A_c+k_{21}A_p $$

The peripheral compartment gains the central-to-peripheral flow and loses the return flow:

$$ \frac{dA_p}{dt}=k_{12}A_c-k_{21}A_p $$

The signs follow the physical directions. For example, \(k_{12}A_c\) is negative in the central balance because that amount leaves central, and positive in the peripheral balance because the same amount enters peripheral. It is not a loss from the whole body: it is an internal transfer.

Units and interpretation: If amounts are measured in mg and time in h, each \(k\) has units \(\mathrm{h}^{-1}\). The product \(k_{12}A_c\), for example, has units \(\mathrm{mg}/\mathrm{h}\) and is an instantaneous amount flow from central to peripheral. A value such as \(k_{12}=0.4\ \mathrm{h}^{-1}\) is not a flow of \(0.4\ \mathrm{mg}/\mathrm{h}\); the flow also depends on \(A_c\).

A Short Amount-Balance Check

Add the two compartment equations. The central-to-peripheral flow appears once with a minus sign and once with a plus sign, as does the return flow. Both cancel because these are transfers within the modeled system. The remaining balance is:

$$ \frac{d(A_c+A_p)}{dt} = R_{\mathrm{in}}(t)-k_{10}A_c $$

This result is a useful check on the meaning and placement of the constants. With no input, total modeled drug amount decreases only through \(k_{10}A_c\). It cannot decrease merely because drug moves from central to peripheral, or increase merely because drug returns from peripheral to central. This balance does not say that \(A_c\) must always decrease: return from the peripheral compartment can temporarily exceed central elimination and distribution.

Worked Example: Calculating the Three Directed Flows

Suppose, hypothetically, that at a particular instant \(A_c=80\ \mathrm{mg}\) and \(A_p=30\ \mathrm{mg}\), with \(k_{10}=0.10\ \mathrm{h}^{-1}\), \(k_{12}=0.25\ \mathrm{h}^{-1}\), and \(k_{21}=0.20\ \mathrm{h}^{-1}\). There is no input at that instant. Calculate the three flows and the instantaneous changes in the two amounts.

The elimination flow is \(k_{10}A_c=(0.10\ \mathrm{h}^{-1})(80\ \mathrm{mg})=8.0\ \mathrm{mg}/\mathrm{h}\), directed from central to outside the modeled system.

The central-to-peripheral flow is \(k_{12}A_c=(0.25\ \mathrm{h}^{-1})(80\ \mathrm{mg})=20\ \mathrm{mg}/\mathrm{h}\). The peripheral-to-central flow is \(k_{21}A_p=(0.20\ \mathrm{h}^{-1})(30\ \mathrm{mg})=6.0\ \mathrm{mg}/\mathrm{h}\).

With no input, the central balance is \(dA_c/dt=-8.0-20+6.0=-22\ \mathrm{mg}/\mathrm{h}\). The peripheral balance is \(dA_p/dt=20-6.0=14\ \mathrm{mg}/\mathrm{h}\). Therefore total amount changes at \(-22+14=-8.0\ \mathrm{mg}/\mathrm{h}\), exactly the negative elimination flow. The peripheral amount is increasing at this instant even though the total amount is decreasing.

Rate Constants Are Not Clearance

A rate constant has units of inverse time, such as \(\mathrm{h}^{-1}\). A clearance has units of volume per time, such as \(\mathrm{L}/\mathrm{h}\), and describes a volume of fluid cleared of drug per unit time. The two quantities are related, but they are not the same parameter and cannot be substituted for one another without a concentration-volume relationship.

For central elimination, let \(C_c\) be the central concentration and \(V_c\) the central volume, with \(A_c=V_cC_c\). The modeled elimination flow can be written either as \(k_{10}A_c\) or as clearance times concentration, \(\mathrm{CL}\,C_c\). Equating these expressions gives:

$$ k_{10}A_c = k_{10}V_cC_c = \mathrm{CL}\,C_c \qquad\Longrightarrow\qquad \mathrm{CL}=k_{10}V_c $$

This relationship applies when the central amount-concentration relationship uses \(V_c\), and \(\mathrm{CL}\) represents the clearance associated with the modeled central elimination process. The units confirm the distinction: \((\mathrm{h}^{-1})(\mathrm{L})=\mathrm{L}/\mathrm{h}\). A numerical value for \(k_{10}\) alone is not enough to give clearance; the relevant volume is also needed.

Likewise, \(k_{12}\) and \(k_{21}\) are directed transfer rate constants, not elimination clearances. They multiply the amount in their respective source compartments to give internal amount flow. Knowing one of these constants without the source amount does not determine the flow, and the constant itself does not have clearance units.

Worked Example: Converting Central Elimination to Clearance

For a hypothetical model, suppose \(k_{10}=0.12\ \mathrm{h}^{-1}\) and \(V_c=25\ \mathrm{L}\). Find the central elimination clearance and check that it agrees with the amount-flow expression when \(C_c=2.0\ \mathrm{mg}/\mathrm{L}\).

First, \(\mathrm{CL}=k_{10}V_c=(0.12\ \mathrm{h}^{-1})(25\ \mathrm{L})=3.0\ \mathrm{L}/\mathrm{h}\). The central amount at the stated concentration is \(A_c=V_cC_c=(25\ \mathrm{L})(2.0\ \mathrm{mg}/\mathrm{L})=50\ \mathrm{mg}\).

Using the rate-constant form, the elimination flow is \(k_{10}A_c=(0.12\ \mathrm{h}^{-1})(50\ \mathrm{mg})=6.0\ \mathrm{mg}/\mathrm{h}\). Using the clearance form, it is \(\mathrm{CL}\,C_c=(3.0\ \mathrm{L}/\mathrm{h})(2.0\ \mathrm{mg}/\mathrm{L})=6.0\ \mathrm{mg}/\mathrm{h}\). The matching flows confirm the conversion and show why \(k_{10}\) and clearance are not numerically interchangeable.

Reading the Subscripts and Checking a Model

The direction is easiest to remember by naming the source first: \(k_{12}\) takes amount from compartment 1 and sends it to compartment 2; \(k_{21}\) takes amount from compartment 2 and sends it to compartment 1. For \(k_{10}\), the source is compartment 1 and the destination is the external sink. In each case, the constant multiplies the amount in the source—not the destination.

Before using the equations, check that the model assumptions match the intended process. These equations describe linear, first-order flows with fixed rate constants and elimination from central only. They do not automatically describe saturable elimination, time-varying rate constants, or elimination from the peripheral compartment. If such a process is present, its amount-flow term must be represented in the model rather than silently assigned to one of these constants.

A practical diagnostic is to inspect the units and signs in each balance. Every term in an amount derivative must have units of amount per time. A transfer term must appear with opposite signs in the two compartment equations. With no input, adding the equations must leave only the external elimination term. Failure of any of these checks often indicates a reversed subscript, a missing flow, or a transfer incorrectly counted as elimination.

A common alternative description uses clearances and volumes rather than microconstants. That parameterization can be useful when the model is specified in terms of volume and clearance, but it does not change the physical direction of any flow. The microconstant form emphasizes the first-order dependence on source amount; a clearance form emphasizes flow proportional to concentration. Conversion requires the appropriate volume relationship, so the labels should not be treated as synonyms.

Common pitfall: It is incorrect to read \(k_{12}\) as a flow from compartment 2 to compartment 1, or to report \(k_{10}\) itself as a clearance. Read the source and destination from the subscripts, multiply the rate constant by the source amount to obtain a flow, and use a specified volume if converting a rate constant to clearance.

Worked Example: Checking Direction and Conservation

Consider a second hypothetical instant with no input, \(A_c=40\ \mathrm{mg}\), \(A_p=60\ \mathrm{mg}\), \(k_{10}=0.05\ \mathrm{h}^{-1}\), \(k_{12}=0.10\ \mathrm{h}^{-1}\), and \(k_{21}=0.15\ \mathrm{h}^{-1}\). Determine whether the peripheral amount is rising or falling and verify the total amount balance.

The central-to-peripheral flow is \((0.10)(40)=4.0\ \mathrm{mg}/\mathrm{h}\). The peripheral-to-central flow is \((0.15)(60)=9.0\ \mathrm{mg}/\mathrm{h}\). Therefore \(dA_p/dt=4.0-9.0=-5.0\ \mathrm{mg}/\mathrm{h}\): the peripheral amount is falling because more drug returns to central than enters peripheral.

Elimination is \((0.05)(40)=2.0\ \mathrm{mg}/\mathrm{h}\). The central balance is \(dA_c/dt=-2.0-4.0+9.0=3.0\ \mathrm{mg}/\mathrm{h}\). Adding the compartment changes gives \(3.0-5.0=-2.0\ \mathrm{mg}/\mathrm{h}\), equal to the negative elimination flow. Central amount is rising at this instant, while total modeled amount is falling.

Check Your Understanding

Answer each question, then compare with the solutions below.

  1. What flow does \(k_{21}\) describe, and which amount does it multiply?
  2. If \(k_{12}=0.30\ \mathrm{h}^{-1}\) and \(A_c=20\ \mathrm{mg}\), what is the central-to-peripheral flow?
  3. With no input, why do the two intercompartmental flows cancel when the compartment amount equations are added?
  4. For \(k_{10}=0.08\ \mathrm{h}^{-1}\) and \(V_c=15\ \mathrm{L}\), what is the associated central elimination clearance under the stated central amount-concentration relationship?

Solutions. (1) \(k_{21}\) describes transfer from peripheral compartment 2 to central compartment 1, and its flow is \(k_{21}A_p\). (2) The flow is \((0.30\ \mathrm{h}^{-1})(20\ \mathrm{mg})=6.0\ \mathrm{mg}/\mathrm{h}\). (3) Each internal transfer leaves one compartment and enters the other at the same rate, so its negative and positive terms sum to zero; only external input and elimination affect total modeled amount. (4) \(\mathrm{CL}=k_{10}V_c=(0.08\ \mathrm{h}^{-1})(15\ \mathrm{L})=1.2\ \mathrm{L}/\mathrm{h}\).