Tutorials › Pharmacokinetic-Pharmacodynamic (PK/PD) Modeling › Deriving the First-Order Effect-Compartment Equation

Effect-Compartment Models and Hysteresis · Tutorial 84 of 598

Deriving the First-Order Effect-Compartment Equation

See how balanced first-order exchange produces the effect-compartment equation and what its assumptions imply for effect-site dynamics.

Advanced 10 min read

What You'll Learn

  • Derive the effect-compartment equation from first-order influx and efflux.
  • Define plasma concentration, effect-site concentration, exchange clearance, and the equilibration rate constant.
  • Check the units and signs of the concentration-change equation.
  • Show mathematically why a fixed plasma concentration is the equilibrium target.
  • Recognize assumptions, diagnostic warning signs, and alternatives to the effect-compartment representation.

From First-Order Exchange to an Effect-Site Equation

A plasma concentration can change before a measurable pharmacologic response does. The effect-site concept represents this delay by allowing the concentration associated with the response to equilibrate with plasma over time, rather than requiring an immediate match. The essential modeling question is how to express that gradual adjustment mathematically.

The answer follows from a simple exchange assumption: movement into and out of a notional effect compartment is proportional to the concentration on the corresponding side. If the two exchange processes are balanced so that equal plasma and effect-site concentrations are an equilibrium, a mass balance gives the first-order effect-compartment equation.

Definition: \(C_p\) is the plasma drug concentration, and \(C_e\) is the concentration in the effect compartment, a latent state used to represent concentration equilibration associated with drug effect. \(K_{e0}\) is the first-order equilibration rate constant, with units of inverse time. The first-order effect-compartment equation is
$$ \frac{dC_e}{dt}=K_{e0}(C_p-C_e), \qquad K_{e0}>0. $$
The equation says that the instantaneous rate of change of \(C_e\) is proportional to the difference between plasma and effect-site concentration.

Concentrations must use the same units, such as mg/L, and time must use one consistent unit, such as hours. Thus, if concentrations are in mg/L and time is in hours, \(dC_e/dt\) is in mg/L per hour and \(K_{e0}\) is in h\(^{-1}\). The rate constant is positive by definition; the direction of change comes from the sign of \(C_p-C_e\).

Deriving the Equation by a Mass Balance

Consider a hypothetical effect compartment with fixed volume \(V_e\). Let \(A_e\) be the amount of drug in it, so \(A_e=V_eC_e\). Suppose drug enters from plasma at a rate proportional to \(C_p\), and leaves the effect compartment at a rate proportional to \(C_e\). A convenient representation uses the same exchange clearance \(Q_e\) for both directions. Since clearance has units of volume per time, each concentration multiplied by \(Q_e\) gives an amount-per-time flux.

$$ \text{influx}=Q_eC_p,\qquad \text{efflux}=Q_eC_e. $$

The amount balance is influx minus efflux. With fixed \(V_e\), dividing the balance by \(V_e\) converts the change in amount to the change in concentration:

$$ \frac{dA_e}{dt}=Q_eC_p-Q_eC_e, \qquad A_e=V_eC_e, \qquad \frac{dA_e}{dt}=V_e\frac{dC_e}{dt}. $$
$$ V_e\frac{dC_e}{dt}=Q_e(C_p-C_e) \quad\Longrightarrow\quad \frac{dC_e}{dt}=\frac{Q_e}{V_e}(C_p-C_e). $$

Defining \(K_{e0}=Q_e/V_e\) yields the stated equation. This ratio has units \((\text{volume}/\text{time})/\text{volume}=1/\text{time}\), as required. The balance also shows why the same coefficient multiplies \(C_p\) and \(C_e\): with equal coefficients, influx and efflux are equal when the two concentrations are equal.

Model assumption: The standard equation assumes first-order exchange, a fixed effect-compartment volume, a positive constant \(K_{e0}\), and equal effective exchange coefficients for influx and efflux. It also treats the effect compartment as having negligible influence on the plasma concentration profile. These are modeling assumptions, not consequences of the mass-balance algebra.

The equal-coefficient condition matters. If the amount balance instead used different clearances, \(Q_{\mathrm{in}}C_p-Q_{\mathrm{out}}C_e\), its zero-rate condition would be \(C_e=(Q_{\mathrm{in}}/Q_{\mathrm{out}})C_p\), not necessarily \(C_e=C_p\). The familiar equation is therefore a particular first-order exchange model chosen to make plasma concentration the equilibrium target.

Direction of Change and the Equilibrium Target

The sign of the derivative follows directly from the concentration difference. When \(C_p>C_e\), the net exchange is into the effect compartment, so \(dC_e/dt>0\). When \(C_p<C_e\), net exchange is outward, so \(dC_e/dt<0\). When \(C_p=C_e\), the two exchange fluxes balance and \(dC_e/dt=0\). These statements describe the instantaneous direction; they do not say that \(C_e\) instantly equals \(C_p\).

The attraction toward plasma concentration can be justified without solving the differential equation. First hold \(C_p\) fixed and define the discrepancy \(e=C_e-C_p\). Because \(C_p\) is constant, \(de/dt=dC_e/dt\). Substitution into the model gives \(de/dt=-K_{e0}e\). Now consider the nonnegative squared discrepancy \(V=e^2/2\):

$$ \frac{dV}{dt} =e\frac{de}{dt} =e(-K_{e0}e) =-K_{e0}e^2\leq 0. $$

Since \(K_{e0}>0\), the squared discrepancy decreases whenever \(e\ne0\), and its derivative is zero at \(e=0\). Therefore the equilibrium \(C_e=C_p\) is stable: a concentration above the fixed plasma value is driven downward, and one below it is driven upward. This establishes the direction and stability of the target without deriving a time-course solution.

If plasma concentration is changing, the same sign rule still applies at each instant, but the target is moving. The effect-site concentration tends to lag behind a rising plasma concentration and may remain above plasma during a decline. Thus, “drives toward plasma” describes the instantaneous exchange tendency; it does not guarantee that \(C_e\) is always moving toward a fixed value or that the two concentrations are equal during changing exposure.

Worked Example: Deriving the Rate Constant From Exchange Parameters

Suppose a hypothetical effect compartment has \(V_e=0.30\) L and an exchange clearance \(Q_e=0.12\) L/h. At a particular instant, \(C_p=5.0\) mg/L and \(C_e=2.0\) mg/L. Find \(K_{e0}\), then calculate the instantaneous rate of change of effect-site concentration.

From the mass-balance derivation, \(K_{e0}=Q_e/V_e\). Substituting the values gives \(K_{e0}=(0.12\ \text{L/h})/(0.30\ \text{L})=0.40\ \text{h}^{-1}\).

The concentration equation gives \(dC_e/dt=K_{e0}(C_p-C_e)\). Therefore, \(dC_e/dt=(0.40\ \text{h}^{-1})(5.0-2.0\ \text{mg/L})=1.2\ \text{mg/L per hour}\).

The positive derivative is consistent with \(C_p>C_e\): effect-site concentration is increasing at that instant. The flux calculation checks the same result. Influx is \((0.12\ \text{L/h})(5.0\ \text{mg/L})=0.60\ \text{mg/h}\), and efflux is \((0.12\ \text{L/h})(2.0\ \text{mg/L})=0.24\ \text{mg/h}\). The net gain is \(0.36\ \text{mg/h}\); dividing by \(0.30\) L gives \(1.2\ \text{mg/L per hour}\).

Worked Example: Checking the Sign During a Concentration Decline

In a hypothetical model, \(K_{e0}=0.25\ \text{h}^{-1}\), plasma concentration at the time of interest is \(3.0\) mg/L, and effect-site concentration is \(4.2\) mg/L. Calculate the instantaneous derivative and interpret its direction.

Substitution gives \(dC_e/dt=(0.25\ \text{h}^{-1})(3.0-4.2\ \text{mg/L})=(0.25)(-1.2)\ \text{mg/L per hour}=-0.30\ \text{mg/L per hour}\).

The negative sign means \(C_e\) is decreasing at that instant. Because \(C_e\) exceeds \(C_p\), the model gives greater outward than inward flux. It does not imply that the effect-site concentration immediately becomes 3.0 mg/L; it specifies only the current rate of change.

Worked Example: Recovering the Rate Constant From an Instantaneous Slope

Suppose a hypothetical observation provides \(C_p=3.2\) mg/L, \(C_e=4.0\) mg/L, and an instantaneous effect-site slope \(dC_e/dt=-0.24\) mg/L per hour. Under the first-order equation, calculate the implied \(K_{e0}\).

Rearrange the model, provided \(C_p\ne C_e\): \(K_{e0}=(dC_e/dt)/(C_p-C_e)\). Substituting gives \(K_{e0}=(-0.24\ \text{mg/L per hour})/(3.2-4.0\ \text{mg/L})=(-0.24)/(-0.8)\ \text{h}^{-1}=0.30\ \text{h}^{-1}\).

The inferred rate constant is positive, and the negative slope is consistent with \(C_e>C_p\). This calculation is an algebraic check under the model, not a general estimation procedure: measured concentration profiles may be noisy, and a slope inferred from sparse observations may not represent an instantaneous derivative.

Assumptions, Diagnostics, and Alternatives

The equation is useful because it reduces a delayed concentration state to one differential equation and one rate parameter. Its simplicity should not be confused with universal validity. It is most defensible when a single first-order equilibration process adequately describes the delay over the observed concentration range and time scale.

Several checks can expose a mismatch between the model and the data or implementation:

  • Check the sign. If \(C_p>C_e\), the computed derivative must be positive; if \(C_p<C_e\), it must be negative. A contrary result often indicates a reversed concentration difference or an inconsistent sign convention.
  • Check units. \(K_{e0}\) must have units of inverse time. Multiplying it by a concentration difference must produce concentration per time.
  • Check the equilibrium implication. Under constant plasma concentration, the model predicts no net change at \(C_e=C_p\). A model intended to target equality but using unequal exchange coefficients is not the same equation.
  • Check whether one rate is adequate. Systematic evidence of multiple delays, changing equilibration behavior, or persistent discrepancies may indicate that a single first-order state is too restrictive.
  • Check whether the plasma input is represented appropriately. Because \(C_e\) responds to the supplied \(C_p\), errors or inappropriate interpolation in plasma concentration can propagate into the effect-site prediction.

The alternatives answer different modeling needs. An instantaneous concentration-effect model sets the effect-driving concentration equal to plasma concentration and has no equilibration delay. An explicit distribution compartment can represent a more mechanistic exchange structure when corresponding concentration or distribution information is available. Transit-compartment models can represent a sequence of delays, while indirect-response models describe delayed changes in response through turnover of the response variable. These alternatives are not interchangeable: the choice depends on whether the observed delay is best represented as concentration equilibration, multiple transit stages, or delayed response dynamics.

The effect compartment is often a latent modeling device rather than a directly sampled anatomical space. Consequently, \(K_{e0}\) summarizes equilibration behavior within the chosen model; it should not automatically be interpreted as a measured tissue transport constant. The equation also assumes linear exchange and does not, by itself, specify how \(C_e\) produces a response. A concentration-effect relationship is a separate model component.

Common pitfall: “Because \(C_e\) moves toward \(C_p\), the two concentrations must be equal at every time.” This is incorrect. Equality is the equilibrium condition when plasma concentration is held fixed. During changing plasma exposure, the effect-site concentration generally responds with a delay, and its instantaneous direction is determined by the current difference \(C_p-C_e\).

Check Your Understanding

Answer each question, then compare your reasoning with the solutions.

  1. If \(K_{e0}=0.50\ \text{h}^{-1}\), \(C_p=6.0\) mg/L, and \(C_e=4.0\) mg/L, what is \(dC_e/dt\), including units?
  2. For the equation in this tutorial, what condition makes \(C_e=C_p\) an equilibrium?
  3. Why does \(K_{e0}\) have units of inverse time when it is derived as \(Q_e/V_e\)?
  4. If plasma concentration is falling and \(C_e>C_p\), what is the sign of the effect-site derivative?
  5. Give one reason the standard first-order effect-compartment equation might not adequately describe a delayed response.

Solutions. (1) \(dC_e/dt=(0.50)(6.0-4.0)=1.0\ \text{mg/L per hour}\), so \(C_e\) is increasing. (2) The equation has zero derivative when \(C_p-C_e=0\), provided \(K_{e0}>0\); equal plasma and effect-site concentrations then have balanced exchange. (3) \(Q_e/V_e\) has units \((\text{volume}/\text{time})/\text{volume}=1/\text{time}\). (4) The derivative is negative because \(K_{e0}>0\) and \(C_p-C_e<0\). (5) A single first-order state may be inadequate if the data indicate multiple delay stages or if response delay reflects turnover rather than concentration equilibration.

Clinical Biostats